IIR滤波器的设计流程

---

前言

前面介绍了IIR滤波器的基本结构，这里将介绍IIR滤波器设计一般流程。

---

一、滤波器技术指标

在滤波器的设计过程中，常常用到滤波器的技术指标。理想滤波器在物理上是不可实现的，其根本原因是从一个频带过渡到另一个频带具有突变。为了物理上可实现，常常在通带和阻带之间设置一个过渡带，且频率响应在通带和阻带应给以较小的容限。一般来说，滤波器的性能要求往往以幅度响应特性的容许误差来表征。

以图中低通滤波器为例，频率响应有通带、过渡带和阻带三个范围。图中 ${\delta }_1$ 为通带容限，${\delta }_2$ 为阻带容限。${\omega }_p$和${\omega }_s$分别为通带截止频率和阻带截止频率。为了逼近理想滤波器，还必须有一个宽度不为零的过渡带 ${\omega }_s-{\omega }_p$。
滤波器的衰减函数为：
$$A(\omega )=-10lg|H(j\omega ){|}^2=-20lg|H(j\omega )|\text{\hspace{1em}(1)}$$
实际中观察滤波器的幅频响应时，常使用增益函数$G(\omega )$：
$$G(\omega )=-A(\omega )=20lg|H(j\omega )|\text{\hspace{1em}(2)}$$
在通带内，幅度响应以最大误差$\pm {\delta }_1$逼近于1：
$$1-{\delta }_1\le H(j\omega )\le 1+{\delta }_1,|\omega |\le {\omega }_p\text{\hspace{1em}(3)}$$
在阻带内，幅度响应以最大误差${\delta }_2$逼近于0：
$$|H(j\omega )|\le {\delta }_2,{\omega }_s\le |\omega |\le \pi \text{\hspace{1em}(4)}$$
实际设计滤波器时，具体指标往往由通带容许的最大衰减 $R_p$ 和阻带应达到的最小衰减 $R_s$ 给出：
$$R_p=20lg\frac{|H(j0)|}{|H(j{\omega }_p)|}=-20lg|H(j{\omega }_p)|=-20lg(1-{\delta }_1)\text{\hspace{1em}(5)}$$
$$R_s=20lg\frac{|H(j0)|}{|H(j{\omega }_s)|}=-20lg|H(j{\omega }_s)|=-20lg{\delta }_2\text{\hspace{1em}(6)}$$

$\omega$ 是用弧度表示的数字角频率，而实际应用时常给出实际频率$f$（单位Hz），因此在数字滤波器设计中还应给出采样频率$f_s$或者采样周期$T$，$\omega$ 和 $f$ 之间的换算关系可表示为 $\omega =\Omega T=2\pi fT$，其中$\Omega$为模拟角频率（$\Omega =2\pi f$，单位为rad/s）。

二、IIR数字滤波器的设计流程

IIR数字滤波器一般通过模拟滤波器进行设计，设计滤波器时需将模拟系统函数$H(s)$转换为数字滤波器系统函数$H(z)$，即将$s$平面映射到$z$平面，该映射必须满足两条基本要求：
（1） $H(z)$的频率响应要能模仿H(s)的频率响应，即$s$平面的虚轴$j\Omega$必须映射到$z$平面的单位圆$e^{j\omega }$上；
（2）$s$平面的左半平面必须映射到$z$平面单位圆的内部（因果稳定的$H(s)$映射成因果稳定的$H(z)$）。
通过模拟滤波器设计数字滤波器主要有以下步骤：
（1）将数字滤波器的设计转换为模拟滤波器的设计（主要是将数字滤波器指标转换为模拟滤波器指标）；
（2）设计满足技术指标的模拟滤波器；
（3）按一定方式将模拟滤波器转换为数字滤波器。

• 模拟滤波器主要包括贝塞尔型、巴特沃斯型、切比雪夫型和椭圆形滤波器等，模拟滤波器的设计一般通过原型低通滤波器为基础，通过复频率转换将模拟低通滤波器转换为模拟高通、带通及带阻滤波器。从模拟滤波器映射为数字滤波器的方法主要有脉冲响应不变法、双线性变换法。
  

---

总结

本篇推文对IIR数字滤波器的设计流程进行了简要介绍，后续将介绍模拟滤波器的具体设计方法。

参考资料

1. 刘纪红、孙宇舸、叶柠 等.数字信号处理原理与实践(修订版):电子信息学科基础课程系列教材[M].北京:清华大学出版社,2014.
2. 【数字信号处理(北京交通大学 陈后金)-哔哩哔哩】 https://b23.tv/cCBgSrV