姜冬竹，扎其勞，2，3

（1.內(nèi)蒙古師范大學(xué) 數(shù)學(xué)科學(xué)學(xué)院，呼和浩特 010022；2.內(nèi)蒙古自治區(qū)應(yīng)用數(shù)學(xué)中心，呼和浩特 010022；3.無(wú)窮維哈密頓系統(tǒng)及其算法應(yīng)用教育部重點(diǎn)實(shí)驗(yàn)室，呼和浩特 010022）

這里Φ=(?1，?2)T。近年來(lái)，對(duì)非線性演化方程周期怪波解的研究得到了學(xué)者們?cè)絹?lái)越多的關(guān)注。在雅可比橢圓函數(shù)［2］的背景下，利用雅可比橢圓函數(shù)展開(kāi)法，結(jié)合達(dá)布變換［3］與譜問(wèn)題的非線性化［4］的方法，許多非線性演化方程的周期怪波解被構(gòu)造。如高階修正Korteweg-de Vries （mKdV）方程［5］、Hirota 方程［6］、NLS 方程［7-8］。因此，本文重點(diǎn)是分析方程（3）的調(diào)制不穩(wěn)定性。在此基礎(chǔ)上，構(gòu)造其在雅可比橢圓函數(shù)dn 和cn背景上的兩類(lèi)不規(guī)則的周期怪波解。

1 調(diào)制不穩(wěn)定性分析

首先給出方程（3）的平面波種子解

這里A為實(shí)常數(shù)，ω表示背景頻率，κ表示波數(shù)。將式（7）代入方程（3）中，可得

在式（7）中加入一個(gè)小的擾動(dòng)參數(shù)ε，則擾動(dòng)解可以表達(dá)成如下形式

這里r1和r2都是小參數(shù)。將式（9）代入方程（3）中，可以得到關(guān)于r1和r2的齊次線性方程組。如果令r1和r2的系數(shù)行列式結(jié)果為零，并求解關(guān)于Ω的方程組，則可求得MI 增益G為

其中

由式（10）可得，在|K|<2A的條件下G會(huì)存在虛部，則擾動(dòng)函數(shù)r（x，t）呈指數(shù)增長(zhǎng)，從而破壞了系統(tǒng)的穩(wěn)定性并產(chǎn)生怪波。MI 波的寬度與參數(shù)A呈正相關(guān)，如圖1 所示。

圖1 方程（3）在α=-0.1，β=0.1，γ=0.1，δ=0.2 時(shí)的調(diào)制不穩(wěn)定性分析Fig.1 The MI analysis of equation （3） with α=-0.1，β=0.1，γ=0.1，δ=0.2

2 雅可比橢圓函數(shù)解

3 dn 周期波和cn 周期波的特征值

4 周期特征函數(shù)與非周期解

4.1 Lax 對(duì)的周期特征函數(shù)

通過(guò)計(jì)算可得

根據(jù)式（18）、（33）、（29）和（36），得

式（38）取正或負(fù)的平方根，則以下關(guān)系可以被推出

根據(jù)式（29）和（31），以下關(guān)系可以被推導(dǎo)出

其滿足

4.2 Lax 對(duì)的非周期解

假設(shè)Φ=(?1，?2)T是Lax 對(duì)式（4）-（5）在λ=λ1時(shí)的周期解，為推導(dǎo)怪波解，需要通過(guò)(ψ1，ψ2)T構(gòu)造非周期解，其表示為L(zhǎng)ax 對(duì)式（4）—（5）的二階線性無(wú)關(guān)解，其形式為

其中θ=θ(x，t)是待定函數(shù)，將式（42）代入式（4），可得

根據(jù)式（12）和式（31），可以將式（43）改寫(xiě)為

根據(jù)式（34）和式（35），式（44）可以退化為

對(duì)式（44）關(guān)于x進(jìn)行積分，得

其中θ0(t)是與t相關(guān)的積分常數(shù)。將式（42）代入式（5），利用雅可比橢圓函數(shù)展開(kāi)法，得

其中

通過(guò)計(jì)算可得Δ1=0，因此

這里η是積分常數(shù)。將式（49）代入式（46）得到

5 周期背景上的怪波解

根據(jù)文獻(xiàn)［3］，方程（3）的一重達(dá)布變換如下

其中(?11，?21)T為L(zhǎng)ax 對(duì)式（4）—（5）在λ=λ1時(shí)的非零解?；谑剑?1）和3.1 節(jié)，應(yīng)用雅可比橢圓函數(shù)dn和cn，方程（3）在周期背景上的怪波解如下

圖2 當(dāng)取 α=1， β=3，γ=3，δ=3，η=0 時(shí)方程（3）在dn 周期背景上的怪波解Fig.2 Rogue wave solution on the dn-periodic background of equation （3） with α=1， β=3，γ=3，δ=3，η=0

分別在k=0.5 和k=0.99 時(shí)描述方程（3）在dn 周期背景上的怪波解。特別地，當(dāng)k=0 或k=1 時(shí)，可以得到dn 周期背景上的退化怪波解（見(jiàn)圖2―5）。分別在k=0.5 和k=0.99 時(shí)描述方程（3）在cn 周期背景上的怪波解。特別地，當(dāng)k=1 時(shí)，可以得到cn 周期背景上的退化怪波解（見(jiàn)圖5―7）。

圖3 當(dāng)取t =0 時(shí)方程（3）在dn 周期背景上沿x 軸傳播的怪波Fig.3 The wave propagation along the x-axis of rogue wave solution on the dn-periodic background of equation （3） at t =0

圖4 當(dāng)取 α=1， β=3，γ=3，δ=3，η=0 時(shí)方程（3）在dn 周期背景上的退化怪波解Fig.4 Degenerate rogue wave solution on the dn-periodic background of equation （3） with α=1， β=3，γ=3，δ=3，η=0

圖5 當(dāng)取t=0 時(shí)方程（3）在dn 周期背景上沿x 軸傳播的退化怪波Fig.5 The wave propagation along the x-axis of degenerate rogue wave solution on the dn-periodic background of equation （3） at t =0

圖6 當(dāng)取α=1， β=3，γ=3，δ=3，η=0 時(shí)方程（3）在cn 周期背景上的怪波解Fig.6 Rogue wave solution on the cn-periodic background of equation （3） with α=1， β=3，γ=3，δ=3，η=0

圖7 當(dāng)取t =0 時(shí)方程（3）在cn 周期背景上沿x 軸傳播的怪波Fig.7 The wave propagation along the x-axis of degenerate rogue wave solution on the cn-periodic background of equation （3） at t =0

6 結(jié)論

本文在分析調(diào)制不穩(wěn)定性的基礎(chǔ)上，通過(guò)引入巴格曼約束，利用譜問(wèn)題的非線性特性和達(dá)布變換方法，構(gòu)造七階非線性薛定諤方程在dn 和cn 背景上的不規(guī)則周期怪波，并通過(guò)圖像顯示了它們的動(dòng)力學(xué)特性。此外，在雅可比橢圓函數(shù)的背景上，雅可比橢圓函數(shù)展開(kāi)法簡(jiǎn)化了推導(dǎo)過(guò)程，顯示了該方法在計(jì)算過(guò)程中的便捷性。