徐建中　莫嘉琪

摘要：利用廣義變分迭代方法討論了一類非線性強迫擾動Klein-Gordon方程.首先，用雙曲函數待定系數法求得了無擾動方程孤子波.其次，利用泛函變分迭代原理得到了強迫擾動Klein-Gordon方程的一個攝動近似解.最后，論述了解的一致有效性.得到的近似解是解析式，它可對近似解進行解析運算，這對用簡單的模擬方法得到的近似解是達不到的.

關鍵詞：攝動解;孤子波;變分迭代

中圖分類號：0157.29

文獻標志碼：A

DOI： 10.3969/j.issn.1000-5641. 2019.06.003

0 引言

非線性孤子波理論在地球物理、光學、力學、理論物理學等學科中有很重要的應用.例如在地理科學、海洋物理、散射波、燃燒理論等方面，許多學者作了多方面的研究，如Salathiel等利用Riccati方程映射方法研究了電晶格孤子行波[1]，Yu等研究了Bose-Einstein冷凝體的孤子波[2]、 Chow等研究了K-折疊松散的Ricci孤子波[3]，等等[4-9].非線性孤子波理論的定量和定性方法已有很多的改進，非線性孤子波解的一種研究方法是擾動理論的近似方法，主要是用擾動理論的近似展開式將非線性孤子波方程轉化為易求解的方程來求解，這樣就擺脫了對于模擬數值方法的依賴.此方法的優(yōu)點在于思路直接明了，計算簡單，可得到解的較高近似度.且求得的近似解具有解析性態(tài)，故不但能進行定量方面的分析，而且還能進行定性方面的分析.此方法使用面廣，具有較廣的研究前景，

作者等人利用廣義變分迭代、同倫映射、微分不等式和不動點理論等方法也研究了一些孤立波理論及有關的非線性問題[10-15].本文就是研究一類非線性高維擾動KG（Klein-Gordon）方程.對于一般較典型的KG方程已有許多研究，它代表的是許多自然現象的簡化情形.但此類方程已經不完全能滿足當前科學發(fā)展的要求，所以需要研究更能反映自然現象的廣義非線性高維強迫擾動KG方程.

1 非線性高維擾動KG方程

考慮如下廣義非線性擾動KG方程：其中a，m，k為正的常數，f為非線性強迫擾動項，它是關于其變量為足夠光滑的函數.

近年來，關于非線性KG方程解的研究大體集中在兩個方面.一方面是利用分析方法求出各種方程的精確解，如He等[16]利用推廣的F展開方法求得了某類發(fā)展方程的精確解，Zhang等[17]利用變形映射方法給出了方程的精確解;另一方面是定性地研究解的性態(tài)，如Zhang等[18]對KG方程給出波函數和能量方程.Teman[19]證明了一類KG方程整體吸引子的存在性.近來，典型的非線性KG方程還有許多研究，例如文獻[20-22]等，但一般的典型KG方程代表的是各種自然現象的精簡和濃縮，它不能滿足當前科學發(fā)展的需要，故有必要來研究更能表示真實自然現象的廣義擾動KG方程.顯然，復雜的非線性方程一般不能求出其精確解.本文提出了一類更一般的帶有非線性高維擾動KG方程（1）的求近似解方法，得到的這種近似解又可以繼續(xù)進行解析分析，進而能得到更深入的物理性態(tài)。

首先作行波變換：

不失一般性.上述所列出的孤子波行波解（II）、（12），在泛函變分迭代方法下，都可作為擾動KG方程的泛函迭代式的初始近似，本文是利用式（9）作為初始近似.即u0（z）=v1（z），并以此來求出對應的各次近似解.

2 強迫擾動KG方程解的變分迭代

為了求得非線性強迫擾動KG方程（3）具有較好精度的近似解，現采用泛函分析變分迭代方法[16].作泛函F：

取ε= 0.5，由式（23）、（24）得到非線性擾動KG方程（21）的孤子波解的一次微擾方程近似解ulapp （w）與模擬精確解的曲線比較，如圖2所示.

4 強迫擾動KG方程近似行波解的意義

本文選取的初始近似u0。是采用典型KG方程的孤子波解.它能較快求出對應于有擾動項情形下的KG方程在所要求精度范圍內的近似解.

由泛函變分迭代方法得到的非線性強迫擾動KG方程的孤子波近似解uapp（t，x，y，z）是近似的解析關系式，還可通過解析運算，譬如進行微分、積分等運算，繼續(xù)對非線性強迫擾動KG方程的孤子波解作進一步研究而得到其他相關的物理性態(tài).例如，可以通過解析式uapp（t，x，y，z）用微分運算計算出非線性擾動KG方程的孤子波解uapp（t，x，y，z）的關于x，y，z或t的變化率的分布情況，以此了解對應孤子波的相關特性.再如，由近似函數u app算出強迫擾動KG方程的孤子波解的波峰值、波谷值、拐點等.還可得到有關非線性強迫擾動KG方程在大氣、海洋物理、散射波、量子力學、燃燒理論等等的其他物理量的預測.而且還可采取適當措施，改變強迫擾動KG方程的非線性擾動項，以得到所要求的大氣、物理、地理等方面的性狀。

