搜索
您的当前位置:首页正文

非线性波动方程的孤波解及余弦周期波解

来源:筏尚旅游网
维普资讯 http://www.cqvip.com 上海理工 大学 学报 第30卷第1期 J.University of Shanghai for Science and Technology Vo1.30 No.1 2008 文章编号:1007—6735(2008)O1—0015—07 非线性波动方程的孤波解及余弦周期波解 张卫国, 刘 刚, 任迎春 (上海理工大学理学院,上海200093) 摘要:利用假设待定法,求出了非线性波动方程的具有双曲正割函数分式形式且渐近值不为零的 精确孤波解和余弦函数周期波解,并分别讨论了它们的有界性,揭示了行波波速改变对钟状孤波解 与余弦函数周期波解波形变化的影响、 关键词:非线性波动方程;孤波解;余弦函数周期波解 中图分类号:O 175.2 文献标识码:A Solitary wave solutions and periodic cosine wave solutions of nonlinear wave equations ZHANG Wei-guo,LlU Gang,REN Ying.Chun (College ofScience,University ofShanghaifor Science and Technology,Shanghai 200093,China) Abstract:The exact solitary wave oslutions in the fractional form of hyperbolic secant function and the periodic wave solutions in the form of cosine function for nonlinear wave equations are obtained by use of undetermined assumption method.The bounded properties of these oslutions rae provided detailedly. The correlative characteristics between the exact solitary wave solutions and the peridoic cosine wave oslutions,as the traveling wave velocity varies,are revealed. Key words:nonlinear让尥勘B equation;solitary 勘B solution;periodic cosine zua72e soluti0 S 1问题的提出 Utt ““一m u+g U3=0 (4) 方程(3)和方程(4)分别是力学和粒子物理中的重要 非线性波动方程 模型. U“+忌“船+b1 U+b2U2+63“3=0 (1) 如果在方程(1)中对系 数适当取值,可以得到 是一维Klein—Gordon方程[1] Sine—Gordon方程[3,4] U +忌“ + (U)=0 (2) ““~Utt sin U (5) 当 (“)=b1 U+b2U +b3U 时的特殊情形. 的近似方程 适当调整方程(1)中的系数,可以得到Klein— 1 Gordon方程[2] Utt U (6) U 一“m+ “+ =0 (3) 以及Sinh—Gordon方程[3] 及Landou—Ginburg—Higgss方程[ ] U“ sinh U (7) 收稿日期:2007一O1—18 基金项目:上海市重点学科建设资助项目(T0502);上海市教委科研资助项目(07ZZ83) 作者简介:张卫国(1957一),男,教授. 维普资讯 http://www.cqvip.com 16 上海理工大学学报 2008年第30卷 的近似方程 Uxx一 = +吉 3 2方程(1)的精确孤波解 易知方程(1)的行波解U(z,t)=U( )=U(z 因此,研究非线性波动方程(1)的解对于很多物理问 —vt),满足 题具有重要意义. ( +k)U ( )+b x乱+b2U +b3U。=0 文献[5]给了方程(3)和方程(4)的型如 或 U(z,t)=atan[b(z+at+C)] U ( )+ ( )+SU ( )+hu。( )=0(9) U(z,t)=atanh[b(z+at+C)] 的解. = 文献[6]给了方程(3)的型如 本文求满足条件 u(x’ 、 —asin[b(z—at+c)] U ( ),U ( )一0,l l—oo COS 0  lz—a 十 l而C  且U( )的渐近值C =liar U( )满足 P±。。 u(x , 、 asi—nh[b(z—at+c)] cos h l D z一= 十 l_a=丽t C  b xz+b2x2+b3z3=0 的解,其中,口、b、c、a为常数. 的孤波解. 本文则给出了非线性波动方程(1)的型如 易知,方程(9)的解必然是方程(1)的解.因此, ( )= Ao+ A xsech2(a ̄)(可从方程(9)出发求非线性波动方程(1)的行波解. 8) 设方程(9)具有型如 的钟状孤波解,同时还获得了余弦周期波解,并讨论 Ao+A1sech2(a ) 了所求解的有界条件,最后还进一步揭示了波速的改 变对孤波解与余弦函数周期波解波形变化的影响. 的孤波解,则 2a (A1Bo—AoB1)[BlSech6(a ̄)一(2B1+3Bo)sech4(a ̄)+2Bosech2(a ̄)] , ———————— ——————一 将 ( ),乱 ( )代人方程(9),并比较sech(a )各次幂的系数,可得Ao,A1,Bo,B1和a 所满足的代数方 程组 B6+ Bo+hA8=0 ] 4a2B( o(ABiB1—o- AoB1)+Bo(AIB0+2BIAo)I+Ao(Ao B1+2AIBo)s+3A ̄Alh=02a 2( B1—A 11 Bo )()(2B2B1+1+ )3Bo )++ 1BI(AoB 2AI1++ )Bo l) AI++ +(AIBo 2Ao+ +B1)s 3Ao+ =A ̄一h:0 }(I 1。) 2a B1(BoA1一AoB1)+ 1B{+hAi+ {B1=0 j 用符号计算软件Maple解方程组(10)得两组解 a.A。:0,B。:一 1,B1:妒1, 2:一 l ) 其中,y是方程 6ly +4sy+3h=0 (12) . 2[hl (18hg一5s )y3一h/s(15h/~4s ) +h(18h l一34hls +8s ) 一3sh (7h/ 2s2)]^ 一A0 ———■ —一 Bo 2h(sy-6h)(y2+sy+h)==一—一—2/(3h/-s2)y3-s(15h/-4s2)yz-(18h——//-7s2)hy+3sh2A1A1 (13) B1= 1 2 2l(3h/一s ) 一s(15h/一4s )y2一(18M一7s ) +3sh 12hv( 一SV一3h) 维普资讯 http://www.cqvip.com

