阻尼衰减曲线用python_简谐运动之阻尼振动

续上一期的文章，我们今天接着来撸一撸带阻尼的简谐振子（damped oscillator）。

我们假设振子在运动过程中受到的阻力仅与它的运动速度有关，并且阻力的大小简单的正比于速度的大小，即

，这里

是一个描述阻力强弱的参数。结合简谐振子原本就受到的和位移成正比的回复力，我们可以写下阻尼振动的运动方程：

注意这里的两个负号：第一个负号是因为回复力的方向和位移相反，第二个负号是因为阻力的作用方向和运动速度方向相反。顺手将加速度

和速度

写成微分形式，我们得到

不妨定义

，可以理解为描述阻力相对强弱的一个指标，这样运动方程就可以改写成：

这依然是一个常系数的微分方程。求解之后，就可以得到振子的位移随时间的变化关系。在下面的讨论中，我们假定振子在

时刻时，初始位移为

，初始速度为

，即拖到偏离平衡位置的某个地方从静止释放。可以证明，满足这两个初始条件的微分方程的解将具有唯一的形式。

我们需要找的解

，其自身、一阶导数、二阶导数一家亲，可以互相抵消搞出个零，所以可以猜想解依然会具有

的形式。将试解（ansatz）

代入要求解的微分方程中，不难得到：

以上关系要对任意