设 计学 院专 班 姓 姓 指 导 数字信号处理 实验报告
题 目： 设计IIR数字滤波器 名 称： 电子与信息工程学院 业： 电子信息工程 级： 电信133班 名： 何炳钢 学 号 13401180318 名： 朱恩良 学 号 13401180314 教 师： 郑悠
一、 实验目的
（1） （2） （3） （4）
加深对脉冲响应不变法和双线性变换法设计IIR数字滤波器基本方法的了解。 掌握使用模拟滤波器原型进行脉冲响应变换的方法。
掌握用双线性变换法设计低通、高通、带通、带阻滤波器的方法。 了解MATLAB有关脉冲响应变换和双线性变换的子函数。
二、 实验原理
1、 脉冲响应不变法的基本知识
脉冲响应不变法又称为冲激响应不变法，是将系统从s平面到z平面的一种映射方法，使数字滤波器的单位脉冲响应系列h(n)模仿模拟滤波器的冲激响应为z?h(n)。其变换关系式
aesT。
由于
esT是一个周期函数，因而s平面虚轴上每一段2?/T的线段都映射到z平面单位
圆上一周。由于重叠映射，因而冲激响应不变法是一种多值映射关系。数字滤波器的频率响应时原模拟滤波器的频率响应的周期延拓。只有当模拟滤波器的频率响应时有限带宽，且频带宽度???T??2s，才能避免数字滤波器的频率响应发生混叠的现象。由于脉冲响应不
变法只适用于限带的模拟滤波器，因此，在高频区幅频特性不等于零的高通和带阻滤波器不能采用脉冲响应不变法。
用MATLAB冲激响应不变法进行IIR数字滤波器设计的步骤如下： （1） 输入给定的数字数字滤波器设计指标；
（2） 根据公式???/T，将数字滤波器指标转换成模拟滤波器设计指标； （3） 确定模拟滤波器的最小阶数和截止频率； （4） 计算模拟低通原型滤波器的系统传递函数；
（5） 利用模拟域频率变换法，求解实际模拟滤波器的传递函数； （6） 用脉冲响应不变法将模拟滤波器转为数字滤波器。 2、 双线性变换法的基本知识
双线性变换法是将整个s平面映射到整个z平面，其映射关系为s?21?zT1?z?1?1或
z?1?sT/2
1?sT/2双线性变换法克服了脉冲响应不变法从s平面到z平面的多值映射的缺点，消除了频谱混叠现象。但其在变换过程中产生了非线性的畸变，在设计IIR数字滤波器的过程中需要进行一定的欲修正。
用MATLAB双线性变换法进行IIR数字滤波器设计的步骤与脉冲响应不变法设计的步骤基本相同：
（1）输入给定的数字数字滤波器设计指标；
（2）根据公式??2?tan()进行欲修正，将数字滤波器指标转换成模拟滤波器设计T2指标；
（3）确定模拟滤波器的最小阶数和截止频率； （4）计算模拟低通原型滤波器的系统传递函数；
（5）利用模拟域频率变换法，求解实际模拟滤波器的传递函数； （6）用脉冲响应不变法将模拟滤波器转为数字滤波器。 3、实验涉及的MATLAB子函数 Impinvar
功能：用脉冲响应不变法实现模拟到数字的滤波器变换。 调用格式：
[bd,ad] = impinvar(b,a,Fs)；将模拟滤波器系数b,a变换成数字的滤波器系数bd,ad，两者的冲激响应不变。
[bd,ad] = impinvar(b,a)；采用Fs的缺省值1Hz。 bilinear
功能：双线性变换——将s域（模拟域）映射到z域（数字域）的标准方法，将模拟滤波器变换成离散等效滤波器。
调用格式：
[numd,dend] = bilinear(num,den,Fs)；将模拟域传递函数变换为数字域传递函数，Fs为取样频率。
[numd,dend] = bilinear(num,den,Fs,Fp)；将模拟域传递函数变换为数字域传递函数，Fs为取样频率，Fp为通带截止频率。
[zd,pd,kd] = bilinear(z,p,k,Fs)；模拟域零极点增益系数变换为数字域，Fs为取样频率。 [zd,pd,kd] = bilinear(z,p,k,Fs,Fp)；模拟域零极点增益系数变换为数字域，Fs为取样频率, Fp为通带截止频率。
[Ad,Bd,Cd,Dd] = bilinear(A,B,C,D,Fs)；将模拟域状态变量系数变换到数字域，Fs为取样频率。
三、 思考题
1、四种低通模拟滤波器之间的比较。 （1）巴特沃斯低通滤波器
要求：通带截止频率fp=2KHz，通带最大衰减最小衰减
Rp?1dB，阻带截止频率fs=5KHz，阻带
A?20dB。
s
由程序运行结果和频率特性曲线可知，该结果在通阻带截止频率处能满足
Rp?1bB、
As?20dB的要求，通带波动小，Rp仅0.273dB。这个巴特沃斯低通滤波器的传递函数为
H(s)?1
(s17689)4?2.6131(s17689)3?3.4142(ss17689)2?2.6131(17689)?1（2）切比雪夫Ⅰ型低通滤波器
要求：通带截止频率fp=2KHz，通带最大衰减Rp?1dB，阻带截止频率fs=5KHz，阻带
最小衰减
As?20dB。