黃 欣

(1. 四川財(cái)經(jīng)職業(yè)學(xué)院 基礎(chǔ)部, 四川 成都 610101; 2. 四川師范大學(xué) 數(shù)學(xué)與軟件科學(xué)學(xué)院, 四川 成都 610066)

非線性波現(xiàn)象有很多,如色散、耗散、反應(yīng)和對(duì)流等,它們出現(xiàn)在離子體物理、固體物理、流體力學(xué)、光學(xué)纖維等各類科學(xué)領(lǐng)域中.非線性現(xiàn)象常常與非線性偏微分方程聯(lián)系緊密,比如在地震研究中,地震波可以用非線性偏微分方程描述,又可以用環(huán)孤粒子或者是環(huán)孤立波來描述;臺(tái)風(fēng)一些非常重要的結(jié)構(gòu)特征可以通過二維的Euler方程來描述;易引發(fā)極端氣候的大氣阻塞現(xiàn)象可以用KdV孤立子來描述.因此在研究非線性現(xiàn)象的過程中,關(guān)于非線性偏微分方程的研究占到了一個(gè)重要的角色.人們把目光集中到了求這些對(duì)應(yīng)非線性發(fā)展方程的精確解[1-8]上,因?yàn)檫@些解正是揭示各類非線性現(xiàn)象的重要工具,為得到這些非線性發(fā)展方程的精確解,許多直接的方法被提出,包括齊次平衡法[9]、雙曲正切法[10]、Exp-函數(shù)法[11]、廣義Riccati方程法[12]等.

首次積分法是求解非線性波動(dòng)方程的有效方法之一,于2002年由Z. S. Feng[13]在求解Burgers-KdV方程中提出.其基本思想是利用行波變換將非線性偏微分方程轉(zhuǎn)化為常微分方程,再構(gòu)造出這個(gè)常微分方程的一個(gè)首次積分,由此就可以把此常微分方程變成一個(gè)一階可積的常微分方程組,通過直接積分的辦法求解此常微分方程組,進(jìn)而得到原方程的精確解.目前這個(gè)方法已經(jīng)被成功的廣泛應(yīng)用于各種非線性發(fā)展方程中,比如Fisher方程[14]、非線性Schr?dinger方程[15]、Gardner方程[16]等.

本文考慮如下(3+1)維修正KdV-Zakharov-Kuznetsev方程

ut+αu2ux+uxxx+uxyy+uxzz=0,

(1)

其中α為實(shí)常數(shù).文獻(xiàn)[17]用(G′/G)-展開法獲得了方程(1)的精確解,文獻(xiàn)[18]用Exp-函數(shù)法研究了此方程.本文將通過首次積分法和符號(hào)運(yùn)算系統(tǒng)Mathematica研究(3+1)維修正KdV-Zakharov-Kuznetsev方程(1),并得出了方程的新精確解.

1 首次積分法

對(duì)于給定的非線性偏微分方程

P(u,ut,ux,uxx,utt,uxt,…)=0,

(2)

假設(shè)方程存在行波解

u=u(ξ),ξ=x-ct,

其中c表示波速,則原方程就變?yōu)?/p>

P(u,u′,u″,…)=0.

(3)

再令

X(ξ)=u(ξ),Y(ξ)=u′(ξ),

(4)

則方程(3)等價(jià)于常微分方程組

(5)

通過常微分方程的定性理論[19]可知,對(duì)于常微分方程組(5),若能找到其2個(gè)首次積分,就能直接獲得其通解,但一般說來,由于沒有系統(tǒng)的構(gòu)造首次積分的理論,即使想要獲得一個(gè)首次積分都很困難.Z. S. Feng[13]提出的首次積分法用到了代數(shù)環(huán)論中的除法定理[20],使得目標(biāo)常微分方程組(5)能夠轉(zhuǎn)化為一階可積的常微分方程,通過轉(zhuǎn)換之后的方程有很多方法可以比較容易處理求解問題.

