振动和波分析

振动

弹簧振子运动方程

$$m\ddot{x}=-kx-h\dot{x}$$

其中 $h$ 为阻尼系数。

令

$${\omega }_0^2=\frac{k}{m},2\beta =\frac{h}{m}$$

于是方程变为

$$\ddot{x}+2\beta \dot{x}+{\omega }_0^{2}x=0$$

用试探解 $x=e^{rt}$ 代入，解得

$$r_1,2=-\beta \pm \sqrt{{\beta }^2-{\omega }_0^2}$$

于是方程的解可以写成

$$x=A_1{e}^{r_{1}t}+A_2{e}^{r_{2}t}$$

分为三种情况：

• $\beta <{\omega }_0$，欠阻尼情况
• $\beta ={\omega }_0$，临界阻尼情况
• $\beta >{\omega }_0$，过阻尼情况

欠阻尼振动

令 ${\omega }_f=\sqrt{w_0^2-{\beta }^2}$，得

$$x=(A_1{e}^{i{\omega }_{f}t}+A_2{e}^{-i{\omega }_{f}t})e^{-\beta t}$$

取实部得

$$x=A_0{e}^{-\beta t}\cos ({\omega }_{f}t+{\phi }_0)$$

虽然不是周期运动，但仍可看做振幅衰减的周期运动，且振幅为 $A=A_0{e}^{-\beta t}$，$t\to +\infty$ 时 $A\to 0$。

能量

速度为

$$v=\frac{dx}{dt}=-A_0{e}^{-\beta t}[\beta \cos ({\omega }_{f}t+{\phi }_0)+\omega \sin ({\omega }_{f}t+{\phi }_0)]$$

动能为

$$E_k=\frac{1}{2}m{v}^2=\frac{1}{2}m{A}_0^2{e}^{-2\beta t}[\beta \cos ({\omega }_{f}t+{\phi }_0)+{\omega }_f\sin ({\omega }_{f}t+{\phi }_0){]}^2$$

势能为

$$V=\frac{1}{2}k{x}^2=\frac{1}{2}m{\omega }_0^2{A}_0^2{e}^{-\beta t}{\cos }^2({\omega }_{f}t+{\phi }_0)=\frac{1}{2}m({\omega }_f^2+{\beta }^2)A_0^2{e}^{-\beta t}{\cos }^2({\omega }_{f}t+{\phi }_0)$$

得机械能为

$$E=E_k+V=\frac{1}{2}m{A}_0^2{e}^{-2\beta t}[{\omega }_f^2+\beta {\omega }_f\sin 2({\omega }_{f}t+{\phi }_0)+2{\beta }^2{\cos }^2({\omega }_{f}t+{\phi }_0)]$$

机械能不守恒。当$\beta <<{\omega }_0$ 时，有 ${\omega }_f\approx {\omega }_0$，于是

$$E\approx \frac{1}{2}m{\omega }_f^2{A}_0^2{e}^{-2\beta t}$$

设原本能量为 $E$，经过一个周期之后能量损失为 $\Delta E$，品质因数为

$$Q=2\pi \frac{E}{\Delta E}=2\pi \frac{1}{1-e^{-2\beta t}}$$

由 $\beta <<{\omega }_0=\frac{2\pi }{T}$，故$\beta T<<1$，故

$$Q\approx \frac{2\pi }{2\beta t}=\frac{{\omega }_0}{2\beta }$$

受迫振动

系统受外力 $F=F_0\cos \omega t$

运动方程

$$m\ddot{x}=-kx-h\dot{x}+F_0\cos \omega t$$

变为

$$\ddot{x}+2\beta \dot{x}+{\omega }_0^{2}x=f_0\cos \omega t$$

齐次方程通解为

$$x=A_0{e}^{-\beta t}\cos ({\omega }_{f}t+{\phi }_0)$$

用试探解得

$$\tilde{B}(r^2+2\beta r+{\omega }_0^2)e^{rt}=f_0{e}^{i\omega t}$$

得

