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

四阶色散非线性薛定谔方程的明暗孤立波和怪波的形成机制

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

摘要:本文研究了四阶色散非线性薛定谔方程的明暗孤立波和怪波的形成机制, 该模型既可以模拟高速光纤传输系统中超短脉冲的非线性传输和相互作用, 又可以描述具有八极与偶极相互作用的一维海森堡铁磁链的非线性自旋激发现象. 本文首先通过对四阶色散非线性薛定谔方程的相平面分析, 发现由其约化得到的二维平面自治系统具有同宿轨道和异宿轨道, 并在相应条件下求得了方程的明孤立波解和暗孤立波解, 从而揭示了同异宿轨道和孤立波解的对应关系; 其次, 基于非零背景平面上的精确一阶呼吸子解, 给出了呼吸子的群速度和相速度的显式表达式, 进而分析得出呼吸子的速度存在跳跃现象. 最后, 为了验证在跳跃点处呼吸子可以转化为怪波, 将呼吸子解在速度跳跃条件下取极限获得了一阶怪波解, 从而证实怪波的产生与呼吸子速度的不连续性有关.
关键词: 同异宿轨道/
孤立波/
怪波/
速度不连续性

English Abstract


--> --> -->
孤子、呼吸子和怪波是自然界中三种典型的非线性波现象, 由于在非线性光学、玻色-爱因斯坦凝聚、大气物理以及等离子体等领域中的潜在应用而受到密切关注[1-20]. 其中, 孤子在群速度色散和非线性效应均衡作用下可以保持不变的波形和速度稳定地传播, 从而在长距离、大容量和高速率光纤通信系统中有着重要应用[9]. 呼吸子的产生源于小扰动引起的平面波调制不稳定性[9]. 研究发现, 经典的非线性薛定谔方程(NLSE)具有两类呼吸子解: Kuznetsov-Ma呼吸子(简称KMB)[10]和Akhmediev呼吸子(简称AB)[11]. 另外, 怪波在数学上可以通过NLSE的有理解(即KMB或AB呼吸子解的一种退化情形[10,11])进行描述, 同时也可以作为在周期边界条件下不稳定Stokes波的一种同宿轨行为[12,13]. 怪波最早被用来描述在海洋中发现的来无影去无踪的大振幅畸形波[14,15], 近年来也在非线性光学和水槽实验中被观察到[16,17]. 通常情况下, ****们认为只有非线性可以解释由能量的大量聚集而形成的高于周围波平均高度的单个水峰现象[18,19], 并且认为调制不稳定性是怪波产生的初始过程[20].
近几年, 人们在实验和理论上开展了大量关于怪波形成机制的研究. 例如, 文献[21]在一个大的定向波池实验中讨论了平面波的调制不稳定性, 该研究为有限水深条件下怪波的产生给出了一定解释; 文献[22]通过光学实验研究了在连续波泵浦区域超连续统的时间特性, 从而指出怪波可由孤子之间的相互碰撞产生; 文献[23]借助数值模拟发现拉曼效应和三阶色散可以激发怪波的形成; 文献[24,25]在理论上指出由速度跳跃可以使得呼吸子演化为怪波, 并且非线性偏微分方程的行波解与对应的常微分方程的轨道存在对应关系, 即, 常微分方程中的同宿轨道和异宿轨道分别与非线性波方程中的孤立波和冲击波相对应, 而对于非线性薛定谔型方程, 其约化后的常微分系统的同宿轨道和异宿轨道则分别对应钟型的明孤立波解、冲击波解或倒钟型的暗孤立波解[26,27].
本文将从结构不连续性角度研究如下四阶色散非线性薛定谔方程[28-30]的孤立波解和怪波解形式机制:
$ \begin{split} & {\rm i}q_t \, +\alpha_1q_{xx}+\alpha_2q|q|^2 \\ & +\varepsilon^2(\alpha_3q_{xxxx}+\alpha_4|q|^2q_{xx}+\alpha_5q^2{q_{xx}}^* \\ & +\alpha_6{q_x}^2q^*+\alpha_7q|q_x|^2+\alpha_8|q|^4q) = 0, \end{split} $
其中$ q(x, t) $代表波的慢变包络振幅, xt是归一化的空间和时间变量, $ \varepsilon^2 $是一个表示高阶线性和非线性强度的无量纲小参数, αj ($ j = 1, 2, \cdots, 8 $)为实参数. 方程(1)主要来源于光纤光学和磁力学: 在光学中可以模拟高速光纤传输系统中超短脉冲的非线性传播及相互作用[31], 在磁力学中亦可描述具有八极与偶极相互作用的一维海森堡铁磁链的非线性自旋激发现象[32]. 特别地, 当参数取值为α1 = 1, α2 = 2, α3 = 1, α4 = 8, α5 = 2, α6 = 6, α7 = 4, α8 = 6时, 该方程是一个可积模型, 具有Lax对和无穷多守恒律[33]. 文献[34]借助双线性方法获得了方程(1)的双线性形式及N孤子解, 并通过渐近分析揭示了孤子的弹性碰撞现象. 文献[30]利用达布变换方法构造了呼吸子解和高阶怪波解的精确表达式.
本文将首先利用相平面分析方法研究方程(1)经约化后得到的常微分方程的动力学性质, 分别给出同宿轨道和异宿轨道对应的参数条件, 并在相应条件下求得了方程(1)的明孤立波解和暗孤立波解, 揭示了同异宿轨道与孤立波解之间的对应关系; 其次, 基于文献[30]获得的呼吸子解表达式, 分析了呼吸子的群速度和相速度随参数的变化关系, 进而确定速度发生跳跃的参数条件; 为了验证在速度跳跃的参数条件下呼吸子可以转化成怪波, 在该条件下对呼吸子解取极限得到一阶怪波解, 从而揭示怪波形成与速度跳跃之间的关联.
为了分析方程(1)对应常微分方程的平衡点及相轨迹, 对方程(1)做如下行波解约化:
$ \begin{split} & q(x,t) = \phi(\xi){\rm e}^{{\rm i}h(x,t)}, \\ & \xi = a(x-ct), \\ & h = Kx-\varOmega t, \end{split} $
其中a, c, K$ \varOmega $都是实数, $ \phi(\xi) $$ \xi $的实函数. 将其代入方程(1)中, 得到关于$ \phi $的常微分方程:
$ \begin{split} & {\rm i}[ -a(c-2K\alpha_1+4K^3\varepsilon^2\alpha_3)\phi'\\ &+2aK\varepsilon^2(\alpha_4-\alpha_5+\alpha_6)\phi^2\phi'+4a^3K\varepsilon^2\alpha_3\phi''']\\ & +\varepsilon^2\alpha_8\phi^5+[\alpha_2-K^2\varepsilon^2(\alpha_4+\alpha_5+\alpha_6-\alpha_7)]\phi^3\\ &+(\varOmega-K^2\alpha_1+K^4\varepsilon^2\alpha_3+a^2\varepsilon^2\alpha_6\phi'^2\\ & +a^2\varepsilon^2\alpha_7\phi'^2)\phi+a^2\varepsilon^2(\alpha_4+\alpha_5)\phi^2\phi''\\ &+a^2[\alpha_1\phi''+\varepsilon^2\alpha_3(-6K^2\phi''+a^2\phi'''')] = 0. \end{split} $
进一步, 分离方程(3)的实部和虚部, 得到
$ \begin{split} & -a(c-2K\alpha_1+4K^3\varepsilon^2\alpha_3)\phi'+2aK\varepsilon^2\\ & \times(\alpha_4-\alpha_5+\alpha_6)\phi^2\phi'+4a^3K\varepsilon^2\alpha_3\phi''' = 0, \end{split} $
$ \begin{split} & \varepsilon^2\alpha_8\phi^5 +[\alpha_2-K^2\varepsilon^2(\alpha_4+\alpha_5+\alpha_6-\alpha_7)]\phi^3\\ &+(\varOmega-K^2\alpha_1+K^4\varepsilon^2\alpha_3+a^2\varepsilon^2\alpha_6\phi'^2\\ & +a^2\varepsilon^2\alpha_7\phi'^2)\phi+a^2\varepsilon^2(\alpha_4+\alpha_5)\phi^2\phi''\\ &+a^2[\alpha_1\phi''+\varepsilon^2\alpha_3(-6K^2\phi''+a^2\phi'''')] = 0. \end{split} $
将(4)式关于$ \xi $积分一次并取积分常数为零, 可得
$ \begin{split} & 12a^{3}K\varepsilon^{2}\alpha_3\phi'' -a(3c-6K\alpha_1+12K^3\varepsilon^2\alpha_3)\phi \\ &+2aK\varepsilon^2(\alpha_4-\alpha_5+\alpha_6)\phi^{3} = 0. \end{split} $
而对(4)式关于$ \xi $求一阶导, 得到$ \phi'''' $如下:
$ \begin{split} \phi'''' = & \frac{1}{4a^2K\varepsilon^2\alpha_3}(-4K\varepsilon^2\alpha_4\phi\phi'^2+4K\varepsilon^2\alpha_5\phi\phi'^2\\ &-4K\varepsilon^2\alpha_6\phi\phi'^2 +c\phi''-2K\alpha_1\phi''\\ &+ 4K^3\varepsilon^2\alpha_3\phi''-2K\varepsilon^2\alpha_4\phi^2\phi''\\ &+2K\varepsilon^2\alpha_5\phi^2\phi'' -2K\varepsilon^2\alpha_6\phi^2\phi''). \end{split} $
结合(6)式和(7)式消掉方程(5)中的$ \phi'' $, $ \phi''' $$ \phi'''' $, 有
$ \begin{split} -\phi'^{2} &+\frac{c^{2}+4K^{2}[4\varepsilon^{2}(-cK+\varOmega)\alpha_3 -(\alpha_1-4K^{2}\varepsilon^{2}\alpha_3)^{2}]}{16a^{2}K^{2}\varepsilon^{4}\alpha_3(\alpha_4-\alpha_5-\alpha_7)}\\ & +\frac{c(\alpha_4+5\alpha_5-2\alpha_6)+2K\alpha_1(-2\alpha_4-4\alpha_5+\alpha_6)+4K\alpha_3[3\alpha_2+K^{2}\varepsilon^{2}(\alpha_4-\alpha_5-2\alpha_6+3\alpha_7)]}{12a^{2}K\varepsilon^{2}\alpha_3(\alpha_4-\alpha_5-\alpha_7)}\phi^{2} \\ & +\frac{(-\alpha_4-3\alpha_5+\alpha_6)(\alpha_4-\alpha_5+\alpha_6)+12\alpha_3\alpha_8}{12a^{2}\alpha_3(\alpha_4-\alpha_5-\alpha_7)}\phi^{4} = 0. \end{split}$
另一方面, 将方程(6)两边同乘以$ 2\phi' $, 得
$12a^{3}K\varepsilon^{2}\alpha_3\phi''\phi'-a(3c-6K\alpha_1+12K^3\varepsilon^2\alpha_3)\phi\phi'+2aK\varepsilon^2(\alpha_4-\alpha_5+\alpha_6)\phi^{3}\phi' = 0. $
同时, 对方程(8)关于$ \xi $求一阶导数, 有
$ \begin{split} &\frac{c(\alpha_4+5\alpha_5-2\alpha_6)+2K\alpha_1(-2\alpha_4-4\alpha_5+\alpha_6)+4K\alpha_3[3\alpha_2+K^{2}\varepsilon^{2}(\alpha_4-\alpha_5-2\alpha_6+3\alpha_7)]}{6a^{2}K\varepsilon^{2}\alpha_3(\alpha_4-\alpha_5-\alpha_7)}\phi\phi'\\ - & 2\phi'\phi''+\frac{(-\alpha_4-3\alpha_5+\alpha_6)(\alpha_4-\alpha_5+\alpha_6)+12\alpha_3\alpha_8}{3a^{2}\alpha_3(\alpha_4-\alpha_5-\alpha_7)}\phi^{3}\phi' = 0. \end{split} $
由于方程(9)和方程(10)含有相同项, 两者需相容, 则参数满足以下条件:
$ c = \frac{2K[6\alpha_2\alpha_3+\alpha_1(\alpha_4-7\alpha_5+\alpha_6-3\alpha_7)-4K^2\varepsilon^2\alpha_3(\alpha_4-\alpha_5+\alpha_6-3\alpha_7)]}{2\alpha_4-8\alpha_5+2\alpha_6-3\alpha_7}, $
$ \begin{split} & 4\alpha_5^2 +\alpha_6^2+\alpha_4(-4\alpha_5+\alpha_6-\alpha_7)-\alpha_6\alpha_7\\ &+\alpha_5(-5\alpha_6+\alpha_7)+12\alpha_3\alpha_8 = 0 . \end{split} $
因此, 在条件(11)式和(12)式下, 方程(9)和方程(10)可同时化为如下常微分方程:
$ \phi''+\beta_1\phi+\beta_2\phi^{3} = 0,$
其中
$ \begin{split} &\beta_1 = -\frac{4K^{2}\varepsilon^{2}\alpha_3-2\alpha_1+\frac{c}{K}}{4a^{2}\varepsilon^{2}\alpha_3},\\ &\beta_2 = \frac{\alpha_4 - \alpha_5 + \alpha_6}{6a^{2}\alpha_3}. \end{split} $
为了借助平面动力系统分岔理论讨论二阶常微分方程(13)的二维平衡点类型及相轨迹, 在变换$ X\equiv\phi $$ Y = \phi_{\xi} $下, 方程(13)可等价于如下二维平面动力系统:
$X' = Y, \tag{15a}$
$Y' = -\beta_1X-\beta_2X^{3} .\tag{15b}$
该系统是一个哈密尔顿系统, 具有如下哈密尔顿函数:
$ H(X,Y) = \frac{1}{2}Y^{2}+\frac{\beta_1}{2}X^{2}+\frac{\beta_2}{4}X^{4} = h^*. $
为了判断系统(15)的平衡点的类型, 我们得出其相应的雅克比矩阵为
${ J} = \left(\!\!\!{\begin{array}{*{20}{c}}0&1\\{ - {\beta _1} - 3{\beta _2}{X^2}}&0\end{array}}\!\!\!\right),$
该雅克比矩阵的特征值为$ \lambda_1, _2 = \pm\sqrt{-\beta_1-3\beta_2 X^{2}} $.
根据特征值与平衡点的对应关系, 我们对系统(15)的平衡点类型及解的稳定性分析如下:
1) 当$ \beta_1/\beta_2\geqslant 0 $时, 系统只有一个平衡点(0, 0)并且在该点处J的特征值为$ \pm\sqrt{-\beta_1} $. 如果$ \beta_1>0 $, 则(0, 0)点是中心点; 如果$ \beta_1<0 $, 则(0, 0)点是不稳定鞍点.
2) 当$ \beta_1/\beta_2<0 $时, 系统有三个平衡点: $ S_0(0, 0) $, $ S_1(\sqrt{-\beta_1/\beta_2}, 0) $$ S_2(-\sqrt{-\beta_1/\beta_2}, 0) $. 对应这些平衡点, J的特征值分别为$ \pm\sqrt{-\beta_1} $, $ \pm\sqrt{2\beta_1} $$ \pm\sqrt{2\beta_1} $. 若$ \beta_1<0 $, 则$ S_0 $是不稳定鞍点, 而$ S_1 $$ S_2 $是中心点, 此时存在围绕中心点$ S_1 $$ S_2 $绕至$ S_0 $的同宿轨道. 若$ \beta_1>0 $, 那么$ S_0 $是中心点, $ S_1 $$ S_2 $是不稳定的鞍点, 此时系统存在异宿轨道.
为了演示系统(15)存在的同宿轨道和异宿轨道, 选取特定参数画出系统(15)在条件$ \beta_1/\beta_2<0 $下的相位图, 如图1所示. 当取定$ \beta_1 = $ –1/10 和$ \beta_2 = 1/18 $, 图1(a)显示有一个鞍点和两个中心点, 随着时间的变化从鞍点出发的轨道最终会返回形成同宿轨道. 而当$ \beta_1 = 1 $$ \beta_2 = -5/9 $时, 图1(b)显示有一个中心点和两个鞍点, 并形成从一个鞍点到另一个鞍点的异宿轨道.
图 1 系统(15)的相位图 (a)同宿轨道(β1 = –1/10, β2 = 1/18); (b) 异宿轨道(β1 = 1, β2 = –5/9)
Figure1. Phase portraits of System (15): (a) Homoclinic orbits (β1 = –1/10, β2 = 1/18); (b) heteroclinic orbits (β1 = 1, β2 = –5/9).

