非线性偏微分方程

1、非线性偏微分方程及其几种解法综述姓名：柏宝红学号：BY1004120目录3.3.7.1、绪论1.1 背景.1.2 现状2、非线性偏微分方程的几种解法1.0逆算符法10齐次平衡法1.1Jacobi椭圆函数方法1.2辅助方程方法1.4F-展开法15双曲正切函数展开法171、绪论以应用为目的，或以物理、力学等其他学科问题为背景的微分方程的研究，不仅是传统应用数学中一个最主要的内容，也是当代数学的一个重要组成部分.它是数学理论与实际应用之间的一座重要桥梁，研究工作一直十分活跃，研究领域日益扩大。目前微分方程研究的主体是非线性微分方程，特别是非线性偏微分方程(NLPDE),很多意义重大的自然科学和工程技

2、术问题都可归结为非线性偏微分方程的研究.现实生活的许多领域内数学模型都可以用NLPDE来描述，很多重要的物理、力学等学科的基本方程本身就是NLPDE,另外，随着研究的深入，有些原先可用线性微分方程近似处理的问题，也必须考虑非线性的影响，所以对NLPDE的研究，特别是NLPDE求解精确解的研究工作就显示出了很重要的理论和应用价值，但是数学研究的结果，在目前还未能提供一种普遍有效的求精确解的方法.20世纪50年代以来，人们对非线性现象的研究中提出了“孤子”的概念，进而使得对NLPDE求解的研究成为非线性科学中的热点。下面介绍一下孤立子理论的研究背景、研究现状。背景孤立子理论己经成为应用数学和数学物

3、理的一个重要组成部分，在流体力学，等离子物理，经典场论，量子论等领域有着广泛的应用。随着近代物理学和数学的发展，早在1834年由英国科学家Russell发现的孤立波现象近二十多年来引起了人们的极大关注，对这一现象的兴趣与日俱增.这是因为一方面孤立子具有粒子和波的许多性能，在自然界中有一定的普遍性，利用孤立子理论也成功地解释了许多物理上长期用经典理论未能解答的现象；另一方面，随着孤立子物理问题的深入研究，孤立子的数学理论也应运而生，并已初步形成比较完善的理论体系。孤立子理论自1965年由Zabusky和Kruskal对孤立子(Soliton,简称孤子)命名后得到了迅速地发展.究其原因是孤波现象无

4、所不在，从天上涡旋星系的密度波，线，超流氨一3,超导JosePhson结，磁学，结构相变，液晶，流体动力学以及基本粒子等，都与孤子有关.其发展大致可分三个阶段：第一阶段，主要是在19世纪.最早讨论孤立子问题的是ScottRussell。1844年英国工程师Russell发现船在运河中快速行驶着，当这条船突然停止时，在船头附近产生了一个光滑的、像小山包一样的水波，然后这个水波离开船头保持它的形状和速度保持不变，接着这个水波的高度逐渐减少，最后在运河的一个拐弯处消失掉，他把这种水波称为孤立波，认为它就是流体运动的一个稳定解.直到1895年,荷兰阿姆斯特丹大学的Korteweg教授和他的学生devr

5、ies才一成功导出了著名KdV方程，求出了与Russell描述一致的即具有形状不变的脉冲状的孤立波解，在理论上证实了孤立波的存在，并对孤立波现象作了较为完整的分析，解释了Russell的浅水波，解决了这个问题。他们的数学模型为ut6uuxuxxx0(1.1)孤立波解为:c2，cu1(x,t)=-sech(x-ct)22后人称u1(x,t)为1一孤立子解，如果令X=x-ct,那么5在平图1.1光滑孤立子U1在上-u平面上的图形1965年美国数学家Kruskal和abusky对KdV方程的孤立波解进行数学模拟，他们发现两个孤立波相撞之后，各自的运动方向和大小形状都保持不变.这种性质与物理中粒子的性

6、质类似，因此他们称这种孤立波为孤立子.在通常情况下，人们把孤立波和孤立子混为一谈，不把它们区别开来。与此同时，在18761882年发现的Backlund变换，成为后来发展孤子理论的重要基础。第二阶段大致可划在19551975年。1955年，Fermi,PastaUlam(FPU)将64个质点用非线性弹簧连成一条非线性振动弦，用计算机计算了一维非线性晶格在各个振动模之间的转换。初始时，这些谐振子的所有能量都集中在一个质点上，其他63个质点的初始能量为零。按照经典的理论，只要非线性效应存在，就会有能量均分，各态历经等现象出现，即任何微弱的非线性相互作用，可导致系统的非平衡状态向平衡状态的过渡。但实

