删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

离散可积系统: 多维相容性

本站小编 Free考研考试/2021-12-29

摘要:对比已有完善而系统理论的微分方程领域, 差分方程理论尚处于发展之中. 近年来离散可积理论的进展, 带来了差分方程理论的革命. 多维相容性是伴随离散可积系统研究出现的新的概念, 作为对离散可积性的一种理解, 提供了构造离散可积系统的B?cklund变换、Lax对和精确解的工具. 本文旨在综述多维相容性的概念及其在离散可积系统研究中的应用.
关键词: 离散可积系统/
多维相容性

English Abstract


--> --> -->
离散系统泛指含有离散自变量的常差分、微分差分、偏差分系统、以及变换和映射等. 由于缺少导数、积分等局部化的数学工具, 对于非线性离散系统的研究, 往往伴随着新的数学概念、理论和方法的出现.
现代可积理论兴起于20世纪60年代中期孤立子的命名[1]、反散射变换方法的建立[2]和Lax对概念的提出[3], 参见文献[4]. 早期对离散可积系统的探索主要是可积离散化, 早在20世纪70年代: Case和Kac[5]对Schr?dinger谱问题的离散以及Ablowitz和Ladik[68]对Ablowitz-Kaup-Newell-Segur(AKNS) 谱问题的离散 建立了基于差分的反散射变换过程; Hirota[9-11]对双线性方程的离散利用双线性B?cklund变换与Lax对的联系获得了一系列离散可积系统. 进入80年代后, 对离散可积系统的研究逐渐向形成独立于连续系统的研究方法, 系统进展有: 京都学派Date等[12-16]和?Ueno等[17]对Sato理论的离散; 荷兰****Nijhoff等[18-22]和Quispel等[23]基于Fokas和Abolwitz[24]的直接线性化格式以及Levi和Benguria[25]的变换与离散的同等性认识发展起来的构造和研究离散可积系统的系统方法. 进入90年代以后, 离散可积系统的显著进展包括: 超离散可积系统的提出及其连续极限的建立[26,27], 基于奇点囿禁的奇点理论与可积性的联系的发现[28,29], 奇点囿禁在可积性判别上的不充分性的发现[30], 以及基于代数熵理论对可积性的判别[31]. 这些进展相继推动了超离散可积系统、离散Painlevé方程、可积性检验等方面的发展.
进入新世纪后, 离散可积系统继续迎来新的发展. Sakai[32]基于有理曲面理论和Blow-up分析对离散Painlevé方程的分类, 揭示了离散Painlevé方程丰富的代数几何结构, Bobenko和Suris[33]、Adler等[34]和Nijhoff等[35,36]****对于“多维相容性”的理解以及对若干离散可积系统的分类, 各种精确求解方法 在离散可积系统中相继实现[3745], 等等, 一系列进展标志着对离散可积系统的研究进入到一个新的阶段. 2009年在著名的英国剑桥牛顿所(The Isaac Newton Institute for Mathematical Sciences)举办的为期半年的离散可积系统主题研讨活动, 是对当时离散可积系统蓬勃发展的一个反映.
可积系统与数学和物理的众多分支都有联系, 已经渗透到数学物理的各个领域[4]. 在数学方面, 离散可积系统的发展伴随着离散几何和离散复分析的发展, 这与连续意义下可积系统丰富的几何背景以及分析的工具在其研究中的重要作用是分不开的. 对与离散可积系统相关的差分算子理论与复分析方面的研究, 如势函数的渐近性、散射理论、Riemann-Hilbert问题、解的长时间渐近分析、初边值问题等等, 建议参考文献[42,46]. 除了提及的上述进展, 还有如离散几何、对称理论、Galois理论、Lagrangian多形式理论等许多方面的重要进展. 值得指出的是, 最近由Hietarinta, Joshi和Nijhoff [47]完成的一部专著《Discrete Systems and Integrability》已经出版. 读者可以参考了解其更多内容. 国际上每两年一届的SIDE会议(Symmetries and Integrability of Difference Equations)是离散可积系统及其相关问题的首要国际会议, 1994年始于加拿大Montreal, 至今已经成功举办13届, 其中SIDE-10于2012年在我国宁波召开.
本文将在第二、三节重点介绍多维相容性的概念及其应用. 希望通过具体的例子让更多读者了解离散可积系统及其中的概念和方法.
离散可积系统从未独立于连续系统. 下面首先来介绍两种引入离散变量的途径.
2
2.1.离散变量的引入
-->离散化是引入离散变量的方式之一. 从熟悉的AKNS谱问题开始:
$ { {\varPhi}}_x ={ {M \varPhi}},\; \; { M} = \left(\!\!\!\begin{array}{cc}\eta & u \\ v & -\eta\end{array}\!\!\!\right),\; \; \; { {\varPhi}} = (\phi_1, \phi_2)^{\rm T}, $
其中$ u = u(x, t), \; v = v(x, t) $, $ \eta $是谱参数, 为了方便, 记$ { M} = { M}(\eta, U) $, $ U = (u, v)^{\rm T} $. 定义
$ { {\varPhi}}_{n+j} = { {\varPhi}}(x+j\epsilon), $
然后利用差分$ ({ {\varPhi}}_{n+1}-{ {\varPhi}}_n)/\epsilon $替代导数$ { {\varPhi}}_x $, 得到
$ \begin{split} & { {\varPhi}}_{n+1} = \left(\!\!\!\begin{array}{cc} 1+\epsilon \eta & \epsilon u \\\epsilon v & 1- \epsilon\eta \!\!\!\end{array} \right) { {\varPhi}}_n,\\ &{ {\varPhi}}_n = (\phi_{1,n}, \phi_{2,n})^{\rm T}. \end{split}$
Ablowitz-Ladik (AL)谱问题指[7]
$ \begin{split} &{ {\varPhi}}_{n+1} = \left(\!\!\!\begin{array}{cc} \lambda & Q_n \\ R_n & 1/\lambda \end{array}\!\!\!\right) { {\varPhi}}_n,\\&{ {\varPhi}}_n = (\phi_{1,n}, \phi_{2,n})^{\rm T}. \end{split}$
在(2)式和
$ (Q_n, R_n) = \epsilon (u,v),\; \; \; \lambda = {\rm e}^{\epsilon \eta} $
之下, 对$ { {\varPhi}}_{n+1} = { {\varPhi}}(x+\epsilon) $$ \epsilon = 0 $展开, 并取$ \epsilon \to 0 $, 则AKNS谱问题(1)式可以作为领头项从AL谱问题(4)式中恢复出来, 同时谱参数从$ \lambda $-平面变到$ \eta $-平面.
(2)式是我们所熟悉的数值(差分)离散: 如图1, 将区间$ [x_0, x] $等分成n份, 步长为$ \epsilon $. 在x点, $ {\varPhi}_n $$ {\varPhi}(x) $仍然表示相同的数值, 但是自变量的空间已经从实数域$ \mathbb{R} $变为整数域$ \mathbb{Z} $. 差分离散是计算连续极限的基础, 但是对于可积系统而言, 差分离散(也称为直接离散)不足以保持原有系统的可积特征. 关于由AL谱问题(4)式引出的半离散AKNS系统, 读者可以参考文献[4851].
图 1 $[x_0, x]$上的数值离散
Figure1. Numerical discretisation on $[x_0, x]$

