非线性偏微分方程

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

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

3、解的研究成为非线性科学 中的热点。下面介绍一下孤立子理论的研究背景、研究现状。1.1 背景孤立子理论己经成为应用数学和数学物理的一个重要组成部分， 在流体力学，等离子物理，经典场论，量子论等领域有着广泛的应用。随着近代物理学和数学的发展，早在 1834 年由英国科学家 Russell 发现的孤立波现象近二十多年来引起了人们的极大关注，对这一现象的兴趣与日俱增.这是因为一方面孤立子具有粒子和波的许 多性能，在自然界中有一定的普遍性，利用孤立子理论也成功地解释 了许多物理上长期用经典理论未能解答的现象 ；另一方面，随着孤立 子物理问题的深入研究，孤立子的数学理论也应运而生，并已初步形 成比较完善的

4、理论体系。孤立子理论自1965年由Zabusky和Kruskal对孤立子（Soliton， 简称孤子）命名后得到了迅速地发展究其原因是孤波现象无所不在， 从天上涡旋星系的密度波，线，超流氦一3，超导JosePhson结，磁学，结构相变，液晶，流体动力学以及基本粒子等，都与孤子有关.其发展大致可分三个阶段：第一阶段，主要是在19世纪.最早讨论孤立子问题的是ScottRussell。1844年英国工程师Russell发现船在运河中快速行驶着， 当这条船突然停止时，在船头附近产生了一个光滑的、像小山包一样 的水波，然后这个水波离开船头保持它的形状和速度保持不变，接着这个水波的高度逐渐减少，最后在运河