7、际计算的结果却与经典理论是背道而驰.实际上，经过相当长时间之后，能量似乎又回到了原来的初始分布，这就是著名的FPU问题。由于FPU问题是在频域空间考察的，未能发现孤波解，因此该问题未能得到正确的解释。后来，人们发现可以把晶体看成具有质量的弹簧拉成的链条，这恰好是Fermi研究的情况。Toda研究了这种模式的非线性振动，得到了孤波解，使FPU问题得到正确的解答，从而进一步激发起人们对孤立波的研究兴趣。1965年，zabusky和Kxusal对等离子体中孤立波的相互碰撞过程进行计算机数值模拟，进一步证实了孤立波在碰撞前后波形和速度保持不变的论断，并且把它命名为孤立子（soliton）,它是指一大类

8、非线性偏微分方程的许多具有特殊性质的解，以及具有相应的物理现象，它的性质具体为：（1）能量比较集中；（2）孤立子相互碰撞时具有弹性散射现象。从此孤立子理论的研究工作得到了迅速发展。第三阶段（1973至今），把孤子概念及理论广泛应用于物理学，生物学，天文学等各个领域，开展了高维孤子的研究.1980年非线性效应专刊PhysicaD问世，与此同时，光纤中的孤子已在实验中产生出来.此后的发展更是突飞猛进。综上所述，孤立子理论的产生和发展是与近代物理密切相关的.孤立子理论不但包括了有关的数学理论，也包括了物理理论，数学的严密性和物理的启发性和实用性两者相互结合，相互依存，相互渗透，相互促进，使孤立子理论

9、显示出强大的生命力，这也是现代自然科学发展的重要特征之一。孤立子一词虽被广泛引用，但无一般性定义数学中，将孤立子理解为非线性偏微分方程的局部行波解，所谓局部是指微分方程的解在空间的无穷远处趋于零或确定常数的情况。换言之，孤立子指的是稳定的孤立波，即与同类孤波碰撞后不会消失，而且波形、波速和幅度不会改变或只有微弱改变的孤立波.在物理中，孤立子被理解为经典场方程的一个稳定的有限能量的不弥散的解，即能量集中在一个狭小的区域内且相互作用后不改变波形和波速。许多非线性发展方程，如KdV方程、SineGordon方程、Boussinesq方程、KP方程，Toda晶格方程等都具有孤立子解.孤立子除常见的钟型

10、和扭型外还有包络孤子、哨孤子、拓扑性孤子和非拓扑性孤子、呼吸子、亮孤子和暗孤子、正孤子和反孤子以及它们叠加而形成的形形色色的孤立子。现状求解微分方程是古老而在理论和实际上又很重要的研究课题，显示解，特别是行波解可以很好的描述各种物理现象，如振动、传播波等.但由于非线性微分方程的复杂性，至今仍有大量的重要方程无法求出精确解，即使己经求出精确解，也各有各的技巧，至今尚无一般的求解方法。所幸的是孤立子理论中蕴涵着一系列构造精确解的有效方法，如反散射法（1ST）、B?cklund变换法、Darboux变换法、Hirota双线性法、Painlev?有限展开法，延拓法及Lie群法等。随着各种求解方法的出现

11、，不但过去难以求解的方程得到解决，而且许多新的，具有重要物理意义的解不断被发现和利用。1967年，Gardrier等人发明了求解KdV方程的逆散射方法（也称为非线性），这一方法利用量子力学中的Schrodinger方程特征值问题（正散射问题）及其反问题（反散射问题）之间的关系，经过求解Gel'fand-Levitan-Marck-enko线性积分方程而给出KdV方程初值问题的解。它不仅对应用技术提供了崭新的方法和概念，而且对数.学自身的发展也有深远影响。随后，Lax将该方法加以综合和推广，使之能够用于求解其他非线性偏微分方程的初值问题，从而逐步形成一种系统的求解方法。1972年，Zak

12、harov和Shabat推广了这一方法，求出高阶KdV方程，立方Sehrodinger方程等的精确解。Ablowitz,Kaup,Newell和Segur则更加一般化反散射方法。李诩神、田畴、屠规章教授等也为发展反散射方法做了很好的工作。1971年，Hirota所引进的双线性变换法（Hirota方法），是构造非线性偏微分方程N一孤立子解及其Backlund变换的一种重要而直接的方法。1975年，Wahlquit和Estabrook提出延拓结构法，以外微分形式为工具，给出寻找与反散射方法相联系的线性特征值问题的系统的方法。1991年，李诩神教授基于对称约束提出一种非线性偏微分方程的直接的变量分离

