【数字信号处理】第六章 IIR滤波器的设计

---

前言

很多大学生觉得上大学就是混文凭的过程, 为什么? 我想大部分学习的目的是为了考试.
我们学习路上很多人告诉我们学习是为了考高分, 但是人生有多少场考试, 要明白考试的目的是什么, 以考试为目的学习在大环境下没有错误, 但至少我们学习完考试完之后, 慢一些, 慢一些地忘去这些知识…
人的一生实在是有限的, 知识的海洋是无限的.
写这篇的时候已经年过二十, 人生可能已经过去三分之一, 太快太快了.
既然我们已经花了这么多时间学习, 就把这些知识除了应试之外的价值再体现出来.
本篇滤波器的设计已经很贴近应用了(但也只是本科阶段的初步涉猎), 希望我和大家都能珍惜自己学到的知识, 学到的知识来之不易, 不要白白看着他溜走!!!:💣

---

模拟滤波器设计

这门课不是叫数字信号处理吗, 为啥有模拟滤波器的设计?? 别慌, 这是因为设计IIR数字滤波器有两种方法, 下一章介绍的就是直接设计法(在时域或者频域直接设计IIR滤波器).

而这节要介绍的是间接设计法---------具体而言就是先根据设计指标设计相应的模拟滤波器, 再将模拟滤波器转换成数字滤波器

模拟滤波器有很多类型和现成的公式\曲线\图表. 我们这里只需要掌握巴特沃思滤波器即可

模拟滤波器的设计过程

1. 知道设计指标
2. 选择滤波器类型--------这里只掌握巴特沃思滤波器
3. 计算滤波器阶数
4. 通过查表或计算确定滤波器的系统函数$H_a(j\Omega )$

模拟滤波器(AF)设计指标

别问我为啥只放图, 问就是懒

Butterworth LP AF的设计步骤

1. 根据${\Omega }_p{\Omega }_s{\alpha }_p{\alpha }_s$计算阶数N和截止频率${\Omega }_C$
   公式有
   
   ${\lambda }_{sp}={\Omega }_s/{\Omega }_p$
   
   $k_{sp}=\sqrt{(10^{0.1\cdot {\alpha }_s}-1)/(10^{0.1\cdot {\alpha }_p}-1)}$
   
   $阶数N=\lg k_{sp}/\lg {\lambda }_{sp}取最小整数$
   
   $截止频率$$\begin{array}{rl}& {\Omega }_c={\Omega }_p{\left(10^{0.1{\alpha }_p}-1\right)}^{-\frac{1}{2N}}\\ & {\Omega }_c={\Omega }_s{\left(10^{0.1{\alpha }_s}-1\right)}^{-\frac{1}{2N}}\end{array}$

2. 求出 $p_k=e^{j\frac{2k+N+1}{2N}\pi }(k=0,1\dots N-1)$ 写出 $G(p)$ 。
   $$G(p)=\frac{1}{{\prod }_{k=0}^{N-1}\left(p-p_k\right)}$$

3. 将 $G(p)$ 去归一化，将 $p=\frac{s}{{\Omega }_c}$ 代入，得 $H_a(s)$ 。

频率变换与HP\BP AF的设计

PS:这里就不掌握带阻BS AF的设计了(￣_,￣ )

设计HP\BP 模拟滤波器的一般过程是

1. 通过频率变换公式, 将目标滤波器转换为相应的归一化低通滤波器指标
2. 设计相应的归一化低通滤波器的系统函数G(p)
3. 对G(p )进行频率变换$$p=\frac{S}{{\Omega }_C}$$得到目标滤波器的系统函数$H_a(s)$

HP AF的设计

1. 把高通滤波器的通带频率和阻带频率取倒数作为LP的参数, 损失函数不用改
   $HP:{\Omega }_p^{\prime }{\Omega }_S^{\prime }{\alpha }_{\mathrm{p}}{\alpha }_{\mathrm{s}}\Rightarrow LP:{\Omega }_p=\frac{1}{{\Omega }_p^{\prime }}{\Omega }_S=\frac{1}{{\Omega }_S^{\prime }}{\alpha }_{\mathrm{p}}{\alpha }_{\mathrm{s}}$
2. 设计LP A F $G(p)$
3. 通过取到数之后的p $$p=\frac{{\Omega }_C}{S}$$得到理想滤波器系统函数H(s)
   $HP$ $AF$ $H(s)={G(p)|}_{p=\frac{{\Omega }_C}{S}}$

BP AF的设计

1. 带通滤波器的设计指标注意多了两个
   ${\Omega }_{pu},{\Omega }_{pl},{\Omega }_{sl},{\Omega }_{su},{\alpha }_p,{\alpha }_s$
   $$\begin{array}{rl}& {\Omega }_{pu},{\Omega }_{pl},{\Omega }_{sl},{\Omega }_{su},{\alpha }_p,{\alpha }_s\\ & \Rightarrow B={\Omega }_{pu}-{\Omega }_{pl},{\Omega }_0^2={\Omega }_{pu}\cdot {\Omega }_{pl}\\ & \Rightarrow {\eta }_{pu},{\eta }_{pl},{\eta }_{sl},{\eta }_{su},{\eta }_0^2,{\alpha }_p,{\alpha }_s\\ & \Rightarrow {\lambda }_p=1,{\lambda }_s=\frac{{\eta }_{su}^2-{\eta }_0^2}{{\eta }_{su}},-{\lambda }_s=\frac{{\eta }_{sl}^2-{\eta }_0^2}{{\eta }_{pl}},{\alpha }_p,{\alpha }_s\end{array}$$
2. 设计 $LPAFG(P)$
   注意: ${\lambda }_s$ 和 $-{\lambda }_s$ 绝对值不等时, 取绝对值较小的 ${\lambda }_s$
3. BP AF通过$$p=\frac{s^2+{\Omega }_{pl}{\Omega }_{pu}}{s\left({\Omega }_{pu}-{\Omega }_{pl}\right)}$$变换得到H(s)
   
   $$H(s)={G(p)|}_{p=\frac{s^2+{\Omega }_{pl}{\Omega }_{pu}}{s\left({\Omega }_{pu}-{\Omega }_{pl}\right)}}$$

妈的这个好难记

---

IIR滤波器设计

设计方法有脉冲响应不变法和双线性变换法. 但是前者有缺陷会引起频域混叠

因此这里仅介绍双线性变换法 😢
其实是老师不要求掌握脉冲响应不变法👽

这里以HP DF的设计为例

1. 确定数字滤波器的技术指标$${\omega }_p{\omega }_s{\alpha }_p{\alpha }_s$$
2. 根据双线性变换法确定对应模拟滤波器的技术指标${\Omega }_p{\Omega }_s{\alpha }_p{\alpha }_s$
   
   $${\Omega }_k=\frac{2}{T_s}\mathrm{tg}\frac{{\omega }_k}{2}$$
3. 求对应滤波器的$H_a(s)$就和上一节间接法设计一样
4. 将模拟滤波器映射$$s=\frac{2}{T_s}\cdot \frac{1-z^{-1}}{1+z^{-1}}$$为数字滤波器
   
   $$H(z)={H_a(s)|}_{s=\frac{2}{T_s}\cdot \frac{1-z^{-1}}{1+z^{-1}}}$$
   如果是高通或者带通, 就使用上节的频率变换公式, 高通为p的倒数, 带通记忆一下嘿嘿