IIR peaking filter 探秘

引言

一阶滤波器虽好，但是只能压一头，如果需要带通/带阻/peaking/notch一类的IIR滤波器，多需要高阶来帮忙了，本文先介绍一下peaking 滤波器的基本概念，peaking 滤波器是audio音频的parametric equalizer的基本单元，了解了它，我们就是可以一窥神奇和parametric equalizer了。

带通变peaking

首先我们还是给出前文得到的带通滤波器传递函数：
$$H(s)=\frac{BWs}{s^2+{\Omega }_m^2+BWs}$$
带宽和通带中频如下定义
$$BW={\Omega }_h-{\Omega }_l$$
$${\Omega }_m=\sqrt{{\Omega }_h\ast {\Omega }_l}$$
实际当中客户会提出品质因数Q来限定${\Omega }_m$与 ${\Omega }_h{\Omega }_l$的关系：
$$Q=\frac{{\Omega }_m}{{\Omega }_h-{\Omega }_l}=\frac{{\Omega }_m}{BW}$$
按照两步走的方法导出一个peaking 滤波器：

第一步：给带通增加一个增益控制旋钮（系数）

我们定义一个系数$B_0$，令带通滤波器的通带增益可调。
$$H(s)=\frac{B_0\ast BWs}{s^2+{\Omega }_m^2+BWs}$$

第二步：并联一个直通信号 y（n）= x（n）

并联直通信号的结果是传递函数,当s=0时, $|H(s)|=1+B_0$;当s=$\infty$，|H(s)|=1。截止频率${\Omega }_c$。下面是传递函数的数学表达：$$H(s)=1+\frac{B_0\ast BWs}{s^2+{\Omega }_m^2+BWs}=\frac{s^2+{\Omega }_m^2+(1+B_0)\ast BWs}{s^2+{\Omega }_m^2+BWs}$$
peaking滤波器也可以理解为带有增益的带通滤波器和直联信号的并联形式。那么notch 滤波器可以通过带阻并联直通来实现吗？答案是否定的，其实notch 滤波器和peaking滤波器可以根据$B_0$调整增益来实现。

双线性变换到数字域

数字滤波器的通带中心频率：$${\omega }_m=\frac{2\pi f_m}{f_s}$$
预设Q值，我们可以得到带宽：
$$BW=\frac{{\Omega }_m}{Q}$$
代入先整理一下我们的传递函数：$$H(s)=\frac{s^2+{\Omega }_m^2+(1+B_0)\ast \frac{{\Omega }_m}{Q}s}{s^2+{\Omega }_m^2+\frac{{\Omega }_m}{Q}s}$$
$$H(s)=\frac{s^2+{\Omega }_m^2+(1+B_0)\ast \frac{{\Omega }_m}{Q}s}{s^2+{\Omega }_m^2+\frac{{\Omega }_m}{Q}s}$$
数字和模拟频率的转换关系：$${\Omega }_m=\frac{{\omega }_m}{T}$$
双线性变换一般要求对频率作预畸变处理，采用如下公式：$${\Omega }_m=\frac{2}{T}\tan \frac{{\omega }_c}{2}$$
带入双线性变换公式：$$s=\frac{2}{T}\frac{1-z^{-1}}{1+z^{-1}}和{\omega }_m=\frac{2\pi f_m}{f_s}$$

直接变换我们得到
$$H(s)=\frac{(1+(\frac{f_m\pi }{f_s}{)}^2+\frac{(1+B_0)\pi f_m}{Q{f}_s})+2((\frac{f_m\pi }{f_s}{)}^2-1)z^{-1}+(1+(\frac{f_m\pi }{f_s}{)}^2-\frac{(1+B_0)\pi f_m}{Q{f}_s})z^{-2}}{(1+(\frac{f_m\pi }{f_s}{)}^2+\frac{\pi f_m}{Q{f}_s})+2((\frac{f_m\pi }{f_s}{)}^2-1)z^{-1}+(1+(\frac{f_m\pi }{f_s}{)}^2-\frac{\pi f_m}{Q{f}_s})z^{-2}}=\frac{Y(Z)}{X(Z)}$$
我们令$a=1+(\frac{f_m\pi }{f_s}{)}^2+\frac{\pi f_m}{Q{f}_s}$，$b=\frac{\pi f_m}{f_s}$，将上式简化如下
$$H(s)=\frac{(1+\frac{B_{0}b}{Qa})+2(\frac{b^2-1}{a})z^{-1}+(1-\frac{(2+B_0)b}{Qa})z^{-2}}{1+2(\frac{b^2-1}{a})z^{-1}+(1-\frac{2b}{Qa})z^{-2}}$$
这个系数计算比之一阶要复杂了非常多，也为频域分析和增益计算增加了难度。这时该方程只需要客户的中心频率和Q值。双线性变换、差分方程略。

频域分析

利用Octave开源的工具，我们设定采样频率$f_s$=48000Hz和频率$f_m$=6000Hz，Q=10， $B_0$=0.414 绘出$H(z)$频域响应曲线:

可以看出${\omega }_m=2\pi \frac{fm}{fs}=\frac{\pi }{4}$的中心频率被加强了，不过该滤波器的$1+B_0=G$的条件不成立了，实际增益到了27dB，还有一个就是高低频也没有完全收敛到0dB。

结论

至此，我们二阶的IIR peaking滤波器设计完成，本文只为理解Peaking滤波器的一些基本原理，实际应用中有很多成熟的peaking filter的设计方法可以参考，有兴趣的同学可以自己去寻找和比较。