← 所有文章

从简谐振子开始:一篇物理笔记

由牛顿方程到能量守恒,用一个熟悉的模型记录推导,也检验数学公式的排版。

  • 经典力学
  • 数学方法
  • 示例笔记

简谐振子是学习物理时最早遇到的模型之一。这篇示例笔记从运动方程出发,整理通解与能量守恒的关系。

本文是网站的测试文章,可自由修改或删除。文中的公式由本地数学排版资源生成。

1. 从运动方程出发

设质量为 mm 的质点受到线性回复力 F=−kxF=-kx,其中 k>0k>0。牛顿第二定律给出

md2xdt2=−kx.m\frac{\mathrm{d}^2x}{\mathrm{d}t^2}=-kx.

定义角频率 ω0=k/m\omega_0=\sqrt{k/m},则运动方程写成

x¨+ω02x=0.\ddot{x}+\omega_0^2x=0.

在给定初始条件 x(0)=x0x(0)=x_0、x˙(0)=v0\dot{x}(0)=v_0 后,解为

x(t)=x0cos⁡(ω0t)+v0ω0sin⁡(ω0t).x(t)=x_0\cos(\omega_0t)+\frac{v_0}{\omega_0}\sin(\omega_0t).

这也可以写成 x(t)=Acos⁡(ω0t+φ)x(t)=A\cos(\omega_0t+\varphi),其中振幅和相位由初始条件确定。

2. 能量为何守恒?

振子的机械能由动能与势能组成:

E=12mx˙ 2+12kx2.E=\frac{1}{2}m\dot{x}^{\,2}+\frac{1}{2}kx^2.

对时间求导,并代入运动方程:

dEdt=mx˙x¨+kxx˙=x˙(mx¨+kx)=0.\begin{aligned} \frac{\mathrm{d}E}{\mathrm{d}t} &=m\dot{x}\ddot{x}+kx\dot{x} \\ &=\dot{x}\left(m\ddot{x}+kx\right) \\ &=0. \end{aligned}

因此,系统的能量在运动过程中保持不变。取 v0=0v_0=0,便有 E=12kx02E=\tfrac12kx_0^2。

简谐振子的位移随时间变化示意图

图 1:归一化位移随时间作余弦变化,图中以角频率为 1 的情形示意。

3. 再检验几种数学排版

积分与求和

高斯积分和几何级数分别为

∫−∞∞e−αx2 dx=πα,α>0,\int_{-\infty}^{\infty}e^{-\alpha x^2}\,\mathrm{d}x =\sqrt{\frac{\pi}{\alpha}},\qquad \alpha>0, ∑n=0∞rn=11−r,∣r∣<1.\sum_{n=0}^{\infty}r^n=\frac{1}{1-r},\qquad |r|<1.

矩阵

令状态向量为 (x,x˙)T(x,\dot{x})^{\mathsf T},二阶方程可写为一阶方程组:

ddt(xx˙)=(01−ω020)(xx˙).\frac{\mathrm{d}}{\mathrm{d}t} \begin{pmatrix}x\\\dot{x}\end{pmatrix} = \begin{pmatrix}0&1\\-\omega_0^2&0\end{pmatrix} \begin{pmatrix}x\\\dot{x}\end{pmatrix}.

多行推导与窄屏阅读

较长的表达式可以用 aligned 换行。若公式仍宽于手机屏幕,可以在公式区域横向滑动:

12mω02A2sin⁡2(ω0t+φ)⏟kinetic energy+12kA2cos⁡2(ω0t+φ)⏟potential energy=12kA2[sin⁡2(ω0t+φ)+cos⁡2(ω0t+φ)]=12kA2.\underbrace{\frac{1}{2}m\omega_0^2A^2\sin^2(\omega_0t+\varphi)}_{\text{kinetic energy}} +\underbrace{\frac{1}{2}kA^2\cos^2(\omega_0t+\varphi)}_{\text{potential energy}} =\frac{1}{2}kA^2\left[\sin^2(\omega_0t+\varphi)+\cos^2(\omega_0t+\varphi)\right] =\frac{1}{2}kA^2.

4. 继续思考

  • 若加入与速度成正比的阻力,能量将如何变化?
  • 当回复力不再与位移成正比时,周期是否依赖振幅?
  • 数值积分时,怎样判断能量误差是否可以接受?

写笔记时,可以从一个明确的问题开始,保留必要的推导,再记录尚未解决的疑问。