线性振动
线性振动是物理学中一种重要的运动形式,广泛存在于各种物理系统中.它通常指系统在平衡位置附近的小幅度振动,其运动方程可以用线性 微分方程 描述.
简谐振动
如果某一物理量在平衡点附近往复振动,且随时间的变化是余弦(或正弦)的函数形式,我们就把这种振动称为简谐振动.
运动方程推导
简谐振动是最基本的线性振动形式,其运动方程为:
𝑑2𝑥𝑑𝑡2+𝜔2𝑥=0d^2xdt^2 + omega^2 x = 0这一方程可以通过牛顿第二定律推导得到.例如,对于一个质量为 𝑚m 的物体,受到的回复力 𝐹F 满足胡克定律:
𝐹=−𝑘𝑥F = -kx根据牛顿第二定律:
𝐹=𝑚𝑑2𝑥𝑑𝑡2F = md^2xdt^2代入回复力表达式:
𝑚𝑑2𝑥𝑑𝑡2=−𝑘𝑥md^2xdt^2 = -kx整理得:
𝑑2𝑥𝑑𝑡2+𝑘𝑚𝑥=0d^2xdt^2 + kmx = 0令 𝜔 =√𝑘𝑚omega = k/m,即得简谐振动的运动方程.
解的形式
该方程的通解为:
𝑥(𝑡)=𝐴cos(𝜔𝑡+𝜙)x(t) = A cos(omega t + phi)其中:
- 𝐴A 为振幅,振子离开平衡位置的最大位移的绝对值,由初始条件决定.是求解上述微分方程时引入的常量;
- 𝜙phi 为初相位,初始时刻振子的相位,由初始条件决定.是求解上述微分方程时引入的常量;
- 𝜔omega 为角频率(圆频率),在简谐振动的几何表示中,表示旋转矢量的角速度;
- 𝑥x 为位移,表示振子相对于平衡位置的偏移量.
例 1
设一个质量为 𝑚m 的物体连接在弹簧上,弹簧劲度系数为 𝑘k.求该系统的角频率 𝜔omega 及其运动方程.
解:
易得运动方程为:
𝑑2𝑥𝑑𝑡2+𝑘𝑚𝑥=0d^2xdt^2 + kmx = 0 角频率:
𝜔=√𝑘𝑚omega = km 因此,运动方程为:
𝑥(𝑡)=𝐴cos(√𝑘𝑚𝑡+𝜙)x(t) = A cosleft(kmt + phiright) 其中 𝐴A 和 𝜙phi 由初始条件决定.
例 2
设一个质量为 𝑚m 的物体连接在摆绳上,摆长为 𝑙l.求该单摆系统的角频率 𝜔omega 及其运动方程.
解:
𝑀=−𝑚𝑔𝑙sin𝜃=𝐼𝛽=𝑚𝑙2𝑑2𝜃𝑑𝑡2M = -mglsintheta = Ibeta = ml^2 d^2dt^2 对于小角度摆动,sin𝜃 ≈𝜃sintheta ~ theta,单摆的运动方程为:
𝑑2𝜃𝑑𝑡2+𝑔𝑙𝜃=0d^2dt^2 + gltheta = 0 角频率:
𝜔=√𝑔𝑙omega = gl 因此,运动方程为:
𝜃(𝑡)=𝐴cos(√𝑔𝑙𝑡+𝜙)theta(t) = A cosleft(glt + phiright) 其中 𝐴A 和 𝜙phi 由初始条件决定.
小角度近似
在单摆的例子中,我们使用了小角度近似 sin𝜃 ≈𝜃sintheta ~ theta.这一近似仅在摆动角度较小时成立(通常小于约 10 度).对于较大的摆动角度,运动方程将变得非线性,不能再用简谐振动的形式描述.
简谐振动的判断
并非所有往复振动都是简谐振动.简谐振动要求回复力与位移成正比且方向相反,即满足 胡克定律的形式.对于非线性系统或大幅度振动,运动方程可能不再是线性的,不能用简谐振动的形式描述.我们往往也可以通过得到相应的 特征微分方程 来判断振动是否为简谐振动.
周期与频率
振动的周期 𝑇T 和频率 𝑓f 分别为:
𝑇=2𝜋𝜔,𝑓=1𝑇=𝜔2𝜋.T = 2omega, quad f = 1T = 2pi.相位
振动的相位 𝜃theta 定义为:
𝜃=𝜔𝑡+𝜙.theta = omega t + phi.相位表示振动在某一时刻所处的位置状态.规定 0 <𝜃 <2𝜋0 < theta < 2pi.
对于两个同频简谐振动,其相位差 Δ𝜃Delta theta 为:
Δ𝜃=𝜃1−𝜃2=𝜙1−𝜙2.Delta theta = theta_1 - theta_2 = phi_1 - phi_2.若 Δ𝜃 =0Delta theta = 0,则两振动 同相,若 Δ𝜃 =𝜋Delta theta = pi,则两振动 反相. 若 0 <Δ𝜃 <𝜋0 < Delta theta < pi,则称 𝑥1x_1 比 𝑥2x_2 超前,若 𝜋 <Δ𝜃 <2𝜋pi < Delta theta < 2pi,则称 𝑥1x_1 比 𝑥2x_2 落后.
x, v, a 三者的相位关系
𝑥(𝑡)=𝐴cos(𝜔𝑡+𝜙)x(t) = A cos(omega t + phi)𝑣(𝑡)=−𝐴𝜔cos(𝜔𝑡+𝜙+𝜋2)=𝐴𝜔sin(𝜔𝑡+𝜙)v(t) = -Aomega cos(omega t + phi + 2) = Aomega sin(omega t + phi)𝑎(𝑡)=−𝐴𝜔2cos(𝜔𝑡+𝜙+𝜋)=−𝐴𝜔2cos(𝜔𝑡+𝜙)a(t) = -Aomega^2 cos(omega t + phi + pi) = -Aomega^2 cos(omega t + phi)由此可见,速度 𝑣v 比位移 𝑥x 超前 𝜋2pi/2,加速度 𝑎a 与位移 𝑥x 反相.你可以想想为什么会有这样的相位关系.
能量分析
动能与势能
简谐振动的能量包括动能和势能,其表达式分别为:
𝐸𝑘=12𝑚𝑣2=12𝑚𝜔2𝐴2sin2(𝜔𝑡+𝜙)E_k = 12mv^2 = 12momega^2A^2sin^2(omega t + phi)在一个周期内的平均值:
―――𝐸𝑘=1𝑇∫𝑇0𝐸𝑘𝑑𝑡=14𝑘𝐴2E_k = 1Tint_0^T E_k dt = 14kA^2𝐸𝑝=12𝑘𝑥2=12𝑘𝐴2cos2(𝜔𝑡+𝜙)E_p = 12kx^2 = 12kA^2cos^2(omega t + phi)在一个周期内的平均值:
―――𝐸𝑝=1𝑇∫𝑇0𝐸𝑝𝑑𝑡=14𝑘𝐴2E_p = 1Tint_0^T E_p dt = 14kA^2总能量守恒
总能量为动能与势能之和:
𝐸=𝐸𝑘+𝐸𝑝=12𝑘𝐴2=12𝑚𝜔2𝐴2E = E_k + E_p = 12kA^2 = 12momega^2A^2总能量在振动过程中保持不变,动能和势能在振动过程中相互转化.
简谐振动的几何表示
简谐振动可以通过旋转矢量的投影来几何表示.设一个长度为 𝐴A 的旋转矢量以角速度 𝜔omega 绕原点逆时针旋转,其在水平轴上的投影即为简谐振动的位移 𝑥(𝑡)x(t):
简正模
在多自由度系统中,线性振动可以分解为若干个独立的简正模(Normal Modes).每个简正模对应一个特定的频率和振动模式,系统的总振动可以看作这些简正模的叠加.
用复数表示简谐振动
简谐振动
𝑥(𝑡)=𝐴cos(𝜔𝑡+𝜑)x(t) = A cos(omega t + varphi)也可以用复数的实部和虚部表示:
̃𝑥(𝑡)=𝐴𝑒𝑖(𝜔𝑡+𝜑) x(t)=Ae^i(omega t+varphi)其中,𝐴A 为振幅,𝜔omega 为角频率,𝜑varphi 为初相位.
上式的右端又可以写为 (𝐴𝑒𝑖𝜑)𝑒𝑖𝜔𝑡 =̃𝐴𝑒𝑖𝜔𝑡(Ae^ivarphi)e^iomega t=Ae^iomega t,其中 𝐴𝑒𝑖𝜑 =̃𝐴Ae^ivarphi=A 是一个复常数,表示振动的初始状态,称为 复振幅.
̃𝑥(𝑡)=̃𝐴𝑒𝑖𝜔𝑡x(t)=Ae^iomega t如果我们取上式的实部,就得到了简谐振动的标准形式:
𝑥(𝑡)=Re(̃𝑥(𝑡))=Re(̃𝐴𝑒𝑖𝜔𝑡)=𝐴cos(𝜔𝑡+𝜑)x(t)=Re(x(t))=Re(Ae^iomega t)=Acos(omega t+varphi)用 ̃𝑥(𝑡)x(t) 表示简谐振动,则速度和加速度为
⎧{ { {⎨{ { {⎩̃𝑣(𝑡)=𝑑̃𝑥(𝑡)𝑑𝑡=𝑖𝜔̃𝐴𝑒𝑖𝜔𝑡=𝑖𝜔̃𝑥(𝑡)̃𝑎(𝑡)=𝑑̃𝑣(𝑡)𝑑𝑡=𝑑2̃𝑥(𝑡)𝑑𝑡2=−𝜔2̃𝐴𝑒𝑖𝜔𝑡=−𝜔2̃𝑥(𝑡)cases v(t)=dx(t)dt=iomega Ae^iomega t=iomega x(t) a(t)=dv(t)dt=d^2x(t)dt^2=-omega^2 Ae^iomega t=-omega^2 x(t) cases例题
已知一个线性三原子分子 𝐴2𝐵A_2B 的纵向振动模型:质量为 𝑚m 的两个 𝐴A 原子位于两端,质量为 𝑀M 的 𝐵B 原子位于中间;相邻原子间由劲度系数均为 𝑘k 的弹簧连接(𝐴A–𝐵B 与 𝐵B–𝐴A).求简正模与对应角频率.
解: 设三原子沿直线位移为 𝑥1,𝑥2,𝑥3x_1,x_2,x_3(分别对应 𝐴,𝐵,𝐴A,B,A),取平衡位置为零.
弹簧形变为 𝑥1 −𝑥2x_1-x_2 与 𝑥3 −𝑥2x_3-x_2,于是运动方程为
⎧{ {⎨{ {⎩𝑚¨𝑥1=−𝑘(𝑥1−𝑥2),𝑀¨𝑥2=−𝑘(𝑥2−𝑥1)−𝑘(𝑥2−𝑥3)=−𝑘(2𝑥2−𝑥1−𝑥3),𝑚¨𝑥3=−𝑘(𝑥3−𝑥2).cases mddot x_1=-k(x_1-x_2), Mddot x_2=-k(x_2-x_1)-k(x_2-x_3)=-k(2x_2-x_1-x_3), mddot x_3=-k(x_3-x_2). cases 取简正模形式 𝑥𝑗(𝑡) =𝑎𝑗𝑒𝑖𝜔𝑡x_j(t)=a_j e^iomega t,得到代数方程
⎧{ {⎨{ {⎩𝑚𝜔2𝑎1=𝑘(𝑎1−𝑎2),𝑀𝜔2𝑎2=𝑘(2𝑎2−𝑎1−𝑎3),𝑚𝜔2𝑎3=𝑘(𝑎3−𝑎2).cases momega^2 a_1=k(a_1-a_2), Momega^2 a_2=k(2a_2-a_1-a_3), momega^2 a_3=k(a_3-a_2). cases 利用对称性分两类:
整体平移模:𝑎1 =𝑎2 =𝑎3a_1=a_2=a_3,此时弹簧不伸长,𝜔0 =0omega_0=0.
反对称伸缩模:𝑎1 = −𝑎3a_1=-a_3,由对称性得 𝑎2 =0a_2=0.
代入第一式得 𝑚𝜔2𝑎1 =𝑘𝑎1momega^2 a_1=k a_1,因此
𝜔−=√𝑘𝑚,(𝑎1,𝑎2,𝑎3)∝(1,0,−1).omega_- = km,quad (a_1,a_2,a_3)propto(1,0,-1). - 对称伸缩模:𝑎1 =𝑎3 =𝑎a_1=a_3=a,𝑎2 =𝑏a_2=b.
由第一式 𝑚𝜔2𝑎 =𝑘(𝑎 −𝑏)momega^2 a=k(a-b),第二式 𝑀𝜔2𝑏 =2𝑘(𝑏 −𝑎)Momega^2 b=2k(b-a).
联立消去 𝑏b 可得
𝜔2+=𝑘𝑚+2𝑘𝑀,(𝑎1,𝑎2,𝑎3)∝(1,−2𝑚𝑀,1).omega_+^2=km+2kM,quad (a_1,a_2,a_3)proptoleft(1,-2mM,1right). 因此该分子有一个零频平移模与两个非零振动模 𝜔−,𝜔+omega_-,omega_+.
例题:受迫振动中「功率共振」的频率
对线性受迫振动
¨𝑥+2𝛾˙𝑥+𝜔20𝑥=𝐹0𝑚cos(Ω𝑡),ddot x+2gammadot x+omega_0^2x=F_0mcos(Omega t), 证明稳态下平均耗散功率 ⟨𝑃⟩langle Prangle 取最大值时的驱动频率为 Ω =𝜔0Omega=omega_0.
解: 稳态响应 𝑥 =𝐴cos(Ω𝑡−𝛿)x=Acos(Omega t-delta).
阻尼耗散功率 𝑃𝑑 =𝑏˙𝑥2P_d=bdot x^2,周期平均
⟨𝑃⟩=𝑏⟨˙𝑥2⟩=𝑏⋅12𝐴2Ω2.langle Prangle=blangledot x^2rangle=b*frac12A^2Omega^2. 用 𝑏 =2𝑚𝛾b=2mgamma 与
𝐴2=(𝐹0/𝑚)2(𝜔20−Ω2)2+(2𝛾Ω)2A^2=(F_0/m)^2(omega_0^2-Omega^2)^2+(2gammaOmega)^2 得
⟨𝑃⟩∝Ω2(𝜔20−Ω2)2+(2𝛾Ω)2.langle Pranglepropto Omega^2(omega_0^2-Omega^2)^2+(2gammaOmega)^2. 对 ΩOmega 求极值可得最大值发生在 Ω =𝜔0Omega=omega_0.
阻尼振动
运动方程
当系统受到阻尼力作用时,振动会逐渐衰减,其运动方程为:
𝑑2𝑥𝑑𝑡2+2𝛽𝑑𝑥𝑑𝑡+𝜔2𝑥=0d^2xdt^2 + 2beta dxdt + omega^2 x = 0其中 𝛽beta 为阻尼系数,描述阻尼的强弱.
解的形式
根据阻尼的大小,解的形式不同:
- 欠阻尼振动(𝛽 <𝜔beta < omega):
𝑥(𝑡)=𝐴𝑒−𝛽𝑡cos(𝜔𝑑𝑡+𝜙)x(t) = A e^-beta t cos(omega_d t + phi)其中 𝜔𝑑 =√𝜔2−𝛽2omega_d = omega^2 - beta^2 为阻尼振动的角频率.
图例

- 临界阻尼(𝛽 =𝜔beta = omega):
𝑥(𝑡)=(𝐶1+𝐶2𝑡)𝑒−𝛽𝑡x(t) = (C_1 + C_2 t)e^-beta t图例

- 过阻尼(𝛽 >𝜔beta > omega):
𝑥(𝑡)=𝐶1𝑒𝑟1𝑡+𝐶2𝑒𝑟2𝑡x(t) = C_1 e^r_1 t + C_2 e^r_2 t其中 𝑟1r_1 和 𝑟2r_2 为两个负实根.
图例

能量衰减
欠阻尼振动中,系统的总能量随时间指数衰减:
𝐸(𝑡)=𝐸0𝑒−2𝛽𝑡E(t) = E_0 e^-2beta t
受迫振动
运动方程
当系统受到周期性外力作用时,其运动方程为:
𝑑2𝑥𝑑𝑡2+2𝛽𝑑𝑥𝑑𝑡+𝜔2𝑥=𝐹0cos(Ω𝑡)d^2xdt^2 + 2beta dxdt + omega^2 x = F_0 cos(Omega t)其中:
- 𝐹0F_0 为外力的振幅;
- ΩOmega 为外力的角频率.
稳态解
稳态解为:
𝑥(𝑡)=𝐴cos(Ω𝑡−𝛿)x(t) = A cos(Omega t - delta)其中:
- 𝐴 =𝐹0√(𝜔2−Ω2)2+(2𝛽Ω)2A = F_0/(omega^2 - Omega^2)^2 + (2betaOmega)^2 为稳态振幅;
- 𝛿 =arctan(2𝛽Ω𝜔2−Ω2)delta = arctanleft(2betaOmega/omega^2 - Omega^2right) 为相位差.
共振现象
当驱动力频率接近系统的固有频率时,即 Ω ≈𝜔Omega ~ omega,系统的振幅达到最大值,称为 共振.
共振条件下的振幅为:
𝐴res=𝐹02𝛽𝜔A_res = F_02betaomega应用
线性振动在工程、物理和日常生活中有广泛的应用,例如:
- 钟摆的运动:简谐振动的经典例子;
- 弹簧 - 质量系统:描述机械振动;
- 电路中的交流电振荡:电感 - 电容回路的振荡;
- 建筑物的抗震设计:利用阻尼减小振幅.
{"differential-equation": {"name": "微分方程", "english": "Differential equation", "summary": "将未知函数及其导数联系起来的方程。", "href": "../../../math/differential-equations/ode-intro/"}}
本页面最近更新:2026/10/3 14:57:39,更新历史
发现错误?想一起完善? 在 GitHub 上编辑此页!
本页面贡献者:Leafuke, 匿名同学
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供,附加条款亦可能应用