一类五阶非线性发展方程的解法.docxVIP

一类五阶非线性发展方程的解法

非线性方程的求解一直是数学科学家的研究重点之一。为了获得尽可能多的非线性方程，已经找到了许多方法，如jacobi椭圆函数法、齐次平衡法、双曲正切函数法、检测函数法和非线性变换法。这些方法广泛应用于非线性发展方程的计算。我们考虑了非线性发展方程的五种类型。

其中α,β,γ为常数,s=±1.方程(1)是许多类似方程的一般形式:当s=1,α=30,β=20,γ=10时,方程(1)是五阶KdV方程112~113;当s=1,α=45,β=15,γ=15时,方程(1)是SawadaKotera(SK)方程125;当s=-1,α=-20,β=25,γ=10时,方程(1)是Kaup-Kupershmidt(KK)方程;当s=-1,α=-4/5,β=2/5,γ=5/2时,方程(1)是Caudrey-Dodd-Gibbon-Sawada-Kotera(CDGSK)方程.本文通过改进文献的方法,借助Maple软件,获得了(1)式的新的精确解.

1的常微分方程

我们以关于u(x,y,t)的非线性方程

为例来说明改进的F-展开法.设方程(2)的左边是关于u和u的各阶导数的多项式,则改进的F-展开法可分为以下的5个步骤:

步骤1寻找方程(2)的行波解,令

其中k,l,p是常数.把(3)代入方程(2)可得到关于u(ξ)的常微分方程(ODE):

步骤2设f(ξ)可以表示为

这里ia和a-i是待定的常数,m通过平衡方程(4)中的最高阶导数项和最高阶非线性项的次数得出,这里φ(ξ)满足下面的一阶非线性常微分方程(Klein-Gordon方程539)

其中c,λ,μ为常数.

步骤3把(5)式代入方程(4),并利用方程(6)使方程(4)左边变为一个关于φ′(ξ)和φi(ξ)的多项式.然后可以根据φi(ξ)是线性独立的,得到这个多项式关于φiφ′的系数为零,则得到一系列关于ia和a-i(i=0,1,,n)的代数方程.

步骤4从步骤3中的代数方程中解出ia和a-i(i=0,1,,n)(这个过程可以借助Maple软件来求解).把ia和a-i代回(5)式中,就可以得到方程(2)的解.

步骤5从方程(6)中,获得相应的关于方程(2)的精确解,而方程(6)的解可以参见文献.

2以11213.53.510.23.5

利用改进的F-展开法求解方程(1),根据步骤1、2,得到方程(1)的解:

其中ξ=kx+pt(k,p是非零常数),ai(i=-2,-1,0,1,2)是待定常数,φ满足方程(6).把方程(7)和方程(6)带入方程(1),并在Maple的帮助下,得到下面的方程:

注意到1,φ,φ2,φ3,φ4,φ5,φ-2,φ-3,φ-4,φ-5,φ-6,φ-7是线性独立的,方程(8)成立当且仅当fi=0(i=0,,11).因此,求解方程(1)就化为求解如下的代数方程组

(注:令考虑a-1=0,a-2=0,则得到了文献中的解,这里略去.)

或

从f11=0可得到a-2的2个值:

或

令a1=0,a-1=0,则可以求得:

a-1,a-2不同时为零的情况.

从f5=0可得到a2的2个值:

这样,便得到了方程(1)的一般形式的解

其中k是任意常数,a2的2个值可以由式(9)和(10)确定,a-2的2个值可以由(11)和(12)式确定,0a和p分别由(13)和(14)式确定.

应用方程(6)的解,求得方程(1)众多的解.

1)若取c0=1,λ=-(1+m)2,μ=2m2,则方程(6)有解φ(ξ)=snξ,从而得到方程(1)相应的解:

其中snξ是模数为m(0m1)的Jacobi椭圆正弦函数.

当m→1时,椭圆函数snξ趋于双曲函数tanhξ,即snξ→tanhξ.于是在(16)式中令m→1,便可得到方程(1)如下的孤波解:

当m→0时,椭圆函数snξ趋于三角函数sinξ,即snξ→sinξ.于是在(16)式中令m→0,便可得到方程(1)如下的三角函数解

2)若取2c=m,2λ=-(1+m),μ=2,则方程(6)有解φ(ξ)=dcξ=dnξ/cnξ.于是得到方程(1)相应的解:

当m→0时,椭圆函数dcξ趋于三角函数secξ,即dcξ→secξ.于是在(16)式中令m→0,便可得到方程(1)如下的解:

3)若取c=1,λ=2-m2,μ=2(1-m2),则方程(6)有解φ(ξ)=scξ=snξ/cnξ.于是得到方程(1)相应的解:

其中csξ=(scξ)-1.

当m→0时,椭圆函数scξ趋于正切函数tanξ,即scξ→tanξ.于是在(16)式中令m→0,便可得到方程(1)如下的解:

4)若取c=(1-m)2/4,λ=(1+m)2/2,μ=(1-m2)/2,则方程(6)有解φ(ξ)=cnξ/(1+snξ).于是得到方程(1)相应的解

5)若取c=(1+m′)2