(2+1)维变系数非线性KP方程新推广解_靳玲花.pdf

I S S N 1 0 0 9-8 9 8 4C N 2 2-1 3 2 3/N长春工程学院学报(自然科学版)2 0 2 2年 第2 3卷 第4期J.C h a n g c h u n I n s t.T e c h.(N a t.S c i.E d i.),2 0 2 2,V o l.2 3,N o.4 2 5/2 51 2 5-1 2 8d o i:1 0.3 9 6 9/j.i s s n.1 0 0 9-8 9 8 4.2 0 2 2.0 4.0 2 5(2+1)维变系数非线性K P方程新推广解收稿日期:2 0 2 2-1 0-2 7基金项目:河套学院科学技术研究项目自然科学一般项目(H Y Z Y 2 0 2 1 0 5)内蒙古自治区高等学校科学研究项目(N J Z Y 2 1 1 7 3)作者简介:靳玲花(1 9 8 4-),女(汉),山西长治人,硕士,讲师主要研究变系数非线性方程孤子解。靳玲花,白 慧,李珊珊(河套学院数学与计算机系,内蒙古 巴彦淖尔 0 1 5 0 0 0)摘 要:为适应非线性发展方程包括变系数非线性发展方程求解的需要,试图探求辅助方程多样化和解的形式的一般化,对王明亮教授提出的(G/G)-展开法进行了更有意义的推广。为验证此推广的可靠性与有效性,将它再次应用到(2+1)维广义圆柱变系数K P方程中以期寻求内涵更为丰富的精确解,最终取得了成功。说明此推广具有可靠性和安全性。关键词:发展方程;精确解;推广的(G/G)-展开法;(2+1)维广义圆柱K P方程中图分类号:O 1 7 5.2 9文献标志码:A 文章编号:1 0 0 9-8 9 8 4(2 0 2 2)0 4-0 1 2 5-0 40 引言非线性科学涉及范围非常广泛,比如力学、数理化学、天文气象学、生物环境学、医学、航天等一些比较重要的领域,它是2 0世纪自然科学继量子力学和相对论两大领域的又一重大发展领域。从数学和物理的视觉来看,通过求解非线性微分方程可以很好地解释自然界中出现的各类非线性现象。为了更好地去理解这种非线性现象,采用构造非线性微分方程的解的方法,而这种做法已经成为当代非线性科学研究领域的重要课题。目前在还没有完全获得系统地处理非线性问题的方法的情况下,不同专业的学者分别总结出了自己特有的研究方法,但一般被认可的范畴包括:孤立子、分形以及混沌。孤立子也叫孤立波,它是指一大类非线性偏微分方程的许多具有相同特殊性质的解,以及它所蕴含的物理现象。孤立子具有一些主要的特性:1)孤立子是波动问题中一种能量有限的局域解;2)能量大多集中在一个很小的区域内(或能在空间给定的区域内稳定存在);3)孤立波相互作用时会出现弹性散射现象。这些性质揭示了孤立子的内涵,同时我们称具有孤立子解的非线性发展方程为孤子方程。孤立子理论包含的内容和研究方法非常丰富,尤其是近十几年来研究队伍不断扩大,研究成果令人瞩目。孤立子理论中的两大主要问题是构造孤立子方程的精确解和研究非线性系统的可积性。随着孤立子理论的发展,已经总结了许多构造非线性方程精确解的方法,如B c k l u n d变换法、D a r b o u x变换法、相似 约 化 法、H i r o t a双 线 性 方 法、反 散 射 方 法(I S T)、延拓结构法、齐次平衡法、经典和非经典李群法、F-展开法、代数几何法、P a i n l e v 分析法和(G/G)-展开法等。本论文采用的就是(G/G)-展开法。由于非线性微分方程的解不再满足线性方程的叠加原理,所以通常很难像线性方程那样用一个统一的方法对其求解。进入2 0世纪以来,对于常系数非线性微分方程,前述诸方法已被大量应用,而变系数模型能够更精确地描述物理、力学问题,特别是高阶的变系数方程。因此,变系数方程有比较广泛的应用性。考虑被数学家和物理学家普遍感兴趣的方程,(2+1)维广义圆柱变系数K a d o m t s e v-P e t v i a s h-v i l l i(K P)方程:(ut+6u ux+ux x x)x+a(t)ux+b(t)uy y=0,(1)式(1)可以化为具有重要物理意义的某些特殊形式的方程,若取a(t)=12t,b(t)=0,则方程就是著名的圆柱K d V方程:ut+6u ux+ux x x+12tu=0,(2)如若取a(t)=12t,b(t)=3t2,则方程化为著名的圆柱K P方程:(ut+6u ux+ux x x)x+12tux3t2uy y=0。(3)K P方程有着广泛的物理背景,在流体力学、等离子体物理和气体动力学等领域有重要作用,可作为描述(2+1)维浅水波和等离子体中的离子声波的模型方程。任何一个模型或系统,尤其是变系数非线性模型或系统,我们可以直接进行系统行为与性能的分析,建立系统的可靠性模型。可靠性理论是正在发展的理论,对许多方面的认识还在不断地深化。对于一些复杂的可靠性问题,都还需要进一步、更广泛深入的研究。本文中对于(2+1)维广义圆柱变系数K P方程的求解,前期已经做了许多工作1。主要方法选自王明亮等2-32 0 0 8年提出的(G/G)-展开法。此方法简单、有效、实用,但作者发现还是有很大的改进空间。接下来就是将(G/G)-展开法进行更一般化的推广后,再次应用在此方程中。1 推广的(G/G)-展开法简介给出如式(4)形式的非线性微分方程:P(u,ut,ux,ut t,ux x,)=0。(4)1)寻求方程(4)如式(5)的行波变换:u(x,t)=u(),=xc t。(5)若是变系数非线性方程,则可以是如式(6)的行波变换:u(x,t)=u(),=k(x)+l(y)+p(t)+q,(6)式中:k(x),l(y),p(t)是待定函数;q,c是待定常数;利用变换式(5)或式(6),可将式(4)化为如式(7)的非线性微分方程(常系数或者变系数):Q(u,u,u,L)=0。(7)2)假设方程(7)有如式(8)形式的解:u()=a0+ni=1aiG()G()()i+biG()G()()-i,(8)式中:an和bn不同时为0,a0,ai,bi(i=1,2,n)都是待定常数。式(8)中的整数n可以通过平衡式(7)中u()的最高阶非线性项和最高阶微分项来确定。G=G()满足如式(9)的二阶常微分方程:G+G=0,(9)或者G+G+G2=0,(1 0)式中和都是任意常数,此辅助方程的论述可见文献1-3。若ai=0,式(8)将化为u()=a0+ni=1bi(G()G()i,bn0。(1 1)但如果bi=0,i=1,2,n,则(8)式将变为如式(1 2)的形式:u()=a0+ni=1ai(G()G()i,an0。3(1 2)3)把式(8)代入式(7)后再利用式(9)或式(1 0)化为(G/G)的幂次多项式,再令各次的系数为0,则得到一个以a0,ai,bi(i=1,2,n),以及k(x),l(y),p(t)为未知量,或者以q,c为未知数的非线性代数方程组。4)求解3)中的代数方程组确定相关的未知量。另外,式(9)或式(1 0)的解都是我们所熟悉的(这里就不再赘述),将其及a0,ai,bi(i=1,2,n),k(x),l(y),p(t),q或者c代入式(8),便可以得到方程(4)的精确行波解。2(2+1)维广义圆柱变系数K P方程推广精确解(2+1)维广义圆柱变系数K P方程为:(ut+6u ux+ux x x)x+a(t)ux+b(t)uy y=0,(1 3)因为该方程是变系数的,所以应该引入行波变换式(6),则该方程变形为:a(t)k(x)+b(t)l(y)+k(4)6(x)u+6k 2(x)u 2+6k(x)u+c(t)k(x)+b(t)l 2(y)+2k(x)+k 2(x)+2k(x)k(x)u,+6k 2(x)u u+(5k 2(x)+k 2(x)k(x)u+k 4(x)u(4)=0。(1 4)由齐次平衡原则得n=2,则应设式(1 4)解的形式为u()=a0+a1G G+a2G G()2+b1G G()-1+b2G G()-2=a0+a1()+a22()+b1-1()+b2-2(),(1 5)这里G()满足二阶线性常微分方程式(9)。由式(1 5)利用式(9),则得到:u=-2a23-a12-2a2+2b2-3+b1-2+2b2-1+b1-a1,u=6a24+2a13+8a22+2a1+6b22-4+2b12-3+8b2-2+2b1-1+2b1+2a22,u=-2 4a25-6a14-4 0a23-8a12-1 6a22+2 4b23-5+6b13-4+4 0b22-3+621长春工程学院学报(自然科学版)2 0 2 2,2 3(4)8a12-2+1 6b2-1+2b1-2a12,u(4)=1 2 0a26+2 4a15+2 4 0a24+4 0a13+1 3 6a222+1 6a12+1 2 0b24-6+2 4b14-5+2 4 0b23-4+4 0b13-3+1 3 6b22-2+1 6b12-1+1 6b2+1 6a23。将前述各式代入式(1 4),合并为关于的各幂次多项式,并令各幂次系数为0,得到包含1 4个方程的一组很庞大的代数方程组,这里由于篇幅有限就不一一呈现。这样的方程组求解过程非常复杂,必须借助M a t h e m a t i c a软件的帮助,并经过认真仔细的分析筛选,最终求出了一组解:a0=-43,a1=2,a2=-2,b1=-2,b2=-22k 2(x),k(x)=x+C1,b(t)=C a(t),l(y)=C2y2+C4y+C5,c(t)=tl 2(y)t0b(t)d t+C3,式 中:C1,C2、C3、C4和C5均 为 积 分 常 数;=-14和C均为常数。将前述解连同式(9)的解一并代入式(1 5)即得到方程(1 3)的解:当0时,u1()=-43+2A1s i n()-A2c o s()A1c o s()+A2s i n()-2A12+A22A1c o s()+A2s i n()2-2A1c o s()+A2s i n()A1s i n()-A2c o s()-22A1c o s()+A2s i n()2A12+A22,(1 6)当0时,u2()=-43+2A1-s i n h(-)+A2-c o s h(-)A1c o sh(-)+A2s i n h(-)-2A2c o s(-)+A1s i n h()2A1c o s()+A2s i n h()2-2A1c o s h(-)+A2s i n h(-)A1-s i n h(-)+A2-c o s h(-)-22A1c o s()+A2s i n h()2A2c o s(-)+A1s i n h()2,(1 7)式中:=x+C2y2+C4ytl 2(y)t0b(t)dt+l;l=C1+C3;k(x),l(y),c(t)满 足 关 系 式c(t)k(x)=1-b(t)l 2(y)。接下来给出当辅助方程为式(1 0)时(2+1)维变系数K P方程的其中一组解:u3()=-26k 2(x)+2 k 2(x)-2c 2s e c h22()1+t a n h2()|-C+(C+)t a n h2()|)-2k 2(x)c 2s e c h22()1+t a n h2()|-C+(C+)t a n h2()|)|2+162c 2s e c h22()1+t a n h2()|-C+(C+)t a n h2()|)|-1-11 23c 2s e c h22()1+t a n h2()|-C+(C+)t a n h2()|)|-2,(1 8)式中:=-C1x+C2y2+C4y+tC1-l 2(y)C1t0b(t)dt+l;l=C0+C3+C5;k(x),l(y),c(t),满足关系式c(t)k(x)=1-b(t)l 2(y)。至此,从前述的求解过程可以看出:在变系数的非线性发展方程中,维数只简单增加一维,其计算量却是增加数倍;由于其计算过程比较复杂,所以在计算的时候需要谨慎分析;与文献1中所求的解作比较,发现本章中的