下载:
-
1955年, Fermi, Pasta和Ulam领导的科学家小组用数值方法计算了用非线性弹簧联结的64个质点组成的谐振子的振动, 其目的是从数值实验上验证统计力学中的能量均分定理, 这在后来被人们称为著名的FPU实验[1]. 1967年, Toda[2]考虑晶体的非线性振动, 提出了著名的Toda晶格方程近似模拟这种情况, 并得到了该模型的孤子解, 从而使FPU试验问题得到合理的解释和正确解答[1]. Toda晶格方程作为一类可积的半离散的非线性微分差分方程可以描述一些物理学中的非线性波的传播现象, 该方程的提出引起了人们对可积的离散孤子方程研究的热潮. 可积系统中的非线性微分差分方程(又称离散的孤子方程)是一类重要的半离散的非线性偏微分方程, 近年来国际上对这类模型的研究有着极大的兴趣, 这些方程与元胞自动机、DNA的研究、辛算法有密切的关系, 在电学、光学、磁性流体、超导、生物和等离子体中有着广泛的应用, 有十分广阔的应用前景, 目前离散问题的研究是国际上研究热点之一[3]. 日本科学家Hirota[4]指出连续的孤子方程可以被离散化, 并且不远的将来将会是差分方程的时代, 这里所说的差分方程就是离散的孤子方程.
近年来, 孤子方程的局域波解受到了数学家和物理学家的广泛关注, 按传播特性, 局域波主要包括孤子、呼吸子、怪波等[5-18], 广泛存在于非线性光学、玻色斯坦凝聚、等离子物理等各种非线性物理系统[5-16]. 孤子又称孤立波, 是一种在传播过程中形状、幅度和速度都维持不变的脉冲状行波, 并且孤子与其他同类孤子相遇后, 能维持其幅度、形状和速度不变[5,6]. 呼吸子是一种特殊的局域周期振荡孤子, 按传播方向有空间周期的Akhmediev呼吸子和时间周期的Kuznetsov-Ma呼吸子[5,6], 从光学和流体力学到玻色-爱因斯坦凝聚体和等离子体等各种物理情况下都能广泛地观察到呼吸子的传播现象, 关于呼吸子碰撞特性的研究最近取得一系列有趣的结果, 包括棋盘干涉时空结构[7], 具有不同传播特征的super-regular呼吸子作用现象[8-11]等. 怪波最初是用来描述海洋中出现的大振幅波, 突然出现然后又很快消失得无影无踪, 怪波也被称为畸形波, 是一种短时间内存在的大振幅的局域波动[5,6,14-18]. 怪波发生物理机理的研究是怪波研究的重点, 调制不稳定性通常被认为是呼吸子和怪波发生的主要机制, 并且一些文献认为它可能是形成怪波和呼吸子的必要条件[10,11]. 目前不同类型的局域波作用解是国际上研究的热点问题, 然而大部分问题和方法仅限于连续的孤子方程的局域波及其作用[5-24]. 例如文献[5,6]研究了高维非线性模型的局域波作用解包括呼吸子和怪波以及Lump解之间的相互作用解; 文献[14]研究了非线性光纤中呼吸子和怪波的相互作用现象. 这些局域波作用对于理解新的物理现象具有重要的理论意义. 然而对于离散孤子方程的局域波的研究, 由于研究的困难性, 据我们所知, 研究还不充分, 不系统[25,26], 因此本文将研究离散的非线性微分差分方程丰富的局域波解. 作为例子, 研究下面的非线性微分差分方程
$\begin{split} r_{n,t} =\, & {\rm i}(1-\sigma |r_{n}|^2)[\sigma r_{n}^*(r_{n-1}^2+r_{n+1}^2)\\ &+\sigma r_{n}(r_{n-1}r_{n+1}^*+r_{n-1}^*r_{n+1})\\ &-r_{n-2}(1-\sigma |r_{n-1}|^2)\\ &-r_{n+2}(1-\sigma |r_{n+1}|^2)+2(r_{n-1}+r_{n+1})],\end{split}$ 其中
$ r_{n} = r(n, t), r_{n, t} = {{\rm d}r_n}/{{\rm d}t} $ , 星号$ * $ 表示共轭, i是虚数单位,$ \sigma = \pm1 $ , 这里加号和减号分别表示聚焦和散焦情况. 方程(1)可以通过AKNS方法得到[3], 这里为节省篇幅不给出具体推导过程. 由于方程(1)和Ablowitz-Ladik方程[3,27] (即${\rm i}r_{n, t} = $ $ r_{n-1}+r_{n+1}-2 r_{n}-\sigma |r_{n}|^2(r_{n-1}+r_{n+1}) $ )属于同一个梯队, 是这个梯队的第二个方程, 我们称方程(1)为高阶的Ablowitz-Ladik方程. 正如五阶KdV方程可以像KdV方程能描述浅水波的运动一样[3], 我们有理由相信方程(1)也可能像它对应的低阶Ablowitz-Ladik方程一样可以描述光纤中光孤子的传播, 因此研究方程(1)具有重要的理论和物理意义. 虽然对离散的Ablowitz-Ladik方程的局域波已经有了一定的研究结果, 特别是应用双线性方法得到了它的怪波[28-30], 然而对于离散的孤子方程(1)还没有系统的方法进行研究, 特别是仍然没有通过我们提出的广义$ (M, N-M) $ 波Darboux变换进行系统的研究. 因此本文将应用广义(M, N–M)波Darboux变换从研究方程(1)的精确解出发, 寻找其新奇的局域波结构, 特别是不同类型的局域波相互作用的新奇局域波结构.本文的主要结构安排如下: 第二节构造出方程(1)的Lax对和广义
$ (M, N-M) $ 波Darboux变换; 第三节将应用广义$ (M, N-M) $ -波Darboux变换给出方程(1)不同类型的局域波解, 并通过数值模拟研究其传播稳定性; 最后是本文的结论. -
$ E{{{\varphi}}_{{{n}}}} = {{{U}}_{{{n}}}}{{{\varphi}} _{{{n}}}},\ \ \ \ {{{\varphi}} _{{{n}},{{t}}}} = {{{V}}_{{{n}}}}{{{\varphi}} _{{{n}}}}, $ 这里
$ {{{U}}_{{n}}} = \left(\!\!\!{\begin{array}{*{20}{c}} {{\lambda ^2}}&{\lambda {r_n}}\\ {\sigma \lambda r_n^*}&1 \end{array}}\!\!\! \right),\;\;\;{{{V}}_{{n}}} = \left(\!\!\!{\begin{array}{*{20}{c}} {{V_{11}}}&{{V_{12}}}\\ {{V_{21}}}&{{V_{22}}} \end{array}}\!\!\!\right), $ 其中
$ \begin{aligned} V_{11} =\; & -\frac{\rm i}{2}\lambda^4+{\rm i}(1+\sigma r_{n}r_{n-1}^*)\lambda^2\\ &+\sigma {\rm i} r_{n+1}r_{n-1}^*(1-\sigma|r_{n}|^2)\\ &+{\rm i}\sigma r_{n} r_{n-2}^*(1-\sigma|r_{n-1}|^2)-{\rm i} r_{n}^2 r_{n-1}^{*2}\\ &-\frac{\rm i}{2\lambda^4}+\frac{{\rm i}(1+\sigma r_{n-1}r_{n}^*)}{\lambda^2}, \\ V_{12} =\, & {\rm i}(\sigma|r_{n}|^2 r_{n+1}^*+\sigma r_{n}^2 r_{n-1}^*-r_{n+1}+2 r_{n})\lambda \\ &-{\rm i} r_{n}\lambda^3+\frac{{\rm i} r_{n-1}}{\lambda^3}\\ &+\frac{{\rm i}(r_{n-2}\!-\!\sigma r_{n-1}^2 r_{n}^*\!-\!\sigma r_{n-1} r_{n-1}^* r_{n-2}\!-\!2 r_{n-1})}{\lambda}, \\ V_{21} =\, & {\rm i}(-\sigma r_{n-2}^*+r_{n-1}^{*2}r_{n}+|r_{n-1}|^2 r_{n-2}^*\\ &+2\sigma r_{n-1}^*)\lambda-\sigma {\rm i} r_{n-1}^* \lambda^3+\frac{{\rm i}\sigma r_{n}^*}{\lambda^3}\\ &+\frac{{\rm i}(\sigma r_{n+1}^*-r_{n}^{*2}r_{n-1}-|r_{n}|^2 r_{n+1}^*-2\sigma r_{n}^*)}{\lambda},\\ V_{22} =\, &\frac{\rm i}{2}\lambda^{4}-{\rm i}(1+\sigma r_{n}r_{n-1}^*)\lambda^2\\ &-{\rm i}\sigma r_{n-1}r_{n+1}^*(1-\sigma|r_{n}|^2)\\ &-{\rm i}\sigma r_{n-2}r_{n}^*(1-\sigma|r_{n-1}|^2)+{\rm i}r_{n}^{*2}r_{n-1}^2\\ &+\frac{\rm i}{2\lambda^4}-\frac{{\rm i}(1+\sigma r_{n-1}r_{n}^*)}{\lambda^2},\end{aligned} $ 这里
$ \varphi_{n} = (\phi_{n}, \psi_{n})^{\rm T} $ 表示Lax对(2)式的解, T表示转置. 借助符号计算Maple, 容易验证方程(2)的相容性条件$ U_{n, t} = V_{n+1}U_{n}-U_{n}V_{n} $ 等价于方程(1). 接下来基于Lax对(2)式, 构造方程(1)的广义$ (M, N-M) $ - 波Darboux变换, 然后通过它构造方程(1)的局域波解. 为此考虑下面的规范变换$ {\widetilde{\varphi}_{n}} = {T_n}{\varphi_{n}}, $ 其中是
$ { T}_{n} $ 是$ 2\times 2 $ 的矩阵, 根据Darboux变换的知识,$ \widetilde{ {\varphi}}_{n} $ 必须满足$ {\widetilde{{{\varphi}}}_{{{n}}+1}} = {{\widetilde{U}_{n}}}{{\widetilde{\varphi}_{n}}},\quad {\widetilde{{{\varphi}}}_{{{n}},{{t}}}} = {\widetilde{{{V}}}_{{{n}}}}{\widetilde{{{\varphi}}}_{{{n}}}}, $ 这里
$ {\widetilde{U}_{n}}, {\widetilde{V}_{n}} $ 和$ {{U}_{n}}, {{V}_{n}} $ 有相同的形式, 是将$ {{U}_{n}}, $ ${{V}_{n}} $ 的$ {r}_{n} $ 替换为$ \widetilde{r}_{n} $ 得到的. 由(3)式和(4)式, 可知$ {T_n} $ ,$ {\widetilde{U}_{n}} $ ,$ {\widetilde{V}_{n}} $ ,$ {U_{n}} $ ,$ {V_{n}} $ 必须满足$\begin{split} & {\widetilde{U}_{n}} = {{{T}}_{{{n}}+1}}{U_{n}}{{T_n}}^{-1},\\ & {\widetilde{V}_{n}} = ({T_{n,t}}+{T_{n}}{V_{n}}){{T_n}}^{-1}, \end{split}$ 为此, 我们构造一个特殊的Darboux阵如下
$ {{{T}}_{{n}}} \!=\! \left(\!\!\!\!\!\begin{array}{*{20}{c}} \lambda^{2N} \!+\! \displaystyle\sum\limits_{j = 1}^{N}a^{(2j-2)}_n \lambda^{2j-2} & \displaystyle\sum\limits_{j = 1}^{N}b^{(2j-1)}_n \lambda^{2j-1}\\ -\sigma\! \displaystyle\sum\limits_{j = 1}^{N} \! b^{(2N-2j+1)*}_n \lambda^{2j-1} &1\!+\! \displaystyle \sum\limits_{j = 1}^{N}\!a^{(2N-2j)*}_n \lambda^{2j} \end{array}\!\!\!\!\!\right), $ 这里N是正整数,
$ a^{(k)}_{n}, b^{(k)}_{n} $ 是$ n, t $ 的未知函数. 当选择N个合适的$ \lambda_{k} (k = 1, 2, ..., N) $ ,$ a^{(k)}_{n}, b^{(k)}_{n} $ 可以通过解下面的方程组得到:$ \left\{ \begin{aligned} &{{{T}}^{{\bf{(0)}}}_{{{n}}}}(\lambda_{i}){{{\varphi}}^{{\bf{(0)}}}_{{{n}}}}(\lambda_{i}) = 0, \\ &{{{T}}^{{\bf{(0)}}}_{{{n}}}}(\lambda_{i}){{{\varphi}}^{{\bf{(1)}}}_{{{n}}}}(\lambda_{i})+{{{T}}^{{\bf{(1)}}}_{{{n}}}}(\lambda_{i}){{{\varphi}}^{{\bf{(0)}}}_{{{n}}}}(\lambda_{i}) = 0, \\ & \qquad \cdots, \\ &\sum\limits_{j = 0}^{m_{i}}{{{T}}^{({{j}})}_{{{n}}}}(\lambda_{i}){{{\varphi}}^{({{m}}_{{{i}}}-{{j}})}_{{{n}}}}(\lambda_{i}) = 0, \end{aligned} \right. $ 这里
$ N = M+ \displaystyle\sum\limits_{i = 1}^{M}m_{i} $ , 其中$ 1\leqslant M\leqslant N $ 是正整数.通过上面的分析并根据文献[7-10]中的步骤, 若
$ {{{\varphi}}({{\lambda}}_{{{k}}})} = (\phi(\lambda_{k}), \psi(\lambda_{k}))^{\rm T}, (k = 1, 2, ..., M) $ 是Lax对的M个解, 且$ r_{n} $ 是方程(1)的种子解, 则可以得到关于方程(1)如下的广义$ (M, N-M) $ -波Darboux变换定理:定理1 当方程(1)的新解
$ \widetilde{r}_{n} $ 和旧解$ r_{n} $ 的变换如下$ \widetilde{r}_n = b^{(1)}_{n+1}+r_n a^{(0)}_{n+1}, $ 则由(4)式确定的矩阵
$ { {\widetilde{U}}}_{n}, { {\widetilde{V}}}_{n} $ 和$ { U}_{n}, { V}_{n} $ 具有相同的形式. 这里在(8)式中$ a^{(0)}_n = {\Delta a^{(0)}_n}/{\Delta_N^{\epsilon}} $ ,$ b^{(1)}_n = {\Delta b^{(1)}_n}/{\Delta_N^{\epsilon}} ,$ $ \begin{split} & \varDelta_N^{\epsilon} = \det([\varDelta^{(1)}, \varDelta^{(2)}, ..., \varDelta^{(n)}]^{ T}) ,\\ &\varDelta^{(i)} = (\varDelta^{(i)}_{j, s})_{2(m_{i}+1)\times2 N},\end{split}$ 其中
$ \varDelta^{(i)}_{j,s} = \begin{cases} \displaystyle\sum\nolimits_{k = 0}^{j-1}C^{k}_{2N-2s} {\lambda_i}^{(2N-2s-k)} {\phi^{(j-1-k)}_{i}}, & 1\leqslant j \leqslant m_{i}+1, ~1 \leqslant s\leqslant N, \\ \displaystyle\sum\nolimits_{k = 0}^{j-1}C^{k}_{4N-2s+1} {\lambda_i}^{(4N-2s+1-k)} {\psi^{(j-1-k)}_{i}}, & 1\leqslant j\leqslant m_{i}+1, ~N+1\leqslant s\leqslant 2N, \\ \displaystyle\sum\nolimits_{k = 0}^{j-(N+1)}C^{k}_{2N-2s} {\lambda_i}^{2N-2s-k^*} {\psi_{i}^{(j-1-k)*}}, & m_{i}+2\leqslant j\leqslant 2(m_{i}+1), ~1\leqslant s\leqslant N, \\ -\sigma \displaystyle \sum\nolimits_{k = 0}^{j-(N+1)}C^{k}_{4N-2s+1} {\lambda_i}^{2s-1+k^*} {\phi_i^{(j-N-1-k)*}}, & m_{i}+2\leqslant j\leqslant 2(m_{i}+1), ~N+1\leqslant s\leqslant 2N, \end{cases}$ $ \Delta a^{(0)}_n $ 和$ \Delta b^{(1)}_n $ 由向量$ (g^{(1)}, g^{(2)}, \cdots, g^{(n)}) $ 替换$ \varDelta_N^{\epsilon} $ 中第N列和$ 2 N $ 列得到, 将$ \Delta a^{(0)}_n $ 和$ \Delta b^{(1)}_n $ 中的n替换为$ n+1 $ 可以得到$ \Delta a^{(0)}_{n+1} $ 和$ \Delta b^{(1)}_{n+1} $ , 这里$ g^{(i)} = (g^{(i)}_{j})_{2(m_i+1)\times1} $ , 这里当$ 1\leqslant j\leqslant (m_{i}+1) $ 时,$ g^{(i)}_j = -\displaystyle\sum\nolimits_{k = 0}^{j-1}C^{k}_{2 N} {\lambda_i}^{2 N-k} {\phi_{i}^{(j-1-k)}} $ , 当$ m_{i}+2\leqslant $ $j\leqslant 2(m_{i}+1) $ 时,$ g^{(i)}_j = - {\psi_{i}^{(j-N-1)*}} $ .这里需要说明的是: 对于定理1, (3)式和(8) 式(即
$ (r_{n}, \varphi_{n})\to (\widetilde r_{n}, \widetilde\varphi_{n}) $ )称为方程(1)的广义$ (M, N-M) $ -波Darboux变换, 这里M表示使用谱参数的个数,$ N-M $ 表示Darboux阵中使用泰勒展开式最高阶导数的和. 当M = N, mi =０$ (i = 1, 2, \cdots, M) $ , 通常的N-波Darboux变换就是定理1的特殊情况, 可以从零种子解(平面波种子解)求得方程(1)的多孤子解(多呼吸子解); 当$ M = 1, m_{1} = N-1 $ , 定理1就约化为$ (1, N-1) $ -波Darboux变换, 可以从平面波种子解求得方程(1)的怪波解; 当$ M = 2, m_{2} = N-2 $ , 定理1就约化为$ (2, N-2) $ -波Darboux变换, 可以从平面波种子解求得方程(1)的怪波与呼吸子的混合作用解; 当M取3到$ N-1 $ 之间的整数时, 可以得到更丰富的局域波作用解. -
在这一节将讨论如何应用定理1中的
$ (M, N-M) $ -波Darboux变换得到方程(1)丰富的局域波解. 下面考虑两类种子解情形下的Lax对的解.情形1 把种子解
$ r_{n} = 0 $ 代入(2)式中得到Lax对的解为$ {\varphi_{n}} = (\phi_{n}, \psi_{n})^{\rm T} = (\lambda^{2 n}{\rm e}^\xi, {\rm e}^{-\xi})^{\rm T} $ , 这里$ \xi = \left(-\dfrac{{\rm i}\lambda^4}{2}+{\rm i}\lambda^2-\dfrac{\rm i}{2\lambda^4}+\dfrac{\rm i}{\lambda^2}\right)t $ .情形2 把种子解
$ r_{n} = c {\rm e}^{{\rm i}(c^2+1)(-6 c^2+2)t} $ 代入(2)中得到Lax对的解为$ \begin{array}{l} {\varphi_{n}} = \left( \begin{array}{*{20}{c}} \phi_{n} \\ \psi_{n} \end{array} \right) = \left( \begin{array}{*{20}{c}} C_1 \tau_+^n {\rm e}^{\rho_+ t+\delta(\varepsilon)}+C_2 \tau_-^n {\rm e}^{\rho_- t-\delta(\varepsilon)} \\ \dfrac{{\rm e}^{2{\rm i}(c^2+1)(3c^2-1)t}}{c}\left[ C_1 (-\lambda^2+\tau_{+}^n)\tau_+^n {\rm e}^{\rho_{+}\tau+\delta(\varepsilon)}+C_2 (-\lambda^2+\tau_{-}^n)\tau_{-}^n{\rm e}^{\rho_{-}t-\delta(\varepsilon)}\right] \end{array} \right), \end{array} $ 其中
$ \begin{aligned} &\tau_{\pm} = \dfrac{\lambda^2}{2}+\dfrac{1}{2}\pm\dfrac{1}{2}\sqrt{\lambda^4-4\lambda^4 c^2-2\lambda^2+1} ,\\ & \rho_{\pm} = -\dfrac{{\rm i}(c^2+1)(6 c^2-2)\lambda^4\pm {\rm i}(\lambda^2-1)( 2\lambda^2 c^2+\lambda^4+1)\sqrt{\lambda^4-4\lambda^2 c^2-2\lambda^2+1}}{2\lambda^4}, \\ & \delta(\varepsilon) = \sqrt{-4\lambda^2 c^2+\lambda^4-2\lambda^2+1}\sum\limits_{\omega = 0}^{N}(e_{\omega}+{\rm i} d_{\omega})\varepsilon^{2 k}, \end{aligned} $ 这里
$ e_{\omega}, d_{\omega} $ 是任意的实常数. 为了得到怪波解, 需要对情形2的解进行泰勒展开, 为此令$ \lambda = -c+\sqrt{1+c^2}+\varepsilon^2, C_1\! =\! -C_2 = {1}/{\varepsilon}, c = {3}/{4} $ , 在$ \varepsilon = 0 $ 处对情形2的解进行泰勒展开, 有$ {{\varphi}}(\varepsilon^2) = $ $\displaystyle\sum\limits_{\omega = 0}^{\infty}{{{\varphi}}^{({{\omega}})}}\varepsilon^{(2 k)} = {{{\varphi}}^{{\bf{(0)}}}}+{{{\varphi}}^{{\bf{(1)}}}}\varepsilon^2+{{{\varphi}}^{{\bf{(2)}}}}\varepsilon^4+... $ , 这里略去$ {{{\varphi}}^{({{i}})}}(i \geqslant 1) $ 的展开式, 仅列出$ {{{\varphi}}^{{\bf{(0)}}}} = (\phi^{(0)}, \psi^{(0)})^{\rm T }$ , 这里$ \phi^{(0)} = \dfrac{ \sqrt{15} }{80} \left(\dfrac{5}{8}\right)^n(64{\rm i} n-645 t){\rm e}^{-\textstyle\frac{275 {\rm i}}{256}t} $ ,$\psi^{(0)} = $ $ \dfrac{ \sqrt{15} }{240} \left(\dfrac{5}{8}\right)^n(192 {\rm i} n+320 {\rm i}-1935 t){\rm e}^{\textstyle\frac{275 {\rm i}}{256}t} $ . 下面将讨论三种特殊情况的Darboux变换得到不同的局域波解. -
根据定理1, 当
$ M = N $ , 广义$ (M, N-M) $ -波Darboux变换可以约化为通常的N-波Darboux变换, 当$ N = 1 $ , 使用情形1中Lax对的解, 一孤子解$ \widetilde{r}_{n} = b^{(1)}_{n+1} $ 可以整理为$ \widetilde{r}_n = -\frac{\lambda_1(|\lambda_1|^4-1){\rm e}^{\eta_1-\eta^*_1}}{2 \lambda^{*2}_1|\lambda_1|} \rm sech(\eta_1+\eta^*_1+\ln|\lambda_1|), $ 这里
$ \eta_1 = n\ln\lambda_1+\left(-\dfrac{{\rm i}\lambda^4_1}{2}+{\rm i}\lambda^2_1-\dfrac{\rm i}{2\lambda^4_1}+\dfrac{\rm i}{\lambda^2_1}\right)t $ . 关于一孤子解(9)式, 我们能计算出一些重要物理量: 振幅为$ \left|\dfrac{(|\lambda_1 |^4-1)}{2 |\lambda_1 |^2}\right| $ , 波数为$ 2\ln|\lambda_1 | $ , 波宽为$ \dfrac{1}{2\ln|\lambda_1 |} $ , 速度为$ \dfrac{ {\rm i}\left(\!-\dfrac{\lambda^4_1}{2}\!+\!\lambda^2_1\!-\!\dfrac{1}{2\lambda^4_1}\!+\!\dfrac{1}{\lambda^2_1}\!+\!\dfrac{\lambda^{*4}_1}{2}\!-\!\lambda^{*2}_1\!+\!\dfrac{1}{2\lambda^{*4}_1}\!-\!\dfrac{1}{\lambda^{*2}_1}\!\right)}{2\ln|\lambda_1 |},$ 能量为
$ \dfrac{(|\lambda_1 |^4-1)^2}{4 |\lambda_1 |^6 |\ln|\lambda_1 ||} $ , 初相为$ \ln|\lambda_1 | $ . 当参数$ \lambda_1 = -1+2 {\rm i} $ 时, 图1(a1)显示了一孤子解(9)式的钟型亮孤子结构; 下面进一步通过向前差分格式[31] 研究一孤子解(9)式的动力学传播稳定性, 图1(b1)显示了一孤子精确解(9)式作为初值时的数值演化, 从图像上看, 数值解几乎重现了精确解, 这也说明(9)式的精确解和数值算法的正确性; 图1(c1)显示了加入了$ 2\% $ 小噪声时一孤子解的数值演化, 可以看出小噪声在较短的时间内对一孤子的稳定传播影响很小. 这里需要说明的是, 在下面的图像中左边第1列表示精确解, 中间第2列表示没有加入噪声的数值解, 右边第3列表示加入$ 2\% $ 小噪声的数值解.
图 1 (a1)−(c1)取参数
$ \lambda_1 = -1+2 {\rm i} $ 时的一孤子解; (a2)−(c2)取参数$ \lambda_1 = 3+{\rm i}, \lambda_2 = 1+2 {\rm i} $ 时的二孤子解. 左列: 精确解; 中列: 数值解; 右列: 加$ 2\% $ 小噪声的数值解. 圈中的数字1和2分别表示参数$ \lambda_1 $ 和$ \lambda_2 $ 对应的孤子Figure 1. (a1)−(c1) One-soliton solution with parameter
$ \lambda_1 = -1+2 {\rm i} $ ; (a2)−(c2) two-soliton solution with parameters$ \lambda_1 = 3+{\rm i}, \lambda_2 = 1+2 {\rm i} $ . Left column: Exact solutions; Middle column: Numerical solutions without noise; Right column: Numerical solutions with$ 2\% $ small noise.当
$ N = 2 $ , 由定理1中(8)式, 我们能够得到二孤子解, 这里为节省篇幅, 不列出具体表达式, 通过渐近分析, 可以知道二孤子之间的碰撞是弹性作用. 图1(a2)显示了二孤子的弹性作用结构, 两孤子在碰撞前后, 形状振幅速度没有发生改变. 图1(b2)显示了二孤子没有加入噪声时的数值演化, 图1(c2)显示了二孤子加入$ 2\% $ 小噪声时的数值演化. 类似于一孤子, 在较短的时间内小噪声对二孤子的碰撞作用影响很小.下面利用N-波Darboux变换和情形2时的Lax对的解研究方程(1)的呼吸子解. 这里取参数
$ C_1 = -C_2 = c = 1, \delta(\varepsilon) = 0 $ . 当$ N = 1, 2 $ , 由(8)式可以得到一呼吸子和二呼吸子解. 图2(a1)–图2(c1)和图(a2)–图2(c2)分别显示了一和二呼吸子的精确解、数值解以及加$ 2\% $ 小噪声的数值解. 在较短的时间内, 一、二呼吸子都有较好的数值稳定性, 随着时间的增加, 不论是否有噪声, 数值演化都显示出较大的振动, 表现出不稳定.
图 2 (a1)−(c1)取参数
$ \lambda_1 ={11}/{4} $ 时的一呼吸子解; (a2)−(c2)取参数$ \lambda_1 = 2, \lambda_2 = 3 $ 时的二呼吸子解. 左列: 精确解; 中列: 数值解; 右列: 加$ 2\% $ 小噪声的数值解. 圈中的数字1和2分别表示参数$ \lambda_1 $ 和$ \lambda_2 $ 对应的呼吸子Figure 2. (a1)−(c1) One-breather solution with parameter
$ \lambda_1 ={11}/{4} $ ; (a2)−(c2) two-breather solution with parameters$ \lambda_1 = 2, \lambda_2 = 3 $ . Left column: Exact solutions; Middle column: Numerical solutions without noise; Right column: Numerical solutions with$ 2\% $ small noise. -
根据定理1, 当
$ M = 1 $ 时, 我们仅使用一个谱参数, 利用$ (1, N-1) $ -波Darboux变换和情形2时Lax对的解结合泰勒展开式可以得到方程(1)的怪波解. 这里取参数$ C_1 = -C_2 = {1}/{\varepsilon} $ ,$ c ={3}/{4} $ , 当$ N = 1, 2 $ , 由(8)式可以得到一阶和二阶怪波解.当
$ N = 1 $ 时, 由$ (1, 0) $ -波Darboux变换可以得到一阶怪波解$ \widetilde{r}_n = b^{(1)}_{n+1}+r_n a^{(0)}_{n+1} $ , 通过整理为$ \widetilde{r}_n = \frac{3715200{\rm i} t+110592n^2+11232675t^2+294912n-61440}{16384(3n+4)^2+14976900t^2+65536}{\rm e}^{-\frac{275{\rm i}}{128}t},$ 该解是非奇性的, 并且在时间和空间上都是局域的. 图3(a1)—图3(c1)显示了一阶怪波(10)式的精确解、数值解以及加
$ 2\% $ 小噪声的数值解. 在较短的时间内, 一阶怪波解具有较好的数值稳定性.
图 3 (a1)—(c1)一阶怪波解; (a2)—(c2)取参数
$ e_1 = d_1 = 0 $ 时具有强作用的二阶怪波解; (a3)—(c3)取参数$ e_1 = 100, d_1 = 0 $ 时具有弱作用的二阶怪波解. 左列: 精确解; 中列: 数值解; 右列: 加$ 2\% $ 小噪声的数值解Figure 3. (a1)−(c1) First-order rogue wave solution; (a2)−(c2) strong interaction second-order rogue wave solution with
$ e_1 = d_1 = 0 $ ; (a3)−(c3) weak interaction second-order rogue wave solution with$ e_1 = 100, d_1 = 0 $ . Left column: Exact solutions; Middle column: Numerical solutions without noise; Right column: Numerical solutions with$ 2\% $ small noise.由定理1中(8)式, 当
$ N = 2 $ 时, 由$ (1, 1) $ -波Darboux变换可以得到二阶怪波解, 这里为节省篇幅, 不列出具体表达式. 图3(a2)—图3(c2)和图(a3)—图3(c3)分别显示了具有强作用和弱作用的二阶怪波解的精确解、数值解以及加$ 2\% $ 小噪声的数值解. 在较短的时间内, 二阶怪波具有较好的数值稳定性. 随着时间的增加, 不论是否有噪声, 数值演化都表现出较强的振动和不稳定. -
根据定理1, 当
$ M = 2 $ 时, 我们将使用两个谱参数, 利用$ (2, N-2) $ -波Darboux变换和情形2时Lax对的解结合泰勒展开式, 可以得到方程(1)的怪波与呼吸子的混合作用解. 由(8)式, 当$ N = 2 $ 时, 需要说明的是, 对于两个谱参数, 一个谱参数使用情形2时Lax对的解结合泰勒展开式, 参数选取与3.2小节中求怪波的参数相同, 对另一个谱参数仍使用情形2时Lax对的解但不使用泰勒展开, 参数$ C_1 = -C_2 = 1, c = {3}/{4}, \delta(\varepsilon) = 0 $ . 由$ (2, 0) $ -波Darboux变换可以得到一阶怪波和一呼吸子的混合作用解, 这里不列出具体表达式. 如果在泰勒展开式中对变量n做一定的平移变换, 就可以控制怪波的位置. 图4(a1)—图4(c1)和图4(a2)—图4(c2)分别显示了一阶怪波和一呼吸子的强作用和弱作用情况下的精确解、数值解以及加$ 2\% $ 小噪声的数值解. 两行图的区别是弱作用时, 对怪波沿n轴正半轴平移了10个单位, 实现了怪波和呼吸子的分离. 数值模拟显示了这些混合解的正确性, 在较短时间内, 不论是否有噪声, 混合作用解均具有较好的数值稳定性.
图 4 (a1)—(c1)一呼吸子和一阶怪波的混合强作用; (a2)—(c2)一呼吸子和一阶怪波的混合弱作用. 左列: 精确解; 中列: 数值解; 右列: 加
$ 2\% $ 小噪声的数值解. 圈中的数字1表示呼吸子, 数字2表示怪波Figure 4. (a1)−(c1) Mixed strong interaction between one-breather and first-order rogue wave; (a2)−(c2) mixed weak interaction between one-breather and first-order rogue wave. Left column: Exact solutions; Middle column: Numerical solutions without noise; Right column: Numerical solutions with
$ 2\% $ small noise. -
本文构造了高阶Ablowitz-Ladik方程的广义
$ (M, N-M) $ -波Darboux变换, 借助符号计算和数值计算, 通过不同特殊情形的Darboux变换, 得到了该模型多孤子解、多呼吸子解、高阶怪波解以及怪波与呼吸子的相互作用解, 并且应用数值模拟研究了它们的稳定性. 需要说明的是当$ 3\leqslant M\leqslant $ $N-1 $ , 利用广义$ (M, N-M) $ -波Darboux变换, 也可以得到更丰富的局域波作用结构, 本文不做讨论. 本文给出的广义$ (M, N-M) $ -波Darboux变换方法为研究Lax可积模型丰富的局域波提供了很好的解决思路, 得到的结果更为全面, 比通常的Darboux变换更具有普遍性, 关于方程(1)的局域波结构的研究丰富了非线性微分差分方程的已知结果, 将为实际应用提供可靠的理论支持, 希望本文得到的结果为解释实际的物理现象提供理论依据.感谢北京信息科技大学理学院孤子与可积系统讨论班成员的讨论.
-
本文构造了一类高阶Ablowitz-Ladik方程的广义
$(M, N-M)$ -波Darboux变换, 借助符号计算从不同背景出发研究了该模型丰富的局域波解, 并利用数值模拟研究了这些解的动力学稳定性.-
关键词:
- 高阶Ablowitz-Ladik方程 /
- 广义(M, N-M)-波Darboux变换 /
- 局域波 /
- 怪波
[1] 尤福财 2008 博士学位论文 (上海: 上海大学)
You F C 2008 Ph. D. Dissertation (Shanghai: Shanghai University) (in Chinese)
[2] Toda M 1967 J. Phys. Soc. Jpn. 23 501
[3] 李翊神 1999 孤子与可积系统 (上海: 上海科技教育出版社) 第25页
Li Y S 1999 Soliton and Integrable System (Shanghai: Shanghai Scientific and Technological Education Publishing House) p25 (in Chinese)
[4] 広田良吾 著 (王红艳, 李春霞, 赵俊霄, 虞国富 译) 2004 孤子理论中的直接方法(北京: 清华大学出版社) 第59页
Hirota R (translated by Wang H Y, Li C X, Zhao J X, Yu G F) 2004 The Direct Method in Soliton Theory (New York: Cambridge University Press) p59
[5] Yue Y F, Huang L L, Chen Y 2018 Comput. Math. Appl. 75 2538
[6] Huang L L, Yue Y F, Chen Y 2018 Comput. Math. Appl. 74 831
[7] Liu C, Yang Z Y, Yang W L, Akhmediev N 2019 J. Opt. Soc. Am. B 36 1294
[8] Ren Y, Wang X, Liu C, Yang Z Y, Yang W L 2018 Commun. Nonlinear Sci. Numer. Simul. 63 161
[9] Liu C, Ren Y, Yang Z Y, Yang W L 2017 Chaos 27 083120
[10] Liu C, Yang Z Y, Yang W L 2017 Chaos 28 083110
[11] Ren Y, Liu C, Yang Z Y, Yang W L 2018 Phys. Rev. E 98 062223
[12] Li H, Lou S Y 2019 Chin. Phys. Lett. 36 050501
[13] Zhao L C, Yang Z Y, Yang W L 2019 Chin. Phys. B 28 010501
[14] Liu X S, Zhao L C, Duan L, Peng G, Yang Z Y, Yang W L 2018 Chin. Phys. Lett. 35 020501
[15] 李淑青, 杨光晔, 李禄 2014 物理学报 63 104215
Li S Q, Yang G Y, Li L 2014 Acta Phys. Sin. 63 104215
[16] 孙俊超, 张宗国, 董焕河, 杨红卫 2014 物理学报 21 210201
Sun J C, Zhang Z G, Dong H H, Yang H W 2014 Acta Phys. Sin. 21 210201
[17] Wen X Y, Yan Z Y 2015 Chaos 25 123115
[18] Wen X Y, Yang Y Q, Yan Z Y 2015 Phys. Rev. E 92 012917
[19] Wen X Y, Yan Z Y, Yang Y Q 2016 Chaos 26 063123
[20] Wen X Y, Yan Z Y 2017 Commun. Nonlinear Sci. Numer. Simul. 43 311
[21] Wang H T, Wen X Y, Wang D S 2019 Wave Motion 91 102396
[22] Huang L L, Qiao Z J, Chen Y 2018 Chin. Phys. B 27 020201
[23] Zhang G, Yan Z Y, Wen X Y, Chen Y 2017 Phys. Rev. E 95 042201
[24] Zhang G, Yan Z Y, Wen X Y 2018 Physica D 366 27
[25] Yang J, Zhu Z N 2018 Chin. Phys. Lett. 35 090201
[26] Wen X Y, Yan Z, Malomed B A 2016 Chaos 26 123110
[27] Ablowitz M J, Clarkson P A 1991 Soliton, Nonlinear Evolution Equations and Inverse Scattering (New York: Cambridge University Press) p55
[28] Ankiewicz A, Akhmediev N, Soto-Crespo J M 2010 Phys. Rev. E 82 026602
[29] Ankiewicz A, Akhmediev N, Lederer F 2011 Phys. Rev. E 83 056602
[30] Ohta Y, Yang J 2014 J. Phys. A: Math. Theor. 47 255201
[31] Lee L, Lyng G, Vankova I 2012 Physica D 241 1767
-
图 1 (a1)−(c1)取参数
$ \lambda_1 = -1+2 {\rm i} $ 时的一孤子解; (a2)−(c2)取参数$ \lambda_1 = 3+{\rm i}, \lambda_2 = 1+2 {\rm i} $ 时的二孤子解. 左列: 精确解; 中列: 数值解; 右列: 加$ 2\% $ 小噪声的数值解. 圈中的数字1和2分别表示参数$ \lambda_1 $ 和$ \lambda_2 $ 对应的孤子Fig. 1. (a1)−(c1) One-soliton solution with parameter
$ \lambda_1 = -1+2 {\rm i} $ ; (a2)−(c2) two-soliton solution with parameters$ \lambda_1 = 3+{\rm i}, \lambda_2 = 1+2 {\rm i} $ . Left column: Exact solutions; Middle column: Numerical solutions without noise; Right column: Numerical solutions with$ 2\% $ small noise.
图 2 (a1)−(c1)取参数
$ \lambda_1 ={11}/{4} $ 时的一呼吸子解; (a2)−(c2)取参数$ \lambda_1 = 2, \lambda_2 = 3 $ 时的二呼吸子解. 左列: 精确解; 中列: 数值解; 右列: 加$ 2\% $ 小噪声的数值解. 圈中的数字1和2分别表示参数$ \lambda_1 $ 和$ \lambda_2 $ 对应的呼吸子Fig. 2. (a1)−(c1) One-breather solution with parameter
$ \lambda_1 ={11}/{4} $ ; (a2)−(c2) two-breather solution with parameters$ \lambda_1 = 2, \lambda_2 = 3 $ . Left column: Exact solutions; Middle column: Numerical solutions without noise; Right column: Numerical solutions with$ 2\% $ small noise.
图 3 (a1)—(c1)一阶怪波解; (a2)—(c2)取参数
$ e_1 = d_1 = 0 $ 时具有强作用的二阶怪波解; (a3)—(c3)取参数$ e_1 = 100, d_1 = 0 $ 时具有弱作用的二阶怪波解. 左列: 精确解; 中列: 数值解; 右列: 加$ 2\% $ 小噪声的数值解Fig. 3. (a1)−(c1) First-order rogue wave solution; (a2)−(c2) strong interaction second-order rogue wave solution with
$ e_1 = d_1 = 0 $ ; (a3)−(c3) weak interaction second-order rogue wave solution with$ e_1 = 100, d_1 = 0 $ . Left column: Exact solutions; Middle column: Numerical solutions without noise; Right column: Numerical solutions with$ 2\% $ small noise.
图 4 (a1)—(c1)一呼吸子和一阶怪波的混合强作用; (a2)—(c2)一呼吸子和一阶怪波的混合弱作用. 左列: 精确解; 中列: 数值解; 右列: 加
$ 2\% $ 小噪声的数值解. 圈中的数字1表示呼吸子, 数字2表示怪波Fig. 4. (a1)−(c1) Mixed strong interaction between one-breather and first-order rogue wave; (a2)−(c2) mixed weak interaction between one-breather and first-order rogue wave. Left column: Exact solutions; Middle column: Numerical solutions without noise; Right column: Numerical solutions with
$ 2\% $ small noise. -
[1] 尤福财 2008 博士学位论文 (上海: 上海大学)
You F C 2008 Ph. D. Dissertation (Shanghai: Shanghai University) (in Chinese)
[2] Toda M 1967 J. Phys. Soc. Jpn. 23 501
[3] 李翊神 1999 孤子与可积系统 (上海: 上海科技教育出版社) 第25页
Li Y S 1999 Soliton and Integrable System (Shanghai: Shanghai Scientific and Technological Education Publishing House) p25 (in Chinese)
[4] 広田良吾 著 (王红艳, 李春霞, 赵俊霄, 虞国富 译) 2004 孤子理论中的直接方法(北京: 清华大学出版社) 第59页
Hirota R (translated by Wang H Y, Li C X, Zhao J X, Yu G F) 2004 The Direct Method in Soliton Theory (New York: Cambridge University Press) p59
[5] Yue Y F, Huang L L, Chen Y 2018 Comput. Math. Appl. 75 2538
[6] Huang L L, Yue Y F, Chen Y 2018 Comput. Math. Appl. 74 831
[7] Liu C, Yang Z Y, Yang W L, Akhmediev N 2019 J. Opt. Soc. Am. B 36 1294
[8] Ren Y, Wang X, Liu C, Yang Z Y, Yang W L 2018 Commun. Nonlinear Sci. Numer. Simul. 63 161
[9] Liu C, Ren Y, Yang Z Y, Yang W L 2017 Chaos 27 083120
[10] Liu C, Yang Z Y, Yang W L 2017 Chaos 28 083110
[11] Ren Y, Liu C, Yang Z Y, Yang W L 2018 Phys. Rev. E 98 062223
[12] Li H, Lou S Y 2019 Chin. Phys. Lett. 36 050501
[13] Zhao L C, Yang Z Y, Yang W L 2019 Chin. Phys. B 28 010501
[14] Liu X S, Zhao L C, Duan L, Peng G, Yang Z Y, Yang W L 2018 Chin. Phys. Lett. 35 020501
[15] 李淑青, 杨光晔, 李禄 2014 物理学报 63 104215
Li S Q, Yang G Y, Li L 2014 Acta Phys. Sin. 63 104215
[16] 孙俊超, 张宗国, 董焕河, 杨红卫 2014 物理学报 21 210201
Sun J C, Zhang Z G, Dong H H, Yang H W 2014 Acta Phys. Sin. 21 210201
[17] Wen X Y, Yan Z Y 2015 Chaos 25 123115
[18] Wen X Y, Yang Y Q, Yan Z Y 2015 Phys. Rev. E 92 012917
[19] Wen X Y, Yan Z Y, Yang Y Q 2016 Chaos 26 063123
[20] Wen X Y, Yan Z Y 2017 Commun. Nonlinear Sci. Numer. Simul. 43 311
[21] Wang H T, Wen X Y, Wang D S 2019 Wave Motion 91 102396
[22] Huang L L, Qiao Z J, Chen Y 2018 Chin. Phys. B 27 020201
[23] Zhang G, Yan Z Y, Wen X Y, Chen Y 2017 Phys. Rev. E 95 042201
[24] Zhang G, Yan Z Y, Wen X Y 2018 Physica D 366 27
[25] Yang J, Zhu Z N 2018 Chin. Phys. Lett. 35 090201
[26] Wen X Y, Yan Z, Malomed B A 2016 Chaos 26 123110
[27] Ablowitz M J, Clarkson P A 1991 Soliton, Nonlinear Evolution Equations and Inverse Scattering (New York: Cambridge University Press) p55
[28] Ankiewicz A, Akhmediev N, Soto-Crespo J M 2010 Phys. Rev. E 82 026602
[29] Ankiewicz A, Akhmediev N, Lederer F 2011 Phys. Rev. E 83 056602
[30] Ohta Y, Yang J 2014 J. Phys. A: Math. Theor. 47 255201
[31] Lee L, Lyng G, Vankova I 2012 Physica D 241 1767
| 引用本文: |
|
| Citation: |
|
计量
- 文章访问数: 1290
- PDF下载量: 68
- 被引次数: 0