除法定理假設(shè)P(w,z)、Q(w,z)都是定義在復(fù)數(shù)域C[w,z]上的多項(xiàng)式,并且多項(xiàng)式P(w,z)在復(fù)數(shù)域C[w,z]上是不可約的．如果多項(xiàng)式Q(w,z)不包含多項(xiàng)式P(w,z)的所有零點(diǎn),則復(fù)數(shù)域C[w,z]上必定存在一個(gè)多項(xiàng)式G(w,z),使得

Q(w,z)=P(w,z)·G(w,z).

2 (3+1)維修正KdV-Zakharov-Kuznetsev方程的精確解

假設(shè)方程(1)有如下形式的行波解

u(x,y,z,t)=u(ξ),ξ=x+y+z-ct,

(6)

其中c為待定參數(shù).將(6)式帶入(1)式得

-cu′+αu2u′+3u?=0,

(7)

將(7)式積分一次,并令積分常數(shù)為0,可得

(8)

令X(ξ)=u(ξ),Y(ξ)=u′(ξ),則方程(8)等價(jià)于

(9)

按照首次積分法的步驟,假定X(ξ)和Y(ξ)是方程組(9)的非平凡解,則方程組(9)的首次積分為

其中ai(X(ξ))(i=0,1,…m)是X(ξ)的多項(xiàng)式,am(X(ξ))≠0.根據(jù)除法定理可知,在復(fù)數(shù)域C上會(huì)有

q(X(ξ),Y(ξ))=g(X(ξ))+h(X(ξ))Y(ξ),

使得

僅考慮(11)式中的2種情況,即假定m=1和m=2.

情形1假定m=1,則有

p(X,Y)=a0(X)+a1(X)Y,

(12)

代入(11)式后比較兩邊Yi(i=2,1,0)的系數(shù)得

(13)

(14)

因?yàn)閍i(X)(i=0,1)為多項(xiàng)式，則從(13)式可知a1(X)為常數(shù),且h(X)=0.為運(yùn)算的簡便,不妨取a1(X)=1,通過平衡a0(X)和g(X)的系數(shù),可以得出deg[g(X)]=1,相應(yīng)的deg[a0(X)]=2.

假定

g(x)=A1X+B0,A1≠0,

則有

(16)

其中A0為積分常數(shù).

李離愀然不樂：“我父親說，龍化為魚，十年為期，意思是長安城不出十年，就會(huì)被叛軍打破，你說到那個(gè)時(shí)候，寫這首詩的吏部郎中王維大人，會(huì)不會(huì)降？”

將(16)式代入(15)式并比較兩邊Xi的系數(shù)可得

將(17)式代入(10)式得

(19)

將(19)式代入(9)式中,即可得(9)式的精確解,故原方程的精確解為

(20)

其中ξ0為任意常數(shù).同理將(18)式代入(10)式得

(21)

將(21)式代入(9)式中,即可得(9)式的解,故原方程的精確解為

(22)

其中ξ0為任意常數(shù).

情形2假定m=2,則有

p(X,Y)=a0(X)+a1(X)Y+a2(X)Y2, (23)

通過比較方程兩邊Yi的系數(shù),可得

(24)

(25)

=g(X)a1(X)+h(X)a0(X),

(26)

(27)

因?yàn)閍i(X)(i=0,1)為多項(xiàng)式,則從(24)式可知a2(X)為常數(shù),且h(X)=0.為運(yùn)算簡便,不妨取a2(X)=1,平衡a0(X)、a1(X)、g(X)的系數(shù),可以得出deg[g(X)]=1.

假定

g(X)=A1X+B0A1≠0,

則可得a1(X)和a0(X)為

(28)

(29)

其中D為積分常數(shù).

將a0(X)、a1(X)、g(X)代入(27)式, 令Xi的同次冪系數(shù)為0,則可以得到一個(gè)系統(tǒng)的非線性代數(shù)方程組,借助符號(hào)計(jì)算系統(tǒng)Mathematica解此方程可得

(30)

將(30)式代入(10)式可得

(31)

將(31)式代入(9)式,可得(9)式的解,故(3+1)維修正KdV-Zakharov-Kuznetsev方程的精確解為

u(x,y,z,t)=

其中ξ0為任意常數(shù).

3 結(jié)論和展望