进一步, 根据相位轨道与非线性偏微分方程的解之间的对应关系[26,27], 可以得出如下结论: 如果$ \beta_1/\beta_2<0 $$ \beta_1<0 $, 则方程(1)具有明孤立波解; 而当$ \beta_1/\beta_2<0 $$ \beta_1>0 $时, 方程(1)有冲击波解或暗孤立波解. 通过(14)式, 可以将条件$ \beta_1/\beta_2<0 $$ \beta_1<0 $具体表示为$ 6 K^2\varepsilon^2\alpha_3-3\alpha_1+ $ 3c/2K < 0, α5α4α6 > 0 且α3 < 0或者$ 6 K^2\varepsilon^2\alpha_3-3\alpha_1+$ 3c/2K < 0, α5α4α6 > 0且α3 < 0. 同样地, 条件$ \beta_1/\beta_2<0 $$ \beta_1>0 $可具体写成$ 6 K^2\varepsilon^2\alpha_3-3\alpha_1+3 c/2 K<0 $, α5α4α6 > 0且α3 < 0或者$ 6 K^2\varepsilon^2\alpha_3-3\alpha_1+ $ 3c/2K < 0, α5α4α6 > 0且α3 < 0.
下面为了验证方程(1)中明孤立波解和暗孤立波解的存在性, 分别在同宿轨和异宿轨条件下利用哈密尔顿函数对系统(15)进行求解. 令$ h^* = h_i = H(S_i)(i = 0, 1, 2) $, 则有
$ h_1 = h_2 = H(S_1) = \frac{-\beta_1^{2}}{4\beta_2}, \;\; h_0 = H(S_0) = 0\,. $
事实上, 由相容性条件(11)式和(12)式可知表达式(8)与哈密尔顿函数(16)式是等价的, 故可得
$ h^{*} = \frac{c^{2}\!+\!4K^{2}[4\varepsilon^{2}(-cK+\varOmega)\alpha_3\!- \! (\alpha_1\!-\! 4K^{2}\varepsilon^{2}\alpha_3)^{2}]} {16a^{2}K^{2}\varepsilon^{4}\alpha_3(\alpha_4-\alpha_5-\alpha_7)}\,. $
1) 当$ \beta_1<0 $$ \beta_2>0 $时, 则平衡点$ S_0 $是不稳定鞍点, 而平衡点$ S_1 $$ S_2 $是中心点, 此时存在围绕中心点$ S_1 $$ S_2 $$ S_0 $的同宿轨道. 当$ h^* = h_0 $时, 由表达式(16)得到
$ \begin{align} \phi^{'2} = -\beta_1\phi^{2}-\frac{\beta_2}{2}\phi^{4}\,. \end{align} $
求解得
$ \phi = \pm\sqrt{\frac{-2\beta_1}{\beta_2}}\rm{sech}(\sqrt{-\beta_1}\xi)\,. $
将解(21)式代入变换(2)式中得到方程(1)的明孤立波解为
$ q(x,t) = \pm {\rm e}^{{\rm i}(kx-\varOmega t)}\sqrt{\frac{-2\beta_1}{\beta_2}}{\rm sech}\sqrt{-\beta_1}a(x-ct)\,, $
其中
$ \begin{aligned} & c = \frac{{2K[6{\alpha _2}{\alpha _3} + {\alpha _1}({\alpha _4} - 7{\alpha _5} + {\alpha _6} - 3{\alpha _7}) - 4{K^2}{\varepsilon ^2}{\alpha _3}({\alpha _4} - {\alpha _5} + {\alpha _6} - 3{\alpha _7})]}}{{2{\alpha _4} - 8{\alpha _5} + 2{\alpha _6} - 3{\alpha _7}}} ,\\ & \varOmega = - \frac{{{c^2} - 4{K^2}[4cK{\varepsilon ^2}{\alpha _3} + {{({\alpha _1} - 4{K^2}{\varepsilon ^2}{\alpha _3})}^2}]}}{{16{K^2}{\varepsilon ^2}{\alpha _3}}} ,\\ & 4\alpha _5^2 + \alpha _6^2 + {\alpha _4}( - 4{\alpha _5} + {\alpha _6} - {\alpha _7}) - {\alpha _6}{\alpha _7} + {\alpha _5}( - 5{\alpha _6} + {\alpha _7}) + 12{\alpha _3}{\alpha _8} = 0 . \end{aligned} $
2) 当$ \beta_1>0 $$ \beta_2<0 $时, 那么平衡点$ S_0 $为中心点, 平衡点$ S_1 $$ S_2 $是不稳定的鞍点, 此时存在异宿轨道. 当$ h^* = h_1 = h_2 $时, 由表达式(16)得到
$ \begin{align} \phi^{'2} = -\beta_1\phi^{2}-\frac{\beta_2}{2}\phi^{4}-\frac{\beta_1^{2}}{2\beta_2}\,. \end{align} $
可求得
$ \phi = \pm\sqrt{\frac{-\beta_1}{\beta_2}}\tanh\left(\sqrt{\frac{\beta_1}{2}}\xi\right)\,. $
将解(24)式代入变换(2)式中得到方程(1)的暗孤立波解为
$ q(x,t) = \pm {\rm e}^{{\rm i}(kx-\varOmega t)}\sqrt{\frac{\beta_1}{-\beta_2}}\tanh\sqrt{\frac{\beta_1}{2}}a(x-ct)\,, $
其中
$ \begin{aligned} c =\, & \frac{{2K[6{\alpha _2}{\alpha _3} + {\alpha _1}({\alpha _4} - 7{\alpha _5} + {\alpha _6} - 3{\alpha _7}) - 4{K^2}{\varepsilon ^2}{\alpha _3}({\alpha _4} - {\alpha _5} + {\alpha _6} - 3{\alpha _7})]}}{{2{\alpha _4} - 8{\alpha _5} + 2{\alpha _6} - 3{\alpha _7}}} ,\\ \varOmega =\, &\, ({\alpha _4} \!-\! {\alpha _5} \!-\! {\alpha _7})\left[ { \!-\! \frac{{3{{(c \!-\! 2K{\alpha _1} \!+\! 4{K^3}{\varepsilon ^2}{\alpha _3})}^2}}}{{16{K^2}{\varepsilon ^2}{\alpha _3}({\alpha _4} \!-\! {\alpha _5} \!+\! {\alpha _6})}} \!+\! \frac{{cK}}{{{\alpha _4} \!-\! {\alpha _5} \!-\! {\alpha _7}}}} \right. \!-\! \frac{{{c^2}}}{{16{K^2}{\varepsilon ^2}{\alpha _3}({\alpha _4} \!-\! {\alpha _5} \!-\! {\alpha _7})}} \left. { \!+\! \frac{{{{({\alpha _1} \!-\! 4{K^2}{\varepsilon^2}{\alpha _3})}^2}}}{{4{\varepsilon^2}{\alpha _3}({\alpha _4} \!-\! {\alpha _5} \!-\! {\alpha _7})}}} \right] ,\\4\alpha _5^2=\, &\, -\alpha _6^2 - {\alpha _4}( - 4{\alpha _5} + {\alpha _6} - {\alpha _7}) + {\alpha _6}{\alpha _7} - {\alpha _5}( - 5{\alpha _6} + {\alpha _7}) - 12{\alpha _3}{\alpha _8}. \end{aligned} $
根据参数条件, 在解(22)式和(25)式中选取恰当的参数值得到明暗孤立波的传播图形, 如图2所示.
图 2 (a)由明孤立波解(22)式描述的明孤立波传输图形, 其中参数选取为 α1 = 1, α2 = 2, α3 = 1, α4 = 8, α5 = 2,α6 = 6, α7 = 4, α8 = 6, c = 1, K = 1, $\varOmega $ = 51/16, ε = 1, a = 1; (b) 由暗孤立波解(25)式描述的暗孤立波传输图形, 其中参数选取为α1 = –1, α2 = 2, α3 = 1, α4 = –8, α5 = –2, α6 = –6, α 7 = –4, α8 = 6, c = –7, K = 1, $\varOmega $ = –123/32, ε = 1, a = 1
Figure2. (a) Propagation of bright solitary wave via Solution (22) with the parameters chosen as α1 = 1, α2 = 2, α3 = 1, α4 = 8, α5 = 2, α6 = 6, α7 = 4, α8 = 6, c = 1, K = 1, $\varOmega $ = 51/16, ε = 1, a = 1; (b) propagation of dark solitary wave via Solution (25) with the parameters chosen as α1 = –1, α2 = 2, α3 = 1, α4 = –8, α5 = –2, α6 = –6, α7 = –4, α8 = 6, c = –7, K = 1, $\varOmega $ = –123/32, ε = 1, a = 1.

目前, 怪波解已经被发现存在于很多非线性发展方程中, 例如NLSE、Hirota方程、导数非线性薛定谔方程、Sasa-Satsuma方程、离散Ablowitz-Ladik方程以及变系数高阶非线性薛定谔方程等[35-37]. 本节将基于方程(1)在可积条件下求得的呼吸子解[30], 分析呼吸子的速度跳跃现象以及怪波产生的条件.
在可积条件下, 分别取方程(1)中的参数为α1 = 1, α2 = 2, α3 = 1, α4 = 8, α5 = 2, α6 = 6, α7 = 4, α8 = 6, 得到如下形式:
$ \begin{split} {\rm i}q_t & +q_{xx}+2q|q|^2+\varepsilon^2(q_{xxxx}+6{q^2_x}q^*+4q|q_x|^2\\ & +8|q|^2q_{xx}+2q^2{q^*_{xx}}+6|q|^4q) = 0 .\\[-10pt] \end{split} $
基于AKNS谱问题, 方程(26)的Lax对表示如下[30]:
${{ {\varPhi}} _x} = { M}{ {\varPhi}} ,\;{ M} = \left( {\begin{array}{*{20}{c}}{ - {\rm i}\lambda }&q\\{ - {q^*}}&{{\rm i}\lambda }\end{array}} \right),\tag{27a}$
${{ {\varPhi}} _t} = { N}{ {\varPhi}} ,\;{ N} = \left( {\begin{aligned}&{8{\rm i}{\varepsilon ^2}{A_4} - 2{\rm i}{\lambda ^2} + {\rm i}|q{|^2}}\\ &\quad\quad\quad\quad\quad {8{\rm i}{\varepsilon ^2}{B_4} + 2\lambda q + {\rm i}{q_x}}\\&{8{\rm i}{\varepsilon ^2}{C_4} - 2\lambda {q^*} + {\rm i}q_x^*}\\ &\quad\quad\quad { - 8{\rm i}{\varepsilon ^2}{A_4} + 2{\rm i}{\lambda ^2} - {\rm i}|q{|^2})}\end{aligned}} \right),\tag{27b}$
其中
$ \begin{split}{A_4} =\,& {\lambda ^4} - \frac{1}{2}|q{|^2}{\lambda ^2} + \frac{\rm i}{4}(qq_x^* - {q_x}{q^*})\lambda \\ & + \frac{1}{8}(3{q^2}{q^{*2}} + {q^*}{q_{xx}} + qq_{xx}^* - {q_x}q_x^*) ,\\{B_4} =\,& {\rm i}q{\lambda ^3} - \frac{1}{2}{q_x}{\lambda ^2} - \frac{\rm i}{4}({q_{x,x}} + 2{q^2}{q^*})\lambda \\ & + \frac{1}{8}({q_{xxx}} + 6q{q^*}{q_x}) ,\\{C_4} =\,& - {\rm i}{q^*}{\lambda ^3} - \frac{1}{2}q_x^*{\lambda ^2} + \frac{\rm i}{4}(q_{xx}^* + 2q{q^{*2}})\lambda \\ & + \frac{1}{8}(q_{xxx}^* + 6q{q^*}q_x^*) .\end{split}$
文献[30]基于Lax对(27)式给出了方程(26)的一阶达布变换, 其中势函数变换关系如下:
$ q^{[1]} = q^{[0]}-\dfrac{2i\Delta_1}{\varDelta_2}\,, $
这里$ \varDelta_1 \!=\! -2 if_{1, 1}f_{2, 1}\mathrm{Im}(\lambda_1) $, $ \varDelta_2 \!=\! f_{1, 1}f_{2, 2} \!-\! f_{1, 2}f_{2, 1},$ 其中$ f_{2, 1} = -f^{*}_{1, 2} $, $ f_{2, 2} = f^{*}_{1, 1} $, $ (f_{i, 1}, f_{i, 2}) $ $ (i = 1, 2) $是Lax对(27)式相应$ \lambda = \lambda_i $($ \lambda_2 = \lambda^{*}_1 $)的两组线性无关解.
以方程(26)的平面波解作为种子解:
$ \begin{align} q^{[0]} = c{\rm e}^{{\rm i}\rho}, \quad \rho = ax+bt, \end{align} $
其中$ b = \varepsilon^2(a^4-12 a^2 c^2+6 c^4)+2 c^2-a^2 $, $ a, \, b, \, c $分别代表波的波数、频率和振幅. 此时, 与$ q^{[0]} $对应的Lax对(27)式的解为
$ f_{1,1} = k_1c{\rm e}^{{\rm i}[\frac{\rho}{2}+\kappa(x,t)]}+{\rm i}k_2\left(\frac{a}{2}+h+\lambda_1\right){\rm e}^{{\rm i}[\frac{\rho}{2}-\kappa(x,t)]},\tag{31a}$
$ f_{1,2} \!=\! k_2c{\rm e}^{-{\rm i}[\frac{\rho}{2}+\kappa(x,t)]}\!+\!{\rm i}k_1\left(\frac{a}{2}\!+\!h\!+\!\lambda_1\right){\rm e}^{-{\rm i}[\frac{\rho}{2}-\kappa(x,t)]},\tag{31b}$
其中
$ \begin{split} \kappa =\, & h(x + dt),\\ {k_1} =\, & {{\rm e}^{{\rm i}h({s_1}\varepsilon + {s_2}{\varepsilon ^2})}},\\ {k_2} =\, & {{\rm e}^{ - {\rm i}h({s_1}\varepsilon + {s_2}{\varepsilon^2})}},\\ h =\, & \sqrt {{c^2} + {{\left( {{\lambda _1} + \frac{a}{2}} \right)}^2}} \equiv {h_R} + {\rm i} {h_{\rm I}} ,\\ d =\, & 2{\lambda _1} - a + {\varepsilon ^2}[a({a^2} - 6{c^2}) - 8\lambda _1^3 \\ & + 4a\lambda _1^2 + (4{c^2} - 2{a^2}){\lambda _1}] \equiv {d_R} + {\rm i} {h_{\rm I}}. \end{split} $
将平面波解(30)式及特征函数(31a)式和(31b)式代入变换(29)式中, 得到方程(26)的一阶呼吸子解为[30]:
$ q^{[1]} = c {\rm e}^{{\rm i}\rho} +\frac{2\eta\{\omega_1\cos(2G)-\omega_2\cosh(2F) - {\rm i}[(\omega_1-2c^2)\sin(2G) -\omega_3\sinh(2F)]\}{\rm e}^{{\rm i}\rho}}{\omega_1\cosh(2F)-\omega_2\cos(2G)}, $
其中
$ \begin{split} & \omega_1 = c^2+(h_{\rm I}+\eta)^2+\left(\xi+h_{\rm R}+\frac{a}{s}\right)^2,\\ & \omega_2 = 2c\left(h_{\rm I}+\eta\right),\\ & \omega_3 = 2c\left(\xi+h_{\rm R}+\frac{a}{2}\right),\\ & F = h_{\rm I} x +(d_{\rm R}h_{\rm I}+d_{\rm I}h_{\rm R})t\equiv k_1 x+w_1 t,\\ & G = h_{\rm R} x +(d_{\rm R}h_{\rm R}-d_{\rm I}h_{\rm I})t\equiv k_2 x+w_2 t, \\ & \lambda_1 = \xi+{\rm i}\eta, \end{split}$
这里下标R和I分别代表相关参数的实部和虚部. 图3展示了由呼吸子解(32)式描述的单个呼吸子的动力学演化, 其在演化过程中呈现出随时间的周期性震荡现象.
图 3 解(32)式描述的一阶呼吸子的动力学演化, 其中参数选取为$\xi=0$, $\eta={1}/{2}$, $c={2}/{5}$$a=0$
Figure3. The propagation of one breather via Solution (32) with the parameters chosen as $\xi=0$, $\eta={1}/{2}$, $c={2}/{5}$ and $a=0$.

