陷波器原理与嵌入式实现技术详解

1. 陷波器基础与核心原理

陷波器（Notch Filter）本质上是一种特殊的带阻滤波器，它的核心功能是选择性滤除信号中特定频率的噪声成分。与普通带阻滤波器不同，陷波器的阻带极窄（通常只有几Hz到几十Hz宽度），就像在频率响应曲线上"挖"了一个深坑，因此得名"陷波"。这种特性使其特别适合处理混叠在有用信号中的单频干扰，比如工频干扰（50Hz/60Hz）、设备谐振噪声等。

1.1 陷波器的数学本质

从传递函数角度看，数字陷波器可以表示为二阶IIR滤波器的标准形式：

code复制H(z) = (1 - 2cos(ω0)z⁻¹ + z⁻²) / (1 - 2rcos(ω0)z⁻¹ + r²z⁻²)

其中ω0=2πf0/fs是被抑制频率的归一化数字频率（f0是陷波频率，fs是采样率），r是极点半径（0<r<1），决定陷波带宽。r越接近1，带宽越窄。这个传递函数在z=e^(jω0)处产生一个零点，恰好抵消掉该频率成分。

注意：极点半径r的选择需要权衡。r=0.99时带宽约2Hz（fs=1kHz），但过渡带会有明显振铃效应；r=0.95带宽约10Hz，但瞬态响应更快。

1.2 模拟vs数字实现对比

传统模拟陷波器通常采用双T型RC网络，其中心频率f0=1/(2πRC)。这种方案在硬件电路中有成熟应用，但存在三个明显缺陷：

1. 频率调整需同步改变多个电阻/电容值
2. 受元件精度和温漂影响大
3. Q值难以做高（通常<50）

数字陷波器则通过算法实现，具有参数可编程、一致性好的优势。以16位ADC系统为例，频率分辨率可达fs/65536，且不受元件老化影响。现代嵌入式处理器（如STM32H743）的单周期乘加指令（MAC）能在1μs内完成二阶IIR滤波运算，满足实时性要求。

2. 陷波器的工程实现方案

2.1 IIR方案设计与优化

IIR陷波器的直接II型实现是最常用的结构，其差分方程为：

code复制y[n] = b0*x[n] + b1*x[n-1] + b2*x[n-2] - a1*y[n-1] - a2*y[n-2]

系数计算公式：

code复制b0 = (1 + r²)/2
b1 = -2cos(ω0)
b2 = (1 + r²)/2
a1 = -2rcos(ω0) 
a2 = r²

在S