[參考文獻]

[1] SALATHIEL Y， AMADOU Y， GARMBO B G. et al. Soliton solutions and traveling wave solutions for adiscrete electrical lattice with nonlinear dispersion through the generalized Riccati equation mapping method[J]. Nonlinear Dynamics， 2017： 87（4）： 2435-2443.

[2] YU F J，LI L.Vector dark and bright soliton wave solutions and collisions for spin-l Bose-Einstein condensate[J]. Nonlinear Dynamics， 2017， 87（4）： 2697-2713.

[3] CHOW B， LU P. On K-noncollapsed complete noncompact shrinking gradient Ricci solitons which split at infinity[J]. Math Annalen， 2016， 366（3/4）：1195-1206.

[4] SAZONOV s V. New type of extremely short vector solitons in a medium of asymmetric molecules [J]. JETPLetters， 2015， 102 （12）： 834-840.

[5] 歐陽成，姚靜蓀，石蘭芳，等一類塵埃等離子體孤子波解[J].物理學報，2014， 63（11）：110203-110208.

[6] 石蘭芳，朱敏，周先春，等.一類非線性發(fā)展方程孤立子行波解[J].物理學報，2014， 63（13）： 130201-130205.

[7] OUYANG C， CHENG L H， MO J Q. Solving a class of burning disturbed problem with shock layers [J]. ChinPhys B： 2012， 21（5）： 15-18

[8] SHI L F，CHEN c s，ZHOU x C. The extended auxiliary equation method for the KdV equation with variablecoefficients [Jl. Chin Phys B， 2011， 20（10）： 803-811.

[9] 石蘭芳，周先春，一類擾動Burgers方程的孤子同倫映射解[J].物理學報，2010. 59（5）：2915-2918.

[10]MO J Q. Homotopic mapping solving method for gain fluency of a laser pulse amplifier [J]. Science in China G.2009， 39（7）： 1007-1010.

[11]MO J Q， LIN s R. The homotopic mapping solution for the solitary wave for a generalized nonlinear evolutionequation [J]. Chin Phys B， 2009， 18（9）： 3628-3631.

[12]MO J Q. Solution of travelling wave for nonlinear disturbed long-wave system [J]. Commun Theor Phys， 2011，55（3）： 387-390.

[13] MO J Q， CHEN x F. Homotopic mapping method of solitary wave solutions for generalized complex Burgersequation [J]. Chin Phys B， 2010， 19（10）： 20-23.

[14] XU J z，ZHOU z F. Existence and uniqueness of anti-periodic solutions for a kind of nonlinear nth-orderdifferential equation with multiple deviating arguments [J]. Ann Diff Eqn： 2012， 28（1）： 105-114.

[15] 徐建中，周宗福.一類四階具有多個偏差變元p-Laplacian中立型微分方程周期解的存在性[J].重慶工商大學學報，2012.29（11）： 9-16.

[16]HE J H， WU x H. Construction of solitary solution and compaction-like solution by variational iteration method[J]. Chaos Solitons& Fractals， 2006， 29（1）： 108-113.

[17] ZHANG x A， CHEN K， DUAN z I.Bound states of KG equation and Dirac for ring-shaped non-sphericaloscillator scalar and vector potential [J]. Chin Phys， 2005： 14（1）： 42-44

[18] ZHANG Q， YUE P，GONG L X. New exact solutions of nonlinear KG equation [J]. Chin Phys： 2006， 15（1）：35-38.

[19]TEMAN R. Infinite-dimensional Dynamical System in Mechanics and Physics [Ml. 2nd ed. New York： Springer，1997.

[20] 張建文，王旦霞，吳潤衡，一類廣義強阻尼Sine-Gordon方程的整體解[J].物理學報，2008， 57（4）：2021-2025.

[21] 韓兆秀.非線性KG方程新的精確解[J].物理學報，2005， 54（4）： 1481-1484.

[22] 盧殿臣，楊立波，洪寶劍，非線性KG方程新的精確解[J].江蘇大學學報，2010， 31（1）：120-124.

[23] CHANG K W. HOWES F A. Nonlinear Singular Perturbation Phenomena： Theory and Applications [M]. NewYork： Springer-Verlag， 1984.

[24]DE JAGER E M， JIANG F R. The Theory of Singular Perturbation [M]. Amsterdam： North-Holland PublishingCo. 1996.

[25] BARBU L，MOROSANU G. Singularly Perturbed Boundary-Value Problems [M]. Basel： Birkhauser， 2007.

（責任編輯：林磊）

收稿日期：2018-08-14

基金項目：國家自然科學基金（41275062）;安徽省高校自然科學研究重點項目（KJ2017A704，KJ2019A1303）;安徽省高校優(yōu)秀青年人才支持計劃項目（gxyq2018116）;安徽省優(yōu)秀教學團隊基金（2016jytd080）;毫州學院自然科學研究重點項目（BY22018803）

第一作者：徐建中，男，副教授，研究方向為應用數學、生態(tài)數學.E-mail： xujianzhongok@163.com.

通信作者：莫嘉琪，男，教授，研究方向為應用數學、生態(tài)數學、數學物理、工程數學.E-mail： mojiaqi@mail.ahnu.edu.cn.