Chapter7 波

Chapter7 波

6.1 机械波的产生和传播

波长、频率和波速(常识)

uu为波速,TT为周期

波速u=λT=λν角频率ω=2πT频率ν=1T=uλ波长λ=uT=uν\text{波速} \, u = \dfrac{\lambda}{T} = \lambda \nu \quad \text{角频率} \, \omega = \dfrac{2\pi}{T} \\ \text{频率} \, \nu =\dfrac{1}{T}=\dfrac{u}{\lambda} \quad 波长\lambda = uT = \dfrac{u}{\nu}

波的几何描述

  • 波(阵)面:振动相位相同的点连成的面。
  • 波前:最前面的波阵面。
  • 波线(波射线):波的传播方向。在各向同性媒质中,波线恒与波面垂直。

6.2 平面简谐波

原点O的振动式

y=Acos(ωt+φ)y = Acos(\omega t + \varphi)

image-20240106145444583

假设波以速度uu向x轴正方向传播

txut-\dfrac{x}{u}时刻点O的相位    \implies t时刻点P的相位

ttxut \rightarrow t-\dfrac{x}{u}

直接就可以得到波函数(行波表达式):

y=Acos[ω(txu)+φ]y = A cos[\omega(t-\dfrac{x}{u})+\varphi]

可以改写成平面简谐波的一般形式

y=Acos(ωtkx+φ)y = Acos(\omega t -kx+\varphi)

其中k=ωu=2πλk = \dfrac{\omega}{u} = \dfrac{2\pi}{\lambda} (角波数)

其他表示形式

y=Acos[2π(tTxλ)+φ]y=Acos[2π(νtxλ)+φ]y = A \cos\left[2\pi(\dfrac{t}{T}-\dfrac{x}{\lambda})+\varphi\right] \\ y = A \cos\left[2\pi (\nu t - \frac{x}{\lambda}) + \varphi\right]

假设波以速度uu向x轴负方向传播

很简单,改个符号就是

y=Acos[ω(t+xu)+φ]y = A cos[\omega(t+\dfrac{x}{u})+\varphi]

任意点行波法

经常使用任意点行波法来计算波函数

给出任意点Q的振动式

yQ=Acos(ωt+φ)y_Q = Acos(\omega t + \varphi)

  1. 若波向x轴正方向传播
image-20240106151120518

ttxxQuy=Acos[ω(txxQu)+φ]t \rightarrow t-\dfrac{x-x_Q}{u} \\ y = A cos[\omega(t-\dfrac{x-x_Q}{u})+\varphi]

  1. 若波向x轴负方向传播

    tt+xxQuy=Acos[ω(t+xxQu)+φ]t \rightarrow t+\dfrac{x-x_Q}{u} \\ y = A cos[\omega(t+\dfrac{x-x_Q}{u})+\varphi]

行波法物理意义

相距Δx\Delta x两点的相位差

Δϕ=ωuΔx=2πTΔxu=2πΔxλ\Delta \phi = -\dfrac{\omega}{u}\Delta x = -\dfrac{2\pi}{T} \dfrac{\Delta x}{u} = -2\pi \dfrac{\Delta x}{\lambda}

Δx=x2x1\Delta x = x_2 - x_1可称为波程差

平面波动方程

具体推导过程见数理方法笔记,此处直接写结论了

对于一个线密度为 μ\mu ,绳中张力为 FF 的绳,容易得到其波动满足

2yt2=Fμ2yx2\frac{\partial^2 y}{\partial t^2} = \frac{F}{\mu} \frac{\partial^2 y}{\partial x^2}

代入平面简谐波表达式

y=Acos[ω(txu)+φ]y = A cos[\omega(t-\dfrac{x}{u})+\varphi]

可以验证:平面简谐波表达式是它的解,同时得到波动方程系数和波速间的一个关系式

u=Fμu = \sqrt{\frac{F}{\mu}}

从而波动方程可以化为 (其中 uu 为波速)

2yt2=u22yx2\frac{\partial^2 y}{\partial t^2} = u^2 \frac{\partial^2 y}{\partial x^2}

注意:这里包含一个常用结论:软绳中波速 uu 满足 u=Fμu = \sqrt{\frac{F}{\mu}}

  • 关于这个方程的解的形式更深刻的理解,参考数理方法的行波法相关内容