递推关系可以很自然地视为离散的方程. 例如Hermite正交多项式$ \{H_n(x)\} $满足递推关系
$ H_{n+1}(x) = 2xH_n(x)-2nH_{n-1}(x). $
第I型Bessel函数
$ {J_\alpha }(x) = {\left( {\frac{x}{2}} \right)^\alpha }\sum\limits_{k = 0}^\infty {\frac{{{{( - {x^2}/4)}^k}}}{{k!{\mkern 1mu} \Gamma (\alpha + k + 1)}}} $
满足方程常微分方程
$ x^2 w''+xw'+(x^2-\alpha^2)w = 0, $
其中x是自变量, $ \alpha $是参数. 由此方程可以建立Bessel函数的递推关系(参考文献[47]的第37页)
$ x(w_{\alpha+1}+w_{\alpha-1})-2\alpha w_{\alpha} = 0, $
此时x$ \alpha $互换了角色, x是参数, $ \alpha $成为自变量. 再如, Painlevé II方程
$ f{''}(t) = 2f^3+tf-\alpha $
的解满足递推关系
$ f_{\alpha+1}(t) = -f_{\alpha}(t)-\dfrac{\alpha+\dfrac{1}{2}}{f^{'}_{\alpha}(t)-f^2_{\alpha}(t)-\dfrac{t}{2}}, $
这里$ f'(t) $表示对t的导数.
在上面Bessel函数和Painlevé II方程的解的例子中, 递推关系分别来自于两个微分方程的解之间的变换, 可以视为通过变换建立起来的解之间的叠加关系. 这是由连续系统到离散系统的常见过程. 尤其, KdV方程的B?cklund变换的非线性叠加公式提供了对KdV方程的完美的离散化.
对于著名的KdV方程,
$ u_t+6uu_x+u_{xxx} = 0, $
Wahlquist和Estabrook[52]在1973年发现形如
$ (\widetilde w+w)_x = 2\lambda-\frac{1}{2}(\widetilde w -w)^2, \tag{7a}$
$ (\widetilde w-w)_t = \frac{1}{2}[(\widetilde w-w)^3]_x-6\lambda (\widetilde w-w)_{x}-(\widetilde w-w)_{xxx}\tag{7b} $
的非线性B?cklund变换, 其中w满足势KdV方程
$ w_t+3(w_x)^2+ w_{xxx} = 0, $
$ u = w_x $满足KdV方程(6).
利用B?cklund变换(7a)式可以建立KdV方程的解的非线性叠加公式. 首先在(7a)式中, 从同样的种子解w出发, 分别记由$ \lambda = \lambda_1 $$ \lambda = \lambda_2 $引出的$ \widetilde w $$ w_1 $$ w_2 $, 即
$ (w_1+w)_x = 2\lambda_1-\frac{1}{2}(w_1 -w)^2,\tag{9a}$
$ (w_2+w)_x = 2\lambda_2-\frac{1}{2}(w_2 -w)^2.\tag{9b}$
接下来, 在(7a)式中取$ w = w_1, \; \lambda = \lambda_2 $, 记$ \widetilde w = w_{12} $, 有
$ (w_{12}+w_1)_x = 2\lambda_2-\frac{1}{2}(w_{12} -w_1)^2;\tag{10a}$
$ w = w_2, \; \lambda = \lambda_1 $, 记$ \widetilde w = w_{21} $, 有
$ (w_{21}+w_2)_x = 2\lambda_1-\frac{1}{2}(w_{21} -w_2)^2.\tag{10b}$
上述过程可以描述为图2.
图 2 B?cklund变换解的交换性质
Figure2. Permutation property of B?cklund transformation

可以证明$ w_{12} $$ w_{21} $能够相等(参考文献[49]). 进一步, 从(9)式和(10)式得到
$4(\lambda_1-\lambda_2) = (w_1-w_2)(w_{12}-w),$
称为(势)KdV方程解的非线性叠加公式, 也称为Bianchi等式, 还称为离散的势KdV方程[47]. 作为离散的方程时, (11)式通常写为
$ (w_{n+1,m}-w_{n,m+1})(w_{n,m}-w_{n+1,m+1}) = p-q, $
其中pq分别是对应于n–方向和m–方向的方向参数.
modified KdV(mKdV)方程
$ u_t = 6u^2u_x -u_{xxx} $
和sine-Gordon方程
$ \theta_{xx}-\theta_{tt} = \sin \theta $
拥有形式相同的非线性叠加公式(文献[53,54]):
$ \tan \frac{\theta_{12}-\theta}{4} = \frac{p+q}{p-q}\tan \frac{\theta_{1}-\theta_2}{4}, $

$ p \sin \frac{\theta_{12}-\theta_1+\theta_2-\theta}{4} = q\sin \frac{\theta_{12}+\theta_{1}-\theta_2-\theta}{4}. $
此方程经过变换$ u = {\rm e}^{\frac{{\rm i}\theta}{2}} $以后, 可以写为
$ p(u_{12} u_2 - u_1 u) = q (u_{12} u_1 - u_2 u), $
此方程也称为离散势mKdV方程. 有意思的是, 离散的sine-Gordon方程形如
$ p \sin \frac{\varphi_{12}+\varphi-\varphi_1-\varphi_2}{4} = q\sin \frac{\varphi_{12}+\varphi+\varphi_{1}+\varphi_2}{4}, $
由Hirota[11]和Orfanidis[55]先后得到, 与(16)式只有部分符号差别, 经过同样的变换$ u = {\rm e}^{\textstyle\frac{{\rm i}\varphi}{2}} $以后, 可以写为[18]
$ p(u_{12} u - u_1 u_2) = q(u u_1u_2 u_{12} -1). $
2阶AKNS方程组
$ u_t = u_{xx}-2u^2v,\; \; v_t = -v_{xx}-2v^2u $
的B?cklund变换的非线性叠加公式形如
$ (u_1-u_2)(uv_{12}+1)+(p-q)u = 0, \tag{21a}$
$ (v_1-v_2)(uv_{12}+1)-(p-q)v_{12} = 0, \tag{21b}$
由Konopelchenko[58]于1982年获得.
B?cklund变换提供了引入离散变量的一种方式. (9a)式中$ w_1 $可视为w沿一个方向上的平移, (9b)式中$ w_2 $可视为w沿另一个方向上的平移. 非线性叠加公式(11)式是这两个方向上平移相容性的结果. 同时, (9a)式作为一个独立方程, 也可视为一个微分-差分方程.
B?cklund变换应用于非线性模型, Darboux变换的目标对象则是在线性层面(Lax对). 对于给定的连续谱问题, 存在与之相应的连续等谱发展方程族; 它的Darboux变换可视为一个离散的谱问题, 与原有连续谱问题之间的相容性引出的微分差分方程 自然给出连续等谱发展方程族的解之间的B?cklund变换. 这一漂亮的联系首先由Levi和Benguria[25,59]发现, 对于研究和理解离散可积系统具有重要意义. AKNS谱问题(1)式存在如下Darboux变换[60]:
$ \widetilde { {\varPhi}} = { T}{ {\varPhi}},\; { T} = { T}(\gamma,{ U},\widetilde { U}) = \left(\!\!\!\begin{array}{cc} 2(\eta-\gamma) + u \widetilde v & u \\ \widetilde v & 1 \end{array}\!\!\!\right), $
式中$ { U} = (u, v)^{\rm T} $. 作为Darboux变换, $ \widetilde { {\varPhi}} $满足AKNS谱问题(1)式, 即$ \widetilde { {\varPhi}} _x = \widetilde { M} { {\varPhi}} $, $ \widetilde { M} = { M}(\eta, \widetilde { U}) $, 这相当于要求(1)式与(22)式相容, 即$ (\widetilde { {\varPhi}} )_x = (\widetilde{{ {\varPhi}}}_x) $, 由此有
$ { T}_x-\widetilde { M}{ T}+{ {TM}} = 0, $
引出的方程
$ u_{x} = \widetilde u- u- u^2 \widetilde v,\; \; v_{x} = v - \underset{\sim }{\mathop{v}}\, + \underset{\sim }{\mathop{u}}\, v^2, $
可视为连续的等谱AKNS方程族的B?cklund变换. 作为离散的谱问题, (22)式引出一个新的半离散AKNS族[61,62], 可视为微分-差分Kadomtsev-Petviashvili (KP) 方程族Lax三重组的对称约束的结果[63]. 若考虑两个不同方向上的Darboux变换:
$ \widetilde { {\varPhi}} = { T}(p/{2},{ U},\widetilde { U}){ {\varPhi}},\; \; \widehat{{ {\varPhi}} } = { T}(q/{2},{ U},{\widehat{{ U} }}){ {\varPhi}}, $
作为离散的Lax对, 其相容性得到方程(21), 且在适当的连续极限下, (21)式引出(20)式(参考文献[43]). 显然, 由Darboux变换作为Lax对引出的4点方程$ \, { Q}({ U}, \widetilde { U}, {\widehat{ U }}, {\widehat{\widetilde{ U }}}; p, q) = 0 $即为4个解之间的非线性叠加公式. 更多的关于Darboux变换作为离散谱问题的例子可见文献[64].
当作为离散的方程时, 变换的符号意义可以理解为
$\begin{split} & u = u_{n,m},\\ & \widetilde u = u_{n+1,m},\\ &{\widehat{u }} = u_{n,m+1}, \\ & {\widehat{\widetilde{u }}} = u_{n+1,m+1}, \end{split}$
变换引入的参数pq分别表示n–方向和m–方向上的链参数. 如图3所示.
图 3 变换与平面网格
Figure3. Map and lattice

2
2.2.多维相容性
-->对于离散可积系统的研究需要引入新的概念. 下面介绍离散可积系统的“多维相容性”的概念.
回到(势)KdV方程解的非线性叠加公式(12)式, 这一叠加关系可以重复下去. 为了方便,我们采用(26)式中的符号, 将(12)式写为
$ (u-{\widehat{\widetilde{u }}})(\widetilde u- {\widehat{u }}) = p-q. $
这是一个定义在平面网格上的四方格方程. 我们引入第三个方向l, 该方向上的平移表示为$ \, {\overline{u}} = u_{n, m, l+1} $, 方向参数为r, 如图4所示.
图 4 相容立方体
Figure4. Consistent cube

观察方程在$ (^{\widetilde{\; \; \;}}, p) $$ (^{\widehat{\; \; \;}}, q) $互换下保持不变, 要求$ (\, \overline{\; \; }, r) $与另外两个方向之间也存在对称性, 则有
$ (u-{\widetilde{\overline{u}}})(\widetilde u- {\overline{u}}) = p-r,\tag{28a}$
$ (u-{\widehat{\overline{u}}})({\overline{u}}- {\widehat{u }}) = r-q,\tag{28b}$
视(27)式定义在图4中立方体的底面, (28a)式和(28b)式则分别定义在左侧和后侧. 如何体现这种维数扩充后的相容性? 首先, 由上述3个方程有
$ {\widehat{\widetilde{u }}} = \frac{p-q}{\widehat{u}-\widetilde{u}} + u,\tag{29a}$
$ {\widehat{\overline{u}}} = \frac{r-q}{{\widehat{u }}-{\overline{u}}}+u,\tag{29b}$
$ {\widetilde{\overline{u}}} = \frac{p-r}{{\overline{u}}-\widetilde u}+u,\tag{29c}$
引入方向l以后的相容性体现在$\overline{ (29{\rm{a}})} = \widetilde{ (29{\rm{b}})} = $$ \widehat{ (29{\rm{c}})} $. 将(29a)式和(29b)式代入到$ \widehat{ (29{\rm{c}})} $算得
$ \widehat{\widetilde{\overline{u}}} = \frac{(q-p)\widetilde u{\widehat{u }}+(r-q){\widehat{u }}{\overline{u}}+(p-r){\overline{u}}\widetilde u}{(q-r)\widetilde u+(p-q){\overline{u}}+(r-p){\widehat{u }}}. $
此式右端在任意交换$ (^{\widetilde{\; \;\; }}, p) $$ (^{\widehat{\; \;\; }}, q) $$ (\, \overline{\phantom{M}}, r) $下都保持不变, 说明从$ \overline{ (29{\rm{a}})} $$ \widetilde{ (29{\rm{b}})} $也将得到同样的$ \widehat{\widetilde{\overline{u}}} $.
方程(27)的这种性质称~为该方程的3维相容性, 也称为立方体相容(consistent-around-the-cube (CAC)), 它体现了(27)式自身的性质. 方程(17)也具有同样的性质. 一般地, 对于一个定义在平面四方格上的方程
${ Q}(u,\widetilde u,{\widehat{u }},\widehat{\widetilde u};p,q) = 0, $
将其嵌入到图4中立方体的6个面中, 得到
$ \begin{split} & { Q}(u,\widetilde u,{\widehat{u }},\widehat{\widetilde u};p,q) = 0, \;\; { Q}({\overline{u}},\overline{\widetilde u},\overline{{\widehat{u }}},\overline{\widehat{\widetilde u}};p,q) = 0,\\ & { Q}(u,\widetilde u,{\overline{u}},\overline{\widetilde u};p,r) = 0, \;\; { Q}({\widehat{u }},\widehat{\widetilde u},\overline{{\widehat{u }}},\overline{\widehat{\widetilde u}};p,r) = 0, \\ & { Q}(u,{\widehat{u }},{\overline{u}},\overline{{\widehat{u }}};q,r) = 0, \;\;{ Q}(\widetilde u,\widehat{\widetilde u},\overline{\widetilde u},\overline{\widehat{\widetilde u}};q,r) = 0. \end{split} $
如果从上式左侧的3个方程中分别解出$ {\widehat{\widetilde{u }}} $, $ {\widetilde{\overline{u}}} $$ {\widehat{\overline{u}}} $, 代入到右侧的3个方程后得到同样的$ \overline{\widehat{\widetilde u}} $, 则称方程(30)具有3维相容性.
对于离散系统而言, 上述多维相容性将提供一系列的可积特征, 后文将进一步介绍. 对于一个d-维的离散方程, 如果嵌入到一个$ (d+1) $-维链后, 所有的d-维子链都是相容的, 则称原d-维离散方程具有$ (d+1) $-维相容性. 对于2维离散方程, 这种多维相容性即为3维相容性, 或CAC性质.
2
2.3.ABS方程
-->本世纪初, 多维相容性逐渐被系统地认识并作为工具应用到离散系统的研究中[65,35,36,33]. 2003年, Adler, Bobenko和Suris (ABS)[34]发表了他们在多维相容性的基础上对四方格方程的分类, 他们得到的方程列表被统称为ABS方程.
ABS在假设方程(30)具有3维相容性的基础上, 进一步要求(30)式满足:
(i) 仿线性(affine linear),
(ii) 在正方形对称群$ D_4 $下不变,
(iii) 四面体性质(tetrahedron property).
在前两条要求下方程(30)的一般形式可以设为
$ \begin{split} \,& k_0 u\widetilde u{\widehat{u }}{\widehat{\widetilde{u }}} + k_1(u\widetilde u {\widehat{u }}+ u \widetilde u {\widehat{\widetilde{u }}} + u{\widehat{u }} {\widehat{\widetilde{u }}}+\widetilde u{\widehat{u }}{\widehat{\widetilde{u }}}) \\ & ~~ + k_2(\widetilde u{\widehat{u }}+u {\widehat{\widetilde{u }}}) + k_3(u\widetilde u+{\widehat{u }}{\widehat{\widetilde{u }}})+ k_4(u{\widehat{u }}+\widetilde u {\widehat{\widetilde{u }}})\\ & ~~+k_5(u+\widetilde u+{\widehat{u }}+{\widehat{\widetilde{u }}})+k_6 = 0,\\[-10pt]\end{split} $
其中$ k_i $都是待定系数; 四面体性质指$ \widehat{\widetilde{{\overline{u}}}} $只与$ {\widehat{\widetilde{u }}} $$ {\widetilde{\overline{u}}} $$ {\widehat{\overline{u}}} $有关, 而与u无关, 即($ \widehat{\widetilde{{\overline{u}}}} $, $ {\widehat{\widetilde{u }}} $, $ {\widetilde{\overline{u}}} $, $ {\widehat{\overline{u}}} $)四个点在立方体(图4)中构成一个四面体. ABS证明了满足上面3条要求的3维相容方程(30)只有9个(允许存在M?bius变换), 并分别命名为H1, H2, H3(δ), A1(δ), A2, Q1(δ), Q2, Q3(δ), Q4:
$ {\rm H}1: \ \ (u-\widehat{\widetilde{u}})(\widetilde{u}-\widehat{u})-p+q = 0,\tag{33a}$
$ {\rm H}2: \ \ (u-\widehat{\widetilde{u}})(\widetilde{u}-\widehat{u})-(p-q)(u+\widetilde{u}+\widehat{u}+\widehat{\widetilde{u}}+p+q) = 0, \tag{33b}$
$ {\rm H}3(\delta): \ \ p(u\widetilde{u}+\widehat{u}\widehat{\widetilde{u}})-q(u\widehat{u}+\widetilde{u}\widehat{\widetilde{u}})+\delta(p^2-q^2) = 0, \tag{33c}$
$ {\rm A}1(\delta): \ \ p(u+\widehat{u})(\widetilde{u}+\widehat{\widetilde{u}})\!-\!q(u+\widetilde{u})(\widehat{u}+\widehat{\widetilde{u}})-\delta^2pq(p-q) \!=\! 0,\tag{33d}$
$ \begin{split}{\rm A}2: \ \ p(1\, & -q^2)(u\widetilde{u}+\widehat{u}\widehat{\widetilde{u}})-q(1-p^2)(u\widehat{u}+\widetilde{u}\widehat{\widetilde{u}})\\ & -(p^2-q^2)(1+u\widetilde{u}\widehat{u}\widehat{\widetilde{u}}) = 0,\\[-10pt]\end{split}\tag{33e} $
$ \begin{split}{\rm Q}1(\delta): \ \ p(u\, &-\widehat{u})(\widetilde{u}-\widehat{\widetilde{u}})-q(u-\widetilde{u})(\widehat{u}-\widehat{\widetilde{u}})\\ &+\delta^2pq(p-q) = 0,\end{split} \tag{33f}$
$ \begin{split} {\rm Q}2: \ \ & p(u-\widehat{u})(\widetilde{u}-\widehat{\widetilde{u}})-q(u-\widetilde{u})(\widehat{u}-\widehat{\widetilde{u}}) \\ +\, &pq(p-q)(u+\widetilde{u}+\widehat{u}+\widehat{\widetilde{u}}-p^2+pq-q^2) = 0,\end{split} \tag{33g}$
$ \begin{split}{\rm Q}3(\delta):\ \ \, & p(1-q^2)(u\widehat{u}+\widetilde{u}\widehat{\widetilde{u}})-q(1-p^2)(u\widetilde{u}+\widehat{u}\widehat{\widetilde{u}}) \\ -\, &(p^2\!-\!q^2)\left(\widetilde{u}\widehat{u}\!+\! u\widehat{\widetilde{u}}\!+\!\frac{\delta^2(1\! - \!p^2)(1\!-\!q^2)}{4pq}\right) \!=\! 0,\end{split}\tag{33h}$
$ \begin{split}{\rm Q}4: \ \ &{\rm{sn}}(p)(u\widetilde{u}+\widehat{u}\widehat{\widetilde{u}})-{\rm{sn}}(q)(u\widehat{u}+\widetilde{u}\widehat{\widetilde{u}})\\-\,&{\rm{sn}}(p-q)(\widetilde{u}\widehat{u}+u\widehat{\widetilde{u}} ) \\ +\,&{\rm{sn}}(p){\rm{sn}}(q){\rm{sn}}(p-q)(1+k^2u\widetilde{u}\widehat{u}\widehat{\widetilde{u}}) = 0, \end{split} \tag{33i}$
其中, H1即为离散势KdV方程(27), H3(δ = 0)为离散势mKdV方程(17) ($ u \to i^{n+m} u $), Q1(δ = 0)为离散的Schwarzian KdV方程, 也称为交比(cross-ratio)方程, A1(δ)在$ u\to (-1)^{n+m} u $下即为Q1(δ), A2在$ u\to u^{(-1)^{n+m}} $下为Q3(δ = 0), Q4是著名的Krichever-Novikov方程的B?cklund变换的非线性叠加公式, 原型中系数用Weierstrass椭圆函数表示[34,66], (33i)中Jacobi椭圆函数的参数化形式来自于Hietarinta[67], 式中k为Jacobi椭圆函数的模(modulus).
2
2.4.离散Boussinesq型方程
-->离散的势Boussinesq方程是平面9点方程, 定义在$ 3\times 3 $格子上:
$ \begin{split} &\frac{P-Q}{x_{11}-x_{20}}-\frac{P-Q}{x_{02}-x_{11}}\\ =\, & (x_{00}-x_{21})(x_{10}-x_{01})\\ & -(x_{01}-x_{22})(x_{21}-x_{12}),\end{split}$
式中$ P, Q $分别表示$ n, m $的方向参数, $ x_{ij} = x_{n+i, m+j} $. 这类方程早期由Nijhoff等[68-70]****构造并研究. 借助于其他(因)变量, (34)式可以改写为(参考文献[65])
$ {\widetilde{y}} - x{\widetilde{x}}+z = 0,\; \; {\widehat{y}} - x\widehat{x}+z = 0, \tag{35a}$
$ y - x\widehat{\widetilde{x}}+\widehat{\widetilde{z}}-\frac{P-Q}{\widetilde{x}-\widehat{x}} = 0. \tag{35b}$
取向量$ { u} = (x, y, z)^{\rm T} $, 方程可以视为向量意义下的四方格系统, 并具有3维相容性.
2011年, Hietarinta[71]在(35)式形式的基础上, 寻找了可能的离散Boussinesq型的3维相容系统, 代表方程为(考虑到对称性和M?bius变换):
$ [{\rm B}-2]:\; \; \widetilde{y} = x \widetilde{x}-z, \; \; \widehat{y} = x \widehat{x}-z, \tag{36a}$
$ y = x \widehat{ \widetilde{x}}- \widehat{ \widetilde{z}} +b_0( \widehat{ \widetilde{x}}-x)+b_1+\frac{P-Q}{ \widetilde{x}- \widehat{x}},\tag{36b}$
$ [{\rm A}-2]:\; \;\widetilde{y} = z \widetilde{x}-x, \; \; \widehat{y} = z \widehat{x}-x, \tag{37a}$
$ y = x \widehat{ \widetilde{z}}-b_0 x+\frac{P \widetilde{x}-Q \widehat{x}}{ \widehat{z}- \widetilde{z}},\tag{37b}$
$ [{\rm C}-3]:\; \; {\widetilde{y}} \,z = {\widetilde{x}} - x,\; \; {\widehat{y}} \,z = {\widehat{x}} - x, \tag{38a}$
$ \widehat{\widetilde{z}}\,y = b_0\,x + b_1 + z \frac{P\,{\widetilde{y}}\,{\widehat{z}}-Q\,{\widehat{y}}\,{\widetilde{z}}}{{\widetilde{z}}-{\widehat{z}}},\tag{38b} $
式中$ P, Q $仍然分别表示$ n, m $的方向参数, $ b_i $都是任意常数, 有的$ b_i $可以通过适当的变换去掉. 显然, [B-2]是(35)式的推广. 文献[71]中还得到[C-4]方程, 现已证明与[C-3]之间存在M?bius变换, 故不再列出.
2
2.5.不同方程间的3维相容性
-->离散的sine-Gordon方程(19)不是3维相容的, 但是它与离散势mKdV方程(17)一起可以实现立方体相容[72]. 在符号(26)式下, 两个方程分别表示为
${\rm B}1(u,\widetilde u,{\widehat{u }},\widehat{\widetilde u};p,q) = p(u {\widehat{\widetilde{u }}} -\widetilde{u} \widehat{u})-q( u\widetilde u{\widehat{u }}{\widehat{\widetilde{u }}} -1) = 0, $

$ {\rm B}2(u,\widetilde u,{\widehat{u }},\widehat{\widetilde u};p,q) = p({\widehat{u }} {\widehat{\widetilde{u }}} -u \widetilde{u})- q( \widetilde u \,{\widehat{\widetilde{u }}}- u{\widehat{u }}) = 0. $
方程(39)在变换$ u\to (-1)^n u $下成为(40)式.
图 5 离散sine-Gordon和势mKdV方程的相容立方体
Figure5. The consistent cube for the discrete sine-Gordon equation and potential mKdV equation

将方程(39)和方程(40)以如下方式放置于立方体的6个面:
$ \begin{split}& {{\text{底}}: \;\;{\rm B}1(u,\widetilde u,{\widehat{u }},\widehat{\widetilde u};p,q) = 0,}\\&{{\text{顶}}:\;\;{\rm B}1({\overline{u}},{{\widetilde{\overline{u}}}}, {{\widehat{\overline{u}}}}, {\widehat{{\widetilde{\overline{u}}}}};p,q) = 0,}\\&{{\text{左}}:\;\; {\rm B}1(u,\widetilde u,{\overline{u}}, {{\widetilde{\overline{u}}}};p,r) = 0,}\\&{{\text{右}}: \;\;{\rm B}1({\widehat{u }},\widehat{\widetilde u},{{\widehat{\overline{u}}}}, {\widehat{{\widetilde{\overline{u}}}}};p,r) = 0,}\\&{{\text{后}}: \;\;{\rm B}2(u,{\widehat{u }},{\overline{u}}, {{\widehat{\overline{u}}}};q,r) = 0,}\\&{{\text{前}}: \;\;{\rm B}2(\widetilde u,\widehat{\widetilde u}, {{\widetilde{\overline{u}}}}, {\widehat{{\widetilde{\overline{u}}}}};q,r) = 0.}\end{split}$
可以检验, 在给定初值$ u, \widetilde u, {\widehat{u }}, {\overline{u}} $以后, $ {\widehat{{\widetilde{\overline{u}}}}} $是唯一的. 这说明两个方程按照如此方式放置于立方体的6个面时, 它们(在一起)具有3维相容性.
ABS方程和离散Boussinesq型方程都是利用同一方程构成相容立方体, 而上面这个例子说明允许用 不同的方程构成相容立方体. 文献[73]给出了更多的例子. Adler等[74]和Boll[75]对这种情况进行了讨论. 详细结果列于Boll[76]的博士学位论文中.
2
2.6.高维的多维相容方程
-->对于定义在6面体上的8点3维(标量)离散方程, 如图6(a), 目前已知的具有四维相容性(见图7)的方程有[34]
图 6 定义3维方程的6面体以及8面体
Figure6. Cube and octahedron for 3D equations

图 7 围绕超立方体的4D相容性
Figure7. 4D consistency around the hyper cube

$ u\widehat{\widetilde{{\overline{u}}}}-\widetilde u{\widehat{\overline{u}}}+ {\widehat{u }}{\widetilde{\overline{u}}} -{\overline{u}} {\widehat{\widetilde{u }}} = 0, \tag{42a}$
$ \frac{(\widetilde u-{\overline{u}})({\widehat{u }}-\widehat{\widetilde{{\overline{u}}}})}{({\overline{u}}-{\widehat{u }})(\widehat{\widetilde{{\overline{u}}}}-\widetilde u)} = \frac{(u-{\widetilde{\overline{u}}})({\widehat{\widetilde{u }}}-{\widehat{\overline{u}}})}{({\widetilde{\overline{u}}}-{\widehat{\widetilde{u }}})({\widehat{\overline{u}}}-u)} . \tag{42b}$
对于定义在8面体上的6点3维离散方程, 如图6(b), 要求这些方程具有4维相容性, Adler等[77]于2011年进行了分类. 满足条件的方程共有5个:
$ {\rm i})\; \; \; \; \; {\widehat{\widetilde{u }}}{\overline{u}} - {\widehat{\overline{u}}} \widetilde u + {\widetilde{\overline{u}}} {\widehat{u }} = 0,\tag{43a}$
$ {\rm ii})\; \; \; \; \frac{({\widehat{\widetilde{u }}}-{\widehat{u }})({\widehat{\overline{u}}}-{\overline{u}})({\widetilde{\overline{u}}}-\widetilde u)}{({\widehat{\widetilde{u }}} - \widetilde u)({\widehat{\overline{u}}}-{\widehat{u }})({\widetilde{\overline{u}}}-{\overline{u}})} = 1, \tag{43b}$
$ {\rm iii})\; \; \; ({\widehat{\widetilde{u }}}-{\widehat{\overline{u}}}){\widehat{u }} + ({\widehat{\overline{u}}}- {\widetilde{\overline{u}}}){\overline{u}} + ({\widetilde{\overline{u}}}-{\widehat{\widetilde{u }}}) \widetilde u = 0, \tag{43c}$
$ {\rm iv}) \; \; \; \frac{({\widehat{\widetilde{u }}}-{\widehat{\overline{u}}})}{{\widehat{u }}} + \frac{({\widehat{\overline{u}}}- {\widetilde{\overline{u}}})}{{\overline{u}}} + \frac{({\widetilde{\overline{u}}}-{\widehat{\widetilde{u }}})}{ \widetilde u} = 0, \tag{43d}$
$ {\rm v}) \; \; \; \; \frac{({\widehat{\overline{u}}}-{\widetilde{\overline{u}}})}{{\overline{u}}} = {\widehat{\widetilde{u }}} \left( \frac{1}{\widetilde u} - \frac{1}{ {\widehat{u }}}\right). \tag{43e}$
这些方程可以存在方向参数, 但是这些参数可以通过适当的变换移除. 这些方程均已经出现于早期的文献中. 如: (42a)式是离散的BKP方程, 又称为Miwa方程[78]; (42b)式被认为是离散Schwarzian BKP, 首先从几何中得到[79]; (43a)式由Hirota[80]在1981年给出, 又称为离散广义Toda方程; (43b)式, (43c)式和(43d)式分别称为离散Schwarzian KP、离散势KP和离散modified KP方程, 最早由Nijhoff等系统构造(见文献[20]中(4)式和(23)式和文献[22]中(4.16)式); (43e)式由Date等于1982年获得(见文献[13]中(N-1)方程).
“多维相容性”不仅作为一类离散系统的可积性的理解, 还提供了研究手段, 可用于构造离散可积系统的B?cklund变换、Lax对、精确解等等.
2
3.1.B?cklund变换
-->B?cklund变换始于对伪球面的构造. Bianchi[81]首先证明了B?cklund给出的含任意参数的变换 可以保持sine-Gordon方程不变. 变换不仅是引入离散变量的一个渠道, 也是联系离散与连续的一座桥梁. 通常, auto-B?cklund变换指同一方程的解之间的变换, nonauto-B?cklund变换指不同方程的解之间的变换.
“多维相容性”提供了非常直观的构造auto-B?cklund变换的方法. 对于任意具有3维相容性的方程(30), 将其置于相容立方体的6个面中, 如图8所示.
图 8 相容立方体
Figure8. The consistent cube

