用jacobi椭圆函数展开法求五阶非线性发展方程3.docxVIP

用jacobi椭圆函数展开法求五阶非线性发展方程3

由于工程、物理和生物领域中存在大量非线性方程，非线性方程的解决方案引起了数学和物理工作者的广泛关注。特别是近年来，一些关于非线性方程解的有效方法，如积分法和jacobi椭圆函数传输法，为非线性发展方程的解提供了强有力的工具。

考虑如下五阶非线性发展方程

ut+αu2ux-βuxuxx-γuuxxx+suxxxx=0(1)ut+αu2ux?βuxuxx?γuuxxx+suxxxx=0(1)

其中α,β,γ是常数,s=±1.

最近,文献用构造辅助方程及把函数u表成一组线性独立的函数的方法求出了方程(1)的新显式解.本文用Jacobi椭圆函数展开法,求出了方程(1)的新周期解,我们的结果推广了文献的相应结论.

为寻求方程(1)的行波解,作变换

u(x?t)=φ(ξ)ξ=k(x-ct)+ξ0(2)

其中k,c,ξ0为常数.将(2)代入方程(1)中,然后关于ξ积分一次,得

sk4φξξξξ-γk2φφξξ+12(γ-β)k2φ2ξ+α3φ3-cφ+d=0(3)

其中d为积分常数.

现在用Jacobi椭圆函数展开法求解方程(3).根据齐次平衡法的思想,为使方程(3)中的非线性项φφξξ与最高阶导数项φξξξξ相平衡,构造方程(3)的解为下述形式

φ=a0+a1snξ+a2sn2ξ(4)

其中a0,a1,a2为待定的常数,snξ是模数为m(0m1)的Jacobi椭圆正弦函数.把(4)代入方程(3),并利用如下公式

ddξsnξ=cnξdnξddξcnξ=-snξdnξddξdnξ=-m2snξcnξcn2ξ=1-sn2ξdn2ξ=1-m2sn2ξ

可得

f6sn6ξ+f5sn5ξ+f4sn4ξ+f3sn3ξ+f2sn2ξ+f1snξ+f0=0(5)

其中

f6=120sa2m4k4-4γa22m2k2-2βa22m2k2+13αa32f5=24sa1m4k4-6γa1a2m2k2-2βa1a2m2k2+αa1a22f4=-120sa2(m2+m4)k4-6γa0a2m2k2-32γa1m2k2-12βa1m2k2+2γa22(1+m2)k2+2βa22(1+m2)k2+αa0a22+αa21a2f3=-20sa1(m2+m4)k4-2γa0a1m2k2+3γa1a2(1+m2)k2+2βa1a2(1+m2)k2+2αa0a1a2+13αa31f2=8sa2(2+13m2+2m4)k4-2βa22k2+4γa0a2(1+m2)k2+12γa21(1+m2)k2+12βa21(1+m2)k2+αa0a21+αa20a2-ca2f1=sa1(1+14m2+m4)k4+γa0a1(1+m2)k2-2βa1a2k2+αa20a1-ca1f0=-8sa2(1+m2)k4-2γa0a2k2+12γa21k2-12βa21k2+13αa30-ca0+d

由于1,snξ,sn2ξ,…,sn6ξ线性独立,因此方程(5)成立当且仅当fi=0(i=1,2,…,6),于是求解方程(1)就转化为求解代数方程组

fi=0i=1?2???6

现在分s=1和s=-1两种情况求解方程(1).

(ⅰ)s=1.

从方程f6=0可求得a2的两个值:

a21=3m2k2α(β+2γ+√(β+2γ)2-40α)(6)a22=3m2k2α(β+2γ-√(β+2γ)2-40α)(7)

若a1=0,则f5=0,f3=0,f1=0.把a1=0代入方程f4=0,f2=0,f0=0,可求得

a0=2(1+m2)k2[60m2k2-βa2-γa2]αa2-6γm2k2(8)c=8(2+13m2+2m4)k4-2βa2k2+4γa0(1+m2)k2+αa20(9)d=8a2(1+m2)k4+2γa0a2k2-13αa30+ca0(10)

于是,得到了方程(1)的周期解

u=a2sn2[k(x-ct)+ξ0]+a0(11)

其中a2的两个值由(6)和(7)确定,a0由(8)确定,c由(9)确定.

由于当m→1时,snξ→tanhξ,故在(11)中令m→1,便可得方程(1)的孤波解

u=a2tanh2[k(x-ct)+ξ0]+a0(12)

其中a2的两个值是在(6)和(7)中令m→1来确定,a1和c分别是在(8)和(9)中令m→1来确定.(12)是文献中得到的解.

在(6)～(12)中,如果取α=30,β=20,γ=10,由(11)和(12)便得五阶KdV方程的新周期解和孤波解

u=6m2k2sn2[k(x-56(1-m2+m4)k4t)+ξ0]-2(1+m2)k2(13)u=6k2tanh2[k(x-65k4t)+ξ0]-4k2

如果取α=45,β=15,γ=15,由(11)和(12)便得S