实验 无限冲激响应(iir)数字滤波器的设计.docVIP

一、 实验目的 1．熟悉用双线性变换法设计IIR数字滤波器的原理与方法； 2．掌握数字滤波器的计算机仿真方法。 二、 实验原理 1、Matlab信号处理工具箱中提供了设计巴特沃思数字滤波器的函数buttord、buttap和butter，使用格式如下： （1） 用于计算巴特沃思低通滤波器的阶N和3dB截止频率的归一化值Wc (即课本符号。参数Wp 和Ws分别是数字滤波器的通带截止频率和阻带截止频率的归一化值（关于归一化），要求和，其中1表示数字频率（对应模拟频率，为采样频率）。和分别是通带最大衰减和阻带最小衰减，单位dB。当时，设计高通滤波器；当和是二元向量时，设计带通（）或带阻（）滤波器，这时返回参数也是二元向量。 （2） 用于计算巴特沃思数字滤波器系统函数中分子和分母多项式系数向量B和A。参数N 和WC分别是巴特沃思数字滤波器的阶和3dB截止频率归一化值（关于归一化）。当ftype缺省时，设计低通滤波器；当ftype=high时，设计高通滤波器；ftype=stop时，设计带阻滤波器，此时为二元向量，和分别是带阻滤波器的通带3dB下，上截止频率的归一化值（关于归一化）；为向量，且ftype缺省时，设计带通滤波器，带通的频率区间为。注意设计出的带通和带阻数字滤波器是2N阶的，这是因为带通滤波器可表示为一个N阶低通滤波器与一个N阶高通滤波器的级联。 由函数的返回向量B和A可写出数字滤波器的系统函数为： 其中，分别是向量B和A的第k个元素。 （2） [H,w]=freqz(B,A,N? ?? ? B和A分别为离散系统的系统函数分子、分母多项式的系数向量，N为正整数，返回量H则包含了离散系统频响在 0~pi范围内N个频率等分点的值，向量w则包含范围内N个频率等分点。调用中若N默认，默认值为512。 三、 实验设备 微型计算机、Matlab7.0教学版。 四、 实验内容 （1）用冲激响应不变法设计IIR数字滤波器。 （2）绘出数字滤波器在频率区间的幅频响应特性曲线； （3）对结果进行分析； 五、 实验数据 程序框图 相频特性 输入波形 输出波形