接下来, 将基于呼吸子解的精确表达式(32)来分析呼吸子的速度跳跃现象. 根据非线性波速度的相关理论[38,39], 从表达式(32)中可以看出, 呼吸子沿时间和空间轴方向发生周期性震荡, 其震荡性体现在三角函数部分$ {\rm{cos}}(2 G) $$ {\rm{sin}}(2 G) $, 且 G中的$ k_2 $$ w_2 $分别代表震荡的空间和时间频率, 震荡速度对应于呼吸子相速度$ V_{\rm p} $. 另外, 呼吸子可以在$ (x, t) $平面上沿直线传播, 其传播轨迹由$ F = 0 $决定, 传播速度对应于呼吸子的群速度$ V_{\rm g} $. 因此, 得到呼吸子解(32)式的群速度$ V_{\rm g} $和相速度$ V_{\rm p} $的表达式如下:
$\begin{split}{V_{\rm g}} =\,& - \frac{{{w_1}}}{{{k_1}}} = - \frac{{{d_{\rm R}}{h_{\rm I}} + {d_{\rm I}}{h_{\rm R}}}}{{{h_{\rm I}}}} =\frac{{ - 8{\eta ^2}(a + 2\xi )(1 - {a^2} + 2{c^2} + 4{\eta ^2} + 4a\xi - 12{\xi ^2})}}{{{a^2} + 4{c^2} - 4{\eta ^2} + 4a\xi + 4{\xi ^2} + \sqrt {{D_1}} }}\\&+ \frac{{[{a^3} - 2{a^2}\xi - a(1 + 6{c^2} + 4{\eta ^2} - 4{\xi ^2}) + 2\xi (1 + 2{c^2} + 12{\eta ^2} - 4{\xi ^2})]({a^2} + 4{c^2} - 4{\eta ^2} + 4a\xi + 4{\xi ^2} + \sqrt {{D_1}} )}}{{{a^2} + 4{c^2} - 4{\eta ^2} + 4a\xi + 4{\xi ^2} + \sqrt {{D_1}} }},\end{split}$
$\begin{split}{V_{\rm p}} =\, & - \frac{{{w_2}}}{{{k_2}}} = - \frac{{{d_{\rm R}}{h_{\rm R}} - {d_{\rm I}}{h_{\rm I}}}}{{{h_{\rm R}}}}= 4\xi ({c^2} + 6{\eta ^2} - 2{\xi ^2}) - [{a^3} + 2\xi - 2{a^2}\xi - a(1 + 6{c^2} + 4{\eta ^2} - 4{\xi ^2})]\\ & + \frac{{( - 1 + {a^2} - 2{c^2} - 4{\eta ^2} - 4a\xi + 12{\xi ^2})(4{c^2} - 4{\eta ^2} + {{(a + 2\xi )}^2}) + \sqrt {{D_2}} }}{{a + 2\xi }},\end{split}$
其中
$\begin{aligned}& {D_1} = 16{\eta ^2}{(a + 2\xi )^2} + {[{a^2} + 4a\xi + 4({c^2} - {\eta ^2} + {\xi ^2})]^2}{\mkern 1mu} ,\\& {D_2} = (4{(c - \eta )^2} + {(a + 2\xi )^2})(4{(c + \eta )^2} + {(a + 2\xi )^2}){\mkern 1mu} .\end{aligned}$
$ \eta\to\pm c $时, 群速度$ V_{\rm g} $和相速度$ V_{\rm p} $化简如下:
$\begin{split}{V_{\rm g}} =\,& \frac{{ - 8{\eta ^2}(a + 2\xi )(1 - {a^2} + 6{c^2} + 4a\xi - 12{\xi ^2})}}{{{a^2} + 4a\xi + 4{\xi ^2} + \sqrt {{{(a + 2\xi )}^2}({a^2} + 16{c^2} + 4a\xi + 4{\xi ^2})} }}\\& + [{a^3} - 2{a^2}\xi - a(1 + 10{c^2} - 4{\xi ^2}) + 2\xi (1 + 14{c^2} - 4{\xi ^2})],\end{split}$
$\begin{split}{V_{\rm p}} =\, & 4\xi (7{c^2} - 2{\xi ^2}) - [{a^3} + 2\xi - 2{a^2}\xi - a(1 + 10{c^2} - 4{\xi ^2})]\\& + \frac{{( - 1 + {a^2} - 6{c^2} - 4a\xi + 12{\xi ^2})({{(a + 2\xi )}^2} + \sqrt {{{(a + 2\xi )}^2}[16{c^2} + {{(a + 2\xi )}^2}]} )}}{{a + 2\xi }}.\end{split}$
从表达式(35)式和(36)式中可以看出, 当$ \xi $沿不同方向趋于$ -a/2 $时, $ \sqrt{(a+2\xi)^2} $的符号是不同的, 从而导致群速度$ V_{\rm g} $和相速度$ V_{\rm p} $趋于不同的值, 即产生跳跃现象. 为了图形演示速度关于参数a的跳跃现象, 我们在条件$ \eta\to\pm c $下, 保留参数a且将其他参数取定, 即$ \varepsilon^2 = 1 $, $ \eta = c $, $ c = 1 $$ \xi = 1 $, 得到,
$ \begin{align} V_{\rm g} = \frac{\delta_1}{\delta_2}, \quad V_{\rm p} = \frac{\delta_3}{\delta_4}, \end{align} $
其中
$\begin{split}{\delta _1} =\, & 2{a^4} + {a^5} + 2{a^2}[\sqrt {{{(2 + a)}^2}(20 + 4a + {a^2})} - 15]\\& - {a^3}[3 + \sqrt {{{(2 + a)}^2}(20 + 4a + {a^2})} ] \\&- a\left\{ {48 + 3[4 + \sqrt {{{(2 + a)}^2}(20 + 4a + {a^2})} ]} \right.\\& \left. {- 2[48 + 5\sqrt {{{(2 + a)}^2}(20 + 4a + {a^2})} ]} \right\}\\& - 2\left\{ {48 - 3[\sqrt {{{(2 + a)}^2}(20 + 4a + {a^2})} - 4]} \right.\\& \left. {+ 2[7\sqrt {{{(2 + a)}^2}(20 + 4a + {a^2})} - 72]} \right\}, \\[-15pt]\end{split}$
$\begin{array}{*{20}{l}}{{\delta _2} = \sqrt {{{(2 + a)}^2}(20 + 4a + {a^2})} - 4 - 4a - {a^2},}\end{array}$
$\begin{split}{\delta _3} =\,& {a^4} - 80 - 7\sqrt {{{(2 + a)}^2}(20 + 4a + {a^2})}\\& + {a^2}[\sqrt {{{(2 + a)}^2}(20 + 4a + {a^2})} - 15]\\& + 4[37 + 3\sqrt {{{(2 + a)}^2}(20 + 4a + {a^2})} ] \\&- 4a[\sqrt {{{(2 + a)}^2}(20 + 4a + {a^2})} - 3],\end{split}$
$ \begin{align} \delta_4 = -2(a+2). \end{align} $
图4给出了群速度$ V_{\rm g} $和相速度$ V_{\rm p} $随参数a的变化关系. 以$ V_{\rm g} $为例, 由图4可清晰地看到群速度在$ a = -2 $处出现突然的跳跃. 为了更加说明该跳跃点, 我们计算了群速度在该点的左右极限, 发现$ \lim\limits_{a \to -2^-}V_{\rm g} = -54 \neq \lim\limits_{a \to -2^+}V_{\rm g} = 14 $, 从而表明群速度在$ a = -2 $处不连续. 类似地, 相速度$ V_{\rm p} $$ a = -2 $处的左右极限分别为14和–54, 也是不连续的, 即$ a = -2 $为跳跃点. 通过对表达式(37)的分析, 速度跳跃点的出现是由于分母中包含$ \sqrt{(2+a)^2} $, 从而导致$ V_{\rm g} $$ V_{\rm P} $$ a = -2 $点处不连续. 根据文献[24]中关于临界频率的定义, 由于速度在波数$ a = -2 $处表现出了速度的不连续性, 因此该不连续点可以定义为方程(26)的临界波数.
图 4 呼吸子的群速度$V_{\rm g}$(红实线)和相速度$V_{\rm p}$(蓝虚线)随参数a的变化关系
Figure4. Group velocity $V_{\rm g}$ (red-solid line) and phase velocity $V_{\rm p}$ (blue-dot line) of the breather