$$\{\begin{array}{l}r=i\omega \\ \tilde{B}=\frac{f_0}{r^2+2\beta r+{\omega }_0^2}=\frac{f_0}{-{\omega }_0^2+2i\beta \omega +{\omega }_0^2}=\frac{f_0[({\omega }_0^2-{\omega }^2)-2i\beta \omega ]}{({\omega }_0^2-{\omega }^2)-4{\beta }^2{\omega }^2}\end{array}$$

$$\tilde{x}=(B_r+i{B}_i)e^{i\omega t}=(B_r\cos \omega t-B_i\sin \omega t)+i(B_i\cos \omega t+B_r\sin \omega t)$$

特解为实部

$$x=Re\tilde{x}=B_r\cos \omega t-B_i\sin \omega t=B\cos (\omega t-\phi )$$

其中 $B=\sqrt{B_r^2+B_i^2}=\frac{f_0}{\sqrt{({\omega }_0^2-{\omega }^2{)}^2-4{\beta }^2{\omega }^2}}$，$\tan \phi =\frac{-B_i}{B_r}=\frac{2\beta \omega }{{\omega }_0^2-{\omega }^2}$

$$x=A_0{e}^{-\beta t}\cos ({\omega }_{f}t+{\phi }_0)+B\cos (\omega t-\phi )$$

稳定状态下

$$x=B\cos (\omega t-\phi )$$

共振

$\frac{dB}{d\omega }=0$ 时，$B$ 最大，振幅最大，此时称为振幅共振。

$$\frac{dB}{d\omega }=\frac{d}{d\omega }\left(\frac{f_0}{\sqrt{({\omega }^2-{\omega }_0^2{)}^2+4{\beta }^2{\omega }^2}}\right)=2f_0[({\omega }^2-{\omega }_0^2{)}^2+4{\beta }^2{\omega }^2{]}^{-3/2}\omega ({\omega }^2-{\omega }_0^2+2{\beta }^2)=0$$

得共振频率 ${\omega }_r=\sqrt{{\omega }_0^2-{\beta }^2}\approx {\omega }_0$，若 $\beta <<{\omega }_0$

共振振幅 $B_r=\frac{f_0}{2\beta \sqrt{{\omega }_0^2-{\beta }^2}}\approx \frac{f_0}{2\beta {\omega }_0}$

共振相位 ${\phi }_r=\arctan \frac{\sqrt{{\omega }_0^2-2{\beta }^2}}{\beta }=\arctan \frac{{\omega }_0}{\beta }\approx \frac{\pi }{2}$

得运动方程为 $x\approx \frac{f_0}{2\beta \omega }\sin (\omega t)$

波

描述波的方程（右行波）

$$y=A\cos [\omega (t-\frac{x}{v})+{\phi }_0]$$

$v$ 为波的相位速度。

波前进一个波长是时间为周期，与波的振动周期相等。

转化为周期变化更直观的表达

$$y=A\cos (2\pi (\frac{t}{T}-\frac{x}{\lambda })+{\phi }_0)$$

令 $k=\frac{2\pi }{\lambda }$，有

$$y=A\cos (\omega t-\lambda x+{\phi }_0)$$

若 $v$ 不随 $\omega$ 的变化而变化，则称波是无色散是。

波动方程

$$\frac{{\partial }^{2}y}{\partial t^2}=v^2\frac{{\partial }^{2}y}{\partial x^2}$$

波动方程

柔软弦的波动方程（原本张力为 $T_0$，波动后可看做张力的水平分力为 $T_0$）

对于弦上的一小段，设两边张力为 $T_1,T_2$，与水平的夹角为 ${\theta }_1,{\theta }_2$，则有

沿 $x$ 方向无位移：$T_1\cos {\theta }_1=T_2\cos {\theta }_2=T_0$

弦受合力 $\Delta F_{ex}=T_2\sin {\theta }_2-T_1\sin {\theta }_1=T_0\tan {\theta }_2-T_0\tan {\theta }_1=T_0\frac{\partial y}{\partial x}{|}_{x+\Delta x}-T_0\frac{\partial y}{\partial x}{|}_x=T_0\frac{{\partial }^{2}y}{\partial x^2}\Delta x$

由牛二得 $\Delta F_{ex}=\rho \Delta x\frac{{\partial }^{2}y}{\partial t^2}$

得 $\frac{{\partial }^{2}y}{\partial t^2}-\frac{T_0}{\rho }\frac{{\partial }^{2}y}{\partial x^2}=0$，联系上面的式子得

$$u=\sqrt{\frac{T_0}{\rho }}$$

$T_0$ 反映介质的弹性，$\rho$ 反映介质的惯性。

纵波波动方程

单位面积和纵向相对形变

$$T=Y\frac{dy}{dx}$$

$Y$ 为弹性模量。

推得

$$u=\frac{Y}{\rho }$$

横波波动方程

有

$$T=G\frac{dy}{dx}$$

$G$ 为剪切模量。

推得

$$u=\sqrt{\frac{G}{\rho }}$$

驻波

$$\begin{gathered} y_1=A\cos (\omega t-kx+{\phi }_1) \\ {y}_2=A\cos (\omega t+kx-{\phi }_2) \end{gathered}$$

两个波合成得

$$y=y_1+y_2=2A\cos (kx+\frac{{\phi }_2-{\phi }_1}{2})\cos (\omega t+\frac{{\phi }_2+{\phi }_1}{2})$$

在依一些点振幅为零。

对于入射波 $y_1=A\cos (\omega t-kx)$，有两种情况：

自由端点，反射波为 $y_2=A\cos [\omega t-k(2l-x)]=A\cos (\omega t+kx-2kl)$

$$y=y_1+y_2=2A\cos (kx-kl)\cos (\omega t-kl)$$

对固定端点，反射波为 $y_2=A\cos (\omega t+kx-2kl-\pi )$

$$y=y_1+y_2=2A\cos (kx-kl-\frac{\pi }{2})\cos (\omega t-kl-\frac{\pi }{2})$$

此时我们称波在端点具有半波损失。

多普勒效应

$\nu =\frac{u}{\lambda }=\frac{1}{T}$

若波源和观察者相对介质运动时，观察者观测到的波速 $u^{\prime }$ 与观测到的波长 ${\lambda }^{\prime }$ 之比称为观测频率 ${\nu }^{\prime }$

$${\nu }^{\prime }=\frac{u^{\prime }}{{\lambda }^{\prime }}$$

$u_R$ 为观察者相对介质的速度，趋近波源为正

$u_S$ 为波源相对介质的速度，趋近观察者为正

$u$ 为介质中的波速

$${\nu }^{\prime }=\frac{u^{\prime }}{{\lambda }^{\prime }}=\frac{u+u_R}{\lambda -u_{S}T}=\frac{u+u_R}{u-u_S}\nu$$