底面和顶面的方程分别为
$ { Q}(u,\widetilde u,{\widehat{u }},\widehat{\widetilde u};p,q) = 0, \tag{44a}$
$ { Q}({\overline{u}},{{\widetilde{\overline{u}}}},{\widehat{\overline{u}}},\widehat{\widetilde{{\overline{u}}}};p,q) = 0.\tag{44b}$
当将$ {\overline{u}} $也视为方程(44a)的解时, 相容立方体侧面的两个方程
$ { Q}(u,\widetilde u,{\overline{u}},{\widetilde{\overline{u}}};p,r) = 0,\tag{45a}$
$ { Q}(u,{\overline{u}}, {\widehat{u }},{\widehat{\overline{u}}};r,q) = 0.\tag{45b}$
自然地构成方程(44a)的一个B?cklund变换. 对于H1方程(33a)来说, 多维相容性为它提供的B?cklund变换是(28)式.
即使是对于两个方程合作构成相容立方体时, 上述思想依然有效. 对于离散sine-Gordon方程(39)和离散势mKdV方程(40)而言, 它们合作构成相容立方体(见(41)式). 此时, 左侧和后侧的两个方程, 即
$ p(u {\widetilde{\overline{u}}} -\widetilde{u} \overline{u})-r( u\widetilde u{\overline{u}}{\widetilde{\overline{u}}} -1) = 0, \tag{46a}$
$ r({\widehat{u }} {\widehat{\overline{u}}} -u \overline{u})- q( {\overline{u}} \,{\widehat{\overline{u}}}- u{\widehat{u }}) = 0, \tag{46b}$
构成离散sine-Gordon方程(39)的B?cklund变换. 不过这个变换没有对称性.
多维相容性也可以用于构造不同方程间的nonauto-B?cklund变换. 下面以Q2和Q1(δ)为例来解释如何利用相容立方体建立nonauto-B?cklund变换. 首先, 将Q2方程(33g)置于相容立方体的6个面. 然后在顶层方程$ {\rm Q}2({\overline{u}}, \widetilde{{\overline{u}}}, \widehat{{\overline{u}}}, \widehat{\widetilde{{\overline{u}}}}; p, q) = $ 0中引入关系
${\overline{u}} = \frac{v}{\epsilon}+\frac{\delta^2}{4\epsilon^2}, $
对方程乘以$ \epsilon^2 $以后取极限$ \epsilon \to 0 $, 顶层方程变为
$ p(v-\widehat{v})(\widetilde{v}-\widehat{\widetilde{v}})-q(v-\widetilde{v})(\widehat{v}-\widehat{\widetilde{v}})+\delta^2 pq(p-q) = 0, $
即Q1(δ). 此时, 立方体的左侧和后侧的两个方程, 经过替换(47)式以及取极限以后, 得到
$ \delta(u-\widetilde u)(v-\widetilde v) = -p[\delta^2(u+\widetilde u)-2 v\widetilde v]+p^2\delta (v+\widetilde v+p\delta)],\tag{48a}$
$ \delta(u-{\widehat{u }})(v-{\widehat{v}}) = -q[\delta^2(u+{\widehat{u }})-2 v{\widehat{v}}]+q^2\delta (v+{\widehat{v}}+q\delta)],\tag{48b}$
这就是Q2与Q1(δ)之间的nonauto-B?cklund变换.
2008年, Atkinson构造了ABS方程之间一系列nonauto-B?cklund变换, 详细结果请参考文献[73]中表3. 除了上述方法, B?cklund变换可以利用Yang-Baxter映射以及ABS方程的分解性质来构造, 相关内容读者可以分别参考文献[82]和文献[83].
2
3.2.Lax对
-->众所周知, B?cklund变换与Lax对存在密切联系, 对于多维相容系统来说, 这种联系更加直接. 下面以H1方程(33a)为例, 利用其3维相容性引出的B?cklund变换构造它的Lax对.
将H1的B?cklund变换(28)式改写为
$ {\widetilde{\overline{u}}} = \frac{p-r}{{\overline{u}}-\widetilde u}+u,\; \; {\widehat{\overline{u}}} = \frac{r-q}{{\widehat{u }}-{\overline{u}}}+u. $
$ {\overline{u}} = g/f $, 得到
$\begin{split} & \frac{\widetilde{g}}{{\widetilde{f}}} = \frac{ug+(p-r-u\widetilde u)f}{g-\widetilde u f},\\& \frac{\widehat{g}}{{\widehat{f}}} = \frac{ug+(q-r-u {\widehat{u }}) f}{g-{\widehat{u }} f}. \end{split}$
引入$ { {\varPhi}} = (g, f)^{\rm T} $, 上式改写为
$ \widetilde{{ {\varPhi}}} = \gamma_1 { M}{ {\varPhi}},\; \; \; \widehat{{\rm {\varPhi}}} = \gamma_2 { {N\varPhi}}, \tag{51a}$
$ { M} = \left(\!\!\!\begin{array}{cc} u & p-r-u\widetilde u\\ 1 & -\widetilde u \end{array}\!\!\!\right),\; \; \; { N} = \left(\!\!\!\begin{array}{cc} u & q-r-u {\widehat{u }}\\ 1 & -{\widehat{u }} \end{array}\!\!\!\right),\tag{51b} $
式中$ \gamma_i $可视为 分式线性形式(50)式写成矩阵形式(51)式后保留的原分子、分母的公因子, 一般取为$ \gamma_1 = \dfrac{1}{\sqrt{|{ M}|}}, \; \gamma_2 = \dfrac{1}{\sqrt{|{ N}|}} $, 对于H1来说, $ \gamma_i $可以取为任意常数. 相容性$ {\widehat{\widetilde{{ {\varPhi}} }}} = \widetilde{\widehat{{ {\varPhi}} }} $引出
$ {\widehat{ M}} { N} = {\widetilde{ N}} { M}, $
计算后可得H1方程(33a).
在H1的Lax对(51)式中, 取$ {\varPhi} = (\phi_1, \phi_2)^{\rm T} $, $ \gamma_1 = \gamma_2 = 1 $, 从$ \widetilde { {\varPhi}} = { M}{ {\varPhi}} $中消去$ \phi_1 $, 得到
$ \widetilde{\widetilde{\phi}}_2 = (u-\widetilde{\widetilde{u}}) \widetilde{\phi}_{2}-(p-r)\phi_{2}, $
此为H1方程谱问题的标量形式. 由此出发 可以构造H1方程的无穷守恒律[84].
上面以H1为例给出了利用B?cklund变换构造Lax对的方法, 事实上, 所有的ABS方程都具有相同结构的Lax对, 且存在一般形式. 具体地, 对于任一ABS方程$ { Q}(u, \widetilde u, {\widehat{u }}, {\widehat{\widetilde{u }}}, p, q) = 0 $, 其Lax对为(参考文献[85])
$ \widetilde{{ {\varPhi}}} = \gamma_1 { M}{ {\varPhi}},\; \; \; \widehat{{ {\varPhi}}} = \gamma_2 { N}{ {\varPhi}}, \tag{54a}$
其中
$\begin{split} & { M} = \left(\begin{array}{cc} -\partial_{{\widehat{u }}}Q& -Q\\ \partial_{{\widehat{\widetilde{u }}}}\partial_{{\widehat{u }}}Q& \partial_{{\widehat{\widetilde{u }}}}Q \end{array} \right)_{{\widehat{u }} = {\widehat{\widetilde{u }}} = 0,\; q = r}, \\&{ N} = \left(\begin{array}{cc} -\partial_{\widetilde u}Q& -Q\\ \partial_{{\widehat{\widetilde{u }}}}\partial_{\widetilde u}Q& \partial_{{\widehat{\widetilde{u }}}}Q \end{array} \right)_{\widetilde u = {\widehat{\widetilde{u }}} = 0,\; p = r},\end{split}\tag{54b}$
$ \gamma_1 = \sqrt{|{ M}|}, \; \gamma_2 = \sqrt{|{ N}|} $. 显然, 将矩阵谱问题改写成标量形式后, 可以发现所有的ABS方程都具有类似于(53)式的二阶离散谱问题.
通常, 对于由一个方程构成相容立方体时, 由于方程的对称性, 其B?cklund变换具有对称性, Lax对也具有对称性, 即在上式M中, 将$ (^{\widetilde{\; \; }}, p) $换为$ (^{\widehat{\; \; }}, q) $, 即得到N. 在文献[86]中搜集了更多的利用3维相容性构造Lax对的例子. 对于由两个方程构成相容立方体时, 例如2.5节中的离散sine-Gordon方程和势mKdV方程, 由它们的相容立方体构造的离散sine-Gordon方程的B?cklund变换具有非对称性, 由此引出的离散sine-Gordon方程的Lax对也是非对称的:
$ \widetilde { {\varPhi}} = \frac{-1}{\sqrt{u\widetilde u}}\left(\begin{array}{lr}-p \widetilde u & r\\ ru\widetilde u & -pu \end{array}\right){ {\varPhi}}, \tag{55a}$
$ \widehat{{ {\varPhi}} } =\frac{1}{\sqrt{u{\widehat{u }}}}\left(\begin{array}{lr} -r u & q u {\widehat{u }}\\ q & -r {\widehat{u }} \end{array}\right){ {\varPhi}}. \tag{55b}$

