本文主要是介绍微分方程(Blanchard Differential Equations 4th)中文版Section4.4,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!
稳态的振幅和相位系统
在本节中,我们回到方程
d 2 y d t 2 + p d y d t + q y = cos ω t \frac{d^2 y}{dt^2} + p \frac{dy}{dt} + qy = \cos \omega t dt2d2y+pdtdy+qy=cosωt
用于周期性强迫的阻尼谐振子。我们的目标是建立解决方案行为与参数之间的定量关系——特别是决定强迫频率的 ω \omega ω 和决定阻尼的 p p p。
正如我们在第4.2节中所见,所有解决方案的长期行为是相同的。它由强迫方程的特解 y p ( t ) y_p(t) yp(t) 决定,形式可以写成
y p ( t ) = A cos ( ω t + ϕ ) . y_p(t) = A \cos(\omega t + \phi). yp(t)=Acos(ωt+ϕ).
因此,对于大时间 t t t,所有解的振荡频率与强迫相同。这是解决方案的稳态行为。振荡的振幅是 A A A,相位角是 ϕ \phi ϕ。在本节中,我们将找到 A A A 和 ϕ \phi ϕ 关于 ω \omega ω、 p p p 和 q q q 的表达式。
稳态行为
为了找到强迫方程
d 2 y d t 2 + p d y d t + q y = cos ω t \frac{d^2 y}{dt^2} + p \frac{dy}{dt} + qy = \cos \omega t dt2d2y+pdtdy+qy=cosωt
的一个特解,我们遵循第4.2节中介绍的技术,并使用应用于复杂化方程的未定系数法
d 2 y d t 2 + p d y d t + q y = e i ω t . \frac{d^2 y}{dt^2} + p \frac{dy}{dt} + qy = e^{i\omega t}. dt2d2y+pdtdy+qy=eiωt.
接下来,我们假设 y c ( t ) = a e i ω t y_c(t) = ae^{i\omega t} yc(t)=aeiωt 并确定使 y c ( t ) y_c(t) yc(t) 成为解的 a a a 的值。将 y c ( t ) y_c(t) yc(t) 代入微分方程的左侧,我们得到
d 2 y c d t 2 + p d y c d t + q y c = − a ω 2 e i ω t + p ( i a ω e i ω t ) + q ( a e i ω t ) = a [ ( q − ω 2 ) + i p ω ] e i ω t . \frac{d^2 y_c}{dt^2} + p \frac{dy_c}{dt} + qy_c = -a\omega^2 e^{i\omega t} + p(i a \omega e^{i\omega t}) + q(a e^{i\omega t}) = a \left[ (q - \omega^2) + i p \omega \right] e^{i\omega t}. dt2d2yc+pdtdyc+qyc=−aω2eiωt+p(iaωeiωt)+q(aeiωt)=a[(q−ω2)+ipω]eiωt.
为了使 y c ( t ) y_c(t) yc(t) 成为解,我们必须取
a = 1 b a = \frac{1}{b} a=b1
其中 b b b 是复数 ( q − ω 2 ) + i p ω (q - \omega^2) + i p \omega (q−ω2)+ipω。
正如我们在第4.2节中所见,我们可以从复数 a a a 的极坐标形式
a = ∣ a ∣ e i ϕ = A e i ϕ a = |a| e^{i\phi} = A e^{i\phi} a=∣a∣eiϕ=Aeiϕ
中确定 y p ( t ) y_p(t) yp(t) 的振幅 A A A 和相位角 ϕ \phi ϕ。(有关复数的极坐标表示的回顾,请参见复数附录。)在这种情况下,由于 a = 1 / b a = 1/b a=1/b,我们有
∣ a ∣ = 1 ∣ b ∣ |a| = \frac{1}{|b|}
这篇关于微分方程(Blanchard Differential Equations 4th)中文版Section4.4的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!