一类非线性偏微分方程的改进的Jacobi椭圆函数精确解
非线性差分-微分方程的Jacobi椭圆函数解
![非线性差分-微分方程的Jacobi椭圆函数解](https://img.taocdn.com/s3/m/6a4b4914c5da50e2524d7f22.png)
c ns r tng t a e i v o uto o n lne r d s r t ys e s i o s d. By usng t e m e ho o t uc i r v lng wa e s l i nst on i a i c e e s t m spr po e i h t d,t he
中 图 分 类 号 :O 15 2 7 . 4
文 献标 识 码 :A
文 章 编 号 :10 —8 X (0 7 0 —0 10 0 19 8 20 ) 40 4 —5
J c b l p i u c in s l t n O t e n n i e r a o i l tcf n to o u i st h o ln a ei o
e a i n; e a ts uton qu to x c ol i s
近年来 ,随 着 离 散 孤 子 在生 物 系统 、原 子 链 、 固态 物理 、光 子结 构等 领域 的发 现 ,离散 非线 性 系
1 方 法 及 步骤
对 于给 定 的差分一 分方 程 微
△( 计p( ,), np ( ,) … , np ( ,) l £ U+ 2 £ , U+I £ ,
l t ie e ua i r bt i d a tc q ton a e o a ne .
Ke r s a o i lp i u ci n d sr t y wo d :J c b l tcf n t ; ic ee mKd ltie e u to ei o V a t q a in; ( c 2+ 1 _ i n in lHy rd ltie ) dme so a b i a t c
( ,), np ( , , U+ 2 z,) … 计 ( ,) … , , ( ,), £ ( £ …, z,), ( ,) 一 0 £ )_ ,
非线性振动、非线性波与Jacobi椭圆函数(续)
![非线性振动、非线性波与Jacobi椭圆函数(续)](https://img.taocdn.com/s3/m/ca4d1c00876fb84ae45c3b3567ec102de2bddf01.png)
教学研究非线性振动、非线性波与Ja cobi 椭圆函数(续)佘守宪(北京交通大学物理系,北京 100044)5 正弦戈登方程、扭结、反扭结波、呼吸孤立波正弦戈登方程(si ng -G ordon 方程,简记为S G方程)φxx -φtt =sin φ(38)是自然科学和应用科学中出现的另一个重要的非线性发展方程(演化方程).在晶体位错的传播、磁旋波在铁磁材料中的传播、约瑟夫森(Josephs on )结中继传输线等乃至生物物理的诸多问题中都应用了这一方程.本文从单摆弹簧链模型这一特例引入SG 方程.考虑如图3所示的N 个相连弹簧-单摆运动方程(称为单摆弹簧链).N 个单摆悬挂在横置的诸弹簧上,每个摆的长度和质量均相同,每个螺旋弹簧的长度和劲度系数也相同,则第i 个单摆在垂直挂弹簧的平面作扭转运动的运动微分方程可写为I d 2φd t 2=K φi -1-2φi +φi +1-F sin φ式中,I 为摆的转动惯量,K 是螺旋弹簧的扭转系数,F sin φ是由重力提供的恢复力矩,F 可视为常量.图3 单摆-弹簧链模型若把相邻两摆之间的距离记作Δx ,在上式两边同乘以1Δx2,则方程就变为1Δx φφ+φ+I/Δx K Δx ·52φ5=F/ΔxK Δx φ在N →∞,Δx ~的连续系统的极限情况下,上式左边第一项取连续近似,就可近似表示为52φ5x2,这样,按(K/F )1/2Δx 和(I/F )1/2的比例,分别取新的位置和时间坐标(x 和t ),就可以得到S G 方程(38),这里,si n φ给出非线性作用(如果没有这一项,SG 方程就变为线性的波动方程).下面,把SG 方程写成一般形式:u tt -c 20u xx +f 20si n u =0(39)用行波法求解,设ξ=x -ct ,u =u (ξ),把它化为c 2-c 2d 2u d ξ2+f 20si n u =0(40)取c 2>c 20c 2<c 20情况与之类似,从略,记作m 2=f 20c 2-c 2,则方程(40)化为与无阻尼单摆运动方程相同的数学形式:d 2u d ξ2+m 2si n u =0(41)仿照单摆的情况,我们先求出它的首次积分:12d u d ξ2+m 2(1-cos u )=H(42)其中,H 为积分常数,然后令:x =si nu2, k 2=H 2m 2(43)则方程(42)化为d x d ξ2=m 21-x 2k 2-x2(44)和单摆情况一样,若k 2<1,则可求得方程的解为椭圆正弦函数解:x =si nu2=±k sn m ξ-ξ0,k(45)其中:k =H 2m 2, m 2=f 20c 2-c 20这就是在c 2>c 20时SG 方程的周期解(椭圆正弦波解),其中ξ为积分常数[与式(5)及式()对应,这里的对应于单摆情况下的ω,而和ξ对应于单摆情况下的θ和]第23卷第2期 大 学 物 理 Vol.23No.22004年2月COLL EGE P HY S ICSFeb.2004©82i -1-2i i 1-t 2sin 00117m 0u t .当k →0时,式(45)化为sinu2=±k si n [m ξ-ξ0], (-π<u <π)这是S G 方程的线性波解.当k →1时(即H →2m 2时),式(45)化为双曲函数型的解:sin u2=±tanh m ξ-ξ0=tanh ±m ξ-ξ0(46)利用双曲正切函数的定义,式(46)很容易化为e ±mξ-ξ0=1+sin(u/2)1-sin (u/2)=1+tan(u/4)1-tan (u/4)=tan u 4+π4于是式(46)可以改写为u =-π+4arctan e±m ξ-ξ0这叫做S G 方程(39)的孤立波解(这是广义的孤立波,它包含冲击波),其中u +=-π+4arctan e mξ-ξ0(47)称为扭结波(kink waves ),也称为拓扑孤立子(topo 2logical s oliton ),见图4(a ),而u -=-π+4arct an e -m ξ-ξ0(48)称为反扭结波(anti 2kink waves ),也称为反孤立子(ant i 2soliton ),见图4(b ).图4 扭结波与反扭结波将式(46)对ξ求导数就得到d u d ξ=2m sech mξ-ξ0这也是一种孤立波的形式.SG 方程除了扭结孤立子解外,还有一种呼吸孤立子解.对S G 方程(39)作变量变换:t 1=f 0t ,x 1=f 0c 0x ,则可以把它化为u x 1x 1-u t 1t 1=sin u以下求方程u xx -u tt =sin u(49)的解这一方程的求解比较复杂所幸的是,受到扭结波、反扭结波的解的形式=±(x )即tanu4=e ±mx ·e �mct的启发,可以把tanu4用分离变量法分离成X (x )和T (t )两者之比来求解.另外,由三角学公式:si n u =2tan u21+t an 2u2, cos u =1-tan 2u21+tan 2u2可导出si n u =4tan u41-t an 2u41+t an2u42(50)这样,如果我们要求特解,可以设tanu4=T (t )X (x )(51)用分离变量法解.用微分法容易求出:u x x =4TX 2+T22X 2+T2X ″-2XX ′2(52)u tt =4XX 2+T22X 2+T2T ″-2T T ′2(53)将式(51)~(53)代入方程(49),就得到X 2+T2X ″X+T ″T -2X ′2+T ′2= T 2-X 2(54)容易看出,如果取X ′2=α1X 2+β1,T ′2=α2T 2+β2其中α1,α2,β1,β2为待定系数,则X ″/X =α1,T ″/T =α2方程(54)就化为 2α1X 2+β1+α2T 2+β2-X 2+T2α1+ α2=X 2-T 2于是,当α1-α2=1,β1+β2=0时,尝试解能满足方程.取α1=m 2,α2=m 2-1,其中m 2<1,并取β1=-n 2,β2=n 2,即得:X (x )=nmcosh m x +c 1T (x )=n1-m2sin1-m 2t +c 0故特解为=+x +4 大 学 物 理 第23卷©8..:u 4arctan em -ct u 4arct anmsi n1-m 2t c 01-m 2cosh m c 1亦即tanh u4=m sin1-m 2t +c 01-m 2c osh mx +c 1(55)这叫做呼吸子解,它描述一种脉动式的扰动,是一种强度作周期T =2π1-m 2的变化的波包,如图5所示,图中取c 0=c 1=0,m =018,它的时间周期是2π/1-m 2,在x 轴上方和下方周期性地变化,很像不断呼吸的样子,同时,又像是孤立子和反孤立子的一对振荡,故称呼吸子(breat her soli ton 或bion ),又称为双孤子(doublet ),它是一种定态波.图5 sine -G or d on 方程的呼吸孤立子 另一方面,如果设特解为u (x ,t )=4arctanX (x)T (t )(56)代入S G 方程,还可以得到描述两个沿相反方向运动的扭结孤立波和反扭结孤立波的迎头碰撞,是孤立波间相互作用的一个实例(事实上,各种孤立波都是由于孤立波之间的相互作用而出现“碰撞”的情况,但分析较为复杂,这个例子是比例简单易懂的一个).将特解式(56)代入后,仿前,就得到2X ′2+T ′2-X 2+T2X ″X +T ″T=X 2-T 2(57)它和式(54)一样,只是X 和T 彼此对换.取X ′2=α1X 2+β1,T ′2=α2T 2+β2,可知,当α1-α2=1,β1+β2=0时,能满足方程.若取α1=m 2,α2=m 2-1,并取β1=n 2,β2=-n 2,就得到:X (x)=nsi n mx ,T (t)=m c osh m 2-1t于是,写作=,就求得特解为u (x ,t )=4arctan c ·sinhx1-c 2c oshct1-c 2(58)这是1962年Perring 和Skyrme 用计算机数值计算实验后猜想出来的解,这里用试探解法证明了它确是精确解.我们注意到,u (x ,t )的渐近形式是:t →-∞时,t anu4∝c expx +ct1-c 2-c exp-x -ct1-c 2t →+∞时,tanu4∝-c exp-x +ct1-c 2+c expx -ct1-c 2这说明,特解式(58)给出了两个沿相反方向运动的扭结孤立波和反扭结孤立波的迎头碰撞,碰撞后形状和速度都保持不变,只有相位(可把指数函数前的因子c 吸收到指数上)改变,这样两个扭结、反扭结波的对撞如图6所示.图6 扭结波与反扭结波的迎头碰撞6 NLS 方程,椭圆函数波与包络型孤立子非线性薛定谔方程(NLS 方程),又称立方薛定谔方程,它是描写非线性波的调制(即非线性波包)方程,其一般形式为:i 5u 5t +α52u 5x 2+β|u |2u =0,i =-1(59)它和薛定谔方程的差别主要在于出现与波函数的立方成比例的非线性项β||,故名我们用在有色散并且折射率与光强有关的光学介质(例如平板光波导和光纤)中波包的传播为例,第2期 佘守宪:非线性振动、非线性波与Jacobi 椭圆函数(续)5 ©8mm 2-1m 11-c2u 2u .引入这一方程.设有一调幅波�E (x ,t )=E (x ,t )exp i k 0x -ω0t式中,k 0和ω0各为载波的波数和角频率,E (x ,t )为包络函数,表示场强E 的频谱在ω0附近有一定的宽度Δω0.可以把包络函数表示为E (x ,t )=∫A (ω)·expi k -k 0x -i ω-ω0t dω当有色散时,k ≈k 0+d kd ω0ω-ω0+12d 2kd ω2ω-ω02上式中,d k d ω0简记为k ′,是群速度v g 的倒数,d 2k d ω20简记为k ″,给出色散(正常色散k ″>0,反常色散k ″<0).若考虑到非线性效应,则介质中的波数k =n 0+n 2|E |2ωc,其中,n 0为线性折射率,n 2为K err 型(自聚焦)介质的非线性系数.因此,当色散效应与非线性效应同时存在时,有:k ≈k 0+k ′ω-ω0+12k ″ω-ω02+g |E |2g =n 2·2πλ容易证明,包络函数E (x ,t )满足的非线性偏微分方程为i55x +k ′55t E -12k ″52E 5t2+g |E |2E =0若取以群速v g =1/k ′随波包运动的动坐标系来描述,并令ε=Δω0ω0为脉冲宽度,则可把(x ,t )换成动坐标(ξ,τ),其中:ξ=ε2x , τ=εt -k ′x于是就得到i 5E 5ξ-12k ″52E 5τ2+g |E |2ε2E =0(60)在反常色散情况(k ″<0)下,可把这方程改写成量纲为一的形式:i 5u 5X +1252u 5T2+|u |2u =0(61)它的数学形式与非线性(立方)薛定谔方程(取�=m =1)5ψ5+5ψ5x+|ψ|ψ=完全相似,我们可把X 换成,T 换成x ,将方程(61)改写为i5u 5t +1252u 5x2+|u |2u =0(62)在正常色散情况(k ″>0)下,仿此,可得方程i 5u 5t -1252u 5x2+|u |2u =0(63)方程(62)和(63)分别对应于NLS 方程(58)中α=±12,β=1的情形.下面求一般形式的NLS 方程(59):i 5u 5t +α52u 5x2+β|u |2u =0的行波解,很有意思的是,这方程通常只有线性方程才具有的形式为u =Aexp i kx -ωt(64)的单色波解,其中A 、ω和k 分别是振幅、角频率和波数.把式(64)代入NLS 方程(59),立刻得到色散(频散)关系式:ω=αk 2-β|A |2(65)它说明,非线性波的色散关系既与波数k 有关,又与振幅A 有关.由上式得到群速度为v g =d ωd k=2αk(66)因为式(65)是根据式(64)得到的,所以式(65)可以看作是NLS 方程(59)的最低阶的近似,相应的式(64)是它的最低阶的近似解(把窄波包近似成单色波,其频率为波包中心的频率).我们来求包络波形式的解,即设解为u (x ,t)=φ(ξ)exp i (kx -ωt ),ξ=x -v g t 这相应于取以群速v g 沿x 正向运动的动坐标描述包络波.代入即得αd 2φd ξ2+i 2αk -v g d φd ξ+ω-αk 2φ+βφ3=0这里,因为v g =2αk ,所以第二项等于零,方程满足φ(ξ)是实函数的通常要求.记作ω-αk 2=-γ,(γ>0)(67)则方程简化为d 2φd ξ2=γαφ-βαφ3(68)下面分三种情况求解.1)α>0,β>0情况这相当于上面提到的在反常色散的K 型自聚焦介质中传播的非线性波的情形用导数公式可求出的二阶导数为6 大 学 物 理 第23卷©8i t 122220t err .dn ud 2d u2dn u =2-k 2dn u -2(dn u)3将此式与式(68)比较,可知我们可以取尝试解φ(ξ)=A dn C ξ-ξ0,k 其中,ξ0是由初始条件确定的积分常数,k 为模,A 、C 为待定常数.代入方程(68),即求出:A =±2γβ2-k 2,C =γα2-k 2故解为φ(ξ)=±2γβ2-k 2dnγα(2-k 2)ξ-ξ0,k(69) 当1>k >0时,给出的振幅函数φ为第三类Jac obi 椭圆函数的调幅波.当k =0时,dn 退化为1,给出的φ=±γβ=常数,这是平面余弦波(线性波).当k =1时,dn 退化为sech ,给出的包络函数是φ=±2γβsechγαξ-ξ这是称为包络型的孤立波(envelop soliton ),它以群速v g 向前传播,如图7所示.图7 包络型孤立波2)α<0,β>0情况由Jacobi 椭圆函数的导数公式,有d 2d u 2sn u =-1+k 2sn u +2k 2(sn u )3与NLS 方程(68)比较,可知我们可以取尝试解φ(ξ)=A sn C ξ-ξ0,k 代入方程(68),即得到:A =±k 2γβ1+k 2,C =-γα1+k 2故解为φ(ξ)=±2γβ+γα+ξξ,当→时,退化为,有φ(ξ)=±k2γβsi n-γαξ-ξ当k →1时,sn 退化为t anh ,有φ(ξ)=±γβt anh-γ2αξ-ξ0它也是孤立波解(但中心处振幅最小,离中心越远,振幅越大,是一种“暗孤波”).3)α>0,β<0情况由Jacobi 椭圆函数的导数公式有d 2d u2cs u =2-k2cs u +2k (cs u )3这里,cs u =cn usn u(与三角函数的余切函数c ot u =cos usi n u相对应).与NL S 方程(68)比较,可知我们可以取尝试解φ(ξ)=A cs C ξ-ξ0,k代入方程(68),即得到:A =±-2γβ2-k 2,C =γα2-k 2即解为φ=±-2γβ2-k 2csγα2-k 2ξ-ξ0,k当k →0时,cn →cos ,sn →sin ,故cs →cot ,包络解退化为φ=±-γβcotγ2αξ-ξ0这是线性调幅波.当k →1时,cn →sech ,sn →tanh ,故cs →csch ,包络解退化为φ=±-2γβcschγαξ-ξ这是孤立波解.总之,在上述三种情形下,当k 由k =0渐增到k =1时,振幅函数φ由线性波的三角函数解变为椭圆函数解,最后变为孤立波解.这里模k 标志非线性作用的强弱程度.参考文献:[1] Na y fe h A H ,Mook D T.Nonlinear Oscillations [M ].New Y or k :John -Wiley ,1979.Chapter 2,2.1;奈弗,穆克.非线性振动(上册)[M ].北京:高等教育出版社,8第章第1节(下转页)第2期 佘守宪:非线性振动、非线性波与Jacobi 椭圆函数(续)7 ©8k 1k 2sn-1k 2-0k k 0sn sin 190.221.24参考文献:[1] 葛隆祺,黄伟祥.抛体运动在体育运动应用[J ].物理教师,1998(4):33.[2] 许耀球,姚天白.田径运动生物力学[M ].北京:北京航空航天大学出版社,1990.250.[3] 漆安慎,杜婵英.力学[M ].北京:高等教育出版社,1997.75.[4] 普朗特L.流体力学概论[M ].北京:科学出版社,1981.582.Theoretic a na lysis an d simulat ive exper iment in t he bestt hro wing a ngle for the shot exerc isesXU Yue 2ming(Depart ment of Physics ,Z hejiang Education Institute ,Hangzhou ,Zhejiang 310012,China )Abstract :The i nfluence factors for t hrowing distance in t he shot exerci ly initial speed ,t hrowing hi gh ,t hrowing angle are analyzed t heoretically and t heir best com bi nation is obtai ned ,and t he result is proved by simulat ive experi ment.K ey w or ds :t hr owi ng ;t he best t hrowi ng angle ;simulative experiment(上接7页)[2] 刘式适,刘式达.特殊函数[M ].第2版.北京:气象出版社,2002.第11章,椭圆函数,620~715.[3] 刘式适,刘式达.物理学中的非线性方程[M ].北京:北京大学出版社,2000.第3章314、317节,第6章613、614、6112节.[4] 佘守宪.孤立子及其相互作用的初等分析[J ].大学物理,1996,15(11):41~44;1996,15(12):35~39.[5] Drazin P G,Johns on R S.Solitons:An Introduction[M ].Ca mbridge :Ca mbridge Univer s ity ,1989.Chap 2ters 1,2,1~38.[6] Re moissenet M.Waves Called Solitons [M ].SecondEdition.Berlin Heidelberg ,New Y ork :S pringer 2Ver 2lag ,1996.[7] 周凌云,王瑞丽,吴敏等.非线性物理理论及应用[M ].北京:科学出版社,2000.第一章、第二章.N onlinear vibra tions ,nonlinear w a ves andJacobian elliptic f unctionsSHE Shou 2xian(De part ment of Physics ,Beijing Jiaotong University ,Beijing 100044,China)Abstract :Succi nct and efficient met hod to obtain analytic solutions of nonli near vi brations and nonlinear waves by Jacobian elliptic functions are i ntroduced.Important t ypical examples are given and explained ,i nclud 2ing si mple pendulum ,Duffing oscillator ,cnoi dal wave and solitary wave solutions of K dV equation ,sine -G or 2don equation ,nonli near Schr �di nger equat ion ,sech 2profile soli tons ,kink and anti -kink solitons ,breat her ,i n 2teraction of a kink and an anti -ki nk ,and envelop solitons.K ey w or ds :nonlinear vibration ;simple pendulum ;Duffi ng oscillator ;solit ary waves ;Jacobian elliptic f unc 2tion24 大 学 物 理 第23卷©8。
一类非线性薛定谔方程的Jacobi椭圆函数解
![一类非线性薛定谔方程的Jacobi椭圆函数解](https://img.taocdn.com/s3/m/ccf2115c2e3f5727a5e962d1.png)
R=2 ,m=0.5,k=1, =12 ,k1=1,c1=2 。
1.0
1.0
令方程(14)中各项的系数为零,得到超定方程组
0.5
0.5
v
0.0
u
0.0
2 Pb1m 2 + 3Ra12b1 - Rb13 = 0 - Pb1 - Pb1m
2
-0.5
-0.5
-1.0 -80
-60
-40
-20
p
p
)的最高阶数分别为 (6) (7)
O[u (d u/d
)]=(q+1)n+p,q=1,2,3…
双曲函
利用非线性波方程(1)中的非线性项和最高阶导数项 平衡,可定出 n,将(4)式代入(3)式,由函数的无关性 可得一代数方程组,通过吴消元法解方程组可得解,而后将 其代入(4)式可得原方程的解。本文就利用此法求解一类 非线性薛定谔方程, 从而得到 Jacobi 椭圆函数解, 当模 m → 1 时,退化为相应的孤立波解。 2 求解非线性薛定谔方程 该方程形式如下 i
25
Qb0 - Rb03 - 3Rb0 a12 = 0
-6 Rb0b1a1 = 0
1.0
图 2 (a) 、 (b)分别是 x=0 时 v1( )和
u1 = v1 ( ) sin( kx - t ) 的图像,步长 t=0.1
1.0
采用消元法[14]求解方程组(13),得到以下几种情况 (1)b0=0, a1 = 0 , b1 = ± 2 Pm / C 此时条件为
2
0.5
0.5
v
0.0
u
0.0
(16)
-0.5
-0.5
Q=P+Pm2,
一种改进的基于Jacobi椭圆函数的随机平均法
![一种改进的基于Jacobi椭圆函数的随机平均法](https://img.taocdn.com/s3/m/bc94d99c48649b6648d7c1c708a1284ac850050c.png)
一种改进的基于Jacobi椭圆函数的随机平均法【摘要】本文介绍了一种基于Jacobi椭圆函数的改进随机平均法,并对其进行了理论分析和模拟实验。
首先对Jacobi椭圆函数进行了概述,然后比较了传统的随机平均法和改进方法的差异。
通过模拟实验结果可以看出,基于Jacobi椭圆函数的方法在一定程度上提高了计算的准确性和效率。
进一步的理论分析支撑了这一结论。
在我们验证了改进方法的有效性,并提出了未来研究方向。
通过本文的研究,我们可以更深入地理解Jacobi椭圆函数在随机平均法中的应用,为相关领域的研究提供了新的思路和方法。
【关键词】Jacobi椭圆函数、随机平均法、改进、模拟实验、理论分析、有效性、未来研究方向、总结。
1. 引言1.1 研究背景在现代科学和工程领域,随机平均法是一种常用的数值计算方法,用于解决复杂问题的数值求解。
传统的随机平均法通常基于传统的数学函数,如正弦、余弦函数等,来近似描述问题的解。
这些传统的数学函数在某些特定问题中可能存在局限性,无法有效地描述问题的复杂特性。
为了克服传统随机平均法的局限性,基于Jacobi椭圆函数的随机平均法应运而生。
Jacobi椭圆函数是一种复杂的特殊函数,具有更加丰富的特性和更广泛的适用性。
通过利用Jacobi椭圆函数的特性,可以更准确地描述问题的解,从而提高求解的精度和效率。
本研究旨在探究一种改进的基于Jacobi椭圆函数的随机平均法,以解决复杂问题的数值求解。
通过引入Jacobi椭圆函数的特性,我们可以提高数值计算的准确性和稳定性,同时拓展随机平均法的适用范围。
这将为科学研究和工程实践提供新的方法和思路,推动数值计算方法的发展和应用。
1.2 研究目的研究目的:本研究旨在通过改进基于Jacobi椭圆函数的随机平均法,提高计算效率和准确性,以解决传统方法在处理复杂计算问题时存在的局限性。
通过深入探讨Jacobi椭圆函数的特性,结合随机平均法的原理和优缺点,我们希望能够从理论和实际应用角度出发,提出一种更加有效的算法。
一类非线性偏微分方程精确解的表达
![一类非线性偏微分方程精确解的表达](https://img.taocdn.com/s3/m/77a1d202bc64783e0912a21614791711cc79797f.png)
一类非线性偏微分方程精确解的表达近年来,随着科学技术的发展,人们越来越重视非线性偏微分方程在解决数学、物理和工程问题等方面的应用。
因此,有关非线性偏微分方程的精确解的研究也受到了广泛的关注。
一般来说,由于非线性偏微分方程本身的复杂性,很难求解出精确解。
幸运的是,研究人员发展了各种推导一类非线性偏微分方程精确解的方法,提出了有效的数学算法,大大简化了解非线性偏微分方程的过程。
首先,要求出一类非线性偏微分方程的精确解,需要建立起正确的方程模型。
可以通过多种方法获得正确的方程模型,其中最常见的就是参数估计法,它能够根据实际观测数据来估计适当的参数值。
一旦确定一类非线性偏微分方程的精确解,就可以采用替代求解方法,如拉格朗日法等,根据不同情况使用不同拉格朗日多项式把方程化为某种可以求解的形式,从而得到精确解。
此外,研究人员还提出了一种新的求解方法展开式求解方法,它能够通过把非线性微分方程表达为一系列无穷级数,来求解出精确解。
这种方法使用级数展开把原始方程表达为展开后的连续微分方程,再对连续微分方程进行复分,最后根据分析关系和条件得出精确解。
最后,可以采用迭代法来求解一类非线性偏微分方程的精确解。
迭代法的基本思想是:设定一个初始猜想值,然后不断迭代,每次迭代后,就可以得到越来越接近正确解的新猜想值,最终可以求得精确的解。
总之,本文从理论上讨论了求解一类非线性偏微分方程精确解的方法,介绍了参数估计法、拉格朗日方法、展开式求解方法以及迭代法作为求解一类非线性偏微分方程精确解的主要方法。
这些方法是由研究人员为了进一步研究非线性偏微分方程而提出的,能够有效应对一类非线性偏微分方程的求解,为更复杂的问题的分析与解决提供了有效的数学方法支持。
综上所述,我们可以看出,求解一类非线性偏微分方程的精确解需要在正确的方程模型基础上,使用参数估计法、拉格朗日方法、展开式求解方法以及迭代法等方法,以提高求解准确性和可靠性,为数学、物理和工程等问题的求解提供关键支持。
用一种扩展的jacobi椭圆函数展开法求解mkdv方程
![用一种扩展的jacobi椭圆函数展开法求解mkdv方程](https://img.taocdn.com/s3/m/e7768342001ca300a6c30c22590102020740f23e.png)
用一种扩展的jacobi椭圆函数展开法求解mkdv方程MKDV方程是一种二阶非线性反应扩散方程,它是由Kadomtsev-Petviashvili方程的非线性化而形成的。
它的物理意义是模拟水平和垂直的海面波动。
近年来,由于不同的物理现象包括流体力学和电磁学,MKDV方程受到了广泛的关注和研究,被用来模拟它们的潮汐现象。
本文将介绍一种新的数值方法,即扩展的Jacobian 椭圆函数展开法,用于求解MKDV方程。
扩展的Jacobian椭圆函数Jacobian椭圆函数是一种全局收敛性很强的近似积分方法,可以有效地解决非线性问题。
在Jacobian椭圆函数展开法中,将解的形式表示为多项式的展开式,其中系数由Jacobian椭圆函数的积分结果计算得出。
但是,此方法在求解MKDV方程时存在一些不足之处。
为了弥补此缺陷,我们提出了一种扩展的Jacobian椭圆函数展开方法。
扩展的Jacobian椭圆函数展开法在这种扩展的Jacobian椭圆函数展开法中,我们先将解的形式进行拆分,用一系列由Jacobian椭圆函数展开而成的子解替换原来的解。
这样,就可以将非线性问题转换为一系列更简单的线性问题,并且可以采用插值的方式来解决在时间和空间维度上的不连续性问题。
此外,我们还可以利用多项式拟合的方法,来求解MKDV方程的近似解。
模拟结果为了验证扩展的Jacobian椭圆函数展开法在求解MKDV方程方面的有效性,我们对此进行了模拟研究。
结果表明,该方法能够很好地求解MKDV方程,并且具有良好的精确性和可靠性。
结论本文提出了一种新的扩展的Jacobian椭圆函数展开方法,用于求解MKDV方程。
该方法利用多项式拟合的方法,能够对问题的解进行有效的拆分,并且具有良好的收敛性。
我们的模拟研究结果表明,扩展的Jacobian椭圆函数展开法能够高效地求解MKDV方程,有效率地提高解的精确度。
总结本文介绍了一种新的数值方法,即扩展的Jacobian椭圆函数展开法,用于求解MKDV方程。
扩展的Jacobi椭圆函数展开法
![扩展的Jacobi椭圆函数展开法](https://img.taocdn.com/s3/m/5ecea725773231126edb6f1aff00bed5b9f3737c.png)
㊀㊀㊀㊀㊀138㊀扩展的Jacobi椭圆函数展开法扩展的Jacobi椭圆函数展开法Һ王小艳㊀(太原科技大学晋城校区,山西㊀晋城㊀048000)㊀㊀ʌ摘要ɔ求非线性偏微分方程的精确解有很多有效的方法,本文介绍了扩展的acobi椭圆函数展开法,并用这种方法求解了BBM方程的精确解,且在极限形式下,这些解退化为方程的孤波解和三角函数解.ʌ关键词ɔJacobi椭圆函数展开法;精确解;BBM方程在本文中,我们将Jacobi椭圆函数展开法及F-展开法相结合,得到了扩展的Jacobi椭圆函数展开法.一㊁方法介绍用扩展的Jacobi椭圆函数展开法求解非线性偏微分方程的主要步骤为:设非线性发展方程的一般形式Pu,ut,ux,uxx,utt,uxt= ()=0(1.1)其中P是关于u,ux,ut,uxt,uxx,uxt, 的多项式.第一步:在行波变换u(x,t)=u(ξ),ξ=αx+kt下,(1.1)式约化为如下的常微分方程:Q(u,uξ,uξξ,uξξξ, )=0(1.2)第二步:假设(1.2)式中的纯量函数u可展开为:u=u0+pu1+p2u2(1.3)则(1.3)式的二次导数为:uξξ=u0ξξ+pu1ξξ+p2u2ξξ(1.4)其中p(0<p<1)是一个比较小的参数u0,u1,u2分别是零阶㊁一阶和二阶拟设解.第三步:把(1.3)和(1.4)式代入(1.2)式,并令p0,p1和p2的系数为0,得到方程组,解此方程组就能得到(1.1)式的零阶㊁一阶和二阶解.第四步:假设零阶拟设解u0的Jacobi椭圆展开函数为:u0=ðni=0aiFi(ξ)其中ai均为待定系数,F(ξ)满足方程Fᶄ2=P0+P2F2+P4F4(1.5)由(1.5)式得:FᶄFᵡ=P2FFᶄ+2P4F3Fᶄ(1.6)Fᵡ=P2F+2P4F3(1.7)Fm=P2Fᶄ+6P4F2Fᶄ(1.8)把(1.6)㊁(1.7)式代入p0的系数=0,利用Matlab将同幂级的F(ξ)合并在一起,提取多项式系数,并令多项式系数为零,得到一些联立的代数方程.解代数方程组,得到ai,从而得到零阶解u0.第五步:同理,假设一阶拟设解u1的Jacobi椭圆正弦函数为:u1=b0FG(1.9)其中bi均为待定系数,F(ξ)G(ξ)均满足(1.5)式,则可设:Gᶄ2=R0+R2G2+R4G4(1.10)且满足下列恒等式:G2=ρF2+n(1.11)则有恒等式P4=ρR43nR4=P2-R2(1.12)把(1.9)㊁(1.12)式代入方程p1的系数=0,利用软件Matlab同幂级的Fi(ξ)Gj(ξ)合并在一起,提取系数多项式,并令多项式系数为零,得到一些联立的代数方程.解代数方程组可得到G(ξ),从而得到零阶解u1.第六步:将u0和u1代入方程p2的系数=0,假设u2=ðni=0ciFi(ξ)(1.13)其中ci均为待定系数.将(1.13)式代入p2的系数=0,利用Matlab将同幂级的F(ξ)合并在一起,提取多项式系数,并令多项式系数为零,得到一些联立的代数方程.解代数方程组,得到ci,从而得到二阶拟设解u2.第七步:将不同的P0,P2,P4及其所对应的F(ξ)值代入所得的u0,u1,u2可以得到多组精确解.这样我们就得到了(1.1)式的扩展的Jacobi椭圆函数展开解.二㊁用扩展的Jacobi椭圆函数展开法解BBM方程(1+1)维Benjamin⁃Bona⁃Mahony(BBM)方程的表达式如下:λut+μux+γuux+wuxxx=0,(2.1)其中下标表示x和t的偏导数,λ,γ,μ和w是常数,且u是一个实纯量函数u(x,t).在这一部分我们应用Jacobi椭圆函数展开法来求解BBM方程,具体步骤如下:方程的移动类波变量ξ为:ξ=αx+kt,其中α和k是常数,则(2.1)式变成以下形式:kλuξ+αμuξ+αγuuξ+α3wuξξξ=0(2.2)在(2.2)式中对ξ求积分并整合得:c+kλu+αμu+αγu22+α3wuξξ=0(2.3)其中c是积分常数.一般地,纯量函数u可展开为:u=u0+pu1+p2u2(2.4)(2.4)式的二次导数为:uξξ=u0ξξ+pu1ξξ+p2u2ξξ(2.5)其中p(0<p<1)是一个比较小的参数u0,u1,u2分别是零阶㊁一阶和二阶拟设解.㊀㊀㊀139㊀㊀把(2.4)和(2.5)式代入(2.3)式,并令p0,p1和p2的系数为0,则有:p0的系数方程为:2c+2(kλ+αμ)u0+αγu20+2α3wu0ξξ=0(2.6)p1的系数方程为:(kλ+αμ)u1+αγu0u1+α3wu1ξξ=0(2.7)p2的系数方程为:2(kλ+αμ)u2+αγu21+2αγu0u2+2α3wu2ξξ=0(2.8)零阶拟设解u0的Jacobi椭圆正弦函数为:u0=a0+a1sn(ξ,m)+a2sn2(ξ,m)(2.9)其中a0,a1和a2都是常数.把(2.9)式代入(2.6)式,同幂级的sn(ξ,m)合并在一起,并令多项式系数为零,得到以下一系列方程:sn0(ξ,m)的系数方程为:μα+kλ+αγa02æèçöø÷a0+2wα3a2+c=0sn1(ξ,m)的系数方程为:(μα+kλ-wα3(1+m2)+αja0)a1=0sn2(ξ,m)的系数方程为:(μα+kλ-wα3(1+m2)+αγa0)a2+αγ2a21=0sn3(ξ,m)的系数方程为:2wα3m2a1+α1γa1a2=0sn4(ξ,m)的系数方程为:6wα3m2a2+αγ2a22=0解方程组,我们可以得到:a0=4wα31+m2()-(kλ+αμ)rαa1=0a2=-12wm2α2r㊀㊀㊀㊀(2.10)零阶解u0为如下形式:u0=a0+a2sn2(ξ2m)(2.11)把(2.11)式代入(2.7)式,得:u1ξξ+4+4m2-12m2sn2(ξ,m)()u1=0(2.12)由Lamé方程的相关性质得u1(ξ)=g0sn(ξ,m)cn(ξ,m)dn(ξ,m)(2.13)其中g0为常数.利用(2.11)和(2.13)式,(2.8)式可写为:wα3u2ξξ+wα2(4+4m2-12m2sn2(ξ,m))u2+αγ2g20sn2(ξ,m)cn2(ξ,m)dn2(ξ,m)=0(2.14)令u2=b0+b1sn2(ξ,m)+b2sn4(ξ,m)(2.15)其中b0,b1和b2是待定系数.将(2.15)式代入(2.14)式,把同幂级的sn(ξ,m)合并在一起,并令多项式系数为零,得到以下一系列方程:sn0(ξ,m)的系数方程为:2wα3b1+(4+4m2)wα3b0=0sn2(ξ,m)的系数方程为:12wα3b2-12wα3m2b0+αγ2g20=0sn4(ξ,m)的系数方程为:6wα3b1+12wα3m2b2+αγ2g201+m2()=0sn6(ξ,m)的系数方程为:8wα3b2+αγ2g20m2=0解代数方程组,可以得到:b0=-γαg2048wα3m2b1=γαg201+m2()24wα3m2b2=-γαg2016wα3综上,(2.1)式的解为:u(x,t)=a0+a2sn2(ξ,m)+pg0sn(ξ,m)cn(ξ,m)dn(ξ,m)+p2[b0+b1sn2(ξ,m)+b2sn4(ξ,m)]当mң1时,u(x,t)=a0+a2tanh2ξ+pg0tanhβcen2ξ+p2[b0+b1tanh2ξ+b2tanh4ξ]当mң0时,u(x,t)=a0+a2sin2ξ+pg0sinξcosξ+p2∣b0+b1sin2ξ+b2sin4ξ]ʌ参考文献ɔ[1]郭玉翠.非线性偏微分方程引论[M],清华大学出版社,2008:238-245.[2]CesarA,GómezS,AlvaroH,et.al.NewperiodicandsolitonsolutionsfortheGeneralizedBBMandBurgers⁃BBMequations[J],AppliedMathematicsandComputation,2010(217):1430-1434.[3]ZhaoXQ,ZhiHY.AnimprovedF⁃expansionmethodanditsapplicationtocoupledDrinfel d⁃Sokolov⁃Wilsonequation[J],CommunicationsinTheoreticalPhysics,2008(50):309-314.[4]刘式括,付遵涛,刘式达等.一类非线性方程的新周期解[A],北京大学物理学报,2002(51):10-14.[5]AbdouMA.TheextendedF⁃expansionmethodanditsapplicationforaclassofnonlinearevolutionequations[J],Chaos,SolitonsandFractals,2007(31):95-104.[6]刘式括,付遵涛,刘式达,赵强.Jacobi椭圆函数展开法及其在求解非线性波动方程中的应用[A],北京大学物理学报,2001,50(11):2068-2071.[7]KumarH,ChandF.ApplicationsofextendedF⁃expansionandprojectiveRicattiequationmethodsto(2+1)⁃dimensionalsolitonequations[J],AipAdvance3,2013,3(032128):1-20.。
一种改进的基于Jacobi椭圆函数的随机平均法
![一种改进的基于Jacobi椭圆函数的随机平均法](https://img.taocdn.com/s3/m/b3ff0d3d53d380eb6294dd88d0d233d4b14e3f95.png)
一种改进的基于Jacobi椭圆函数的随机平均法
随机平均法是解决椭圆偏微分方程的一种有效方法。
Jacobi椭圆函数是随机平均法所使用的重要工具。
本文将介绍一种改进的Jacobi椭圆函数随机平均法。
首先,我们来回顾一下随机平均法的基本思想。
随机平均法可以将椭圆偏微分方程表
示为各个随机序列的平均值,而这些随机序列是在一个特定的概率密度函数下产生的。
其中,Jacobi椭圆函数在随机平均法中发挥了重要作用。
Jacobi椭圆函数代表将一个椭圆变形为一个矩形的变换。
这个变换可以简化椭圆偏微分方程的解析解。
然而,传统的随机平均法在计算时间和精度方面存在一些问题。
解决这些问题的一种
方法是通过引入一些新的随机序列来改进Jacobi椭圆函数随机平均法。
具体来说,我们可以使用更高阶的Jacobi椭圆函数,这可以提高计算精度。
此外,我们还可以使用多项式随机序列来产生更好的概率密度函数。
这些改进可以显著提高计算效率和精度。
值得注意的是,我们还可以使用并行计算来加速Jacobi椭圆函数随机平均法。
这可以通过将计算分成几个独立的部分并同时计算这些部分来实现。
这样可以大大提高计算速度,同时还可以降低计算误差。
改进的Jacobi椭圆函数随机平均法在解决椭圆偏微分方程方面表现出色。
它具有高精度和高效率的优点,同时还可以使用并行计算来更好地处理大型问题。
因此,它是一个非
常有前途的方法,在数学和工程领域都可以得到广泛的应用。
构造非线性发展方程Jacobi椭圆函数精确解的一种方法
![构造非线性发展方程Jacobi椭圆函数精确解的一种方法](https://img.taocdn.com/s3/m/cdb6e387a0116c175f0e4873.png)
在 非线 性孤 立 子理论 中 , 构造非 线性 发 展方 程 ( ) 组 的精 确 解有 许 多 方 法 , 双 曲正 切 函数 法 , 次平 如 齐
衡 [J。 椭 函 展 法 4辅 方 法5 等 文 [ 7用 助 程( )一c 6 ( + 法2 ab 圆 数 开 l 助 程 C 3. 献 6 ] 辅 方 + z n 3 ci , -  ̄ , )
收 稿 日期 :2 0 —0 —0 07 6 1
( 4 )
基 金项 目 :国家 自然科 学基 金 资助 项 目(0 6 0 6 I内蒙 古 自然 科 学 基 金 资 助 项 目( 0 4 8 2 13 I内蒙 古 高 等 学 校 科 学 研 究 资 助 项 目 14 10 ) 200000 )
Jcb 椭圆函数精确解. ao i 本文受文献 [ — 9 5 ]的启发, 引入一 个双曲函数 型辅助方程 , 2 种情况讨论 了 分 0
Jc b 椭 圆 函数 形 式 的 s hz ao i i ()和 c s () 借 助 符 号 计 算 系 统 Mah mai n oh , te t a构 造 了 正 规 长 波 方 程 和 c mKd 方 程 的 J c b 椭 圆 函数 精确 解 . V ao i
( )一c 6nz s ) a s ) + Sh( h + h i 。 c o
所确 定 . 经过 求 解方 程 ( )得到下 列形 式 的 s h 4 i ()和 c s ( ) n o hz :
当 a= 2 1 k ) ( + z ,b一 2 - 1 k ) C=一 2 1 k )时 , ( + , (+
摘
要 :在 双 曲 正 切 函 数 法 、 次 平衡 法 和 辅 助 方 程 法 的基 础 上 . 出 一 种 双 曲 函数 型辅 助 方 程 和 函 数 变 换 齐 给
离散的非线性薛定谔方程的Jacobi椭圆函数解
![离散的非线性薛定谔方程的Jacobi椭圆函数解](https://img.taocdn.com/s3/m/b5d51109ba68a98271fe910ef12d2af90242a87c.png)
离散的非线性薛定谔方程的Jacobi椭圆函数解一、Jacobi椭圆函数展开法近期提出并发展的Jacobi椭圆函数展开法可用来求解非线性数学物理方程的周期波解,Jacobi椭圆函数展开法又可看成是F-展开法的具体情形。
Jacobi椭圆函数展开法解非线性微分-差分方程分为以下4个步骤。
第一步:设n(t)可表示成Jacobi椭圆函数或tan函数展开法J的有限幂级数形式第二步:由非线性项与最高阶偏导数项的齐次平衡来确定M的值;第三步:将(1)代入所需要求解的方程中,得到关于J的多项式,置各J的系数为零,得到(可能是k)的代数方程组;第四步:解第三步得到的方程组得到代入(1)中得到方程的周期波解,取极限得极限情形下的方程的孤立子解。
二、方法应用给出离散的非线性薛定谔方程如下:1.sn椭圆函数展开。
首先用Jacobi椭圆sin函数展开法解方程(4),设方程(4)的解为将(5)、(6)、(7)代入方程(4)得:化简(8),并令同类项的系数为零,得如下一个代数方程组:2.cn椭圆函数展开。
用Jacobi椭圆cosin函数代替sin函数解方程(4),类似(3.1)中的方法得方程(2)的解为:3.dn椭圆函数展开。
用Jacobi椭圆-函数代替sin函数解方程(4),类似(3.1)中的方法得方程(2)的解为参考文献[1]M·J·AblowiteandP·A·Clarkson.Solitons,NonlinearEvolut ionEquationnsandInverseScattering[M].CambridgeUniversityPress,NewYork(1991)[2]朱加民.非线性离散薛定谔方程的显示精确解[J].江西科学.2022,23(4):402~404基金项目:本文系甘肃政法学院科研资助项目(青年项目)基金(GZF2022XQNLW25)。
Jacobi椭圆函数展开法在两个非线性偏微分方程解中的应用
![Jacobi椭圆函数展开法在两个非线性偏微分方程解中的应用](https://img.taocdn.com/s3/m/a0f8e0cb6137ee06eef91801.png)
2 0 1 3年 9 月
V o 1 . 3 4 N o . 5
S e p . 2 0 1 3
井 冈 山大 学学报 ( 自然科学 版) J o u r n a l o f J i n g g a n g s h a n U n i v e r s i t y( Na t u r a l S c i e n c e ) 4
T W O NoNLI NE AR RTI AL DI F F E RE NTI AL E QUA T I oNS
RUAN Chu a n— t o n g, ZJ ANG Rui . 1 i
( De p a r t me n t o f Ma t h e ma t i c s , Z h o u k o u No r ma l U n i v e r s i t y , Z h o u k o u , H e n a n 4 6 6 0 0 1 , C h i n a )
得 的非线性 偏微 分方程 的精 确解类 型丰 富 ,深受 学
者们 喜爱 。 本文利 用文献 [ 5 ] 提出了 “ 秩 ”的概念 ,分析 了
e x a c t s o l u t i o n
自然 学 科 中 的很 多现 象 是 非 线 性 学科 研 究 的 重 要 内容 ,我们 可 以把 非线 性 的许 多 问题 归结为求 解 非线 性偏微 分 方程 的精确解 的问题 。可 由于非 线 性方程 的复杂性 ,对 于非线 性方 程还 没有通 用 的解 法 。不 过经过 长 时间不 断努力 ,学者们 发现 了很 多
文章编号:1 6 7 4 — 8 0 8 5 ( 2 0 1 3 ) 0 5 — 0 0 0 4 — 0 5
用扩展的Jacobi椭圆函数展开法求解一类耦合方程组
![用扩展的Jacobi椭圆函数展开法求解一类耦合方程组](https://img.taocdn.com/s3/m/4ff3eb3c10661ed9ad51f3ed.png)
k u" p v +qc u +r ” c + k l l u2  ̄ k 一
【 c’ ’ km + k 0 一 v “ + s hv : k + v v
。
’
.
各 方程每 项的 秩分 别为 l ,,:1 , 1 是奇数 ,满足 秩 同类 的 ,1 1 3 ,1 都 1,
《 F +F J ” 2 ② =p F
v=2 。 , =l 时,F( =C r 一1 ~r ) n也 时满 足① ,② 。当
』 =-1 =2 . :, 一l F( = n同时满足① , , 一, , : 时 ) d ② ’
做 行 波 变 换 ( 2 ) , 将 方 程 ( 4 ) 化 为
( 5)
们得到:( ’ 1 ) 一, ) , 一(。 ) +l ) =(一 2 ( . 1 。 =。 , +1 . ,
令: =, + =一, 十1, =f n:F 4 所 以 S 分 . : ( ) .S () n 别 满 足 一阶 常 微 分 方 程 ① 和 二 阶 常 微 分 方 程 ② :
代人 ( ) 1 得到关于 F 4 的多项式。令其同次幂的系数为 0 () ,即可
得到 关 于待 定常 数 a( =l … .以及 k ,r 代数 方程 组 ,我们 i i , ) 2 ,C 的
1 .扩 展 的 Jcb椭 圆 函数展 开法 aoi
J c b 椭圆函数展 开法是将偏微分方程的解展开成 J c b ao i ao i 椭 圆 函 数 的有 穷级 数 形 式 。 扩 展 的 Jc b椭 圆 函数 展 开 法 是 将 而 ao i 这些 J c b 椭圆函数用统一的一个 函数 F来代替 。 ao i
时满足① , 所以我们可以 ②。 看到, 对不同的 J oi a b椭圆函数 F 4 c ()
修正jacobi椭圆函数展开法及其应用
![修正jacobi椭圆函数展开法及其应用](https://img.taocdn.com/s3/m/5e26892fcd1755270722192e453610661fd95a52.png)
修正jacobi椭圆函数展开法及其应用
Jacobi椭圆函数展开法是一种用于解决椭圆型偏微分方程的数值解法,它是由德国数学家Jacobi于1829年提出的。
Jacobi椭圆函数展开法的基本思想是将椭圆型偏微分方程分解
为一系列的椭圆函数,然后用椭圆函数的级数展开式来求解。
Jacobi椭圆函数展开法的优点是可以解决复杂的椭圆型偏微分方程,而且计算结果精度高,可以得到更精确的解。
另外,Jacobi椭圆函数展开法的计算量较小,可以有效地减少计算时间。
Jacobi椭圆函数展开法的应用非常广泛,它可以用于解决复杂的椭圆型偏微分方程,如椭圆型偏微分方程的热传导问题、椭圆型偏微分方程的电磁场问题等。
此外,Jacobi椭圆函数展开法还可以用于解决椭圆型偏微分方程的统计物理问题,如椭圆型偏微分方程的热力
学问题等。
总之,Jacobi椭圆函数展开法是一种非常有效的数值解法,它可以用于解决复杂的椭圆型偏微分方程,具有计算精度高、计算量小、应用广泛等优点,是一种非常有效的数值解法。
一种改进的基于Jacobi椭圆函数的随机平均法
![一种改进的基于Jacobi椭圆函数的随机平均法](https://img.taocdn.com/s3/m/cd6bcfb8360cba1aa811da8d.png)
一种改进的基于Jacobi椭圆函数的随机平均法作者:徐文俊郑丽文马品奎来源:《振动工程学报》2019年第03期摘要:建立了改进的基于Jacobi椭圆函数的随机平均法,用于预测有界噪声激励作用下硬弹簧和软弹簧系统的随机响应。
通过引入基于Jacobi椭圆函数的变换,导出关于响应幅值和激励与响应之间相位差的随机微分方程,应用随机平均原理,将响应幅值近似为一个Markov 扩散过程,建立其平均的It随机微分方程。
响应幅值的稳态概率密度由相应的简化FokkerPlanckKolmogorov方程解出;进而得到系统位移和速度的稳态概率密度。
以DuffingVan der Pol 振子为例,研究了硬刚度及软刚度情形下的随机响应,通过与Monte Carlo数值模拟结果比较证实了此方法的可行性及精度。
由于广义调和函数是基于线性系统的精确解,Jacobi椭圆函数是基于非线性系统的精确解,研究结果表明基于Jacobi椭圆函数的随机平均法得到的结果与Monte Carlo模拟方法更接近。
因此与基于广义调和函数的随机平均相比,基于Jacobi椭圆函数更加精确,因为它是基于保守的非线性系统。
关键词:随机振动; 随机平均; 有界噪声; 硬刚度; 软刚度中图分类号: O324; O322; 文献标志码: A; 文章编号: 1004-4523(2019)03.0444.08引言随机平均法是非线性随机系统响应分析的有效方法之一。
该方法在保留系统本质非线性特性的同时降低了系统维数,应用平均原理后,系统的慢变过程近似为扩散Markov过程,通过求解相应的FokkerPlanckKolmogorov (FPK)方程得到响应概率密度,随机平均技术基于Khasminskii[12]提出的一些定理,迄今的研究可归为以下5类:标准随机平均法[3]、能量包线随机平均法[47]、拟Hamilton系统随机平均法[811]、基于广义谐和函数的随机平均法[12]、基于椭圆函数的随机平均法[13]。
jacobi 椭圆函数
![jacobi 椭圆函数](https://img.taocdn.com/s3/m/e4f902d380c758f5f61fb7360b4c2e3f57272534.png)
Jacobi椭圆函数是一类特殊的椭圆函数,通常用符号$sn(u|m)$、$cn(u|m)$和$dn(u|m)$来表示。
它们是Jacobi椭圆积分的逆函数。
这些函数是由Carl Gustav Jacobi在19世纪开发的。
Jacobi椭圆函数经常在数学和物理领域中出现,特别是在椭圆曲线理论、非线性振动问题、薛定谔方程的椭圆函数解等方面的研究中。
这些函数在实际应用中也具有广泛的用途,例如计算数学表格、信号处理、图像处理等。
Jacobi椭圆函数的定义涉及两个参数,$u$是椭圆函数的自变量,$m$是椭圆模量。
模量$m$的取值范围在0到1之间,可以调整椭圆函数的性质。
具体而言,Jacobi椭圆函数有以下三个基本形式:
1. $sn(u|m)$:Jacobi椭圆正弦函数,它描述了椭圆曲线上的点的纵坐标。
2. $cn(u|m)$:Jacobi椭圆余弦函数,它描述了椭圆曲线上的点的横坐标。
3. $dn(u|m)$:Jacobi椭圆Delta函数,它描述了椭圆曲线上的点的切线斜率的倒数。
这些函数的性质和计算方法在数学文献中有详细的介绍,包括级数展开、递归关系和数值近似等方法。
对于实际的计算,可以借助计算机软件或数学库来进行计算。
需要注意的是,Jacobi椭圆函数可能有多种定义形式和符号约定,具体的使用方式和参数解释应根据具体的数学文献或应用场景而定。
一类非线性发展方程的改进的Jacobi椭圆函数精确解
![一类非线性发展方程的改进的Jacobi椭圆函数精确解](https://img.taocdn.com/s3/m/7bf8fcfe51e2524de518964bcf84b9d529ea2c5b.png)
一类非线性发展方程的改进的Jacobi椭圆函数精确解肖亚峰;薛海丽;张鸿庆【期刊名称】《中北大学学报(自然科学版)》【年(卷),期】2012(033)002【摘要】Expansion method for extended elliptic function is improved and applied to construct the explicit exact solutions of a class of ordinary differential equation. A series of explicit exact solutions are obtained by this method. In some extrem cases, the corresponding solitary wave solutions and the single triangle function solutions are obtained. The exact solutions of many nonlinear evolution equations, such as Modified Improved Boussinesq(MIB) equation, nonlinear Schrodinger equation, MKdV equation etc. , can be obtained by transforming into the ordinary differential equation.%对扩展的Jacobi椭圆函数展开法进行了改进,并将其应用到一类常微分方程中,比较方便地得到了该方程的一系列新的精确解,在极限情况下可得到相应的孤立波解和单周期波解.许多非线性发展方程(如Modified Improved Boussinesq(MIB)方程,非线性薛定谔方程,MKdV方程等)都可借助此方程得到其相应的新的精确解.【总页数】5页(P93-97)【作者】肖亚峰;薛海丽;张鸿庆【作者单位】中北大学理学院,山西太原030051;中北大学软件学院,山西太原030051;大连理工大学数学科学学院,辽宁大连116024【正文语种】中文【中图分类】O175.2【相关文献】1.一类非线性发展方程的显式精确解 [J], 陈自高;蔡洪涛2.一类非线性偏微分方程的改进的Jacobi椭圆函数精确解 [J], 肖亚峰;薛海丽;张鸿庆3.一类非线性发展方程的显式精确解 [J], 陈自高;张愿章4.一类非线性发展方程的精确解 [J], 王振立;刘勇5.构造非线性发展方程Jacobi椭圆函数精确解的一种方法 [J], 套格图桑;斯仁道尔吉因版权原因,仅展示原文概要,查看原文内容请购买。
一类非线性偏微分方程精确解的表达
![一类非线性偏微分方程精确解的表达](https://img.taocdn.com/s3/m/e9888e39eef9aef8941ea76e58fafab069dc44af.png)
一类非线性偏微分方程精确解的表达非线性偏微分方程是一类具有非线性特征的偏微分方程,在许多科学和工程领域中起着重要的作用。
与线性偏微分方程不同,非线性偏微分方程的解决方法更为困难,并且往往只能通过数值方法或近似解方法来求解。
然而,仍然有一些非线性偏微分方程存在着一些特殊的解,称为精确解。
精确解是指满足非线性偏微分方程的全部解析表达式。
这些精确解通常具有简洁的形式和重要的物理意义,因此对于理论研究和实际应用都具有重要的价值。
以下将介绍一些常见的非线性偏微分方程及其精确解。
1. 汉密尔顿-雅可比方程(Hamilton-Jacobi equation)汉密尔顿-雅可比方程在经典力学和量子力学中广泛应用,它的一般形式为:∂S/∂t+H(∇S,x)=0其中,S是哈密顿函数的特征函数,H是哈密顿量。
通常情况下,这个方程只能通过近似解或数值方法来求解。
但是,在一些特殊情况下,汉密尔顿-雅可比方程存在一些精确解。
例如,当哈密顿量H满足一定条件时,可以通过分离变量法或特殊变量变换得到精确解。
2. 伪线性方程(Pseudo-Linear Equation)伪线性方程是一类介于线性和非线性之间的方程,它具有其中一种线性性质但是包含了非线性项。
伪线性方程的精确解可以通过多种方法来求解,如分离变量法、变换法、叠加法等。
3. 密立根方程(Burgers' equation)密立根方程是一种具有非线性性质的守恒型方程,广泛应用于流体动力学和量子场论等领域。
它的一般形式为:∂u/∂t+u∂u/∂x=ν∂^2u/∂x^2其中,u是速度场,ν是粘性系数。
密立根方程的精确解可以通过特殊变量变换、相似变量法、分析解法等多种方法来求解。
4. 黏弹性流体方程(Viscoelastic fluid equation)黏弹性流体方程是描述黏弹性流体动力学行为的方程,具有非线性特点。
黏弹性流体方程的精确解多数情况下较为困难,通常需要通过数值方法来求解。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1 方 法 简述
为了 方便应 用 本文 中提 出的改进 方 法 , 面列 出主要求解 步 骤 : 下
步骤 1 对 给定 的非 线性 偏微 分方 程 ( 不妨 假设 仅含有 两 个变 量 z, £ )
F( t ,‘, 矗, , , )= 0 ,‘ z U 口 扛 … f , () 1
方程( E) NE s的精确 解 尤其是 孤 立子 解起 的重 要作 用 引起 了越来 越 多 学 者 的兴 趣 . 近年来 , 多 数学 家 许 和 物理 学家 经过 不 懈 努 力 , 现 了孤 立 子 理 论 中蕴 藏 着 一 系 列 构 造 精 确 解 的有 效 方 法 , 反 散 射 方 发 如 法 Biku d变 换 、 ro x变 换[ 、 rt 换 。齐 次 平 衡 法[ P ilv 析[ T n u、 ̄ ln c Dab u 3 Hi a变 ‘、 ] o 、 a ee分 、 ah方 法 及 n 推 广 的 Ta h方法 等[ 1 . n 7 0 随着各 种方 法 的出 现 , 多新 的 、 有重 要 物理 意义 的解 不断被 发现 和应用 . -] 许 具 最 近 , i hk o等人 提 出 了 Jcb 椭 圆函数展 开法 , 功地 求 出 了一大 类非线 性 演化方 程的一 系 LuS iu ao i 成
方程 ( ) 如下 形式 的解 : 5有
~ +l^十 s Jl n ‰.
所求 得 的 A、 、。 a a 和 。的值代 入 ( ) , 得如 下精 确解 : 6中 可
第 一组解 :
㈤
根据前 述 步骤 3 5将 () 入 () , - , 6代 5 中 得到 一个 关 于 D‘A (=O 1 .) .( , ,. 的代 数方 程 ( ) 收集 关 S ・ 组 . 于 n A (=O1 .) 同幂 次项 , s( f , ’. 的 ・ 并且 令 它们 的 系数 为零 , 得到 一 个关 于 A、la / 、。和 a a 的超 定 非 线 性代 数方 程组 . 使用 Ma l软 件包 “ h res求 解上 面 的超 定代 数方 程组 , pe C ast” 可得 A、。a 、。和 a 。的值. 将
3 大连 理工 大学 数学科 学学 院 , . 辽宁 大连 1 6 2 ) 10 4 摘 要 : Jcb 椭圆 函数展开法进行 了改进 , 对 ao i 并将其应用 到一类常微 分方 程 中, 比较方便 的得到其 Jcb 椭 ao i
圆 函数解 . 许多非线性发 展方 程都可借助该方程得到其相应 的精 确解 , MK V方 程 、2 ) MKP方 程及 如 d ( +1 维 非线性波方 程等方程的一 系列新 的精 确解 .
文 章 编 号 :1 7 -9X(0 1 0 -0 10 6 26 1 2 1 )40 0-5
一
类 非 线性 偏 微 分 方程 的 改进 的 J cb 椭 圆函数 精 确 解 o i a
肖亚峰 薛海 丽 张 鸿 庆 。 , ,
(. 1 中北大学 数学 系, 山西 太原 0 0 5 ;. 30 12 中北 大学 软件学院 , 山西 太原 0 0 5 ; 30 1
首先考 虑 如下 一类 常微 分 方程
=
h “4 h , i -3
() 5
其 中 , h h , 为 常数 .
由于在行波变换下 , d MK V方程 、2 ) MK ( +1维 P方程及非线性波动方程等方程均可化为方程() 5. 故本 文首先 考 虑方 程 ( ) 5 的形 如 ( ) 4 的解 , 然后 借 助 ( ) 5 的解 可 得到上 述三 个 方程 的解. 根据前 述 的求解 步 骤 2 通 过 平 衡 方程 ( ) 的最 高 阶线 性 项 和非 线 性项 h 。 可 知 = 1 故 设 , 5中 。 , .
精确解. 下面我们 以 MK V方程 、2 ) MK d (+1维 P方程及非线性波方程为例来构造相应的精确解.
( )MKd 方程 1 V
=
一
“ ( 一 ± 0
+ 1
,
+
m
千s i n
1
+
1
一
) + ) +
.
1
一
8
J l ( = ±
2 m+ 1 (
l+ 1 (
+
qs = i n
2 m+ 1 (
! ] , ! ] ,
一
(3 1)
( 4 1)
( 一 土 9
ns
( 8 1)
) +
1 3
( = ±
卜
卜 厂
L
1 5
2 m一 1 (
一-
一 ns
f
h i
! ] ,
( 9 1)
2一 c
l 4
( = 士
ns
( 2 - h ) i
s
.
(O 2)
+
当 一 1时 , 述解 退化 为 : 上
些方 法求 出了一 大类 非线 性 演化 方程 的 精 确 解及 其 相 应 的退 化 解 . 文 受 上 述 文 献 的启 发 , 出基 于 本 提 Jcb 椭 圆函数 的新 的 Jcb 函数 展 开方 法 , ao i ao i 可构造 出非 线性 波 动方 程新 的 Jcb椭 圆函数 解. 该方 ao i 将 法 应用 到 一类 常微 分 方 程 , 到 了该 方 程 丰 富 的 精 确 解 , 可 进 一 步 得 到 了 MKd 方 程 、 2 1 维 得 并 V (+ ) MKP方 程及 非线 性 波动 方程 等方 程新 的精 确解 .
注 3 如果 将式 ( ) 4 中的 Jcb 椭 圆 函数替换 为 另外 l 类 Jcb 椭 圆 函数 , ao i 1 ao i 则可 得到其 它新 的 J— a
4
甘 肃联合 大学学报 ( 自然科 学版)
第2 5卷
3 三类非线性演化方程 的精确解
如果 非线 性 发展方 程 能在行 波 变换 下化 为常微 分 方程 ()则 可 利用 () 5, 5 的解 () 2 ) 7 一(2获得 对应 的
解) .
注 1 Jc b椭 圆函数 n 具 有性 质 : ao i s
( ) s = 1n (
S n
, 2 l s = c t g ,i s  ̄ = c c . ( )i mn ( o h( ) l mn ( ) s (
肿 _1 . 胛_O .
2 一类 常微分方 程 的改进 的 Jc b 椭 圆函数解 ao i
2
甘 肃 联 合 大 学 学报 ( 自然 科 学 版)
第2 5卷
一
蔷
,
㈤
其中, 的值可通过平衡式() 3 中最高阶线性项和非 线性项来确定( : 注 如果得到 7的值不是正整数 , / " 则 可 以作变 换 ( 一 ()代入 ()然 后再 平衡 其 中的最 高 阶线 性项 和非 线性 项 ) , 3, ; 步骤 3 将() 4代入() 可得到一个关于 n A )i- ,, 的代数方程 ( . 3, s (- 1 …) ( _0 组) 收集关 于 n A ) s(
步骤 2 设 ( ) 3 的解 为
收 稿 日期 :0 10 -5 2 1—42 。 基 金 项 目 : 家重 点基 础 研 究 专 项 基 金 (0 4 B 10 0 l 国 20 C 3 80 ) 国家 自然 科 学 基 金 青 年 基 金 (0 0 1 5 . 1 9 14 )
作者 简介 : 肖亚峰( 9 6) 男 。 1 7 一, 山西芮城人 , 中北大学讲师 , 博士 , 主要 从事 符号计算和非线性物理研究。
第2 5卷 第 4期
2 1 年 7月 01
甘 肃联合大 学学报( 自然科 学版 )
J u n lo n u Lin eUnv ri Na u a ce c s o r a fGa s a h iest y( t r l in e ) S
Vo. 5No 4 12 .
J 12 1 u. 0 1
“ = 、 c( ̄ , 土/ s±/ ( 一 c 二
l =√ c千 . l 土一 s  ̄ s ( 。 c/ (
() 1 2
第 4期
肖亚 峰 等 : 类 非 线 性 偏 微 分 方 程 的 改进 的 Jcb 椭 圆 函数 精 确 解 一 ao i
3
第二 组解
一 一 _ 2 m + 1 ( 一 一 :
(=O 1… ) 同幂次 项 , ,, 的 并且 令 它们 的 系数 为零 , 到 一个 关 于 A、。a(=0 1 2 …) 得 、 f ,, , 的超 定 非线 性
代数方程组; 步骤 4 使用王东 明教授基 于吴代数消元法[] " 开发的 Ma l件包“ h r t”求解步骤 3中得 到 pe C as s , e
J
( = 士
+
幅 讪 )
 ̄c ) :t ih 十 o
( 1 2)
l
I
U6 =士 h 1( 一
3
一
幅
)
‘
(2 2)
l
汕 ) +
注 2 根据 所查 文献 , 第二 组 解均 为新 解. 当式 ( ) 4 中 一O时 , 则文 中方 法退化 为 LuS K. 提 出 i . 等 的 J c b 椭 圆 函数 法 引、 n aoi Ta h方 法 和扩展 的 Ta h方 法U t . n -o ] e b 椭 圆 函数 解 , oi 限于篇 幅 , 此不 再一 一列 出. 在
的超 定代 数方 程组 , 可得 A、ia(=0 12 … ) / , f ,, , 的值 I  ̄ 步骤 5 将 步骤 4中求得 的 A、 、 ,, , 代入 () , a(=O 1 2 …) 4 中 即可 得到方 程 () 1 的精确解 . 因此 , 根据 上 面的 步 骤 , 以得 到 方 程 ( ) 可 1 的一 系 列 新 的 Jcb 椭 圆 函 数 解 ( 括 孤 波 ( £ x,)= ( ), = . 一 c , 2 7 t () 2