郭婷婷

(山西工程科技職業(yè)大學(xué) 數(shù)學(xué)教學(xué)研究部,山西 太原 030619)

0 引言

繼相對(duì)論和量子力學(xué)之后,對(duì)非線性科學(xué)的研究在自然科學(xué)中異軍突起。而孤立子理論作為該領(lǐng)域內(nèi)的一個(gè)重要分支,為數(shù)學(xué)物理及相關(guān)學(xué)科開辟出一片新天地,例如非線性光學(xué)中涉及的空間光孤子,電荷密度波,生物科學(xué)中研究的達(dá)維多夫孤子等。數(shù)學(xué)的啟發(fā)性與相關(guān)學(xué)科的應(yīng)用性相得益彰,孤子理論研究領(lǐng)域生機(jī)盎然。

孤子是一種非線性現(xiàn)象,為合理描述其不斷演化的力學(xué)性態(tài),對(duì)非線性孤波模型的求解顯得尤為重要。隨著孤子問題的深入研究,許多有效的求解方法應(yīng)運(yùn)而生,例如:Lie代數(shù)方法[1-2]、Pfaffian化技巧[3]、黎曼Theta函數(shù)方法[4]、Wronskian技巧[5]、Hirota方法[6]、Bell多項(xiàng)式方法[7]等。在直接對(duì)非線性偏微分方程求解較困難的情況下,可以借助Bell多項(xiàng)式、雙線性算子等工具將其雙線性化,進(jìn)而為方程求解拓寬道路。

對(duì)于(3+1)維非線性偏微分方程

(1)

1 雙線性算子與雙線性方程

一般地,對(duì)于給定的自然數(shù)q,引入雙線性Dq微分算子:

(2)

這里k≥1,δj=(-1)r(j),j≥0,r(j)表示被q整除所得的余數(shù)。當(dāng)q取自然數(shù)1,2,3時(shí),δj的取值情況如下：

q=1:+,+,+,+,+,+,…,j=0,1,2,…

q=2:+,-,+,-,+,-,…,j=0,1,2,…

q=3:+,-,+,+,-,+,…,j=0,1,2,…

可以看出,當(dāng)q=1時(shí),Dq算子為一般的導(dǎo)數(shù);當(dāng)q=2時(shí),Dq算子為Hirota雙線性算子;當(dāng)q=3時(shí),Dq算子的表達(dá)式將不同于Hirota算子。

對(duì)于非線性方程(1),通過做應(yīng)變量變換φ(t,z,y,x)=-3[lng(t,z,y,x)]xx,可轉(zhuǎn)化為

(3gxz-2gyt-g3xy)g+3gx(g2x,y-gz)-3g2xgxy+gy(2gt+g3x)=0

(3)

結(jié)合雙線性Dq算子(2)的定義,我們考慮q=2時(shí)的情形,將其轉(zhuǎn)化為雙線性型方程

(4)

下面我們來構(gòu)造方程(1)的三波解。

2 一個(gè)(3+1)維非線性方程的三波解

對(duì)于一個(gè)可積的非線性方程,我們可以運(yùn)用擬設(shè)解的方法來構(gòu)造方程的解,通過平衡方程,確定待定系數(shù),進(jìn)而得到給定方程更多的精確解,來豐富非線性波動(dòng)模型的解集。根據(jù)推廣的三波方法,首先假設(shè)雙線性方程(4)存在解：

g(t,z,y,x)=p1(eα1+e-α1)+p2cosα2+p3sinα3

(5)

其中αi=dit+ciz+biy+aix,i=1,2,3,這里ai,bi,ci,di為待定系數(shù)。將表達(dá)式(5)代入雙線性方程(4),并平衡方程中包含e-α1,eα1,cosα2,sinα2,cosα3,sinα3項(xiàng)的系數(shù),將得到一組關(guān)于ai,bi,ci,di,pi(i=1,2,3)的代數(shù)方程,通過符號(hào)計(jì)算軟件Maple求解出該代數(shù)系統(tǒng)。

2.1 第一種情況

(6)

這里a1,b3,c1,c2,c3,p1,p2,p3為非零的常數(shù),ε=±1,與之相應(yīng)可以確定

(7)

為雙線性方程(4)的解,

(8)

為(3+1)維非線性方程(1)的精確解。為分析該孤波不斷演化的性態(tài),圖1給出t=0和t=4兩幅波形圖。隨著時(shí)間的推移,孤波向著y軸負(fù)半軸方向運(yùn)動(dòng)。

(a) t=0

(b) t=4

2.2 第二種情況

(9)

這里a2,b2,b3,c1,p1,p2,p3為非零的常數(shù),ε=±1,將以上復(fù)系數(shù)解代入(5)式,有以下表達(dá)式成立:

(10)

因而,(3+1)維方程(1)有以下解：

(11)

為描述孤子解的傳播態(tài)勢(shì),圖2給出兩幅將y,z確定后,孤波隨x和t變化的波形圖,可以看出,孤波位于xot坐標(biāo)面的下方。

(a) y=-3

(b) y=3

2.3 第三種情況

(12)

其中ε=±1,a3,b1,b2,b3,c1為非零的常數(shù),結(jié)合以上系數(shù)表達(dá)式,雙線性方程(4)存在解：

(13)

相應(yīng)地,(3+1)維偏微分方程(1)存在解：

(14)

為了解孤子動(dòng)態(tài),我們給出一種x=0,y=1時(shí)的波動(dòng)圖,如圖3所示。

圖3 當(dāng)p1=b1=c1=y=1,p2=2,p3=-2i,a3=i,b2=-i,b3=i,x=0時(shí)的三波解(式(14))Fig.3 The three-wave solution(14)with p1=b1=c1=y=1,p2=2,p3=-2i,a3=i,b2=-i,b3=i,x=0

2.4 第四種情況

(15)

其中,a2,a3,b1,b2,c1,p1,p2,p3為可以自由選取的非零實(shí)常數(shù),將以上實(shí)系數(shù)代入(5)式,得雙線性方程(4)的解：

(16)

于是,我們得到非線性方程(1)的三波解：

(17)

特別地,當(dāng)x=z=2時(shí),孤波(17)的力學(xué)性態(tài)如圖4所示。

圖4 當(dāng)p1=p2=1,b2=-1,p3=a2=b1=z=x=2,a3=c1=3時(shí)的三波解 (式(17))Fig.4 The three-wave solution(17)p1=p2=1,b2=-1,p3=a2=b1=z=x=2,a3=c1=3

3 結(jié)論

對(duì)于高維孤子方程(1),為方便尋求其孤子解,可引入雙線性Dq微分算子,將其轉(zhuǎn)化為相應(yīng)的雙線性方程(4),依據(jù)推廣的三波方法,提前擬設(shè)雙線性型方程(4)具有形式解(5),其中的待定參數(shù)將通過符號(hào)計(jì)算軟件Maple獲取,構(gòu)造出雙線性方程(4)的四種解,其中解(7)、(10)、(13)是復(fù)系數(shù)解,(16)式為實(shí)系數(shù)解。相應(yīng)地,(3+1)維孤子方程(1)的三波解為(8)、(11)、(14)、(17)式,并分別給出四種情況下的六幅波形圖,這種求三波解的方法可以推廣到許多高維孤子方程上。