参数
$$f_{obs} = f_0\,\frac{c \pm v_{o}}{c \mp v_{s}}\qquad \Delta f = f_{obs} - f_0$$
$$Ma = \dfrac{v_s}{c},\qquad \sin\theta = \dfrac{1}{Ma}\;(Ma>1)$$
实时绘制移动音源发出的圆形波面。音源前方波面密集(音调升高),后方波面拉伸(音调降低)。当音源速度超过波速时绘制马赫圆锥(冲击波),并实时计算观测频率 f′。
$$f_{obs} = f_0\,\frac{c \pm v_{o}}{c \mp v_{s}}\qquad \Delta f = f_{obs} - f_0$$
$$Ma = \dfrac{v_s}{c},\qquad \sin\theta = \dfrac{1}{Ma}\;(Ma>1)$$
多普勒效应模拟器的物理模型基于波动基本原理。静止媒质中,声以音速 $v$ 传播。当速度为 $v_s$ 的音源发声时,观测者静止的观测频率 $f'$ 由原频率 $f$ 表示为 $f' = f \frac{v}{v - v_s}$。音源越接近观测者,分母越小,频率越高;越远离,频率越低。反之,观测者以速度 $v_o$ 运动时,音源静止则 $f' = f \frac{v + v_o}{v}$,观测者越接近音源频率越高。当音源速度超过音速时,波面重叠形成马赫圆锥,产生冲击波。圆锥半顶角 $\theta$ 由 $\sin\theta = \frac{v}{v_s}$ 给出,这是超音速领域的特征现象。本工具实时计算这些公式,通过圆形波面动画、频率图表和冲击波可视化,帮助直观理解。
当声源或观测者运动时,观测到的频率会发生变化(多普勒效应)。对声速 $v$、声源频率 $f$,观测频率 $f'$ 由下式给出。
$f' = f\,\dfrac{v \pm v_o}{v \mp v_s}$
其中 $v_o$ 为观测者速度,$v_s$ 为声源速度。符号取法使得靠近时频率升高(音调变高),远离时降低(观测者靠近时分子取 $+$,声源靠近时分母取 $-$)。救护车警笛靠近时听起来音调高、驶过后变低,正是这一效应。
当声源速度超过声速($v_s > v$)时,波前相互叠加,形成圆锥状的冲击波(音爆)。以马赫数 $M=v_s/v$ 表示时,冲击波的半顶角由 $\sin\theta = 1/M$ 给出。
多普勒效应不仅发生于声波,也发生于光(电磁波);天体远离时波长被拉长,表现为红移,是宇宙膨胀的证据。其应用十分广泛,用于速度测量(多普勒雷达、测速枪)、医学超声血流测量(彩色多普勒)、气象雷达等。本模拟器可改变声源与观测者的速度,观察频率的变化。
工业实际应用例
汽车工业中,被动声纳和多普勒速度计用于车辆接近警报系统开发。具体来说,丰田和博世采用的毫米波雷达通过多普勒效应精确检测前方车辆的相对速度,实现碰撞回避支援和自适应巡航控制(ACC)。风电行业中,装在风机舱顶的激光雷达(如 Vaisala 产品)通过测量风速的多普勒频移,进行叶片桨距控制优化。
研究和教育应用
大学物理实验和声学工程课程中,本模拟器作为教材可视化波面的压缩和膨胀。特别是在航空航天学院,学生可以实时观察超音速流中马赫圆锥的形成过程,直观理解冲击波产生条件。医学院超声诊断实习前期,学生通过本工具先学习血流速度测量(多普勒法)原理,提高临床实习效率。
与CAE分析的协联及实务定位
本模拟器在CAE工具(如 ANSYS Fluent、COMSOL Multiphysics)执行的声学分析之前的阶段使用。设计阶段通过快速参数扫描了解频率移动趋势后,再转入详细的有限元分析,从而减少车辆风切音降低和飞机发动机噪声对策等开发工时。在实务中,还可用于实验测量的前期检讨和 CAE 结果的有效性确认(基准检验)。
很多人认为"观测者远离音源时观测频率必然降低",但实际上,即使观测者远离,如果音源以足够高的速度接近,频率仍可能上升。相对速度的方向和大小都会影响结果,有时会出现反直觉的情况。
有人以为"冲击波(马赫圆锥)仅在音源超音速时瞬间产生",但实际上超音速状态持续期间,圆锥形波面一直存在。本模拟器中保持音源速度在音速以上时移动,你可以看到马赫圆锥角度随速度变化实时改变。
有人说"观测周频率公式仅对音源和观测者直线运动有效",实际上存在包含角度分量的通用公式。本模拟器为简化起见仅考虑一维运动,讨论斜向多普勒效应时需要额外补正。
救护车以时速126 km/h(35 m/s)速度接近,鸣笛频率2000 Hz,观测者静止时接收频率为 f'=2000×(340+0)/(340-35)=2230 Hz。同样速度远离时 f'=2000×340/375=1813 Hz,频率差约417 Hz。飞机以马赫1.2(408 m/s)超音速飞行时,冲击波圆锥角 θ=arcsin(1/1.2)≈56°