13、方法；随后，楼森岳教授等提出另一种更有效的直接变量分离法得到了许多的（2+l）维非线性发展方程的精确解。精确求解非线性发展方程的工作具有重复性、固定的套路和规律、计算量大的特点，计算机代数的出现使人们摆脱了刻板、大量而重复的计算，提高了速度保证了准确率.1996年，Parkes和Duffy给出了求非线性发展方程孤立波解的双曲正切函数法的Mathematiea程序包。王明亮教授等基于非齐次项与高阶导数项平衡的原则，将非线性方程齐次化、代数化，提出了齐次平衡法。近年来提出并发展起来的齐次平衡方法，实际上是求非线性偏微分方程精确解的一种指导原则，故也称为齐次平衡原则。依据该原则，可事先判定某类非线性

14、偏微分方程是否有一定形式的精确解存在，如果回答是肯定的，则可按一定的步骤求出它来，并同时得到其满足某些条件的Backlund变换。因而齐次平衡原则具有直接、简洁、步骤分明的特点，再者，还适用于计算机的符号计算系统进行计算，且得到的是精确的结果.至今，齐次平衡原则在非线性数学物理中已得到广泛的应用，且其应用范围正在不断的扩展，己成为处理非线性数学物理相关问题的有效工具之一。所以，近年来在齐次平衡原则下又发展了多种求解非线性偏微分方程精确解的方法：像Tanh一函数法，SineCosine方法，Jacobi椭圆函数展开法，Riccati方程方法及F一展开法等。这些方法一般都借助于计算机代数系统（Ma

15、thematica或Maple）,求解方便、直接，而且可以对解进行数值模拟以便于直观分析解的性质。2、非线性偏微分方程的几种解法2.1逆算符法据逆算符方法的基本思想，把偏微分方程Au=A(u,ut,ux,uxx,.)=0改写为：Lu+Ru+Nu=0(2.1)其中L和R是线性微分算子，Nu是非线性项。算子L是可逆的，作用逆算子L,于上式两边得到u=f-L"(Ru)十L-1(Nu)(2.2)其中f满足(2.1)及初始条件，根据逆算符方法u可以分解为一系列分量之和-beu=£4(2.3)n-0利用回归关系可以得到11u。=f(x),u=-L(Rq)+L(Nuk)(2.4)非线性项

16、F(u)=Nu可以表示为无限级数之和-HeF(u)-A(2.5)n-0其中A是Adomian多项式，定义为1dn二iA=XL【FQ九ui)?x，n=0,12.(2.6)n!di=0利用(2.3)和(2.4)可以依次解出山，5，!,.，从而得到方程的解u=u0*,u2+.(2.7)业已证明Adomian分解法是收敛的，而且收敛速度相当快，能够得到精确解。齐次平衡法齐次平衡法是一种求解非线性偏微分方程非常重要的方法，它将非线性发展方程的求解问题转化为纯代数运算。利用这种方法不仅可以得到方程的Backlund变换，而且能得到非线性偏微分方程的新解.该方法的大致步骤如下：对于给定一个非线性偏微分方程(

17、2.8)这里P一般是其变元的多项式，其中含有非线性项及线性出现的最高阶偏导数项。一个函数w=w（x,t）称为是方程（2.8）的拟解，如果存在单变元函数f=f（w）,使得f（w）关于x,t的一些偏导数的适当的线性组合，即(2.9)式（犬/）=3+我（工工）（人）关于苏山的低于加+门阶的偏导数的适当线性组合）精确的满足（2.8）和（2.9）中的非负整数m,n,单位元函数f=f（w）以及函数w=w（x,t）都是待定的，将（2.9）代入（2.8）中可以通过以下步骤确定它们：首先，使高阶偏导数项中包含的w=w（x,t）的偏导数的最高哥次和非线性项中包含的关于w=w(x,t)的偏导数的最高哥次相等，来决定

18、非负整数m,n是否存在。其次，集合w=w(x,t)的偏导数的最高哥次的全部项，使其系数为零，而得f(w)满足的ODE,解之可得f=f(w),一般是对数函数。第三，将f(w)的各阶导数的非线性项，用f(w)的较高阶的导数来代替，再将f(w)的各阶导数项分别合并在一起，并令其系数为零，而得w=w(x,t)的各次齐次型的PDE组，可适当选择(2.9)中线性组合的系数，使PDE组有解。最后，若前三步的解答使肯定的，将这些结果代入(2.9),经过一些计算就得(2.8)的精确解。从(2.9)中可以看出，如果v(x,t)是方程(2.8)的一个解，则通过上述步骤就可以得到方程的Backlund变换。Jacob

