韓冰冰
(盤錦職業(yè)技術(shù)學(xué)院 基礎(chǔ)部,遼寧 盤錦 124000)
近年來(lái),非線性方程在自然科學(xué)與工程應(yīng)用方面占據(jù)非常重要的地位。其中求解非線性方程的精確解是一項(xiàng)具有重要現(xiàn)實(shí)意義的科研工作[1]。同時(shí),非線性方程解的唯一性問(wèn)題一直都是科研工作者的關(guān)注熱點(diǎn)[2]。尤其是求解非線性方程的方法有很多,包括雙曲正切函數(shù)展開(kāi)法[3]、sine-cosine方法[4]、齊次平衡法[5]以及試探函數(shù)法[6]等。盡管這些方法能夠求解非線性方程的孤波解或者沖擊解,但對(duì)于周期解的解析并不是都適宜[3-6]。針對(duì)這一問(wèn)題科研工作者提出了許多新的求解方法。例如,楊娟等人通過(guò)Riccati展開(kāi)法和復(fù)變換獲得非線性分?jǐn)?shù)階Sharma-Tasso-Olever方程和時(shí)空分?jǐn)?shù)階耦合Burgers方程的精確解[7]。熊淑雪等人利用復(fù)變換和整合分?jǐn)?shù)階導(dǎo)數(shù)方法求出了分?jǐn)?shù)階JM方程的新精確解[8]。王輝采用tanh函數(shù)方法對(duì)耦合Kaup-Kupershmidt方程求解,通過(guò)行波約化和Riccati方程進(jìn)行方程轉(zhuǎn)化最終求得顯式行波解[9]。目前,尋找非線性方程新形式的精確解仍然是一件很有意義的工作。為尋求非線性Klein-Gordon方程的更多的周期解,將Jacobi橢圓函數(shù)展開(kāi)法作進(jìn)一步推廣,得到了該方程許多豐富的行波解,然后通過(guò)參數(shù)取值得到了該方程的一些特殊的精確解,包括三角函數(shù)解、雙曲函數(shù)解及它們的混合解,使以前的一些結(jié)論得到了有效推廣。
Jacobi橢圓函數(shù)的定義和性質(zhì)可以參看文獻(xiàn)[10],其中sn(ξ,k)稱為Jacobi橢圓正弦函數(shù),cn(ξ,k)稱為Jacobi橢圓余弦函數(shù),dn(ξ,k)稱為第三類Jacobi橢圓函數(shù),用到Jacobi橢圓函數(shù)的以下性質(zhì):
1.1.1 轉(zhuǎn)化關(guān)系
1.1.2導(dǎo)數(shù)關(guān)系
1.1.3 平方關(guān)系
1.1.4 極限關(guān)系
當(dāng)k→0時(shí):
sn(ξ,k)=sin(ξ),cn(ξ,k)=cos(ξ),dn(ξ,k)=1,
ns(ξ,k)=csc(ξ),cs(ξ,k)=cot(ξ),ds(ξ,k)=csc(ξ),
sc(ξ,k)=tan(ξ),nc(ξ,k)=sec(ξ),dc(ξ,k)=sec(ξ),
sd(ξ,k)=sin(ξ),cd(ξ,k)=cos(ξ),nd(ξ,k)=1
當(dāng)k→1時(shí):
sn(ξ,k)=tanh(ξ),cn(ξ,k)=sech(ξ),dn(ξ,k)=sech(ξ),
ns(ξ,k)=coth(ξ),cs(ξ,k)=csch(ξ),ds(ξ,k)=csch(ξ),
sc(ξ,k)=sinh(ξ),nc(ξ,k)=cosh(ξ),dc(ξ,k)=1,
sd(ξ,k)=sinh(ξ),cd(ξ,k)=1,nd(ξ,k)=cosh(ξ)
其中0≤k≤1是Jacobi橢圓函數(shù)的模數(shù)。
非線性發(fā)展方程的一般形式可寫(xiě)為
式中的F是關(guān)于變?cè)猽,ut,ux,utt,uxt,uxx,...的多項(xiàng)式。引入行波變換
其中k1和c是非零的待定常數(shù),ξ0是任意常數(shù)。將(2)代入(1)得到關(guān)于u(ξ)的常微分方程
設(shè)方程式(3)具有如下形式的行波解
考慮非線性Klein-Gordon方程
將式(2)代入式(5),整理得到
將式(4)代入式(6)。當(dāng)Hi(ξ,k)選取不同的Jacobi橢圓函數(shù)時(shí),計(jì)算得到不一樣的非線性Klein-Gordon方程行波解。本研究根據(jù)以下4種情形進(jìn)行計(jì)算:
2.1.1 情形一
令H1(ξ,k)=sn(ξ,k),H2(ξ,k)=cn(ξ,k),H3(ξ,k)=dn(ξ,k)時(shí),得到
對(duì)式(7)求二階導(dǎo)數(shù)得到
對(duì)式(7)求立方得到
將式(7)、式(8)、式(9)同時(shí)代入式(6)整理得到
在式(10)中令sn(ξ,k)lcn(ξ,k)mdn(ξ,k)n(l,m,n=0,1,2,3)系數(shù)均為零,解得式(5)的以下六個(gè)解:
2.1.2 情形二
令H1(ξ,k)=sd(ξ,k),H2(ξ,k)=cd(ξ,k),H3(ξ,k)=nd(ξ,k)時(shí),得到
u(ξ)=a0+a1sd(ξ,k)+a2cd(ξ,k)+a3nd(ξ,k)
按照情形一的計(jì)算步驟,計(jì)算得到以下六個(gè)式(5)的解:
2.1.3 情形三
令H1(ξ,k)=ns(ξ,k),H2(ξ,k)=cs(ξ,k),H3(ξ,k)=ds(ξ,k)時(shí),得到
u(ξ)=a0+a1ns(ξ,k)+a2cs(ξ,k)+a3ds(ξ,k)
按照情形一的計(jì)算步驟,計(jì)算得到以下六個(gè)式(5)的解:
2.1.4 情形四
令H1(ξ,k)=sc(ξ,k),H2(ξ,k)=nc(ξ,k),H3(ξ,k)=dc(ξ,k)時(shí),得到
u(ξ)=a0+a1sc(ξ,k)+a2nc(ξ,k)+a3dc(ξ,k)
按照情形一的計(jì)算步驟,計(jì)算得到以下六個(gè)式(5)的解:
根據(jù)所求得的行波解,利用Maple 7軟件繪制解的波形圖。以u(píng)1(x,t)圖形為例,u1(x,t)的波形圖見(jiàn)圖1。
圖1 解u1(x,t)的波形圖:ξ0=0,c0=2,α=1,β=1,c=1,k=,x=1
上述得到的非線性Klein-Gordon方程的24個(gè)行波解,如果考慮它們的極限情況,即當(dāng)k→1時(shí),Jacobi橢圓函數(shù)退化為雙曲函數(shù);當(dāng)k→0時(shí),Jacobi橢圓函數(shù)退化為三角函數(shù)。
2.2.1情形一
當(dāng)k→0時(shí),非線性Klein-Gordon方程,即式(5)的24個(gè)行波解退化為以下三角函數(shù)形式的解:
2.2.2 情形二
當(dāng)k→1時(shí),得到的非線性Klein-Gordon方程(5)的24個(gè)行波解退化為以下雙曲函數(shù)形式的解:
本研究以Klein-Gordon方程為例,介紹了Jacobi橢圓函數(shù)的性質(zhì)和展開(kāi)法,利用該方法求得非線性Klein-Gordon方程一系列新的精確周期解,并且充分考慮到極限情況,得到的三角函數(shù)及雙曲函數(shù)的解,從而揭示了求解非線性Klein-Gordon方程精確行波解的理論和技巧。此外,該方法具有一定的普遍性,可以用來(lái)求解更多的非線性發(fā)展方程,例如非線性Schrodinger方程、KP方程、Ginzburg-Landau方程等。