跳转至

振动的合成与分解

本页讨论「线性系统」中振动的叠加原理:多个振动可以直接相加;也讨论把复杂振动分解成简单振动(尤其是简谐分量)的基本方法.

注意

本页默认讨论线性范围(小振幅、线性回复力、线性阻尼等).当系统显著非线性时,叠加原理一般不再成立,参见 非线性振动.

1. 线性叠加原理

若系统满足线性微分方程

L[x](t)=f(t)L[x](t)=f(t)

其中 LL 是线性算子(例如 mx¨+bx˙+kxmddot x+bdot x+kx),则有:

  • 若 x1x_1 对应外力 f1f_1,x2x_2 对应外力 f2f_2,则 x=x1+x2x=x_1+x_2 对应外力 f=f1+f2f=f_1+f_2.
  • 特别地,齐次方程(f=0f=0)的解满足任意线性组合仍为解.

这就是叠加原理在振动中的数学基础.

2. 同方向简谐振动的合成(相同频率)

考虑两个同频简谐振动

x1(t)=A1cos⁡(ωt+φ1),x2(t)=A2cos⁡(ωt+φ2).aligned x_1(t) &= A_1cos(omega t+varphi_1), x_2(t) &= A_2cos(omega t+varphi_2). aligned

合成振动 x=x1+x2x=x_1+x_2 仍为同频简谐振动:

x(t)=Acos⁡(ωt+φ).x(t)=Acos(omega t+varphi).

2.1 相量法(矢量合成)

令复相量

x~=Aeiφ,x(t)=Re⁡{x~eiωt}.tilde x = A e^ivarphi,quad x(t)=Re tilde x e^iomega t .

则

x~=x~1+x~2=A1eiφ1+A2eiφ2.tilde x = tilde x_1+tilde x_2 = A_1e^ivarphi_1+A_2e^ivarphi_2.

由复数模与辐角得:

A2=A12+A22+2A1A2cos⁡(Δφ),Δφ=φ2−φ1.A^2=A_1^2+A_2^2+2A_1A_2cos(Deltavarphi),quad Deltavarphi=varphi_2-varphi_1.

相位满足

tan⁡φ=A1sin⁡φ1+A2sin⁡φ2A1cos⁡φ1+A2cos⁡φ2.tanvarphi=A_1sinvarphi_1+A_2sinvarphi_2A_1cosvarphi_1+A_2cosvarphi_2.

旋转矢量法

如果你觉得计算比较复杂,那么不妨通过矢量作图来理解相量法.将两个振动的振幅作为矢量的长度,相位作为矢量与水平轴的夹角,然后将两个矢量首尾相接,得到的合矢量的长度和方向即为合成振动的振幅和相位.

superposition_vector

你可以在上图中作合适的辅助线来计算合成振幅 AA 和相位 φvarphi.

2.2 特殊情形

  • 同相:Δφ=0⇒A=A1+A2Deltavarphi=0Rightarrow A=A_1+A_2(最大增强).
  • 反相:Δφ=π⇒A=|A1−A2|Deltavarphi=piRightarrow A=|A_1-A_2|(相消).
  • 正交相位:Δφ=π/2⇒A=A12+A22Deltavarphi=pi/2Rightarrow A=A_1^2+A_2^2.
例题:两简谐同频合成

已知 x1=3cos⁡(ωt)x_1=3cos(omega t),x2=4cos⁡(ωt+π/3)x_2=4cos(omega t+pi/3),求合成振幅.

解:

A2=32+42+2⋅3⋅4cos⁡π3=9+16+24⋅12=37A^2=3^2+4^2+2* 3* 4cos3=9+16+24*frac12=37

所以 A=37A=37.

3. 拍(Beats):频率接近的简谐振动合成

设

x1=Acos⁡(ω1t),x2=Acos⁡(ω2t),ω1≈ω2.aligned x_1 &= Acos(omega_1 t), x_2 &= Acos(omega_2 t),quad omega_1~omega_2. aligned

用和差化积公式 cos⁡a+cos⁡b=2cos⁡a−b2cos⁡a+b2cos a+cos b=2cosa-b/2cosa+b/2 得

x(t)=2Acos⁡(Δω2t)cos⁡(ω¯t),x(t)=2Acosleft(Delta2tright)cosleft(baromega tright),

其中 Δω=ω1−ω2Deltaomega=omega_1-omega_2,ω¯=(ω1+ω2)/2baromega=(omega_1+omega_2)/2.

  • 快速振动角频率:ω¯baromega.
  • 包络(振幅调制)角频率:|Δω|/2|Deltaomega|/2.
  • 拍频(每秒振幅最大次数):