已有研究表明呼吸子在速度跳跃点会转变为怪波[24,25]. 为了验证方程(26)中怪波的形成与速度跳跃的关系, 我们将呼吸子解(32)式在$ \xi\to-{a}/{2} $$ \eta\to c $条件下取极限, 即 $ \lambda_1 = \xi+{\rm i} \eta$ $ \to \lambda_0 = $a/2 + ic, 得到方程(26)的一阶怪波解如下:
$ q_{\rm{rw}} = \left[c-\frac{4c\theta_1(\lambda_0-{\lambda_0}^2)(\theta_4+\alpha^*_1t)} {4c^2\theta_3\theta_4+\theta_1\theta_2+4c^2[\theta_4\alpha_1+\alpha^*_1(\theta_3+\alpha_1t)]t}\right] {\rm e}^{{\rm i}\rho(x,t)}, $
其中
$\begin{split} & \rho = ax+bt\,, b = (a^4-12a^2c^2+6c^4)\varepsilon^2+2c^2-a^2\,,\\ & \theta_1 = \dfrac{4c^2\alpha_1t}{a+2\lambda_0}+\dfrac{(ax+2\lambda_0 x+2i)\alpha_1+ \alpha_2 [a^4\varepsilon^2-a^2(1+6c^2\varepsilon^2)-8ac^2\varepsilon^2\lambda_0+4{\lambda_0}^2(1+2c^2\varepsilon^2-4\varepsilon^2{\lambda^2_0})]}{a-a^3\varepsilon^2+6ac^2\varepsilon^2-2\lambda_0 +2a^2\varepsilon^2\lambda_0-4c^2\varepsilon^2\lambda_0-4a\varepsilon^2{\lambda^2_0}+8\varepsilon^2{\lambda^3_0}}\,,\\ & \theta_2 = \dfrac{4c^2\alpha_1^*t}{a+2{\lambda^*_0}}+\frac{\alpha_1^*(ax+2x{\lambda_0}^*-2i)+\alpha_2^*[a^2(a^2\varepsilon^2-6c^2\varepsilon^2-1) -8ac^2\varepsilon^2{\lambda^*_0}+{\lambda^{*2}_0}(4+8c^2\varepsilon^2)-16\varepsilon^2{\lambda^{*4}_0}]}{a-a^3\varepsilon^2 +6ac^2\varepsilon^2+2{\lambda^*_0}(a^2\varepsilon^2-2c^2\varepsilon^2-1)-4a\varepsilon^2{\lambda^{*2}_0}+8\varepsilon^2{\lambda^{*3}_0}}\,,\\ & \theta_3 = \dfrac{\alpha_1x +\alpha_2[a^3\varepsilon^2-2a^2\varepsilon^2{\lambda_0} + a(4\varepsilon^2{\lambda^2_0}-1-6c^2\varepsilon^2) +2({\lambda_0}+2c^2\varepsilon^2{\lambda_0}-4\varepsilon^2{\lambda^3_0})]} {a^3\varepsilon^2-2a^2\varepsilon^2{\lambda_0}+a(4\varepsilon^2{\lambda^2_0}-1-6c^2\varepsilon^2) +2({\lambda_0}+2c^2\varepsilon^2{\lambda_0}-4\varepsilon^2{\lambda^3_0})}\,,\\ & \theta_4 = \dfrac{\alpha^*_1 x + \alpha^*_2[a(a^2\varepsilon^2-6c^2\varepsilon^2-1) + {\lambda^*_0}(2-2a^2\varepsilon^2+4c^2\varepsilon^2) + 4a\varepsilon^2{\lambda^{*2}_0} -8\varepsilon^2{\lambda^{*3}_0}]}{a(a^2\varepsilon^2-6c^2\varepsilon^2-1) + {\lambda^*_0}(2-2a^2\varepsilon^2+4c^2\varepsilon^2) + 4a\varepsilon^2{\lambda^{*2}_0}-8\varepsilon^2{\lambda^{*3}_0}}. \end{split} $
这里$ \alpha_1 $$ \alpha_2 $是两个非零的实参数. 图5给出了由解(42)式描绘的一阶怪波的动力学演化, 其中参数选取为$ c = 1 $, $ \eta = 1 $, $ \varepsilon^2 = 1 $$ \alpha_1 = \alpha_2 = 1 $. 通过极限计算, 当$ x\to \infty $$ t\to \infty $时, 一阶怪波$ |q_{\rm{rw}}|^2 $在无穷远处高度为1, 从而说明了怪波的空间和时间局域性. 进一步, 借助数值近似计算可以验证$ |q_{\rm{rw}}|^2 $分别在点(0.3660257, –1)和(–1.3660251, –1)附近取得最小值$ 5.32932\times10^{-13} $$ 5.3456\times10^{-13} $, 在点(–0.4999952, –1)附近取得最大值9, 该值即为一阶怪波的振幅. 因此, 解(42)式符合一阶怪波解具有一个极大值点和两个极小值点的特点, 且最高点振幅为背景波三倍以上. 以上分析表明方程(26)的呼吸子解(42)式在速度跳跃条件$ \xi\to-{a}/{2} $$ \eta\to c $下转化成了怪波, 进而证明速度的不连续性可以产生怪波. 值得注意的是, 该速度跳跃条件与达布变换方法求怪波解[24]时的谱参数取值相一致, 而该谱参数条件对应方程的调制不稳定区域[40], 因此调制不稳定性与速度跳跃从不同角度揭示了怪波产生的条件.
图 5 解(42)式描述的一阶怪波的动力学演化, 其中参数选取为$\xi=1$, $\eta=1$, $c=1$, $\alpha_1=1$, $\alpha_2=1$$a=-2$
Figure5. The propagation of first-order rogue wave via Solution (42) with the parameters chosen as $\xi=1$, $\eta=1$, $c=1$, $\alpha_1=1$, $\alpha_2=1$ and $a=-2$.

