实验三IIR数字滤波器设计..docxVIP

实验三 IIR数字滤波器设计一、实验目的掌握利用脉冲响应不变法和双线性变换法设计IIR数字滤波器的原理及具体方法。2、掌握脉冲响应不变法和双线性变换法设计IIR数字滤波器的优缺点及使用范围。二、实验设备与环境计算机、MATLAB软件环境三、实验基础理论数字滤波器和模拟滤波器的一些指标为通带响应中的容限，为阻带的容限，为通带波动，为阻带衰减，且为通带波动系数，是单位为rad/s的带通截止频率，A为以dB为单位的阻带衰减参数，是单位为rad/s的阻带截止频率。参数和A与滤波器的带通波动及阻带衰减之间有如下关系模拟原型滤波器（1）巴特沃斯滤波器设计方法一：在一般设计中，先把选择为1rad/s，使频率得到归一化，归一化后的巴特沃斯滤波器系统函数的极点分布及分母多项式系数等都有现成的表格可查。然后将归一化巴特沃斯滤波器系统函数中的变量s用s/替换以后，即可得到任意的非归一化的巴特沃斯滤波器。由于巴特沃斯滤波器系统函数不存在零点，则将归一化崩塌沃斯滤波器系统函数的极点倍乘，增益倍乘以，即可实现上述由归一化巴特沃斯滤波器到非归一化巴特沃斯滤波器的变换。方法二：在MATLAB中还可以直接利用butter函数设计巴特沃斯模拟原型滤波器，butter函数既可以用于设计模拟巴特沃斯滤波器也可用于设计数字巴特沃斯滤波器，用于设计模拟巴特沃斯滤波器时，其调用格式为[b,a]=butter(N,Wc,’s’)，N表示滤波器的阶数，Wc表示截止频率，‘s’表示设计模拟滤波器。设计方程由于N比实际需要的要大，所以在或处将满足或超过指标，一般选择下式结果作为的值，（2）切比雪夫低通滤波器切比雪夫I型：方法一：MATLAB提供了函数[z,p,k]=cheb1ap(N,)来设计一个阶数为N,通带波动为的归一化（截止频率为1）切比雪夫I型模拟原型滤波器，它返回数组z和零极点p以及增益k。切比雪夫I型滤波器系统函数不存在零点，将归一化切比雪夫I型滤波器系统函数的极点倍乘以，增益倍乘以s=0时非归一化滤波器系统函数分母多项式和归一化滤波器系统函数多项式的壁纸，即可实现由归一化切比雪夫I型滤波器到非归一化切比雪夫I型滤波器的变换。方法二：在MATLAB中还可以直接利用cheby1函数设计切比雪夫I型模拟原型滤波器。cheby1函数既可用于设计模拟切比雪夫I型滤波器也可用于设计数字切比雪夫I型滤波器，用于设计模拟切比雪夫I型滤波器。设计方程：切比雪夫II型：方法一：MATLAB提供了函数[z,p,k]=cheb2ap(N,)来设计一个阶数为N，阻带衰减为的归一化（截止频率为1）切比雪夫II型模拟原型滤波器，它返回数组z和零极点p以及增益k。切比雪夫II型滤波器系统函数存在零点，可以将归一化切比雪夫II型滤波器系统函数的极点和零点倍乘，增益倍乘s=0时非归一化的有理函数和归一化的有理函数的比值，即可实现由归一化切比雪夫II型滤波器到非归一化切比雪夫II型滤波器的变换。方法二：在MATLAB中还可以直接利用cheby2函数设计切比雪夫II型模拟原型滤波器。cheby2函数既可用于设计模拟切比雪夫II型滤波器也可用于设计数字切比雪夫II型滤波器，用于设计模拟切比雪夫II型滤波器时。（3）椭圆滤波器方法一：MATLAB提供函数[z,p,k]=ellipap(N,,)来设计一个阶数为N，通带波动为，阻带衰减为的归一化（截止频率为）椭圆模拟原型滤波器，它返回数组z和零极点p以及增益k，将归一化椭圆滤波器系统函数的极点和零点倍乘，增益倍乘s=0时非归一化的有理函数和归一化的有理函数的比值，即可实现由归一化切比雪夫II型滤波器到非归一化椭圆滤波器的变换。方法二：在MATLAB中还可以直接利用ellip函数设计椭圆模拟原型滤波器。ellip函数既可用于设计模拟椭圆滤波器也可用于设计数字椭圆滤波器，用于设计模拟椭圆滤波器时。设计公式：其中，, , 模拟滤波器到数字滤波器的变换（1）脉冲响应不变法：基本原理：从时域响应出发，使数字滤波器的单位脉冲响应h(n)模仿模拟滤波器的单位冲激响应, h(n)等于的取样值。变换方法：将数字滤波器性能指标转换成响应的模拟滤波器频率指标将H(s)转换为方法一：利用residue函数和residuez函数实现脉冲响应不变法[r,p,k]=residue(b,a),[b,a]=residue(r,p,k)实现H(s)的分解[r,p,k]=residuez(b,a),[b,a]=residuez(r,p,k)实现H(z)的分解方法二：MATLAB提供了impinvar函数采用脉冲响应不变法实现模拟滤波器到数字滤波器的变换，其只用方法如下：[bz,az]=impinvar(b,a,fs) 采用脉冲响应不变法将模拟滤波器系统函数的系数向量b和a变换成为数字滤波器系统函