5、的一个拐弯处消失掉，他把这 种水波称为孤立波，认为它就是流体运动的一个稳定解.直到1895年, 荷兰阿姆斯特丹大学的 Korteweg教授和他的学生devries才一成功 导出了著名KdV方程，求出了与Russell描述一致的即具有形状不变 的脉冲状的孤立波解，在理论上证实了孤立波的存在， 并对孤立波现 象作了较为完整的分析，解释了 Russell的浅水波，解决了这个问题。 他们的数学模型为(1.1)q 6uux uxxx孤立波解为:c2 vcu1 (x,t)sech ( (x-ct)2 2后人称1J1 (x,t)为1 一孤立子解，如果令二X - Ct，那么q在平面上的图为图1.1所示图1.1

6、光滑孤立子U1在-U平面上的图形1965年美国数学家Kruskal和abusky对KdV方程的孤立波解进行数学模拟，他们发现两个孤立波相撞之后，各自的运动方向和大小 形状都保持不变这种性质与物理中粒子的性质类似，因此他们称这 种孤立波为孤立子在通常情况下，人们把孤立波和孤立子混为一谈， 不把它们区别开来。与此同时，在1876 一 1882年发现的Backlund变换，成为后来发展孤子理论的重要基础。第二阶段大致可划在 1955一 1975年。1955年，Fermi，PastaUlam(FPU)将64个质点用非线性弹簧连成一条非线性振动弦，用计 算机计算了一维非线性晶格在各个振动模之间的转换。初

7、始时，这些谐振子的所有能量都集中在一个质点上，其他63个质点的初始能量 为零。按照经典的理论，只要非线性效应存在，就会有能量均分，各 态历经等现象出现， 即任何微弱的非线性相互作用， 可导致系统的非 平衡状态向平衡状态的过渡。 但实际计算的结果却与经典理论是背道 而驰.实际上，经过相当长时间之后，能量似乎又回到了原来的初始 分布，这就是著名的FPU问题。由于FPU问题是在频域空间考察的， 未能发现孤波解，因此该问题未能得到正确的解释。后来，人们发现 可以把晶体看成具有质量的弹簧拉成的链条， 这恰好是 Fermi 研究的 情况。 Toda 研究了这种模式的非线性振动，得到了孤波解，使 FPU 问

8、题得到正确的解答，从而进一步激发起人们对孤立波的研究兴趣。 1965年，zabusky和Kxusal对等离子体中孤立波的相互碰撞过程进行 计算机数值模拟， 进一步证实了孤立波在碰撞前后波形和速度保持不 变的论断，并且把它命名为孤立子(soliton)，它是指一大类非线性偏 微分方程的许多具有特殊性质的解， 以及具有相应的物理现象， 它的 性质具体为： (1)能量比较集中； (2)孤立子相互碰撞时具有弹性散射 现象。从此孤立子理论的研究工作得到了迅速发展。第三阶段 (1973 至今)，把孤子概念及理论广泛应用于物理学， 生 物学，天文学等各个领域，开展了高维孤子的研究 .1980 年非线性效 应

9、专刊 PhysicaD 问世，与此同时，光纤中的孤子已在实验中产生出 来.此后的发展更是突飞猛进。综上所述，孤立子理论的产生和发展是与近代物理密切相关的 . 孤立子理论不但包括了有关的数学理论， 也包括了物理理论， 数学的 严密性和物理的启发性和实用性两者相互结合， 相互依存，相互渗透，相互促进， 使孤立子理论显示出强大的生命力， 这也是现代自然科学 发展的重要特征之一。孤立子一词虽被广泛引用， 但无一般性定义数学中， 将孤立子理 解为非线性偏微分方程的局部行波解， 所谓局部是指微分方程的解在 空间的无穷远处趋于零或确定常数的情况。 换言之，孤立子指的是稳 定的孤立波，即与同类孤波碰撞后不会消

10、失，而且波形、波速和幅度 不会改变或只有微弱改变的孤立波 .在物理中，孤立子被理解为经典 场方程的一个稳定的有限能量的不弥散的解， 即能量集中在一个狭小 的区域内且相互作用后不改变波形和波速。 许多非线性发展方程， 如 KdV 方程、Sine Gordon 方程、Boussinesq 方程、KP 方程，Toda 晶格方程等都具有孤立子解 .孤立子除常见的钟型和扭型外还有包络 孤子、哨孤子、拓扑性孤子和非拓扑性孤子、呼吸子、亮孤子和暗孤 子、正孤子和反孤子以及它们叠加而形成的形形色色的孤立子。1.2 现状求解微分方程是古老而在理论和实际上又很重要的研究课题，显 示解，特别是行波解可以很好的描述各

11、种物理现象，如振动、传播波 等.但由于非线性微分方程的复杂性，至今仍有大量的重要方程无法 求出精确解，即使己经求出精确解，也各有各的技巧，至今尚无一般 的求解方法。所幸的是孤立子理论中蕴涵着一系列构造精确解的有效 方法，如反散射法 (IST)、 B?cklund 变换法、 Darboux 变换法、 Hirota 双线性法、Painlev?有限展开法，延拓法及Lie群法等。随着各种求 解方法的出现，不但过去难以求解的方程得到解决，而且许多新的， 具有重要物理意义的解不断被发现和利用。1967年，Gardrier等人发明了求解KdV方程的逆散射方法（也称 为非线性），这一方法利用量子力学中的Sch

12、rodinger方程特征值问题 （正散射问题 ）及其反问题 （反散射问题 ）之间的关系，经过求解 Gel' fand Levitan 一 Marck 一 enko线性积分方程而给出 KdV方程初值 问题的解。它不仅对应用技术提供了崭新的方法和概念，而且对数 . 学自身的发展也有深远影响。随后， Lax 将该方法加以综合和推广， 使之能够用于求解其他非线性偏微分方程的初值问题， 从而逐步形成 一种系统的求解方法。1972年，Zakharov和Shabat推广了这一方法， 求出高阶KdV方程，立方Sehrodinger方程等的精确解。Ablowitz， Kaup, Newell和Segur

13、则更加一般化反散射方法。李诩神、田畴、屠 规章教授等也为发展反散射方法做了很好的工作。1971 年, Hirota 所引进的双线性变换法 （Hirota 方法）,是构造非 线性偏微分方程 N 一孤立子解及其 Backlund 变换的一种重要而直接 的方法。1975年，Wahlquit和Estabrook提出延拓结构法，以外微分形式 为工具,给出寻找与反散射方法相联系的线性特征值问题的系统的方 法。1991 年,李诩神教授基于对称约束提出一种非线性偏微分方程 的直接的变量分离方法； 随后,楼森岳教授等提出另一种更有效的直 接变量分离法得到了许多的（2+1）维非线性发展方程的精确解。精确求解非线性

14、发展方程的工作具有重复性、固定的套路和规律、计算量大的特点，计算机代数的出现使人们摆脱了刻板、大量而 重复的计算，提高了速度保证了准确率.1996年，Parkes和Duffy给 出了求非线性发展方程孤立波解的双曲正切函数法的Mathematiea程序包。王明亮教授等基于非齐次项与高阶导数项平衡的原则， 将非线 性方程齐次化、代数化，提出了齐次平衡法。近年来提出并发展起来的齐次平衡方法， 实际上是求非线性偏微 分方程精确解的一种指导原则， 故也称为齐次平衡原则。 依据该原则， 可事先判定某类非线性偏微分方程是否有一定形式的精确解存在， 如 果回答是肯定的， 则可按一定的步骤求出它来， 并同时得到

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

16、可以对解进行数值模拟以便于直观分析解的性质。2、非线性偏微分方程的几种解法2.1逆算符法据逆算符方法的基本思想，把偏微分方程Au二A（u,q,Ux,Uxx,）二0 改写为：Lu Ru Nu = 0（2.1）其中L和R是线性微分算子，Nu是非线性项。算子L是可逆的，作 用逆算子L，于上式两边得到u = f - L1（Ru） L1（Nu）（2.2）其中f满足（2.1）及初始条件，根据逆算符方法 u可以分解为一系 列分量之和-bou 八 un（2.3）n =0利用回归关系可以得到u f（x）,uk-L_1（Ruk） L（Nuk）（2.4）非线性项F（u）= Nu可以表示为无限级数之和F（u）八 A（

17、2.5）n=0其中A是Adomian多项式，定义为1 dn人nFC 1）亠 n = 0,1,2（2.6）n! di =0利用（2.3）和（2.4）可以依次解出,口小2,.,从而得到方程的解u 二 u。比比.（2.7）业已证明Adomian分解法是收敛的，而且收敛速度相当快，能够得 到精确解。2.2齐次平衡法齐次平衡法是一种求解非线性偏微分方程非常重要的方法，它将非线性发展方程的求解问题转化为纯代数运算。利用这种方法不仅可以得到方程的Backlund变换，而且能得到非线性偏微分方程的新解. 该方法的大致步骤如下：对于给定一个非线性偏微分方程(2.8)这里P一般是其变元的多项式，其中含有非线性项及

18、线性出现的 最高阶偏导数项。一个函数w二w（x,t）称为是方程（2.8）的拟解，如果存在单变元 函数f = f（w），使得f （w）关于x,t的一些偏导数的适当的线性组合，即心沪嘿久畑）（如关于湎的低(2.9)于用+幵阶的偏导数的适当线性组合）精确的满足（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)的偏导数的最高幕次相等,来决定非负整数m,n是否存在其次，集合w = w(x,t)

