话题
搜索

达芬微分方程


Duffing方程最普遍的强迫形式是

 x^+deltax(代扣所得税)^+(β^3+/-ω0^2x)=γ(ω+phi)。
(1)

根据所选参数,方程式可以采用多种特殊形式。例如在没有阻尼和强制的情况下,δ=伽马=0取加号,方程变成

 x^+ω0^2x+β3=0
(2)

(Bender和Orszag 1978年,第547页;Zwillinger 1997年,第122页)。该方程可以显示混沌行为。对于β>0,方程式表示“硬弹簧”和用于β<0,它代表“软弹簧”。如果β<0,的相图曲线是闭合的。

如果我们取而代之β=1,ω0=1,重置时钟,以便φ=0,并使用减号,则方程式为

 x^+deltax(代扣所得税)^+(x^3-x)=gammacos(欧米伽特)。
(3)

这可以写成一阶常微分方程组

x ^。=是的,
(4)
年^。=x-x^3-deltay+gammacos(欧米伽特)
(5)

(威金斯1990年,第5页),在非强制情况下,这一点可以简化为

x ^。=年
(6)
年^。=x-x^3-位移
(7)

(威金斯1990年,第6页;奥特1993年,第3页)。

这组耦合微分方程的不动点由下式给出

 x ^=y=0,
(8)

所以y=0,

年^。=x-x^3-位移
(9)
=x(1-x^2)-0
(10)

x=0,+/-1.因此,固定点是(-1,0),(0,0)、和(1,0).

不动点的稳定性分析可以通过线性化方程来进行。差异化提供

x^。。=年^。
(11)
=x-x^3-位移
(12)
y^。。=(1-3x^2)x^-代尔泰^。,
(13)

可以写为矩阵方程

 [x^..;y^..]=[0 1;1-3x^2-增量][x^.;y^.]。
(14)

检查点(0,0)的稳定性:

 |0-λ1;1-δ-lambda |=λ(λ+δ)-1=λ^2+λdadelta-1=0
(15)
 λ_+/-^((0,0))=1/2(-增量+/-平方(增量^2+4))。
(16)

但是增量^2>=0,所以λ_+/-^((0,0))是真实的。sqrt(增量^2+4)>|增量|,总会有一个积极的 ,所以这个不动点是不稳定的。现在看看(+/-1, 0). 特征方程为

 |0-λ1-2-δ-lambda |=λ(λ+δ)+2=λ^2+λdelta+2=0,
(17)

它有根

 λ_+/-^((+/-1,0))=1/2(-增量+/-平方(增量^2-8))。
(18)

对于增量>0,R[λ_+/-^((+/-1.0))]<0,所以这个点是渐近稳定的。如果增量=0,λ_+/-^((+/-1,0))=+/-isqrt(2),所以这个点是线性的稳定(Wiggins 1990,第10页)。然而,如果增量输入(-2sqrt(2),0),根式给出想像的部分实部>0,因此该点不稳定。如果δ=-2sqrt(2),λ_+/-^((+/-1,0))=平方(2),具有积极的 真实的 ,因此该点不稳定。如果增量<-2sqrt(2),然后|增量|<sqrt(增量^2-8),所以两者都是积极的这一点是不稳定的。

Duffing振荡器相图

有趣的是,特殊情况增量=0没有强迫,

x ^。=年
(19)
年^。=x-x^3,
(20)

可以通过求积积分。差异化(19)和堵塞英寸(20)给予

 x^=年^=x-x^3。
(21)

将两边乘以x ^。给予

 x^。。x ^-x^.x+x^.x^3=0。
(22)

但这是可以写的

 d/(dt)(1/2x^.^2-1/2x^2+1/4x^4)=0,
(23)

所以我们有一个运动不变量小时,

 h=1/2x^^2-1/2x^2+1/4x^4。
(24)

解决x ^^2给予

 x ^^2=((dx)/(dt))^2=2小时+x ^2-1/2x ^4
(25)
 (dx)/(dt)=平方(2h+x^2-1/2x^4),
(26)

所以

 t=intdt=int(dx)/(平方(2h+x^2-1/2x^4))
(27)

(威金斯1990年,第29页)。

注意,运动不变量小时满足

 x ^=(部分)/(部分)=(部分)
(28)
 (部分h)/(部分x)=-x+x^3=-y^。,
(29)

因此,Duffing振荡器的方程由哈密顿量系统

x ^。=(部分)/(部分)
(30)
年^。=-(部分)/(部分)
(31)

(威金斯1990年,第31页)。


与Wolfram一起探索| Alpha

工具书类

C.M.本德。和S.A.Orszag。科学家和工程师的高级数学方法。纽约:McGraw-Hill,第547页,1978年。奥特,E。混乱在动力系统中。纽约:剑桥大学出版社,1993年。塔博,M。混乱非线性动力学中的可积性:导论。纽约:威利,第35页,1989年。Trott,M.“数学指南附加材料:Duffing振荡器的Wigner函数。"http://www.mathematicaguidebooks.org/addressions.shtml#N_1_08.威金斯,“阻尼、强迫Duffing振子动力学的应用”§1.2E英寸引言应用非线性动力系统和混沌。纽约:Springer-Verlag,第5-6、10、23、26-32、44-45、50-51和153-175页,1990年。兹威林格,D.(编辑)。CRC公司标准数学表和公式。佛罗里达州博卡拉顿:CRC出版社,第413页,1995D.茨威林格。手册《微分方程》,第三版。马萨诸塞州波士顿:学术出版社,第122页,1997

参考Wolfram | Alpha

达芬微分方程

引用如下:

埃里克·魏斯坦(Eric W.Weisstein)。“Duffing微分方程。”发件人数学世界--Wolfram Web资源。https://mathworld.wolfram.com/DuffingDifferentialEquation.html

主题分类