本文研究了一个在光纤光学和磁力学提出的四阶色散非线性薛定谔方程(1), 分别从动力学分析和速度不连续性两方面探讨了孤立波和怪波的产生条件.
首先, 利用行波解约化将该四阶色散非线性薛定谔方程变为常微分方程, 然后通过分离常微分方程的实虚部得到二维平面自治系统, 即系统(15a)和(15b). 进而, 借助相平面分析方法讨论了系统(15a)和(15b)的平衡点类型及解的稳定性, 同时根据平衡点类型与相轨道的关系发现二维平面自治系统具有同宿轨道和异宿轨道, 并在相应条件下求解得到了方程(1)的明孤立波解和暗孤立波解, 从而验证了同异宿轨道分别对应非线性方程的明暗孤立波解.
其次, 为了从速度不连续性角度研究怪波的产生条件, 我们基于非零背景平面之上的精确一阶呼吸子解, 推导出了呼吸子的群速度和相速度的显式表达式. 通过对速度表达式的理论分析发现, 当$ \xi\to-{a}/{2} $$ \eta\to\pm c $时, 群速度和相速度均关于波数存在不连续点, 即跳跃点或临界波数(见图3所示). 同时, 从数学上指出了该跳跃点的出现是由于速度表达式中含有因子$ \sqrt{(a+2\xi)^2} $. 为了验证呼吸子在速度跳跃点可以转化成怪波, 对一阶呼吸子解在速度不连续点处取极限获得了方程(26)的一阶怪波解, 从而表明怪波的形成与呼吸子速度的不连续性有关. 另外, 如果取定速度表达式(35)式和(36)式中其他参数而只保留振幅参数c, 则群速度和相速度均不会出现不连续点, 如图6所示. 在图6中, 参数选取为$ \xi = {1}/{2}(-a+\sqrt{-4 c^2}) $, $ \eta = 0 $$ a = 1 $, 此时群速度和相速度都是关于振幅参数c的连续函数.
图 6 (a) 群速度$V_{\rm g}$随振幅参数c的变化(红实线)和(b)相速度$V_{\rm p}$随振幅参数c的变化(蓝点线), 其中参数选取为$\xi={1}/{2}(\sqrt{-4 c^2}-a)$, $\eta=0$$a=1$
Figure6. Variation of the group velocity $V_{\rm g}$ (red-solid line) and phase velocity $V_{\rm p}$ (blue-dot line) about the amplitude parameter c with the parameters chosen as $\xi={1}/{2}(\sqrt{-4 c^2}-a)$, $\eta=0$ and $a=1$.