6.3 简谐波的能量传输

质点的能量

image-20250111212715086

质点的机械能由动能和势能组成,表示为:

dE=dEk+dEpdE = dE_k + dE_p

动能部分

由于弦作小横振动,在x轴方向上,该质点的位移很小,速度可以忽略不计,因此只考虑y轴方向上有速度,动能的微分可以表示为:

dEk=12μdx(yt)2dE_k = \frac{1}{2} \mu dx \left( \frac{\partial y}{\partial t} \right)^2

其中:

  • μ\mu 为单位长度的质量
  • yt\frac{\partial y}{\partial t} 为振动速度

势能部分

image-20250111212730855image-20250111212746548

弦的某段长度 dxdx 经过扰动后长度发生变化,扰动后的长度为:

dx[(dx)2+(dy)2]1/2=dx[1+(yx)2]12dx[1+12(yx)2]dx \to \left[ (dx)^2 + (dy)^2 \right]^{1/2} = dx \left[ 1 + \left( \frac{\partial y}{\partial x} \right)^2 \right]^{\frac{1}{2}} \approx dx \left[ 1 + \frac{1}{2} \left( \frac{\partial y}{\partial x} \right)^2 \right]

张力 FF 在弦上各点保持不变,故可以写出张力所作的功

dA=F(dx[1+12(yx)2]dx)=12Fdx(yx)2dA = F \left( dx \left[ 1 + \frac{1}{2} \left( \frac{\partial y}{\partial x} \right)^2 \right] - dx \right) = \frac{1}{2} F dx \left( \frac{\partial y}{\partial x} \right)^2

因此张力做的微小功 dAdA 即为质点的势能增量:

dEp=dA=12Fdx(yx)2dE_p = dA = \frac{1}{2} F dx \left( \frac{\partial y}{\partial x} \right)^2

机械能的总和

平面简谐波的波函数为

y=Acos(ωtkx+φ)y = A \cos(\omega t - kx + \varphi)

由此可以分别求出:

  1. 动能微分:

    dEk=12μdxω2A2sin2(ωtkx+φ)dE_k = \frac{1}{2} \mu dx \omega^2 A^2 \sin^2 (\omega t - kx + \varphi)

  2. 势能微分:

    dEp=12Fdxk2A2sin2(ωtkx+φ)dE_p = \frac{1}{2} F dx k^2 A^2 \sin^2 (\omega t - kx + \varphi)

结合波速公式 u=ωk=Fμu = \frac{\omega}{k} = \sqrt{\frac{F}{\mu}},可得动能和势能相等(并且是同步变化的)

dEk=dEpdE_k = dE_p

因此,机械能总微分为:

dE=2dEk=2dEp=dxμω2A2sin2(ωtkx+φ)dE = 2dE_k = 2dE_p = dx \mu \omega^2 A^2 \sin^2 (\omega t - kx + \varphi)

进一步化简为:

dE=dxμω2(A2y2)dE = dx \mu \omega^2 \left( A^2 - y^2 \right)

结论

  1. 质点的动能和势能大小相等,同步变化。

    image-20250111213949363
  2. 当质点处于平衡位置时(y=0)(y = 0), 其动能、势能和总机械能最大;当质点处于最大位移位置时(y=A)(|y| = A), 其动能 、势能和总机械能为零,为最小。

  3. 质点的机械能不守恒

    dEconstdE \neq \text{const}

    这是由于质点和其相邻质点的相互作用,有外力对其作功,导致其机械能不守恒。

机械波的能量

能量密度

机械波的能量密度,即机械波在单位体积中的机械能

ϵ=dESdx=2ϵk=2ϵp=ρω2A2sin2(ωtkx+φ)=ρω2(A2y2)\epsilon = \dfrac{dE}{Sdx} = 2 \epsilon_k = 2\epsilon_p = \rho \omega^2 A^2 sin^2(\omega t -kx+\varphi) = \rho\omega^2(A^2-y^2)

其中ρ=μS\rho = \dfrac{\mu}{S},为传播介质的体密度

平均能量密度

ϵ=2ϵk=2ϵp=1T0Tϵdt=12ρω2A2\overline{\epsilon} = 2 \overline{\epsilon_k} = 2 \overline{\epsilon_p} = \frac{1}{T} \int_0^T \epsilon \, dt = \dfrac{1}{2} \rho \omega^2 A^2

