本文的pdf文件：link

       设计计算机程序，根据滤波器的主要技术指标设计IIR数字巴特沃斯和切比雪夫低通，高通，带通和带阻滤波器；绘制滤波器的幅频特性和相频特性曲线，验证滤波器的设计结果是否达到设计指标要求；画出数字滤波器的直接型，级联型，并联型结构信号流图。

       指输入输出均为数字信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。因为数字滤波器通过数值运算实现滤波，所以数字滤波器具有如下优点：精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配以及可实现模拟滤波器无法实现的特殊功能等。        从滤波特性上来说可分为低通、高通、带通、带阻和全通。其各自的幅频响应如下图所示：

       N阶巴特沃斯低通滤波器的幅度平方函数为： ∣ H a ( j Ω ) ∣ 2 = 1 / ( 1 + ( Ω / Ω c ) 2 N ) , Ω c 为 3 d B 截 止 频 率 |H_a (jΩ)|^2=1/(1+(Ω/Ω_c )^2N ) ,Ω_c 为3dB截止频率 ∣Ha​(jΩ)∣2=1/(1+(Ω/Ωc​)2N),Ωc​为3dB截止频率        其幅频响应与阶数N的关系如下图所示：

       1.根据技术指标 Ω p 、 α p 、 Ω s 和 α s Ω_p 、α_p 、Ω_s和α_s Ωp​、αp​、Ωs​和αs​，用(1)式求出滤波器的阶数N。

        2.按照(2)式，求出归一化极点 p k p_k pk​，将 p k p_k pk​代入(3)式，得到归一化低通原型系统函数 G a ( p ) G_a (p) Ga​(p)。也可以根据阶数N直接查表（教材P_172）得到 p k p_k pk​和 G a ( p ) G_a (p) Ga​(p)。

       3.将 G a ( p ) G_a (p) Ga​(p)去归一化。将 p = s / Ω c p=s/Ω_c p=s/Ωc​代入 G a ( p ) G_a (p) Ga​(p)，得到实际的滤波器系统函数 H a ( s ) = G ( p ) ∣ ( p = s / Ω c ) H_a (s)=G(p)|(p=s/Ω_c ) Ha​(s)=G(p)∣(p=s/Ωc​)

       切比雪夫I型滤波器的平方幅度响应为：

       切比雪夫II型滤波器的平方幅度响应为：        各自的幅频响应曲线如下图所示：

       1）确定技术指标参数 Ω p 、 α p 、 Ω s 和 α s Ω_p 、α_p 、Ω_s和α_s Ωp​、αp​、Ωs​和αs​。        2）求滤波器阶数N和参数ε。归一化边界频率为 λ p = 1 ， λ s = Ω s ⁄ Ω p λ_p=1，λ_s=Ω_s⁄Ω_p λp​=1，λs​=Ωs​⁄Ωp​        其中

       3）求归一化系统函数 G a ( p ) G_a (p) Ga​(p)。为求 G a ( p ) G_a (p) Ga​(p)，先按照(4)式求出归一化极点 p k ， k = 1 , 2 , … , N . p_k，k=1,2,…,N. pk​，k=1,2,…,N.

       将极点 p k p_k pk​代入式(5)，得到：        4）将 G a ( p ) G_a (p) Ga​(p)去归一化，得到实际的 H a ( s ) H_a (s) Ha​(s) H a ( s ) = G a ( p ) ∣ ( p = s ⁄ Ω p ) H_a (s)=G_a (p)|(p=s⁄Ω_p ) Ha​(s)=Ga​(p)∣(p=s⁄Ωp​)

       根据所要设计滤波器的响应指标就可以求得响应的滤波器参数，然后在进行频域转换得到数字滤波器。以下便是设计IIR滤波器的流程：