2
3.3.孤子解
-->下面以离散势KdV (即H1)方程(33a)为例, 演示如何通过其B?cklund变换构造单孤子解(1-soliton solution (1SS)). 首先需要一个种子解. 取
$ a^2 = -p,\; \; \; b^2 = -q, $
不难发现
$ u_0 = a n +b m +\lambda, $
是方程(33a)的一个解, 其中$ \gamma $是任意常数.
方程(33a)的B?cklund变换是(28)式, 即
$ (u-\widetilde{\overline{u}})(\widetilde{u}-\overline{u}) = -a^2+k^2,\tag{58a}$
$ (u-{\widehat{\overline{u}}})(\overline{u}-\widehat{u}) = -k^2+b^2. \tag{58b}$
相应于参数化(56)式, 这里取$ r = -k^2 $. 取u为(57)作为种子解, 设新解$ \overline{u} $形如
$ \overline{u} = \overline{u}_{0}+v, $
其中$ \overline{u}_{} $是种子解(57)沿l方向的一次平移:
$ \overline{u}_{0} = a n +b m +k+\lambda, $
v是待定函数. 将(59)式和(57)式代入B?cklund变换(58)式得到
$ \widetilde{v} = \frac{Ev}{v+F},\; \; \; \widehat{v} = \frac{Gv}{v+H}, $
其中
$\begin{split} & E = -(a+k),\; \; F = -(a-k),\\ & G = -(b +k),\; \; H = -(b-k).\end{split}$
引入$ v = f/g $$ { {\varPhi}} = (f, g)^{\rm T} $, 将(61)改写为
$ \widetilde{ {\varPhi}} = \left(\begin{array}{cc} E & 0\\ 1 & F \end{array} \right){ {\varPhi}},\; \; \; \widehat{{ {\varPhi}}} = \left(\begin{array}{cc} G & 0\\ 1 & H \end{array} \right){ {\varPhi}}. $
验算可见成立相容性$ \widehat{\widetilde{{ {\varPhi}}}} = \widetilde{\widehat{{ {\varPhi}}}} $, 否则需要在上式的矩阵前乘上平衡因子(类似于(51)式中的$ \gamma_i $). 进一步求解, 先后得到
$ f_{n,m} = E^n f_{0,m},\; \; \; g_{n,m} = \frac{E^n-F^n}{-2k}f_{0,m}+F^n g_{0,m}, $
以及
$ f_{n,m} = G^m f_{n,0},\; \; \; g_{n,m} = \frac{G^m-H^m}{-2k}f_{n,0}+ H^m g_{n,0}. $
从而有
$ {\varPhi}_{n,m} = \left(\begin{array}{cc} E^nG^m & 0\\ \dfrac{E^nG^m-F^nH^m}{-2k} & F^nH^m \end{array} \right){ {\varPhi}}_{0,0}. $