第1期 张卫国,等:非线性波动方程的孤波解及余弦周期波解 17 其中,Y是代数方程 l(9hl一2s ) ~2s(9hl一2s ) 。+ h(18hl一5s ) +6sh Y+9h。=0 的根. 4 一s2+s,/7,2—4hl 。4 ——— —一 (14) A阱: 3√S 一4hl A 记方程(12)的解为Y1,Y2,分别代人(11);记方 程(14)的解为Y3~Y6,分别代入方程(13),化简可 知方程组(10)有6组解. (18) Bo4:一 + ̄,),-4hlAo4 情形1当 I ̄y1=--一2s + ̄/4s2--18hl 时 B14=y4A14 情形5当Y取 一o, = (一十 j 4h/一s2一s 4—S2--—4h/ 口; ) —— B11= 钆,a 一丢 情形2当 :一一2s - ̄/4s2-18hl时 Ao5=—— 8h 2(4 s(s一3 ̄s2-4h/) +s) 15 Ao2=O, 一 z B12= 2A12,a =一 l 情形3当Y取 J l(16) ) 3、 (19) B05= ̄/7s2- 4h/-sA B15=ysA15 情形6当Y取 一s+,//TfsZ-4hl,(一十 一4hl 一 A 。 (・一 4 一s2一s ̄/s2—4 a ——— —一 ) 一、§ ———— ~4hl—s 七s : 5 S‘一4hL (17) Ao6= ——_A16 S) B。3:一S q-4S2,】,3√s2—4M -(20) 4hlA。3 Bo6= B1 =N3A1 4 S2--4hl-sA2l ∞∞ 情形4当Y取 B16= 6A16 (1_2 ̄/ s(s +3.. ̄s2..- 4 h1))…假定式(15)~(20)qu ̄a 2>o(j=1,2,…6)成 立,将其代人到式(8),可得非线性波动方程(1)的行 波解 34t-2bl1(z,t)= 2一——√—=2二6磊二sig二~n=(二9v二6 2二3+二6_ 。1 k)sec’hz 焉_=-.V/ +一— —一  k ) vz— Jt l【 √一十4_ ̄ b2 2 6s;gn 一 (9v623+61k))J1 。 hz 2 /r二 三 (+志 z一 ) 34  ̄blsgn(v2+k)一(21) 2(z,t)= 一一√2—=二6二i二一二二9二6二36二1_ 。 sech’z .丢焉V/ 十一— —一  k ) vz— t 2+[【 ・一,/ 2-^/b2s 6g;n一(v926+36k1 )jJ sech2 12 (22) 维普资讯 http://www.cqvip.com