能流密度

image-20250111220643360

单位时间内通过垂直于波传播方向的某一面积的能量称为通过该面积的能流

定义能流PP

P=ϵSudtdt=ϵSuP = \dfrac{\epsilon Sudt}{dt} = \epsilon S u

能流密度:通过垂直于波传播方向的单位面积的能流

J=PS=ϵuJ = \dfrac{P}{S} = \epsilon u

能流密度为矢量,方向与波速方向相同

J=ϵu\vec{J} = \epsilon \vec{u}

强度(平均能流密度)

定义波的强度II为能量密度在一个周期上的平均值,单位 W/m2W/m^2

I=J=1T0TJdt=uT0Tϵdt=uϵI = \overline{J} = \frac{1}{T} \int_0^T J \, dt = \frac{u}{T} \int_0^T \epsilon \, dt = u\overline{\epsilon}

对于平面简谐波来说

I=12ρω2A2uI = \dfrac{1}{2} \rho \omega^2 A^2 u

不难发现,强度即为平均能流密度

强度的矢量投影

强度为II的波,传播方向与所看平面的夹角为α\alpha,则穿过该平面的平均能流密度(强度)为

I=IcosαI' = Icos\alpha

image-20240106161816804

拓展:从能量角度来看球面简谐波的表达式

image-20250111221530143

球面简谐波的强度(单位时间通过单位面积的能量流量)为:

I1=12ρA12ω2u,I2=12ρA22ω2uI_1 = \frac{1}{2} \rho A_1^2 \omega^2 u, \quad I_2 = \frac{1}{2} \rho A_2^2 \omega^2 u

其中:

  • A1A_1A2A_2 分别为两个球面上的振幅;
  • ρ\rho 为介质密度;
  • ω\omega 为角频率;
  • uu 为波速。
平均能流相等

在介质不吸收能量的假设下,两个球面 S1S_1S2S_2 上的平均能流相等,即:

p1=I14πr12=p2=I24πr22\overline{p_1} = I_1 4\pi r_1^2 = \overline{p_2} = I_2 4\pi r_2^2

由此可得:

A1A2=r2r1\frac{A_1}{A_2} = \frac{r_2}{r_1}

即振幅 AA 随距离 rr 的增大成反比。

球面简谐波的表达式

振动传播的球面波表达式为:

y=Arcos[ω(tru)+φ]y = \frac{A}{r} \cos\left[\omega \left(t - \frac{r}{u}\right) + \varphi \right]

其中:

  • AA 为离波源单位距离处波的振幅;
  • rr 为球面波传播距离;
  • φ\varphi 为初相位。

声波的声强级

规定

I0=1012W/m2I_0 = 10^{-12} W/m^2

声强级定义(贝尔Bel,分贝dB)为

LI=lgII0(Bel)LI=10lgII0(dB)1Bel=10dBL_I = lg\dfrac{I}{I_0} (Bel) \quad L_I = 10 lg\dfrac{I}{I_0} (dB) \\ 1Bel = 10dB

6.4 波的叠加和干涉

惠更斯原理

惠更斯原理:波阵面上每一点都可以看作是发射子波的波源,在其后的任一时刻,这些子波波阵面的包络面就是该时刻的波阵面—惠更斯原理

用惠更斯原理解释波的衍射

衍射:波传播过程中当遇到障碍物时,能绕过障碍物的边缘而传播(偏离直线传播)的现象

image-20240106163040084

用惠更斯原理解释波的反射

image-20240106163209232 $$ A'D = AA'' \implies r=i $$ 即反射角与入射角相等

用惠更斯原理解释波的折射

image-20240106163858723

可推出折射定律:

AD=u1(t2t1)=ADsiniAA=u2(t2t1)=ADsinr    sinisinr=u1u2A'D = u_1(t_2-t_1) = ADsini \\ AA'' = u_2(t_2-t_1) = ADsinr \\ \implies \dfrac{sini}{sinr} = \dfrac{u_1}{u_2}

波的干涉

相干条件:频率相同,振动方向相同,相位差恒定

波的干涉和振动的叠加类似,在相量图上也能看作振动矢量相加

image-20240106181929632