$ \rho_{n,m}\! =\! \biggl(\frac{E}{F})^n \! \biggl(\frac{G}{H})^m \rho_{0,0} \!=\! \biggl(\frac{a +k}{a -k})^n \!\biggl(\frac{b + k}{b-k})^m \rho_{0,0}, $
其中$ \rho_{0, 0} $为常数, 从(63)式可得
$ v_{n,m} = \frac{v_{0,0}\rho_{n,m}/\rho_{0,0} } {1+\dfrac{v_{0,0}}{2k} - \dfrac{v_{0,0}}{2k}\rho_{n,m}/\rho_{0,0}}. $
重新定义常数$ \rho_{0, 0} $后, 有
$ v_{n,m} = \frac{-2k \rho_{n,m}}{1+\rho_{n,m}}. $
代入到(59)式, 离散势KdV方程(33a)的1-孤子解为
$ u_{n,m}^{1SS} = \overline{u}_{0}+v_{n,m} = an + bm + \lambda +\frac{k(1- \rho_{n,m})}{1+\rho_{n,m}}. $
依据上述过程, 若利用B?cklund 变换获得2-孤子解显然不容易. 但是1-孤子解的结构往往可以“暗示”一些2-孤子解的结构信息, 有助于发现合适的变换公式将离散方程双线性化, 并进一步获得N-孤子解. 相关内容可以参考文献[38, 87, 88].
2
3.4.0-孤子解与不动点方法
-->对于(33)式中所列的ABS方程, $ u = 0 $一般都不是解, 而且也不易明显看出一些简单的解. “不动点方法”是在B?cklund 变换的基础上求解0-孤子解的方法. 从上一节的求解过程来看, 对于方程$ { Q}(u, \widetilde u, {\widehat{u }}, {\widehat{\widetilde{u }}};p, q) = 0 $的B?cklund变换(3.1), 其中的参数r是孤子参数, 将作为离散谱在$ {\overline{u}} $中引出一个孤子. 现在, 如果孤子参数r在变换中“失效”, 即$ {\overline{u}} = u $, 则有
$ { Q}(u,\widetilde u, u,\widetilde u;p,r) = 0,\tag{68a}$
$ { Q}(u, u, {\widehat{u }},{\widehat{u }};r,q) = 0. \tag{68b}$
由此变换引出的原方程的解应当是最简单的, 可以作为0-孤子解.
对于H1方程(33a), 有
$ (\widetilde u-u)^2 = r-p,\quad ({\widehat{u }}-u)^2 = r-q. $
由此很容易得到
$ u = an+bm+\gamma,\; \; a^2 = r-p, \; \; b^2 = r-q, $
即(57)式.
如果方程$ { Q}(u, \widetilde u, {\widehat{u }}, {\widehat{\widetilde{u }}};p, q) = 0 $在变换$u\to T(u) $下不变, 可在变换(3.1)中取$ {\overline{u}} = T(u) $, 从而有
$ { Q}(u,\widetilde u, T(u), T(\widetilde u);p,r) = 0, \tag{69a}$
$ { Q}(u, T(u), {\widehat{u }},T({\widehat{u }});r,q) = 0. \tag{69b}$
对于Q1(δ)方程(33f), 在$ T(u) = u+c $下不变, 对于的B?cklund变换(69)式有
$\begin{split} & r(u-\widetilde u)^2 = p(c^2+\delta^2r(p-r)),\\ & r(u-{\widehat{u }})^2 = q(c^2+\delta^2r(q-r)). \end{split}$
经过参数化$ (p, q)\to(\alpha, \beta) $:
$\begin{split} & p = \frac{c^2/r-\delta^2 r}{a^2-\delta^2},\quad q = \frac{c^2/r-\delta^2 r}{b^2-\delta^2},\\ & \alpha: = pa,\quad \beta: = qb, \end{split}$
可得Q1(δ)的线性0-孤子解:
$ u_0 = \alpha n+\beta m+\lambda, $
其中$ \lambda $是任意常数. 此外, Q1(δ)方程(33f), 还在$ T(u) = -u $下不变. 此时的B?cklund变换(69)式给出
$\begin{split} & 4pu\widetilde u +r(u-\widetilde u)^2 = \delta^2pr(p-r),\\ & 4q u{\widehat{u }}+ r(u-{\widehat{u }})^2 = \delta^2r(q-r). \end{split}$
经过参数化$ (p, q)\to(\alpha, \beta) $:
$ p = \dfrac{1}{2} r(1- \cosh(\alpha')) = -\dfrac{1}{4}r(1-\alpha)^2/\alpha, $
$ q = \dfrac{1}{2} r(1- \cosh(\beta')) = -\dfrac{1}{4}r(1-\beta)^2/\beta, $
得到解
$ u_0 = A \alpha^n \beta^m +B \alpha^{-n}\beta^{-m},\quad AB = \delta^2 r^2/16. $
上述所描述的过程称为“不动点方法”. 首先用于构造Q4方程的种子解[89], 其后又系统地应用于其他ABS方程0-孤子解的构造[37,38]. 当然, 也可以利用其他的方法构造0-孤子解, 例如, 利用方程间的B?cklund变换(如(48)式), 从一个方程的0-孤子解得到另一个方程的0-孤子解.
通过简单的描述和例子介绍了多维相容性的概念及其在离散可积系统中的应用. 对于多维相容性, 存在一定的几何背景, 换言之, 经典的初等几何中的点线之间的关系蕴含着若干离散可积系统的多维相容性(如文献[79,90]). 此外, 多维相容性也可以从平面波因子的对称性来理解. Miwa映射提供的离散的平面波因子[78]
$ \psi =\prod\limits_{i}{{{(1-{{p}_{i}}k)}^{{{n}_{i}}}}} $
在各个方向上具有对称性; 考虑到不少离散可积系统可以从平面波因子出发通过Cauchy矩阵方法[37,45] 等途径来构造, 不难理解由此获得的离散可积系统具有多方向上的相容性.
对比连续系统, 目前, 离散可积系统在代数结构和工具、几何背景、离散的分析工具等方面都仍待发展. 例如, 基于谱问题、零曲率方程和Kac-Moody代数的连续可积系统 理论并不适用于离散系统, 许多连续意义下的几何体系仍未实现离散化, 离散的复分析也尚未成熟.
对于离散可积系统的研究也是一个在不断认识离散可积性、发展研究方法和工具、与新兴数学工具结合(如Tropical几何、Cluster代数等)的过程. 关于离散可积系统中值得关注的发展与研究方向, 建议读者关注SIDE会议的网站.
相关话题/系统 文献 方程 过程 数学

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • Boussinesq方程的Lax对、B?cklund变换、对称群变换和Riccati展开相容性
    摘要:Boussinesq方程是流体力学等领域一个非常重要的方程.本文推导了Boussinesq方程的Lax对.借助于截断Painlevé展开,得到了Boussinesq方程的自B?cklund变换,以及Boussinesq方程和Schwarzian形式的Boussinesq方程之间的B?cklun ...
    本站小编 Free考研考试 2021-12-29
  • 高阶Ablowitz-Ladik方程的局域波解及稳定性分析
    摘要:本文构造了一类高阶Ablowitz-Ladik方程的广义$(M,N-M)$-波Darboux变换,借助符号计算从不同背景出发研究了该模型丰富的局域波解,并利用数值模拟研究了这些解的动力学稳定性.关键词:高阶Ablowitz-Ladik方程/M,N-M)-波Darboux变换');"href=' ...
    本站小编 Free考研考试 2021-12-29
  • 一个可积的逆空时非局部Sasa-Satsuma方程
    摘要:本文给出了一个可积的逆空时(逆空间-逆时间)非局部Sasa-Satsuma方程.建立了这个方程的Darboux变换,并且构造了这个逆空时非局部方程在零背景条件下的孤子解.关键词:可积的逆空时非局部Sasa-Satsuma方程/达布变换/孤子解EnglishAbstractAnintegrabl ...
    本站小编 Free考研考试 2021-12-29
  • 四阶色散非线性薛定谔方程的明暗孤立波和怪波的形成机制
    摘要:本文研究了四阶色散非线性薛定谔方程的明暗孤立波和怪波的形成机制,该模型既可以模拟高速光纤传输系统中超短脉冲的非线性传输和相互作用,又可以描述具有八极与偶极相互作用的一维海森堡铁磁链的非线性自旋激发现象.本文首先通过对四阶色散非线性薛定谔方程的相平面分析,发现由其约化得到的二维平面自治系统具有同 ...
    本站小编 Free考研考试 2021-12-29
  • 可积谐振系统中的极端波事件研究进展
    摘要:从微观角度上讲,单个极端异常波事件可视为可积模型方程的时空局域有理函数解.本文主要讨论了三类典型的可积谐振相互作用模型(即长波短波谐振方程,三波谐振相互作用方程,非线性薛定谔和麦克斯韦-布洛赫方程)的基阶Peregrine异常波解及其相关研究进展;明确指出了这些基阶异常波解形式具有普适性,可推 ...
    本站小编 Free考研考试 2021-12-29
  • 可积系统多孤子解的全反演对称表达式
    摘要:多孤子解是非线性数学物理系统的基本激发模式.文献中存在各种类型的表达式,如广田(Hirota)形式,朗斯基(Wronskian)或双朗斯基形式和法夫(Phaffian)形式.最近在多地系统的研究中,我们发现使用一种全新但等价的形式具有极为简洁和方便的优点.本文主要综述多种类型可积非线性系统的多 ...
    本站小编 Free考研考试 2021-12-29
  • 准Λ型四能级系统选择反射光谱
    摘要:运用密度矩阵理论研究了气固界面准$\Lambda$型四能级原子系统的非线性选择反射光谱.基于刘维尔方程给出了一阶近似条件下探测光场的解析式.在探测场为弱场时,分析了信号场拉比频率、失谐量和耦合场失谐量对反射光谱线型的影响.数值模拟表明:信号场参与产生的选择反射峰线宽可以利用信号场拉比频率进行调 ...
    本站小编 Free考研考试 2021-12-29
  • 隔板对流系统的热流特性及热量输入与传递特性
    摘要:采用DNS方法对隔板对流装置进行模拟计算,研究系统中热流特性以及热量输入与传递特性.讨论了热流的纵向和横向输运特性,在此基础上对传热通道和狭缝区域的热通量以及对应底板外界输入热通量进行了定量化分析.研究结果表明,通道中低温流体向下冲击底板而后转入水平运动,流入狭缝区域并不断被加热,只进行水平的 ...
    本站小编 Free考研考试 2021-12-29
  • 钨空位捕获氢及其解离过程的分子动力学
    摘要:氢(H)同位素滞留问题是聚变堆第一壁材料设计的关键,而深入理解H在缺陷(如空位)处的非均匀形核长大过程有助于揭示H起泡及滞留的机制.针对第一壁材料钨(W)中空位捕获H及其解离的动力学过程开展研究,通过耦合捕获和解离两过程,构建新物理模型,避免了原单一过程的物理模型需准确记录相应事件首次发生时间 ...
    本站小编 Free考研考试 2021-12-29
  • Shimizu-Morioka系统与Finance系统生成Lorenz混沌的微分几何策略
    摘要:从一种受控混沌系统生成另一混沌系统可增强保密通信的安全性,具备潜在应用前景.研究了如何通过状态变换以及单输入反馈,驱使受控Shimizu-Morioka系统与受控Finance系统生成Lorenz混沌动态.主要方法是运用微分几何理论,将上述三种系统等价转换为下三角形式,并尽量简化和一致化其方程 ...
    本站小编 Free考研考试 2021-12-29