18 上海理工大学学报 2008年第30卷 b2- ̄/ b2- 4b3blsgn(v2+k)一—. 62(3√b;一4b3b1sgn( +k)十b2) secha( 3(z—vt)) 1+ 26 (23) _『1+ 3二b2= (=3二 二厂二 二二=二 s==二二一二bg=n2二s(g= n二(二 +二2二忌+=)二忌+=)]二I62)-J sgn(v +忌) 2b3 6 secha(Ot3(z—vt)) 62一 Us4(z,t)=一 sech2(口4(z—vt)) 1— 2一√b2(3 Sgn( +忌)十b2) sech2(a4(z—vt)) (24) 『 一 3—b—2 (——3—√— =6=;==一二=4=6=3==6一=1—s—bg2—n—s(—g —n—( —+—2——k+—)—是—+—)]—b一z) JI sgn(v +忌) 2b3 6 其中 1 a3,4 厶 62+ Us5(z,t)=一 sech2(a5(z—vt)) (25) 1— 2一『 一 by/bz(2—sbgzn-( —3—2— ̄—+/。b。忌‘。 2。 —)-‘—+—4—3—b —‘3。b厂。 ’l。’。s。g=。n。 (。 v‘。 2。。+。。 k。。) —) 263 sech2(a5(z—vt)) Us6(z,t) b2+√b;一4b3blsgn(v +k) 6√b;一4b3b1  ̄/b2(b2-—3v—/b—2-—— ——4———b———3——。b。。。l。。s。g‘。。n。。‘。(‘ ‘v‘。。2。 。。+。。。。 k。。。)。——) b2sgn( 2+忌)+3 厂 _=_ sech2(a6(z—vt)) (26) sech2(a6(z—vt)) 1+ 2_f1+ 1 __-—— 其中 口5.6 2 在式(21)~(26)中 +k)≤0,2b;一9b1b3>0满足,则U 1( ,t)为有界 钟状孤波解;否则,Usl(z,t)非孤波; (b)若b2( +k)<0,b3( +k)>0,或62 ('02 q-k tI v1 Z+十k忌息><<U: 定理1 a.设b1( +忌)<0. (v2+k)≥0,2b;一9b1b3>0满足,则U 2(z,t)为 有界钟状孤波解;否则,Us2(z,t)非孤波. (a)若b2( +k)>0,b3( +k)>0,或b2( 维普资讯 http://www.cqvip.com 第1期 张卫国,等:非线性波动方程的孤波解及余弦周期波解 19 b.设b;一4b1b3>0,2b;一9b1b3≠0. (C)若b1( +k)>0,b2( +k)>0,b3( + √62(3√bi一4b3b1+b2)>3√bi一4b3b1—62 由于,当3√6i~4b3b1—62≤0时,上式恒成立. k)>0满足,则Us3(z,t)为有界钟状孤波解;否则, 3(z,t)非孤波; 当3 ̄/bi一4b3b1一b2>0时,上式等价于b2一 >0. (d)若b1( +k)>0,b2( +k)>0,b3≠0, 或b2( +k)<0,b3( +k)<0,2bi一9b1b3<0 由此即可推知,当b2>0,且b3b1>0时, 满足,则 4(z,t)为有界钟状孤波解;否则, 4(z, t)非孤波; (e)若b1( +k)>0,b2( +k)<0,b3≠0, 或b2( +k)>0,b3( +k)<0,2bi一9b1b3<0 满足,则Us5(z,t)为有界钟状孤波解;否则, 5(z, t)非孤波; (f)若b1( +k)>0,b2( +k)<0,b3( + k)>0满足,则Us6(z,t)为有界钟状孤波解;否则, Us6(z,t)非孤波. 证明现以对 3(.17,t)结论的证明为例来论 证说明,并且首先在 +k>0条件下考虑. a.为使Us3(z,t)有意义,参数必须满足下列两 条件: bi一4blb3>0,b2>0,或 62<0,且2bi一9blb3<0 b.考察a;>0的条件,因 2 4blb3一b{+b2、/bl一4b1 b3 —— ~ b3√b;一4b1b3(b2一√bi一4b1b3) 86;( +k) a;>0甘b3(b2一√b;一4b1b3)>0 所以,a;>0等价于下列3个条件之一成立: (a)63<0,b1>0; (b)b2>0,b3>0,b1>0; (c)b2<0,b3<0. c.考察Us3(z,t)分母不为零的条件. 由于seth2(a )∈[0,1],则Bo+B1seth2(a )取 值介于Bo和Bo+B1之间.即只需满足B0(B0+ B1)>0,就可以保证Bo+B1sech (a )不等于零.对 / ̄s3(X, )而言,只需研究满足Bo3(Bo3+B13)>0的 参数条件.由Bo3(Bo3+B13)表达式化简可得,B03 (Bo3+B13)>0等价于 2一lL  √b2(3^+ //bi一4b3b1+b2)j} >0 即等价于 3(z,t)分母不为零. 综合a~c的讨论,可得 当b1>0,b2>0,b3>0,bi一4b1b3>0时, Us3(z,t)有意义、有界. 又由于当l l=l z—vt l一十oo时 ㈦一一 厶c,3 所以, 3( )=Us3(z—vt)是非线性波动方程(1)的 有界钟状孤波解. 当 +k<0时,类似于上面的证明过程,可以 得到:当b1<0,b2<0,b3<0,b;一4b1b3>0时, Us3(z,t)有意义、有界,即 3( )= 3(.17一wt)是 非线性波动方程(1)的有界钟状孤波解. 3非线性波动方程(1)余弦函数周期 波解及其有界性证明 情形1~6中假定a >0( =1,2,…,6)时,得 到了钟状孤波解,若情形1~6中的a;<0( =1,2, …,6),则aj(j=1,2,…,6)为复数. a.从情形1、2可知,当l>0时,有解 A,B;0・,z A。・,z, B,0・,z 。・,z’ (I,2=Bll27) ,2, a1,2 z a 1,2 j 其中,a ±√寺. b.从情形3、4可知,当 h(4hl一 2+5~ )<0 时,有解 B 13,。 ∞4=B13 ,B,o3 ,4,a3,4= ,3.∞4 j }( 一 28)其中, =±√一一4hl-52+5 ̄s2-4h1. c.从情形5,6可知,当 h(4hl— 2—5~ )<0 时,有解 B s,6 Ao5,6, B,o5,6 o5,6 (5,6=Bls,6,a5,6= ,5,6 J 29) 维普资讯 http://www.cqvip.com 20 上海理工大学学报 2008年第30卷 其中, =±√一4一hl-s2-s ̄52-4hl 将式(27)~(29)代入式(8),有 u/t/( )= A +Aljsech2(c ̄¥) B +Bljsech2( ) A +AIjsech2( ) B +Bljsech2( ) Aoj+AIjsec2(口 ) Boj+BIjse 2 / ) Alj+Aojcos2(a ̄ ̄) B1J+Boj∞s2(口 ) 于是,非线性波动方程(1)有余弦函数周期波解 Upl,2(X,£)= 2∞S2(a ( 一 £))一(±3 e2+1) ) 其中 b asgn( +忌) 1一 ——?——————一 2bz一9b3b a b2sgn( +k) ●, 。‘—‘’:—’ 。’’ 。 。。。。。一 2bz一9b3b a ± ( ,£)=一—b2-y/b2- 4b 3blsgn(v2+k)—————一. ,4[1± ](31)4其中 . √6 一4b3b a ^1 —==二=========二二二二二二=二=二二二二= √62(62+3 sgn(732+ )) b2sgn(v + )一3 , √62(62+3√6;一46361sgn( + )) 一 一 ,  Ib;一4b3b1一b2√b;一4b3b1 ( +k)l L 863( + ) _J , 、 b2+ ̄/b;一4b3b1sgn( +k) 6( , ) 一————— ■——一‘ (32) 1 一— 二二二二=二二二二二=二===二二=二二==二二= √b2(b2—3 ̄/6;一4b3b1sgn( +k)) b2sgn( + )+3 2 /62(62—3 sgn( 2+ )) I b;一4b3b1+b2√b;一4b3b1 ( +k)l L 863( +k) _J 考虑使得上述 ( )(i=1,2,…,6)成立的参 数条件,有定理2. 定理2 a.假设b1( +k)>0,2b;一9b3b1> 0,b3( +k)>0,则 p1( ), p2( )为方程(1)的有 界余弦周期解; b.假设6;一4b3b1>0,b3( +k)>0,若条件 2b2—9b3b1<0,b2( +k)<0或b1( +k)<0, b2( +k)>0满足,则 p3( ), p4( )为方程(1) 的有界余弦周期解; c.假设6;一4b3b1>0,b3( +k)>0,若条件 2b2—9b3b1<0,b2( +k)>0或b1( +k)<0, b2( +k)<0满足,则 p5( ), p6( )为方程(1) 的有界余弦周期解. 证明此处证明类似于定理1的证明(略). 4讨论 为了方便起见,以方程(1)的系数满足条件 b1<0,b2≤0,2b;一9b1b3>0 (33) 为例来讨论. .1孤波解与周期波解 a.b3>0时,当 +k>0时,与定理1的b,定 理2的c相比较,发现有孤波“ 同时存在周 期波 p5, p6;当 +k<0,与定理1的d比较,存 在孤波 结论:方程(1)的系数满足式(33)及b3>0时, 有大速度( >一k)的孤波 1、Us2与大速度的周期 波 5、Up6共存; +k<0,即 <一k时,有小速 度的孤波 此时周期波消失. b.b <0时,同上讨论,得出结论:方程(1)的 系数满足式(33)及b3<0时,有大速度的孤波 1; 当 <一k时,有小速度的孤波Us3与小速度的周 期波 1、Up2共存.显然,此条件下当速度变小时出 现周期波. 维普资讯 http://www.cqvip.com

第1期 张卫国,等:非线性波动方程的孤波解及余弦周期波解 21 与 1关于 有如下连续演变关系: 4.2孤波与周期波关于’,的演变 当 2>一k时,行波表现为有界钟状孤波 若方程(1)满足 ( );当.0-2=一k时,本文解释为行波归于平静;当 b1<0,b2≤0,b3<0,2b2—9b3b1>0(34) z<一k时,行波表现为有界余弦函数周期波 由定理1的a,若 +k>0,Usl为有界钟状孤波解; Up1(). 由定理2的a,若 +k<0, l、Up2为有界周期波 图1为 1( )和 1( )波形随 变化的图形, 解.经分析得到方程(1)的系数满足式(34)的行波解 取k=一1,bl=一0.5,b2=一2,b3=~0.5. 一3.5 一5 一一2.5 I4 —3 一1.5 一l -0.5 O (a)当v=1.001时的U ({)图像 (b)当v=0.999时的Upl({)图像 -3.5 一0-2.5 6 一—4 0 一1.5 一2 0.5 f -f 、J (c)当 1.2时的“d({)图像 (d)当v:O 8时的 。({)图像 图1 Usl( )和Up】( )波形随波速变化的图形 Fig.1 Relation between Usl( )and Upl( )aS the wave velocity v ̄ying 根据定理1和定理2,在方程(1)的系数满足其  IM 1.London:Academic Prec Inc Ltd,1982. 他条件下,可得孤波解与周期波解共存的情况,并且 [3] D0oD R K,EILBECK J C,GIB =B0N J D,et a1.Solitons 可以得出波速改变对行波波形的影响. nad Nonlinear Wave Equations[M].London:Academic Pree Ine Ltd,1982. 参考文献: 『4] ABL0vITZ M J,SEGI爪H.Solitons and The Inverse [1]ELEONSKY V M.KULAGIN N E.Methods of asymp一 … Scattering Transform[M].Philiadelphia:SIAM,1981. n)tic expansi。n and。f finite dimensional mode1 in the the— L 3 J 王明亮,周宇斌.一个非线性波动方程的精确解[J].兰 ory of nonlinear waves[J].Advances in Nonlinear Waves, … 州大学学报(自然科学版),1996,32(1):1—5. 1985,2(1):21—29.L6 J 蔡红颖,郭冠平.一类非线性波动方程新的精确孤立波 [2]DOOD R K.Solitons and Nonlinear Wavas Equations 解[J].吉首大学学报(自然科学版),2002,23(2):56— 59. 

因篇幅问题不能全部显示,请点此查看更多更全内容

Top