波的叠加模拟器 返回
高中物理 / 音响工程

波的叠加模拟器

操纵两个正弦波的频率·振幅·位相,实时观察合成波·拍频·完全干涉。通过频谱和李萨如图形标签页多角度可视化波的关系。

参数

波 1(蓝)
波 2(橙)
预设
计算结果
拍频频率 f_beat
拍频周期 T_beat
最大振幅(强干涉)
最小振幅(弱干涉)
波形动画
频谱
李萨如图形
波形
李萨如
深化理解的对话
🙋
我在音乐课上听过"拍频"这个词,但为什么两个接近的音频会产生那种"汪汪汪"的声音呢?
🎓
可以用加法定理计算。f₁ = 440Hz、f₂ = 441Hz 这两个波叠加时,合成波的振幅会以 A·cos(2π·1·t) 的方式变化——也就是说每秒"变大然后变小"只有1次。这就是拍频频率 |f₁-f₂| = 1Hz 的本质。点击预设中的"拍频",你会看到绿色合成波的振幅波动。
🙋
我把位相差设成180°后,合成波变成了零!这就是降噪耳机吗?
🎓
正是!降噪耳机通过麦克风捕捉外部噪音,瞬时生成反向位相(180°错开)的声波来重叠。相同频率·振幅且反向位相时振幅为零,即"声音消失"。不过实时处理有毫秒级延迟——无法完美消除,但对飞机低频发动机噪音这类稳定噪音效果很好。
🙋
我看"频谱"标签页里f₁和f₂各有一条柱子。设f₁=5、f₂=10算不算"倍音"?
🎓
是的,f₂ = 2f₁ 叫"第2倍音(第2高次谐波)"。弦乐器音色之所以因乐器而异,是因为基音(f₁)和各倍音的振幅比不同。纯正弦波没有倍音。小提琴和长笛的"拉"音虽然频率相同,但音色不同就是因为倍音频谱组成差异。在倍音预设中,你会发现f₁和f₂不再"贴在一起"而是相隔一个八度。
🙋
听说李萨如图形在振动频率比为整数时会形成闭合曲线,这在实际中怎么用呢?
🎓
过去用示波器把两信号接入X-Y模式,通过目测李萨如图形判断频率比。现在控制系统测试中,还是用这方法从椭圆的倾斜度读取基准信号和测量信号的位相差。CAE中可以用来可视化2自由度振动系统的频率响应。f₁:f₂=1:2是"8字形",1:3是"扭M形"——比值复杂的话图形也变得复杂。在模拟器上逐个尝试看看。
🙋
波的叠加在结构振动分析中也会用到吗?
🎓
这是核心概念。"模态叠加法"就是把复杂结构的振动响应用各固有模态的叠加来计算的手法——是有限元法求解器的标准算法。每个模态视为"单纯正弦波",其总和就是实际位移。地震响应分析中,用"反应谱法"按模态积分输入能量来求最大响应,也是这个原理。
常见问题
理论·主要公式
有不满足叠加原理的情况吗?
非线性媒质中不适用。例如高强度激光通过非线性光学晶体时会发生"第二谐波产生(SHG)",频率变成原来的两倍。声波振幅过大也会产生非线性效应,形成冲击波。CAE中对应地盘非线性响应、大变形下的材料非线性。通常工程条件下(小位移、线性材料)都能用叠加原理。
从拍频频率的测量能得到什么?
能直接算出两个振动源的频率差。乐器调音时,用基准音叉(440Hz)和乐器音重叠,调到拍频消失(零频率)。超声波流量计中,流体内传播速度的差异(多普勒位移)表现为拍频,可转换为流量。旋转机械监测中,轴承损伤引起的微弱频率变化也能作为拍频检测出来。
傅里叶变换与波的叠加有什么关系?
傅里叶变换是"任意波形分解为正弦波叠加的操作"。这个模拟器反过来演示"两个正弦波的叠加"。频谱标签页显示的柱状图就是合成波经FFT变换后的结果。CAE中对振动数据(时间序列)做FFT可以分解成频率分量,进而分析特定固有振动频率的成分。
定常波(驻波)是怎么产生的?
相同频率·振幅的波反向传播时产生驻波。y = A·sin(kx-ωt) + A·sin(kx+ωt) = 2A·sin(kx)·cos(ωt),形成位置固定的"腹点"和"节点"图案。弦乐器弦的振动、管乐器气柱振动就是这样。两端固定的弦中,只有波长是弦长整数分之一时才能形成稳定驻波,这决定了倍音结构。
位相差φ=90°时,合成波会怎样?
当f₁=f₂、A₁=A₂时,合成振幅是 A·√2≈1.414A,位相偏移45°的正弦波。这介于完全强干涉(φ=0°,振幅2A)和完全弱干涉(φ=180°,振幅0)之间。一般地,当A₁=A₂=A时合成振幅由 2A·cos(φ/2) 给出。拖动位相滑块从0→180°,可看到合成波振幅平滑衰减。
什么是模态叠加法?
用各固有模态的叠加计算结构动力响应的分析手法。从特征值分析得到n个模态形状(φ₁,φ₂,...,φₙ)作为基底,位移展开为 u(t) = Σ qᵢ(t)·φᵢ。各模态响应 qᵢ 遵循单自由度振动方程,大幅降低计算复杂度。地震响应分析、旋转机械不平衡响应计算广泛应用,相比直接求解全自由度可显著减少计算成本。