fbeat=|f1−f2|=|Δω|2π.f_beat=|f_1-f_2|=|Deltaomega|2pi.

值得注意的是,这里的拍频并不是 |Δω|4π|Deltaomega|/4pi,从图中我们也可以看出,振幅(恒正)变化的周期是我们平常所说的振动变化周期的二分之一.

拍

例题:拍频

两音叉频率分别为 440Hz440 ,Hz 与 442Hz442 ,Hz,一起敲响,拍频是多少?

解: fbeat=|442−440|=2Hzf_beat=|442-440|=2 ,Hz.

4. 正交方向简谐振动合成:李萨如图形(Lissajous)

当两个相互垂直的简谐振动频率不同时,合成的轨迹与频率之比和两者的相位都有关系,并且图形一般比较复杂,难以用简单的几何图形描述.但是,当两者的频率成整数比时,轨迹是闭合的,这种图形称为李萨如图形.

设

{x=Acos⁡(ωxt+α),y=Bcos⁡(ωyt+β).cases x=Acos(omega_x t+alpha), y=Bcos(omega_y t+beta). cases
  • 若 ωx=ωy=ωomega_x=omega_y=omega,可消去 tt 得轨迹为椭圆;相位差决定椭圆形状.

以 α=0alpha=0 为例:

(xA)2+(yB)2−2xAyBcos⁡δ=sin2⁡δ,δ=β−α.left(xAright)^2+left(yBright)^2-2xAyBcosdelta=sin^2delta, quad delta=beta-alpha.
  • 若 ωx:ωyomega_x:omega_y 为有理比,则轨迹闭合,常用于示波器测频与测相.

李萨如图形

李萨如图形的频率比判断方式

观察李萨如图形时,可以通过数图形在水平方向和垂直方向上各自的最大交点数来判断频率比.例如,如果图形在水平方向上有 3 个最大交点,而在垂直方向上有 2 个最大交点,那么频率比就是 3:2.

同步锁频现象

当两个频率接近的振动相互作用时,可能会出现同步锁频现象,即两个振动最终以相同的频率振动.

例题:用李萨如图形测相位差(同频)

示波器 XX、YY 两通道分别输入

x=Acos⁡(ωt),y=Acos⁡(ωt+δ),0≤δ≤π.x=Acos(omega t),qquad y=Acos(omega t+delta),quad 0<= delta<= pi.

在屏幕上观察到闭合椭圆.若测得椭圆与 yy 轴交点为 y=±y0y=+/- y_0(即 x=0x=0 时的 |y||y|),求 δdelta 与 y0/Ay_0/A 的关系.

解: 令 x=0⇒cos⁡(ωt)=0⇒ωt=π/2x=0Rightarrow cos(omega t)=0Rightarrow omega t=pi/2 或 3π/23pi/2.
代入 yy:

y=Acos⁡(π2+δ)=−Asin⁡δy=Acosleft(2+deltaright)=-Asindelta

或

y=Acos⁡(3π2+δ)=Asin⁡δ.y=Acosleft(32+deltaright)=Asindelta.

因此 y0=A|sin⁡δ|y_0=A|sindelta|.在题设 0≤δ≤π0<=delta<=pi 下 sin⁡δ≥0sindelta>= 0,得到

 sin⁡δ=y0A ,δ=arcsin⁡(y0A).sindelta=y_0A ,qquad delta=arcsinleft(y_0Aright).

5. 振动的分解:傅里叶观点(最小必要版)

任意「足够好」的周期函数 x(t)x(t)(周期 TT)可展开为

x(t)=a02+∑n=1∞[ancos⁡(nω0t)+bnsin⁡(nω0t)],ω0=2πT.x(t)=a_02+sum_n=1^infleft[a_ncos(nomega_0 t)+b_nsin(nomega_0 t)right],quad omega_0=2T.

系数由正交性给出:

an=2T∫0Tx(t)cos⁡(nω0t)dt,bn=2T∫0Tx(t)sin⁡(nω0t)dt.aligned a_n &= 2Tint_0^T x(t)cos(nomega_0 t) ,dt, b_n &= 2Tint_0^T x(t)sin(nomega_0 t) ,dt. aligned

物理意义:线性系统对不同频率分量「各自响应,再相加」,这使得傅里叶分解在振动与波中极其重要.

练习

用傅里叶级数写出周期方波(幅值 ±A+/- A)的展开,并观察其只含奇次谐波.



评论