微分方程 · 简谐振动

简谐振动

简谐振动是回复力与位移成正比的运动,x″ = −ω²x。它的状态(x, v/ω)以恒定的角速度ω转动,所以不管拉开多远,每次来回的时间都相同。

01

简谐振动与本页画面

简谐振动(也叫简谐运动)是物体在与位移成正比的回复力作用下的运动:mx′′=−kxm x'' = -k x,即x′′=−ω2xx'' = -\omega^2 x,ω=k/m\omega = \sqrt{k/m}。它的解都是x=Acos⁡(ωt−φ)x = A\cos(\omega t - \varphi),周期T=2π/ωT = 2\pi/\omega与振幅A无关。也常写成x=Asin⁡(ωt+φ)x = A\sin(\omega t + \varphi),相位取法不同。本页用弹簧振子说明为什么,并和单摆对照。

状态在转:(x, v/ω) 的速度垂直于半径,大小是半径的 ω 倍
ddt(xv/ω)=(va/ω)=ω(v/ω−x),A=x02+(v0/ω)2,E=12k (x2+(v/ω)2)=12kA2\frac{d}{dt}\begin{pmatrix} x \\ v/\omega \end{pmatrix} = \begin{pmatrix} v \\ a/\omega \end{pmatrix} = \omega\begin{pmatrix} v/\omega \\ -x \end{pmatrix}, \qquad A = \sqrt{x_0^2 + (v_0/\omega)^2}, \qquad E = \tfrac12 k\,(x^2 + (v/\omega)^2) = \tfrac12 kA^2

三层共用一根x轴。在相平面上,状态是点(x, v/ω)。它在相平面上的速度有两部分:横向是v,竖向是a/ω = −ωx;合起来处处垂直于半径,长度是半径的ω倍,所以整张平面上的点都以同一个角速度ω绕中心转,像一张唱片。导轨上的物块是这个点在x轴上的影子,所以x是余弦;点到中心的距离是振幅,方向是相位,距离的平方乘½k是能量——弹簧里的是½kx²,运动的是½k(v/ω)²。纵轴取v/ω才使轨道成圆:直接用v是椭圆,点也不再匀速转。相平面下面的纸带记下x随时间的变化,每0.1 s一道刻度,每次过零留一个记号;不论振幅多大,过零都相隔T/2,对照振子的记号和主振子的对齐。

单摆满足s′′=−gsin⁡(s/L)s'' = -g\sin(s/L),弧长s=Lθs = L\theta:回复力不与s成正比,所以不是简谐振动。它在相平面(s, ṡ/ω₀)上的能量线是卵形而不是圆,周期随摆幅变长,T=T0/AGM(1,cos⁡(θmax⁡/2))T = T_0/\mathrm{AGM}(1, \cos(\theta_{\max}/2));摆幅θmax⁡\theta_{\max}较小时约为T0(1+θmax⁡2/16)T_0(1 + \theta_{\max}^2/16)(θmax⁡\theta_{\max}用弧度),其中T0=2πL/gT_0 = 2\pi\sqrt{L/g}。这里的θmax⁡\theta_{\max}是由能量决定的摆幅,只有从静止放开时才等于起始摆角。单摆的纸带是数值解,并用上面的精确周期核对过;弹簧的纸带是精确解。越过虚线所示的分界轨道,摆会翻过顶点整圈转动;本页把摆角限制在170°以内。

02

值得注意的

  • 回复力与位移成正比,每次来回的时间就相同;单摆是对照,在那里小角近似失效。其他力也可能产生等时运动——惠更斯1673年的摆线摆就是一例,它的回复力与弧长成正比——所以本页不说反过来的结论。
  • 相流箭头是点在相平面上的速度,按同一个比例画(比例写在旁边):离中心越远越长,转得却一样快。换成单摆,箭头不再垂直于半径,绕中心转的快慢沿卵形变化。
  • 改k、m或L时,相平面上的起点不动(振幅和相位不变),所以以m/s计的初速度随ω换算;改之前的记录留作淡影,方便对照。
03

操作

  • 状态点可以拖到最大的圆内任何地方(贴近两条轴时会吸住:静止,或在中心),也可以拖物块或摆锤,从静止拉开。播放按真实时间走——1秒就是1秒——再按一次停下;读数每秒刷新四次。
  • 演示模式打开全屏视图:空格播放 / 停下,← → 从静止拉近或拉远,H隐藏周期和过零时刻。
04

相关概念