19、的偏导数的最高幕次的全部项，使其系 数为零，而得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变换。2.3 Jacobi椭圆函数方法考虑非线

20、性偏微分方程(2.8)，寻求它的行波解为U = u(),二 k(x- Ct)(2.10)其中k和c分别为波数和波速将u()展开为下列Jacobi椭圆正弦函数sn的级数:nU 八 ajSnj(2.11)j =0它的最高阶数为0(u( ) = n (2.12)因为duTn=E jajSnJcdntj =0(2.13)其中cn和dn分别为Jacobi椭圆余弦函数和第三种 Jacobi椭圆函数，且cn2 =1-sn2 ,dn2 = 1 - m2sn2(2.14)m（ m 1）为模数，且dsnddd"cn=-sn dn ,dd_dn二-m2sn cn(2.15)由（2.13）式，可以认为 屯的

21、最高阶数是d -O（乎"n 1（ 2.16）类似地，有O(u汁2n 1,O%“2,O(兽)“3(2.17)在（2.11）式中选择n,使得非线性偏微分方程（2.8）中的非线性项 和最高阶数项平衡，将（2.11）代入非线性偏微分方程（2.8）中，并利用（2.14）和（2.15），可将方程（2.8）变成关于sn（>0,1,.,N）的多项式.置sH©的各次幕次的系数为零，得关于 a°,a1,.,aN,k,c的代数方程 组。解上述方程组，将结果代入（2.11 ）中，得（2.8） Jacobi椭圆函 数解。应该指出的是，因为 m，1时，sn > tani ，（2。

22、11）式 就退化为nu 八 aj tanhj（ 2.18）jq所以此方法包含了双曲正切函数展开法。2.4辅助方程方法考虑非线性偏微分方程（2.8）,寻求它的行波解为U 二 U（ ）二 k（X - Ct）（2.19）其中k和c分别为波数和波速。将式（2.19）代入方程（2.8）中，则（2.8）化为uC ）的非线性常 微分方程（NODE）:G（u,u ,u ,.） = 0（2.20）设u（）可表示为z（）的有限幕级数：Nu（）八 aiZi（ ）（2.21）i =0这里的ai是待定参数，n为一常数，由非线性偏微分方程（2.8）中具有 支配地位的非线性项和最高阶数项平衡得到，z（）满足如下新的辅助常微

23、分方程（乎）2 =Az（ ）2 Bz（ ）4 Cz（ ）6（2.22）d其中A,B,C为待定参数。将（2.21）代入NODE（ 2.20）中，利用（2.22）可将方程（2.20） 左边变成z（）的多项式。令z（）的各幕次的系数为零，可得关于 兔,耳,.8" ,k,c的代数方程组。解上述方程组，可解得兔，和，a”，k,c，将结果代入（2.21 ）中，得（2.8）的行波解的一般形式利用表2.1,适当选取A, B, C, 的值，可得方程（2.8）的一些特殊解表2.14 Br G d与方程（2.22）»解z（）之间的关系表=±1 5般为任意实数ACAxaA>QC为任

24、意.实数为任意实数sec h(風X j一C为任意实数A为任意实数«c h( Xj=声费-川(7(1 + £ cothW X)1亡为任意实数A>0（_ u强A<0C为任意实数A>0(高 0e4Kcos(2A-8AX)C为衽意实数A<Q(财半tV-A $inh(2Vf)- RQ为枉意实数A>0（驾=）片皿沁2匸亦）-/?00心为任意实数sec h( 77)(j:半B + 2eAC Tarth(Vf)A<000A为任童实数曲压）（严Wa00&为任意实数csch(V7f)(*一)再B + 2£y/-ACAOO0心为任意实数沁压)

25、(j=尸l)%B 十 2eV-AC cot(Vi)aC为任意实数A-0（亠凱 5曲（¥幻）怡QO（?为任意实数A=O（-令1北3仪¥幻必JX)C为任意实数也为任意实融I(e2t -4B)2-64AC2.5 F-展开法考虑非线性偏微分方程（2.8），寻求它的行波解为U 二 u（ ）二 k（x - Ct） （2.23）其中k和c分别为波数和波速。将式（2.23）代入方程（2.8）中，则（2.8）化为u）的非线性常微分方程（NODE）:P（u,u ,U ,.）二 0（2.24）设u（）可表示为F）有限幕级数：N u（）二 a。' ajF'（）何 7）（2.25）i

26、討这里Q是待定常数，F）满足下列一阶常微分方程：F 2（ pF4 QF2 R （2.26）这里p, Q，R是待定常数，正整数N是由具有支配地位的非线性项与最高阶偏导数项平衡确定。将（2.25）代入NODE（ 2.24）中，利用（2.26）可将方程（2.24）左边变成的多项式，置F）的各次幕次的系数为零，得关于ag,a1,.,aN ,k,c的代数方程组。求解上述方程组，可解得a°,a.,aN,k,c，将结果代入（2.25）中，得到（2.8）的行波解的一般形式。利I用表2.2,适当选取p，Q，R的值，可得方程（2.8）的由Jacobi 函数表示的周期波解。表2.2 P, 0J?与方程= pFQF1 +R的解F©之间的关系表QRF2()= pF* +歼 + RF©d m1铀占-m22m<lm2CT1$"-12-ffi2mL-严© = (1-戸赫-Ddnf1-(1+ m2)严 ) = (1-尸nsiHn 打 l1-IFt2 (f) = (1 一 尸)3? _ 1) F J 聊'！ncfcnO1