y1=A1cos(ωtkr1+φ1)y2=A2cos(ωtkr2+φ2)y=y1+y2=Acos(ωt+φ)A=A12+A22+2A1A2cosΔΦΔΦ=φ2φ1k(r2r1)φ=arctanA1sin(φ1kr1)+A2sin(φ2kr2)A1cos(φ1kr1)+A2cos(φ2kr2)y_1 = A_1cos(\omega t -kr_1+\varphi_1) \\ y_2 = A_2cos(\omega t -kr_2+\varphi_2) \\ y = y_1 + y_2 = Acos(\omega t+\varphi) \\ A = \sqrt{A_1^2+A_2^2+2A_1A_2cos\Delta \Phi} \quad \Delta \Phi = \varphi_2 - \varphi_1 - k(r_2 - r_1) \\ \varphi = arctan\dfrac{A_1sin(\varphi_1 -kr_1)+A_2sin(\varphi_2 -kr_2)}{A_1cos(\varphi_1 -kr_1)+A_2cos(\varphi_2 -kr_2)}

其中r2r1r_2 - r_1为波程差

φ1=φ2\varphi_1 = \varphi_2

  1. r1r2=±nλr_1-r_2=\pm n\lambda时,干涉相长:A=A1+A2A = A_1 + A_2
  2. r1r2=±(2n+1)λ2r_1-r_2=\pm (2n+1)\dfrac{\lambda}{2}时,干涉相消:A=A1A2A = \lvert A_1 - A_2 \rvert
  3. 其他一般情况:A1A2<A<A1+A2\lvert A_1 - A_2 \rvert <A< A_1 + A_2

提醒:善用向量叠加的几何意义,减少计算

驻波

当两列振幅相同,频率相同,振动方向相同的波以相反方向传播时,叠加形成驻波。

y1=Acos(ωtkx)y2=Acos(ωt+kx)y=y1+y2=2Acoskxcosωty_1 = Acos(\omega t -kx) \\ y_2 = Acos(\omega t +kx) \\ y = y_1 + y_2 = 2Acoskxcos\omega t

驻波各质点均做简谐振动,但振幅随位置做周期性变化,其中振幅部分即为2Acoskx2Acoskx

有些情况,正向传播和反向传播的波有相位差,算一算就行了

当然,记住三角公式

sinαcosβ=12[sin(α+β)+sin(αβ)]cosαsinβ=12[sin(α+β)sin(αβ)]cosαcosβ=12[cos(α+β)+cos(αβ)]sinαsinβ=12[cos(α+β)cos(αβ)]\begin{align*} \sin \alpha \cos \beta &= \frac{1}{2} [\sin(\alpha + \beta) + \sin(\alpha - \beta)] \\ \cos \alpha \sin \beta &= \frac{1}{2} [\sin(\alpha + \beta) - \sin(\alpha - \beta)] \\ \cos \alpha \cos \beta &= \frac{1}{2} [\cos(\alpha + \beta) + \cos(\alpha - \beta)] \\ \sin \alpha \sin \beta &= -\frac{1}{2} [\cos(\alpha + \beta) - \cos(\alpha - \beta)] \end{align*}

sinα+sinβ=2sin(α+β2)cos(αβ2)sinαsinβ=2cos(α+β2)sin(αβ2)cosα+cosβ=2cos(α+β2)cos(αβ2)cosαcosβ=2sin(α+β2)sin(αβ2)\begin{align*} \sin \alpha + \sin \beta &= 2 \sin \left(\frac{\alpha + \beta}{2}\right) \cos \left(\frac{\alpha - \beta}{2}\right) \\ \sin \alpha - \sin \beta &= 2 \cos \left(\frac{\alpha + \beta}{2}\right) \sin \left(\frac{\alpha - \beta}{2}\right) \\ \cos \alpha + \cos \beta &= 2 \cos \left(\frac{\alpha + \beta}{2}\right) \cos \left(\frac{\alpha - \beta}{2}\right) \\ \cos \alpha - \cos \beta &= -2 \sin \left(\frac{\alpha + \beta}{2}\right) \sin \left(\frac{\alpha - \beta}{2}\right) \end{align*}

驻波的特征

波腹与波节

