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

高阶Ablowitz-Ladik方程的局域波解及稳定性分析

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

摘要:本文构造了一类高阶Ablowitz-Ladik方程的广义$(M, N-M)$-波Darboux变换, 借助符号计算从不同背景出发研究了该模型丰富的局域波解, 并利用数值模拟研究了这些解的动力学稳定性.
关键词: 高阶Ablowitz-Ladik方程/
广义(M, N-M)-波Darboux变换/
局域波/
怪波

English Abstract


--> --> -->
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, NM)波Darboux变换从研究方程(1)的精确解出发, 寻找其新奇的局域波结构, 特别是不同类型的局域波相互作用的新奇局域波结构.
本文的主要结构安排如下: 第二节构造出方程(1)的Lax对和广义$ (M, N-M) $波Darboux变换; 第三节将应用广义$ (M, N-M) $-波Darboux变换给出方程(1)不同类型的局域波解, 并通过数值模拟研究其传播稳定性; 最后是本文的结论.
由AKNS方法[3], 可以构造方程(1)的Lax对如下
$ 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 =0 $ (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变换得到不同的局域波解.
2
3.1.应用N-波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 $对应的孤子
Figure1. (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 $对应的呼吸子
Figure2. (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.

2
3.2.应用(1, N–1)-波Darboux变换得到怪波解
-->根据定理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\% $小噪声的数值解
Figure3. (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\% $小噪声的数值解. 在较短的时间内, 二阶怪波具有较好的数值稳定性. 随着时间的增加, 不论是否有噪声, 数值演化都表现出较强的振动和不稳定.
2
3.3.应用(2, N–2)-波Darboux变换得到怪波与呼吸子的混合作用解
-->根据定理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表示怪波
Figure4. (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)的局域波结构的研究丰富了非线性微分差分方程的已知结果, 将为实际应用提供可靠的理论支持, 希望本文得到的结果为解释实际的物理现象提供理论依据.
感谢北京信息科技大学理学院孤子与可积系统讨论班成员的讨论.
相关话题/传播 结构 种子 物理 系统

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 可积谐振系统中的极端波事件研究进展
    摘要:从微观角度上讲,单个极端异常波事件可视为可积模型方程的时空局域有理函数解.本文主要讨论了三类典型的可积谐振相互作用模型(即长波短波谐振方程,三波谐振相互作用方程,非线性薛定谔和麦克斯韦-布洛赫方程)的基阶Peregrine异常波解及其相关研究进展;明确指出了这些基阶异常波解形式具有普适性,可推 ...
    本站小编 Free考研考试 2021-12-29
  • 可积系统多孤子解的全反演对称表达式
    摘要:多孤子解是非线性数学物理系统的基本激发模式.文献中存在各种类型的表达式,如广田(Hirota)形式,朗斯基(Wronskian)或双朗斯基形式和法夫(Phaffian)形式.最近在多地系统的研究中,我们发现使用一种全新但等价的形式具有极为简洁和方便的优点.本文主要综述多种类型可积非线性系统的多 ...
    本站小编 Free考研考试 2021-12-29
  • 二维材料/铁电异质结构的研究进展
    摘要:二维材料是一类具有原子层厚度的层状材料,拥有独特的电学、磁学、光学和力学性能.以石墨烯和过渡金属硫族化合物为代表的二维材料展现出迁移率高、能带可调、可见光透过率高等特点,是近年来微纳科学领域的前沿热点.将二维材料与各种功能材料,如SiO2绝缘体、半导体、金属、有机化合物等结合,可以深化和拓宽二 ...
    本站小编 Free考研考试 2021-12-29
  • 基于密度泛函理论的外电场下C<sub>5</sub>F<sub>10</sub>O的结构及其激发特性
    摘要:采用密度泛函(DFT)B3LYP/6-311g(d)对C5F10O分子进行几何结构优化,研究外加电场(0-0.03a.u.,1a.u.=5.142×1011V/m)对分子的几何结构、能量、前线轨道能级、红外光谱的影响.在相同基组下,采用TD-DFT方法计算和分析C5F10O的轨道成分和激发特性 ...
    本站小编 Free考研考试 2021-12-29
  • 具有内参考热补偿功能的三层膜结构微球腔折射率传感器
    摘要:光学微腔在高灵敏度传感中有着重要的应用前景,而在传感中热漂移是制约其走向实用的重要因素.本文提出了一种镀有三层膜结构的微球腔,可以在实现高灵敏度折射率传感的同时,具备内参考热补偿功能.该结构由内到外分别涂覆折射率为高、低、高的薄膜,内外两高折射率层可以分别支持各自的回音壁模式,称之为内层模式和 ...
    本站小编 Free考研考试 2021-12-29
  • 准Λ型四能级系统选择反射光谱
    摘要:运用密度矩阵理论研究了气固界面准$\Lambda$型四能级原子系统的非线性选择反射光谱.基于刘维尔方程给出了一阶近似条件下探测光场的解析式.在探测场为弱场时,分析了信号场拉比频率、失谐量和耦合场失谐量对反射光谱线型的影响.数值模拟表明:信号场参与产生的选择反射峰线宽可以利用信号场拉比频率进行调 ...
    本站小编 Free考研考试 2021-12-29
  • He离子辐照对石墨烯微观结构及电学性能的影响
    摘要:本文采用5.4keV不同剂量的He离子辐照单层石墨烯,利用X射线光电子能谱(XPS)、拉曼光谱(Raman)和半导体参数分析仪表征辐照前后石墨烯微观结构和电学性能变化.研究结果表明:随着辐照剂量增大,单层石墨烯的缺陷密度逐渐增加,当辐照剂量增至1.6×1013He+/cm2,石墨烯开始由纳米晶 ...
    本站小编 Free考研考试 2021-12-29
  • 隔板对流系统的热流特性及热量输入与传递特性
    摘要:采用DNS方法对隔板对流装置进行模拟计算,研究系统中热流特性以及热量输入与传递特性.讨论了热流的纵向和横向输运特性,在此基础上对传热通道和狭缝区域的热通量以及对应底板外界输入热通量进行了定量化分析.研究结果表明,通道中低温流体向下冲击底板而后转入水平运动,流入狭缝区域并不断被加热,只进行水平的 ...
    本站小编 Free考研考试 2021-12-29
  • Shimizu-Morioka系统与Finance系统生成Lorenz混沌的微分几何策略
    摘要:从一种受控混沌系统生成另一混沌系统可增强保密通信的安全性,具备潜在应用前景.研究了如何通过状态变换以及单输入反馈,驱使受控Shimizu-Morioka系统与受控Finance系统生成Lorenz混沌动态.主要方法是运用微分几何理论,将上述三种系统等价转换为下三角形式,并尽量简化和一致化其方程 ...
    本站小编 Free考研考试 2021-12-29
  • TiAl电子态结构的<i>ab initio</i>计算
    摘要:应用完全活动基自洽场方法,结合N电子价态微扰近似(NEVPT2),对TiAl金属二聚体的基态和若干最低电子激发态的势能曲线进行了计算.完全活动空间由Al的3个价电子(3s23p1)轨道和Ti的4个价电子(3d24s2)轨道构成,计算基组选用Karlsruhegroup的价分裂全电子基组def2 ...
    本站小编 Free考研考试 2021-12-29