利用符號(hào)計(jì)算求解高階非線性演化方程的新方法

史振華， 夏鐵成

(上海大學(xué)理學(xué)院，上海200444)

眾所周知，非線性科學(xué)是現(xiàn)代科學(xué)的核心，尋找非線性演化方程的精確解是數(shù)學(xué)物理研究領(lǐng)域的重要工作之一.到目前為止已出現(xiàn)許多方法，如反散射法［1］、Hirota雙線性方法［2］、Backlund變換法［3］、齊次平衡法［4］、tanh函數(shù)法［5-6］、Exp函數(shù)法［7］、Jacobi橢圓函數(shù)展開法［8-9］、F擴(kuò)展法［10-12］等.

本研究將運(yùn)用推廣的方法求解如下(2+1)維Bogoyavlenskii破裂孤子方程：

和(3+1)維Kadomtsev-Petviashvili(K-P)方程

對(duì)于一個(gè)給定的非線性演化方程，自變量X= (x，y，z，…，t)，因變量為u，即

接下來(lái)，按以下步驟來(lái)求解u：

步驟1 假設(shè)u(x，y，z，…，t)=u(ξ)行波變化，并將式(3)轉(zhuǎn)變?yōu)槌Ｎ⒎址匠?/p>

式中，G=G(ξ)，ξ=c1x+c2y+c3z+…+dt，且滿足方程

式中，am，…，a0，c1，c2，c3，…，d，λ，μ為待定常數(shù)，正整數(shù)m可以通過齊次平衡法來(lái)確定;

步驟4 通過解方程組，求出常數(shù)am，…，a0，c1，c2，c3，…，d，λ和μ.又由于方程(6)的解是已知的，故可以將am，…，a0，c1，c2，c3，…，d，λ，μ和方程(6)的通解代入式(5)，就可以得到式(3)的行波解.

下面采用該方法求解2個(gè)著名的方程.

首先，求解如下(2+1)維Bogoyavlenskii破裂孤子方程：

其等價(jià)形式為

將式(7)和(8)進(jìn)行行波變換，得

因此，可將式(7)和(8)轉(zhuǎn)化為如下常微分方程：

式中，G=G(ξ)，ξ=c1x+c2y+dt滿足方程

利用齊次平衡法，可以將u(ξ)和v(ξ)表示為

式中，c1，b0，b2，a0，λ和μ為任意常數(shù).

利用式(17)，式(15)和(16)可以寫為

將方程(14)的通解代入式(18)和(19)，可以得到方程(7)和方程(8)的3種形式的精確解.

當(dāng)λ2-4μ＞0時(shí)，得到如下的雙曲函數(shù)解：

當(dāng)λ2-4μ＜0時(shí)，得到如下三角函數(shù)解：

當(dāng)λ2-4μ=0時(shí)，可得如下有理解：

對(duì)(3+1)維Kadomtsev-Petviashvili(K-P)方程

同樣可以進(jìn)行行波變換，得

將方程(26)轉(zhuǎn)換為

式中，G=G(ξ)，ξ=c1x+c2y+c3z+dt滿足方程

利用齊次平衡法比較式(28)中的u″″和uu″，得到m=2.可以將u(ξ)表示為

式中，c1，c2，c3，a0，λ和μ為任意常數(shù).

利用式(32)，式(31)可以表示為

將方程(30)的通解代入式(33)，可得方程(26)的3種形式的精確解.

當(dāng)λ2-4μ＞0時(shí)，可得到如下雙曲函數(shù)解：

當(dāng)λ2-4μ＜0時(shí)，可得到如下三角函數(shù)解：

當(dāng)λ2-4μ=0時(shí)，可得到如下有理解：

3 結(jié)束語(yǔ)

［1］ ABLOWITZ M J，CLARKSONP A.Soliton nonlinear evolution equations and inverse scattering［M］.New York：Cambridge University Press，1991.

［2］ HIROTAR.Exact solutions of the Korteweg-de-Vries equation for multiple collisions of solitons［J］.Phys Rev Lett，1971，27：1192-1194.

［3］ 陳登遠(yuǎn).Backlund變換與n孤子解［J］.數(shù)學(xué)研究與評(píng)論，2005，25(3)：479-488.

［4］ 王明亮，李志斌，周宇斌.齊次平衡原則及其應(yīng)用［J］.蘭州大學(xué)學(xué)報(bào)：自然科學(xué)版，1999，35(3)：8-16.

［5］ FANE G.Extended tanh-function method and its applications to nonlinear equations［J］.Phys Lett A，2000，277：212-218.

［6］ 李德生，張鴻慶.改進(jìn)的tanh函數(shù)方法與廣義變系數(shù)KdV和MKdV方程新的精確解［J］.物理學(xué)報(bào)，2003，52(7)：1569-1573.

［7］ HEJ H，WUX H.Exp-function method for nonlinear wave equations［J］.Chaos，Solitons and Fractals，2006，30(3)：700-708.

［8］ LIUS K，F(xiàn)UZ T，LIUS D，et al.Jacobi elliptic function expansion method and periodic wave solutions of nonlinear wave equations［J］.Phys Lett A，2001，289：69-74.

［9］ FUZ T，LIUS K，LIUS D，et al.New Jacobi elliptic function expansion and new periodic solutions of nonlinear wave equations［J］.Phys Lett A，2001，290：72-76.

［10］ ZHOUY B，WANGM L，WANGY M.Periodic wave solutions to a coupled KdV equation with variable coefficients［J］.Phys Lett A，2003，308：31-36.

［11］ LIUJ B，YANGK Q.The extended F-expansion method and exact solutions of nonlinear PDEs［J］.Chaos，Solitons and Fractals，2004，22(1)：111-121.

［12］ WANGD S，ZHANGH Q.Further improved F-expansion method and new exact solutions ofKonopelchenko-Dubrovsky equation［J］.Chaos，Solitons and Fractals，2005，25：601-610.