19、i椭圆函数方法考虑非线性偏微分方程(2.8),寻求它的行波解为u=u(-),-=k(x-Ct)(2.10)其中k和c分别为波数和波速将u(D展开为下列Jacobi椭圆正弦函数sn-的级数:nu=£包水(2.11)j=0它的最高阶数为O(u(D)=n(2.12)因为du/.j_1.”=£jajsnJEcrtdn(2.13)dj=0其中cn-和dn分别为Jacobi椭圆余弦函数和第三种Jacobi椭圆函数，且cn2=1-sn2,dn2=1-m2sn2(2.14)m（0<m<1）为模数，且(2.15)sn=cndn,-d-cn：=一sndn,-dn=一m2sncndd

20、d由（2.13）式，可以认为业的最高阶数是dO(%)=n+1(2.16)类似地，有。喈)=2n1,O(喧2,O(2)=n3(2.17)在（2.11）式中选择n,使得非线性偏微分方程（2.8）中的非线性项和最高阶数项平衡，将（2.11）代入非线性偏微分方程（2.8）中，并利用（2.14）和（2.15）,可将方程（2.8）变成关于snY（i=0,1,.,N）的多项式.置sn化的各次届次的系数为零，得关于a0,a1,.,aN,k,c的代数方程组。解上述方程组，将结果代入（2.11）中，得（2.8）Jacobi椭圆函数解。应该指出的是，因为mT1时，snttanh：（2。11）式就退化为nu=

21、3;ajanh比(2.18)i=o所以此方法包含了双曲正切函数展开法。2.4辅助方程方法考虑非线性偏微分方程(2.8),寻求它的行波解为U=U(-),*=k(X-Ct)(2.19)其中k和c分别为波数和波速。将式(2.19)代入方程(2.8)中，则(2.8)化为U6)的非线性常微分方程(NODE):G(u,u,u化，.)=0(2.20)设u(d可表示为z)的有限哥级数：NU(-)aiZi(-)(2.21)=0这里的ai是待定参数，n为一常数，由非线性偏微分方程(2.8)中具有支配地位的非线性项和最高阶数项平衡得到，z代)满足如下新的辅助常微分方程(*)2=Az(D2+Bz(D4+CzR(2.2

22、2)d其中A,B,C为待定参数。将(2.21)代入NODE(2.20)中，利用(2.22)可将方程(2.20)左边变成z6)的多项式。令zK)的各哥次的系数为零，可得关于a0,a,.,aN,k,c的代数方程组。解上述方程组，可解得3o,a1,.,aN,k,c,将结果代入（2.21）中，得（2.8）的行波解的一般形式。利用表2.1,适当选取A,B,C,的值，可得方程（2.8）的一些特殊解。表2.14B,C/与方程（2.22）的解z（）之间的关系表讨15一般为任总实数）AcAxaA>Qc为任意实数A为任意实数sech(V7fX-;c，为任意实数A为任意实数C5Ch(A-7=jB14(7(1+

23、蠹coth(Vj?f)/>0c为任意实数A>0("_声A<0C为任意实数(-0e4Kcos(2A)-8AX)C为任意实数("声tV-A$inh(2Vf)-RC为任意实数A>0（-法E而皿2后用/?月>000A为任意实数sech(77)(j-：一)%B+2eACtanhJT。A<000A为任意实数5gc()(7：B+2c-JACtan(V-J)介。00A为任意实数csch(V7f)(7AO00A为任意实数csc(47f(1l)*B+2eJ-，Ccot(WIg)QO。为任意实数A=0（4。+£也由（¥§）声QOC

24、为任意实数A=0口却“8山1率勖成JX)C为任意实数为任意实数;M44.总部|(e2t-4B)2-64AC2.5F-展开法考虑非线性偏微分方程（2.8）,寻求它的行波解为u=u（1）,l=k（xCt）（2.23）其中k和c分别为波数和波速。将式（2.23）代入方程（2.8）中，则（2.8）化为u6）的非线性常微分方程（NODE）:P（u,u；u”,.）=0（2.24）设ud）可表示为F）有限哥级数：Nu（?）=a°+£aiF）（aN#0）（2.25）i1这里ai是待定常数，FG）满足下列一阶常微分方程：FU）=pF4+QF2+R（2.26）这里p,Q,R是待定常数，正整数N

25、是由具有支配地位的非线性项与最高阶偏导数项平衡确定。将（2.25）代入NODE（2.24）中，利用（2.26）可将方程（2.24）左边变成的多项式，置F（C）的各次哥次的系数为零，得关于a0,a,.,aN,k,c的代数方程组。求解上述方程组，可解得a0,a1,.,aN,k,c,将结果代入（2.25）中，得到（2.8）的行波解的一般形式。利用表2.2,适当选取p,Q,R的值，可得方程（2.8）的由Jacobi函数表示的周期波解。表2.2P,Q,J?与方程产气J）=pF'。尸+火的解F©之间的关系表PQR尸"(j)=pF+R£m-（】十方）F气g)=Q-F、l-/F')5时-m22m<lm2尸2=（1-尸）（苏尸+1-一）cn4-1mL-9穴6）=（1_尸乂产十疗-1）dnf1-(1+m2)mi926)=(1-7乂加一)蔗/3n0”1-2m<I-n?尸理g)=(1产*/_DF*_ncfcnO&