主题
Search

Duffing 微分方程


Duffing 方程最常见的受迫形式是

 x^..+deltax^.+(betax^3+/-omega_0^2x)=gammacos(omegat+phi).
(1)

根据所选参数,该方程可以呈现多种特殊形式。例如,在没有阻尼和没有外力的情况下,delta=gamma=0 并取加号,则方程变为

 x^..+omega_0^2x+betax^3=0
(2)

(Bender 和 Orszag 1978, p. 547; Zwillinger 1997, p. 122)。此方程可以表现出混沌行为。对于 beta>0,该方程表示“硬弹簧”,对于 beta<0,它表示“软弹簧”。如果 beta<0,则相图曲线是闭合的。

如果改为取 beta=1omega_0=1,重置时钟使得 phi=0,并使用减号,则方程变为

 x^..+deltax^.+(x^3-x)=gammacos(omegat).
(3)

这可以写成一阶常微分方程组,如下所示

x^.=y,
(4)
y^.=x-x^3-deltay+gammacos(omegat)
(5)

(Wiggins 1990, p. 5),在无外力的情况下,简化为

x^.=y
(6)
y^.=x-x^3-deltay
(7)

(Wiggins 1990, p. 6; Ott 1993, p. 3)。

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

 x^.=y=0,
(8)

因此 y=0,且

y^.=x-x^3-deltay
(9)
=x(1-x^2)-0
(10)

给出 x=0,+/-1。因此,定点为 (-1,0)(0,0),和 (1,0)

可以通过线性化方程来分析定点的稳定性。微分得到

x^..=y^.
(11)
=x-x^3-deltay
(12)
y^..=(1-3x^2)x^.-deltay^.,
(13)

这可以写成矩阵方程

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

考察点 (0,0) 的稳定性

 |0-lambda 1; 1 -delta-lambda|=lambda(lambda+delta)-1=lambda^2+lambdadelta-1=0
(15)
 lambda_+/-^((0,0))=1/2(-delta+/-sqrt(delta^2+4)).
(16)

但是 delta^2>=0,所以 lambda_+/-^((0,0)) 是实数。由于 sqrt(delta^2+4)>|delta|,总会有一个,因此该定点是不稳定的。现在看 (+/-1, 0)。特征方程为

 |0-lambda 1; -2 -delta-lambda|=lambda(lambda+delta)+2=lambda^2+lambdadelta+2=0,
(17)

其根为

 lambda_+/-^((+/-1,0))=1/2(-delta+/-sqrt(delta^2-8)).
(18)

对于 delta>0R[lambda_+/-^((+/-1,0))]<0,因此该点是渐近稳定的。如果 delta=0lambda_+/-^((+/-1,0))=+/-isqrt(2),因此该点是线性稳定的 (Wiggins 1990, p. 10)。然而,如果 delta in (-2sqrt(2),0),根式给出虚部实部>0,因此该点是不稳定的。如果 delta=-2sqrt(2)lambda_+/-^((+/-1,0))=sqrt(2),它有一个,因此该点是不稳定的。如果 delta<-2sqrt(2),那么 |delta|<sqrt(delta^2-8),所以两个都是的,并且该点是不稳定的。

DuffingOscillatorPhasePortrait

有趣的是,特殊情况 delta=0 且没有外力,

x^.=y
(19)
y^.=x-x^3,
(20)

可以通过求积分来求解。对 (19) 求导并代入 (20) 得到

 x^..=y^.=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,

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

求解 x^.^2 得到

 x^.^2=((dx)/(dt))^2=2h+x^2-1/2x^4
(25)
 (dx)/(dt)=sqrt(2h+x^2-1/2x^4),
(26)

因此

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

(Wiggins 1990, p. 29)。

请注意,运动不变量 h 满足

 x^.=(partialh)/(partialx^.)=(partialh)/(partialy)
(28)
 (partialh)/(partialx)=-x+x^3=-y^.,
(29)

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

x^.=(partialh)/(partialy)
(30)
y^.=-(partialh)/(partialx)
(31)

(Wiggins 1990, p. 31)。


使用 Wolfram|Alpha 探索

参考文献

Bender, C. M. 和 Orszag, S. A. 科学家和工程师高级数学方法。 纽约:McGraw-Hill, p. 547, 1978。Ott, E. 动力系统中的混沌。 纽约:剑桥大学出版社,1993。Tabor, M. 非线性动力学中的混沌和可积性:导论。 纽约:Wiley, p. 35, 1989。Trott, M. "Mathematica 指南 补充材料:Duffing 振荡器的 Wigner 函数。" http://www.mathematicaguidebooks.org/additions.shtml#N_1_08Wiggins, S. "应用于阻尼、受迫 Duffing 振荡器的动力学。" §1.2E 在 应用非线性动力系统和混沌导论。 纽约:Springer-Verlag, pp. 5-6, 10, 23, 26-32, 44-45, 50-51, 和 153-175, 1990。Zwillinger, D. (Ed.). CRC 标准数学表格和公式。 Boca Raton, FL: CRC Press, p. 413, 1995。Zwillinger, D. 微分方程手册,第 3 版。 Boston, MA: Academic Press, p. 122, 1997。

在 Wolfram|Alpha 中引用

Duffing 微分方程

请引用为

Weisstein, Eric W. "Duffing 微分方程。" 来自 MathWorld--Wolfram Web 资源。 https://mathworld.net.cn/DuffingDifferentialEquation.html

主题分类