主题
Search

欠阻尼简谐运动


SHOUnderdamped

欠阻尼简谐运动是 阻尼简谐运动 的一个特例

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

其中

 beta^2-4omega_0^2<0.
(2)

由于我们有

 D=beta^2-4omega_0^2<0,
(3)

因此,数量

gamma=1/2sqrt(-D)
(4)
=1/2sqrt(4omega_0^2-beta^2)
(5)

是正数。将试探解 x=e^(rt) 代入微分方程,得到满足以下条件的解

 r_+/-=-1/2beta+/-igamma,
(6)

即,解的形式为 如下形式

 x=e^(-(beta/2+/-igamma)t).
(7)

使用 欧拉公式

 e^(ix)=cosx+isinx,
(8)

这可以被重写为

 x=e^(-(beta/2)t)[cos(gammat)+/-isin(gammat)].
(9)

我们对实数解感兴趣。由于我们在这里处理的是线性齐次常微分方程,线性无关 解的线性组合也是解。由于我们在 (9) 中有这样的解的和,因此 虚部实部 分别满足常微分方程,因此是我们要求的解。正弦项前面的常数是任意的,所以我们可以将解识别为

x_1=e^(-(beta/2)t)cos(gammat)
(10)
x_2=e^(-(beta/2)t)sin(gammat),
(11)

所以通解是

 x=e^(-(beta/2)t)[Acos(gammat)+Bsin(gammat)].
(12)

初始值为

x(0)=A
(13)
x^.(0)=-1/2betaA+B,gamma
(14)

所以 AB 可以用初始条件表示为

A=x(0)
(15)
B=(betax(0))/(2gamma)+(x^.(0))/gamma.
(16)

上面的图显示了一个欠阻尼简谐振荡器,参数为 omega=0.3beta=0.4,针对各种初始条件 (A,B)

对于一个具有余弦强制函数的欠阻尼振荡器 g(t)=Ccos(omegat),因此

 x^..+betax^.+omega_0^2x=Ccos(omegat),
(17)

定义

gamma=1/2sqrt(4omega_0^2-beta^2)
(18)
alpha=1/2beta
(19)

为了方便起见,然后注意

4omega_0^2-beta^2=4gamma^2
(20)
omega_0^2=gamma^2+1/4beta^2
(21)
=gamma^2+alpha^2
(22)
beta=2alpha.
(23)

我们现在可以使用参数变分法来获得特解,如下所示

 x^*=x_1v_1+x_2v_2,
(24)

其中

v_1=-int(x_1(t)g(t))/(W(t))
(25)
v_2=int(x_2(t)g(t))/(W(t))
(26)

朗斯基行列式为 朗斯基行列式

W(t)=x_1x^._2-x^._1x_2
(27)
=gammae^(-2alphat).
(28)

这些可以直接积分得到

v_1=-C/gammainte^(alphat)sin(gammat)cos(omegat)dt
(29)
v_2=C/gammainte^(alphat)cos(gammat)cos(omegat)dt.
(30)

因此,

x^*(t)=C((alpha^2+gamma^2-omega^2)cos(omegat)+2alphaomegasin(omegat))/([alpha^2+(gamma-omega)^2][alpha^2+(gamma+omega)^2])
(31)
=C/(sqrt((omega_0^2-omega^2)^2+omega^2beta^2))cos(omegat+delta),
(32)

其中使用了 和角公式 以及

 delta=tan^(-1)((betaomega)/(omega^2-omega_0^2)).
(33)

另请参阅

临界阻尼简谐运动, 阻尼简谐运动, 过阻尼简谐运动, 简谐运动

使用 Wolfram|Alpha 探索

参考文献

Papoulis, A. Probability, Random Variables, and Stochastic Processes, 2nd ed. 纽约:McGraw-Hill,第 525-527 页,1984 年。

请引用本文为

Weisstein, Eric W. "欠阻尼简谐运动。" 来自 MathWorld--Wolfram Web 资源。 https://mathworld.net.cn/UnderdampedSimpleHarmonicMotion.html

主题分类