y=y1+y2=2Acos2πxλcosωty = y_1 + y_2 = 2Acos\dfrac{2\pi x}{\lambda}cos\omega t

  • 波腹:

    cos2πxλ=1    x=±n2λ\lvert cos\dfrac{2\pi x}{\lambda} \rvert = 1 \implies x = \pm \dfrac{n}{2} \lambda

  • 波节:

    cos2πxλ=0    x=±(2n+1)4λcos\dfrac{2\pi x}{\lambda} = 0 \implies x = \pm \dfrac{(2n+1)}{4} \lambda

    image-20240106191210545

相位特征

y=y1+y2=2Acos2πxλcosωt={2Acos2πxλcosωtcos2πxλ02Acos2πxλcos(ωt+π)cos2πxλ<0y = y_1 + y_2 = 2Acos\dfrac{2\pi x}{\lambda}cos\omega t = \begin{cases} 2A \lvert cos\dfrac{2\pi x}{\lambda} \rvert cos\omega t & cos\dfrac{2\pi x}{\lambda} \geq 0 \\ 2A \lvert cos\dfrac{2\pi x}{\lambda} \rvert cos(\omega t+ \pi) & cos\dfrac{2\pi x}{\lambda}<0 \end{cases}

  • 两相邻节点间各质点的振幅部分cos2πxλcos\dfrac{2\pi x}{\lambda}的正负号相同,节点两侧各质点的振幅部分cos2πxλcos\dfrac{2\pi x}{\lambda}正负号相反
  • 两相邻节点间各质点相位相同,节点两侧各质点相位相反

驻波的能量

先看单个质点,直接上结论:单个质点的平均机械能守恒,不随时间改变,只和位置有关

dEk=12μdx(yt)2=12μdx(2Aωcos2πxλsinωt)2dEp=12Fdx(yx)2=12Fdx(4πAλsin2πxλcosωt)2dEp=μdx(Aωsin2πxλ)2dEk=μdx(Aωcos2πxλ)2dE=dEp+dEk=μdx(Aω)2=constdE_k= \frac{1}{2} \mu dx \left( \frac{\partial y}{\partial t} \right)^2 = \dfrac{1}{2}\mu dx(2A\omega cos\dfrac{2\pi x}{\lambda} sin \omega t)^2 \\ dE_p= \frac{1}{2} F dx \left( \frac{\partial y}{\partial x} \right)^2 = \dfrac{1}{2}F dx(\dfrac{4 \pi A}{\lambda} \sin\dfrac{2\pi x}{\lambda} \cos \omega t)^2 \\ \overline{dE_p} = \mu dx(A\omega sin\dfrac{2\pi x}{\lambda})^2 \\ \overline{dE_k} = \mu dx(A\omega cos\dfrac{2\pi x}{\lambda})^2 \\ \overline{dE} = \overline{dE_p} + \overline{dE_k} = \mu dx (A\omega)^2 = const

由上面的式子可以得出以下结论:

  1. 位置与能量的关系

    • 波腹:波腹只有动能

      dEp=0dEk=maxdE_p = 0 \quad dE_k = max

    • 波节:波节只有势能

      dEk=0dEp=maxdE_k = 0 \quad dE_p = max

  2. 时间与能量的关系

    • 当所有各点达到最大位移时cosωt=±1cos \omega t = \pm 1,全部能量为势能

      dEk=0dEp=maxdE_k = 0 \quad dE_p = max

    • 当所有各点达到平衡位置时cosωt=0cos \omega t = 0,全部能量为动能

      dEp=0dEk=maxdE_p = 0 \quad dE_k = max

  3. 驻波不传递能量,能量只在波腹与波节之间振荡,转移

简正模式

两端固定的张紧弦中产生的驻波

image-20240106195241588

y=y1+y2=2Asin2πxλcosωtyx=0=yx=L=0    sin2πλL=0    2πλL=nπ    λ=2Lnn=1,2,3y = y_1 + y_2 = 2Asin\dfrac{2\pi x}{\lambda}cos\omega t \\ y|_{x= 0} = y|_{x=L} = 0 \implies sin \dfrac{2\pi}{\lambda}L = 0 \\ \implies \dfrac{2\pi}{\lambda}L = n\pi \implies \lambda = \dfrac{2L}{n} \quad n =1,2,3\cdots

P.S. 上面的式子里的sin是因为半波反射,具体看下一节

从而可以得到简正频率

ν=uλ=nu2Ln=1,2,3\nu = \dfrac{u}{\lambda} = \dfrac{nu}{2L} \quad n =1,2,3\cdots

对应的驻波称为弦的简正模或固有振动

半波反射和全波反射

入射波:y1=Acos(ωt2πλx+φ1)=A1cosΦ1反射波:y2=Acos(ωt+2πλx+φ2)=A2cosΦ2入射波:y_1 = Acos(\omega t - \dfrac{2\pi}{\lambda}x +\varphi_1) = A_1 cos \Phi_1 \\ 反射波:y_2 = Acos(\omega t +\dfrac{2\pi}{\lambda}x +\varphi_2) = A_2 cos \Phi_2 \\

入射波在固定端反射时相位发生突变π\pi的现象称为半波损失

image-20240106200836046

设反射地点x=Lx=L,对于固定端反射,有半波损失

(y1+y2)x=L=0    (cosΦ1+cosΦ2)x=L=0    (Φ1Φ2)x=L=±π(y_1 + y_2)_{x=L} = 0 \\ \implies (cos \Phi_1+cos \Phi_2)_{x=L} = 0 \\ \implies (\Phi_1 - \Phi_2)_{x=L} = \pm \pi

故可以得出结论:

固定端有半波损失

(Φ1Φ2)x=L=±π(\Phi_1 - \Phi_2)_{x=L} = \pm \pi

自由端无半波损失

(Φ1Φ2)x=L=0(\Phi_1 - \Phi_2)_{x=L} = 0

显然,是否有半波损失和是固定端还是自由端反射有关:

  1. 固定端反射有半波损失,反射面处为波节。故固定端反射为半波反射
  2. 自由端反射无半波损失,反射面处为波腹。故自由端反射为全波反射
image-20240106201706834

产生驻波的常见方式——反射

一类常见题目的思路:

  1. 先写出入射波函数
  2. 看反射端是固定端还是自由端:固定端有半波反射,自由端无半波反射
  3. 确定反射端位置x=Lx = L,根据2的判断,写出反射波函数

6.5 多普勒效应

波源和接收器都静止(废话)

νs=νr\nu_s = \nu_r

波源不动,接收器动

设接收器以vrv_r(相对于介质)靠近波源

u=u+vrλ=λ    νr=uλ=u+vrλ=u+vruνsu' = u + v_r \\ \lambda' = \lambda \\ \implies \nu_r = \dfrac{u'}{\lambda'} = \dfrac{u+v_r}{\lambda} = \dfrac{u+v_r}{u} \nu_s

结论:

νr=u+vruνs\nu_r = \dfrac{u+v_r}{u} \nu_s

波源动,接收器不动

设波源以vsv_s(相对于介质)靠近接收器

λ=λvsT=(uvs)T=uvsνsu=u    νr=uλ=uuvsνs\lambda' = \lambda - v_sT = (u-v_s)T = \dfrac{u-v_s}{\nu_s} \\ u' = u \\ \implies \nu_r = \dfrac{u'}{\lambda'} = \dfrac{u}{u-v_s}\nu_s

结论:

νr=uuvsνs\nu_r = \dfrac{u}{u-v_s}\nu_s

波源和接收器都动

设接收器以vrv_r(相对于介质)靠近波源,波源以vsv_s(相对于介质)靠近接收器

λ=λvsT=(uvs)T=uvsνsu=u+vr\lambda' = \lambda - v_sT = (u-v_s)T = \dfrac{u-v_s}{\nu_s} \quad u' = u + v_r

从而得到结论

νr=u+vruvsνs\nu_r = \dfrac{u+v_r}{u-v_s} \nu_s

vsv_svrv_r不在两者连线上

νr=u+vrcosθruvscosθsνs\nu_r = \dfrac{u+v_r cos\theta_r}{u-v_s cos\theta_s} \nu_s

解题基本思路

  1. 确定波源和接收器:实际题目比较复杂,可能有多个波源和接收器,也有可能存在反射,使得一个物体既是波源又是接收器

  2. 规定速度正方向,确定vsvrv_s v_r正负号

    image-20240106210018010
  3. 确定速度vsv_svrv_r的代数值

  4. 确定波源频率vsv_s

  5. 代入公式,计算接收器频率

0%