趙小山，靳文娟，王 璟

（天津職業(yè)技術(shù)師范大學(xué) 理學(xué)院，天津 300222）

隨著科學(xué)技術(shù)的發(fā)展，非線性科學(xué)在自然科學(xué)、社會(huì)科學(xué)領(lǐng)域中的應(yīng)用越來越廣泛，特別是尋找非線性發(fā)展方程的精確解在非線性的問題研究中占有重要地位。傳統(tǒng)的求解非線性波方程的方法主要有逆散射法[1]、Bucklund 法[2]、Darboux 變換法[3]、Hirota 雙線性法[4]和Painlevé展開法[5]等。近年來，涌現(xiàn)出一系列新的求解方法，如齊次平衡法[6]、雙曲正切函數(shù)展開法[7]、EXP函數(shù)法[8]和變分迭代法[9-11]等。本文利用何吉?dú)g教授[9]提出的變分迭代理論，對一類非線性偏微分方程組進(jìn)行求解，得到了方程的近似解，同時(shí)證明了該解法的快速收斂性。

1 方法簡述

本節(jié)主要簡述變分迭代法[8]的主要內(nèi)容。這一理論是對Inokuti提出的廣義拉氏乘子的改進(jìn)。在變分迭代法中考慮微分方程：

其中：L為線性算子；N為非線性算子；g（t）為非齊次項(xiàng)。用變分迭代法得式（1）的校正泛函：

其中：λ為廣義拉氏乘子；un為第n次近似解為限制變分，即n=0。

在該方法中，首先要確定拉氏乘子λ，λ可由變分理論識別。例如：所選乘子滿足校正泛函取駐值，即δun+1（t）=0；再通過任意初始函數(shù)u0及計(jì)算所得的拉氏乘子λ得到連續(xù)逼近解un，n≥0。若連續(xù)近似解序列收斂，則可以得到精確解。

2 方法的應(yīng)用

考慮如下方程組：

滿足：u（x，0）=ex，v（x，0）=e-x。將式（3）轉(zhuǎn)化為關(guān)于 u、v的變分形式，通過變分迭代法及式（2）對微分方程組（3）構(gòu)建t方向的校正泛函，得：

用變分理論中的限制變分的概念，將校正泛函式（4）重寫：

令δun+1（x，t）=0，δvn+1（x，t）=0，得：

由式（8）、式（9）求得拉氏乘子：

現(xiàn)將式（10）、式（11）分別代入校正泛函式（5），得下列迭代式：

令u0（x，t）=u（x，0）=ex，v0（x，t）=v（x，0）=e-x，并將其帶入式（12），得：

取定 t=0.1，x=0.5，得系統(tǒng)（3）精確解為：

表1給出了各階近似解的近似精度。

表1 系統(tǒng)（3）的近似解比較

由表1可以看出，每迭代一次，在確定其自變量取值后其精度不斷提高，即近似值不斷接近精確值，當(dāng)?shù)降谖宕螘r(shí)已經(jīng)很接近精確解。

3 結(jié)束語

本文通過采用變分迭代法，對非線性偏微分方程組進(jìn)行求解，得到了二元非線性偏微分方程組的近似解，證明了變分迭代法求解非線性系統(tǒng)的有效性。

［1］ABLOW ITZ M J，CLARKSON P A.Soliton，Nonlinear Evolution Equations and Inverse Scattering[M].Cambridge：Cambridge Univ Press，1991.

［2］谷超豪.孤立子理論及其應(yīng)用[M].杭州：浙江科技出版社，1990.

［3］MATVEEV V B，SALLEM A.Daroux Transformations and Solitons[M].Berlin：Springer，1991.

［4］HIROTA R.Exact solution of the Korteweg-de Vries equation for multiple collisions of solitons[J].Phys Rev Lett，1971，27：1192-1194.

［5］樓森岳.推廣的Painlevé展開及KdV方程的非標(biāo)準(zhǔn)截?cái)嘟鈁J].物理學(xué)報(bào)，1998，47：1739-1745.

［6］WANG M L，ZHOU Y B.Application of homogeneous balance method to exact solutions of nonlinear equations in mathematical physics[J].Phys Lett A，1996，213：67-75.

［7］李志斌，張善卿.非線性波方程準(zhǔn)確孤立波解的符號計(jì)算[J].數(shù)學(xué)物理學(xué)報(bào)，1997，17（1）：81-89.

［8］HE J H.New periodic solutions for nonlinear evolution equations using Exp-function method[J].Chaos Soliton Fract，2007（34）：1421-1429.

［9］HE J H.Some asymptotic methods for strongly nonlinear equations[J].Internat J Modern Phys B，2006，20（10）：1141-1199.

［10］ZHU J M，LU Z M，LIU Y U.Codoubly periodic wave solutions of Jaulient-Miodek equations using variational iteration method combined with Jacobianfunction method[J].Comun Theor Phys，2008（49）：1403-1406.

［11］CHEN X M，WANG L J.The variational iteration method for solving a neutral functional-differential equation with proportional delays[J].Computers and Mathematics with Applications，2010，59：2696-2702.