波的叠加模拟器简介

以两个正弦波 \( y_1(t) = A_1 \sin(2\pi f_1 t + \phi_1) \) 和 \( y_2(t) = A_2 \sin(2\pi f_2 t + \phi_2) \) 的叠加为基础,实时计算合成波 \( y(t) = y_1(t) + y_2(t) \)。当频率 \( f_1 \) 和 \( f_2 \) 接近时,合成波振幅会周期变化产生拍频,其周期由 \( T_{\text{beat}} = 1/|f_1 - f_2| \) 给出。位相差 \( \Delta\phi = \phi_2 - \phi_1 \) 为 \( 0 \) 或 \( 2n\pi \) 时发生完全强干涉,为 \( \pi \) 时发生完全弱干涉。该模型允许独立操纵各波的振幅 \( A_1, A_2 \) 和频率 \( f_1, f_2 \)、位相 \( \phi_1, \phi_2 \),通过频谱显示观察频率成分,通过李萨如图形观察波的位相关系,从而直观理解波的叠加原理和干涉现象。

实际应用

工业实际应用例
汽车业中,丰田混动车用本模拟器分析发动机与电机频率干涉产生的"拍频",并在主动噪音控制(ANC)系统中优化反向位相的调整。索尼降噪耳机设计时,在逆位相合成波降噪电路优化中验证低频噪音消除效率,实现业界领先的降噪性能。

研究与教育应用
大学物理实验和音响工程讲座中作为波的叠加原理的视觉教材使用。东京大学基础物理实习将李萨如图形用于学生直观学习位相差与干涉图案的关系,频谱显示作为傅里叶变换概念导入工具而获得好评。

与CAE分析的连接和实务定位
本工具在本格CAE(ANSYS、COMSOL)之前,作为"概念设计阶段"波的基本参数(频率·振幅·位相)影响的直观理解工具使用。建筑音响设计中,在隔音墙形状最优化之前用本模拟器预先验证多重反射波的干涉条件,可降低CAE计算负荷。实务中定位为CAE结果验证和参数研究效率化的简易分析工具。

常见误解与注意事项

"频率不同的两波叠加必然产生可观测的拍频"这是常见误解,实际上要产生明显可听(或可见)的拍频,需要两波振幅相近且频率差足够小(通常数Hz以内)。频率差过大时听不到拍频,只是两个不同音高同时发声的状态。

"完全干涉(强干涉·弱干涉)只在振幅相等的两波间发生"是误解。实际上即使振幅不等,位相差为0度(或180度)也会产生干涉。但完全抵消(振幅为零)需要振幅相等,振幅不等时弱干涉不完全。

"李萨如图形的形状仅由频率比决定"是误解。初始位相差也会大幅影响图形的倾斜度和开度。仅改变位相就能让圆变成椭圆或直线,因此在解释时需同时确认频率比和位相差。

使用指南

  1. 第1波的频率(vf1)设为100Hz、振幅(va1)设为5mm,第2波的频率(vf2)设为105Hz、振幅(va2)设为5mm
  2. 在0°至180°范围内改变位相差(sa1和sa2),观察屏幕上部合成波形如何变化
  3. 在右侧面板的频谱显示中确认100Hz和105Hz处出现峰值,验证拍频频率为5Hz

具体计算例

扬声器音响试验中同时输出基准信号1000Hz和测试信号1008Hz时,拍频频率为f_beat = |1008 - 1000| = 8Hz,拍频周期为T_beat = 1/8 = 0.125秒。两波振幅均为2Pa且位相差0°(强干涉)时最大振幅为4Pa,位相差180°(弱干涉)时振幅接近0。时间轴观察会看到拍音以约125毫秒周期周期性大小变化。

实务中的注意事项

  1. 建筑物制振设计中,若固有周波数195Hz和200Hz接近,拍频现象会以5Hz周期激化共振,必须在速度反应谱中确认峰值分离
  2. 超声波清洗机以40kHz驱动时,若频率误差超过500Hz,拍频的非线性效应变得显著,空化效率下降,需要定期频率校准
  3. 手动调节位相差时,机械系中实际位相延迟随时间变化,需注意动态模拟结果与静态分析结果的差异