20、谐波振荡相关知识解析

谐波振荡相关知识解析

1. 临界阻尼运动的Mathematica代码示例

1.1 问题描述

使用Mathematica求解方程(6.4.3)的解 (x(t))，给定参数 (m = 1) kg，(k = 1) N/m，(x(0) = 1) m，(v(0) = 0) m/s 以及 (b = 2.0) Ns/m。

1.2 求解过程

对于给定的参数值，计算可得 (\omega^2 = \gamma^2 - \omega_0^2 = (b/2m)^2 - (k/m) = 0)。由于该值为零，我们预期会出现临界阻尼运动。下面使用Mathematica给出解析解，它是临界阻尼振荡的两种可能解（即 (e^{-\omega_0t}) 和 (te^{-\omega_0t}) 项）的线性组合，这与临界阻尼运动的解析解的一般形式(6.4.16)一致。注意，在这种情况下 ((b/2m)^2 = (k/m))，在 NDSolve 命令中，常数 (b) 被替换为 (b = 2m\sqrt{k/m} = \sqrt{4km})。

sol1 = DSolve[{m*x''[t]+k*x[t]+Sqrt[4*k*m]*x'[t]==0,x[0]==xo,x'[0]==0},x[t],t]//FullSimplify;
x1 = x[t]/.sol1[[1]]
m = 1;k = 1;xo = 1;

1.3 输出结果

输出结果为：(e^{-\frac{kt}{\sqrt{km}}}(m + \sqrt{kmt