在交換代數(shù)環(huán)論的基礎(chǔ)上,運(yùn)用首次積分法,借助符號(hào)計(jì)算系統(tǒng)Mathematica,本文對(duì)(3+1)維修正KdV-Zakharov-Kuznetsev方程進(jìn)行了研究,獲得了該方程的精確行波解,這些解均為新解.下一步, 打算對(duì)此方法進(jìn)行改進(jìn),使其可以往高階的、高維的非線性波動(dòng)方程上推廣研究,此外還準(zhǔn)備結(jié)合分岔理論,結(jié)合相圖,分析行波解的特性.

[1]Guo B L, Ling L M. Rogue wave, Breathers and Bright-Dark-Rogue solutions for the coupled Schr?dinger equations[J]. Chinese Phys Lett,2011,28:110202.

[2]Guo B L, Ling L M, Liu Q P. Nonlinear Schr?dinger equation: generalized Darboux transformation and rogue wave solutions[J]. Phys Rev,2012,E85:026607.

[3]Guo B L, Ling L M, Liu Q P. High-order solutions and generalized Darboux transformations of derivative nonlinear Schr?dinger equations[J]. Stud Appl Math,2013,130:317-344.

[4]鄧艷萍,周鈺謙. 非線性Schr?dinger方程組初值問題的駐波解[J]. 四川師范大學(xué)學(xué)報(bào):自然科學(xué)版,2013,36(2):243-246.

[5]劉倩,周鈺謙,劉合春. 廣義Hirota-Satsuma耦合KdV方程的精確解[J]. 四川師范大學(xué)學(xué)報(bào):自然科學(xué)版,2011,34(3):335-339.

[6]Ling L M, Liu Q P. A long waves-short waves model: Darboux transformation and soliton solutions[J]. J Math Phys,2011,52:053513.

[7]Huang X. The investigation of solutions to the coupled Schr?dinger-Boussinesq equations[J]. Abstr Appl Anal,2013,2013:170372.

[8]Huang X, Guo B L, Ling L M. Darboux transformation and novel solutions for the long wave-short wave model[J]. J Nonl Math Phys,2013,20:514-528.

[9]Wang M L, Zhou Y B, Li Z B. Application of a homogeneous balance method to exact solutions of nonlinear equations in mathematical physics[J]. Phys Lett,1996,A216:67-75.

[10]Fan E G. Extended tanh-function method and its applications to nonlinear equations[J]. Phys Lett,2000,A277:212-218.

[11]Shin B C, Darvishi M T, Barati A. Some exact and new solutions of the Nizhnik-Novikov-Vesselov equation using the exp-function method[J]. Comput Math Appl,2009,58:2147-2151.

[12]Yan Z Y, Zhang H Q. New explicit solitary wave solutions and periodic wave solutions for Whitham-Broer-Kaupequation in shallow water[J]. Phys Lett,2001,A285:355-362.

[13]Feng Z S. The first integral method to study the Burgers-Korteweg-de Vries equation[J]. Phys Lett,2002,A35:343-349.

[14] Raslan K R. The first integral method for solving some important nonlinear partial differential equations[J]. Nonlinear Dynamics,2008,53:281-286.

[15] Taghizadeh N, Mirzazadeh M, Farahrooz F. Exact solutions of the nonlinear Schr?dinger equation by the first integral method[J]. J Math Anal Appl,2011,374:549-553.

[16]Taghizadeh N, Mirzazadeh M, Farahrooz F. New exact solutions of some nonlinear partial differential equations by the first integral method[J]. Appl Appl Math,2010,5(10):446-456.

[17]Zayed E M E. Traveling wave solutions for higher dimensional nonlinear evolution equations using the (G′/G)-expansion method[J]. Appl Math Informatics,2010,28:383-395.

[18]Naher H, Abdullah F A, Akbar M A. New traveling wave solutions of the higher dimensional nonlinear partial differential equation by the exp-function method[J]. J Appl Math,2012,2012:575387.

[19]丁同仁,李承治.常微分方程教程[M]. 2版. 北京:高等教育出版社,2004.

[20]Bourbaki N. Commutative Algebra[M]. Paris:Herman,Publishers in Arts and Science,1972.