感谢深圳大学高等研究院的贺劲松教授在中国石油大学(北京)所作关于怪波形成机制的学术报告.
相关话题/系统 文献 光学 传播 方程

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 可积谐振系统中的极端波事件研究进展
    摘要:从微观角度上讲,单个极端异常波事件可视为可积模型方程的时空局域有理函数解.本文主要讨论了三类典型的可积谐振相互作用模型(即长波短波谐振方程,三波谐振相互作用方程,非线性薛定谔和麦克斯韦-布洛赫方程)的基阶Peregrine异常波解及其相关研究进展;明确指出了这些基阶异常波解形式具有普适性,可推 ...
    本站小编 Free考研考试 2021-12-29
  • 可积系统多孤子解的全反演对称表达式
    摘要:多孤子解是非线性数学物理系统的基本激发模式.文献中存在各种类型的表达式,如广田(Hirota)形式,朗斯基(Wronskian)或双朗斯基形式和法夫(Phaffian)形式.最近在多地系统的研究中,我们发现使用一种全新但等价的形式具有极为简洁和方便的优点.本文主要综述多种类型可积非线性系统的多 ...
    本站小编 Free考研考试 2021-12-29
  • 非线性光学中的暗孤子分子
    摘要:孤子分子是当前非线性光学中的重要课题.本文首先研究具有高阶色散和高阶非线性效应非线性光学模型中各种周期波(孤子晶格)的严格解,及各种可能的单孤子解.然后在一个可积的情况下,利用推广的双线性形式,给出多孤子解,并从多孤子解的速度共振条件给出暗孤子分子的严格解析表达式.对于本文给出模型的多暗孤子分 ...
    本站小编 Free考研考试 2021-12-29
  • 基于Pancharatnam-Berry相位超表面的二维光学边缘检测
    摘要:提出并设计一种基于Pancharatnam-Berry(P-B)相位超表面的二维光学微分器,并实现对光学图像的二维光学边缘检测.在环形光栅相位的作用下,该P-B相位超表面可将光束的左右旋分量在径向进行分离,在滤除中间重叠部分的线偏振光后,保留下来的光学信息即为二维光学微分结果.同时,通过调节该 ...
    本站小编 Free考研考试 2021-12-29
  • 准Λ型四能级系统选择反射光谱
    摘要:运用密度矩阵理论研究了气固界面准$\Lambda$型四能级原子系统的非线性选择反射光谱.基于刘维尔方程给出了一阶近似条件下探测光场的解析式.在探测场为弱场时,分析了信号场拉比频率、失谐量和耦合场失谐量对反射光谱线型的影响.数值模拟表明:信号场参与产生的选择反射峰线宽可以利用信号场拉比频率进行调 ...
    本站小编 Free考研考试 2021-12-29
  • 隔板对流系统的热流特性及热量输入与传递特性
    摘要:采用DNS方法对隔板对流装置进行模拟计算,研究系统中热流特性以及热量输入与传递特性.讨论了热流的纵向和横向输运特性,在此基础上对传热通道和狭缝区域的热通量以及对应底板外界输入热通量进行了定量化分析.研究结果表明,通道中低温流体向下冲击底板而后转入水平运动,流入狭缝区域并不断被加热,只进行水平的 ...
    本站小编 Free考研考试 2021-12-29
  • 基于聚焦离子束纳米剪纸/折纸形变的三维微纳制造技术及其光学应用
    摘要:高精度的三维微纳制造技术是现代光电子学和微纳光子学发展的重要基础之一,是实现下一代微纳光子集成器件的重要前提.纳米尺度的剪纸和折纸技术由于能够实现丰富的三维形变,正发展成为一门新兴的研究领域.本文系统地介绍了一种新型的片上三维微纳加工方法—基于聚焦离子束的纳米剪纸/折纸技术.该技术利用聚焦离子 ...
    本站小编 Free考研考试 2021-12-29
  • 基于空间角度复用和双随机相位的多图像光学加密方法
    摘要:提出了基于空间角度复用和双随机相位的多图像光学加密新方法.加密过程中,首先将原始图像进行随机相位调制和不同距离的菲涅耳衍射;其次,将携带调制后图像的参考光与携带随机相位且具有不同立体角的参考光相干叠加,产生干涉条纹;最后,将不同方向的干涉条纹叠加形成复合加密图像.解密为加密的逆过程,将复合加密 ...
    本站小编 Free考研考试 2021-12-29
  • Shimizu-Morioka系统与Finance系统生成Lorenz混沌的微分几何策略
    摘要:从一种受控混沌系统生成另一混沌系统可增强保密通信的安全性,具备潜在应用前景.研究了如何通过状态变换以及单输入反馈,驱使受控Shimizu-Morioka系统与受控Finance系统生成Lorenz混沌动态.主要方法是运用微分几何理论,将上述三种系统等价转换为下三角形式,并尽量简化和一致化其方程 ...
    本站小编 Free考研考试 2021-12-29
  • 含石墨烯分界面有耗分层介质的传播矩阵
    摘要:给出一种适用于含导电界面的有耗分层介质的传播矩阵方法.利用相位匹配原理给出斜入射时有耗介质波矢量的实部和虚部,二者方向不同使得在介质中传播非均匀平面波.根据边界条件,推导了跨越石墨烯界面的传播矩阵,以及“无限薄”石墨烯层的反透射系数解析式.最终将传播矩阵方法推广应用于含石墨烯界面的有耗分层介质 ...
    本站小编 Free考研考试 2021-12-29