实验五IIR数字滤波器设计

1、实验四实验四 IIRIIR数字滤波器设计数字滤波器设计 一、实验目的一、实验目的 掌握滤波器的作用、分类。掌握滤波器的作用、分类。 掌握掌握IIRIIR数字滤波器的原理。数字滤波器的原理。 学习数字滤波器的设计方法。学习数字滤波器的设计方法。 二、实验原理及方法二、实验原理及方法 滤波滤波 是信号处理的一种最基本而重要的技术。利用滤波是信号处理的一种最基本而重要的技术。利用滤波从复杂的信号中提取所需要的信号，抑制不需要的部从复杂的信号中提取所需要的信号，抑制不需要的部分。分。滤波器滤波器 是具有一定传输特性的信号处理装置。是具有一定传输特性的信号处理装置。 数字滤波器的工作原理数字滤波器的工作

2、原理数字滤波器是具有一定传输特性的数字信号处理装置。数字滤波器是具有一定传输特性的数字信号处理装置。它的输入和输出均为离散的数字信号，借助数字器件它的输入和输出均为离散的数字信号，借助数字器件或一定的数值计算方法，对输入信号进行处理，改变或一定的数值计算方法，对输入信号进行处理，改变输入信号的波形或频谱，达到保留信号中有用成分去输入信号的波形或频谱，达到保留信号中有用成分去除无用成分的目的。如果加上除无用成分的目的。如果加上A/DA/D、D/AD/A转换，则可以转换，则可以用于处理模拟信号。用于处理模拟信号。滤波器的分类滤波器的分类滤波器的种类很多，分类方法也不同。滤波器的种类很多，分类方法也

3、不同。 按处理的信号划分：模拟滤波器、数字滤波器按处理的信号划分：模拟滤波器、数字滤波器 按频域特性划分；低通、高通、带通、带阻。按频域特性划分；低通、高通、带通、带阻。 按时域特性划分：按时域特性划分：FIRFIR、IIRIIR 滤波器设计滤波器设计按任务要求确定按任务要求确定FilterFilter的性能指标；的性能指标；用用IIRIIR或或FIRFIR系统函数去逼近这一性能要求；系统函数去逼近这一性能要求；选择适当的运算结构实现这个系统函数；选择适当的运算结构实现这个系统函数；用软件还是用硬件实现。用软件还是用硬件实现。IIRIIR滤波器设计滤波器设计由于它的脉冲响应序列由于它的脉冲响应

4、序列h(n)是无限长的，称为无限长是无限长的，称为无限长冲激响应滤波器。冲激响应滤波器。IIR滤波器的设计根据滤波器某些性能指标要求，设计滤波器的设计根据滤波器某些性能指标要求，设计滤波器的分子和分母多项式。滤波器的分子和分母多项式。设计方法：设计方法： 模拟滤波器变换法（经典设计法）模拟滤波器变换法（经典设计法） 直接设计法直接设计法 参数模型设计法参数模型设计法 最大平滑滤波器设计法最大平滑滤波器设计法IIR设计方法比较设计方法比较借助模拟借助模拟filterfilter的设计方法（经典设计法）的设计方法（经典设计法）1 1、对设计性能指标中的频率指标进行转换使其满足模、对设计性能指标中的

5、频率指标进行转换使其满足模拟滤波器原型设计性能指标；拟滤波器原型设计性能指标；2 2、估计模拟滤波器最小阶数和边界频率。、估计模拟滤波器最小阶数和边界频率。MatlabMatlab提供提供的函数（的函数（buttord,cheb1ord,cheb2ord,ellipordbuttord,cheb1ord,cheb2ord,ellipord）；）；3 3、设计模拟低通滤波器原型，、设计模拟低通滤波器原型， MatlabMatlab提供的函数提供的函数（buttap,cheb1ap,cheb2ap,ellipapbuttap,cheb1ap,cheb2ap,ellipap）；）；4 4、由模拟低通

6、原型经频率变换获得模拟滤波器，、由模拟低通原型经频率变换获得模拟滤波器， MatlabMatlab提供的函数（提供的函数（lp2lp,lp2hp,lp2bp,lp2bslp2lp,lp2hp,lp2bp,lp2bs）5 5、将模拟滤波器离散化获得、将模拟滤波器离散化获得IIRIIR数字滤波器，数字滤波器， MatlabMatlab提供的函数（提供的函数（bilinear,impinvarbilinear,impinvar）例：用冲激响应不变法设计例：用冲激响应不变法设计ButterworthButterworth低通数字滤波低通数字滤波器，要求通带频率为器，要求通带频率为 ，通带波纹小于，通带

7、波纹小于1dB1dB，阻带在阻带在 内，幅度衰减大于内，幅度衰减大于15dB,15dB,采样周期采样周期为为Ts=0.01sTs=0.01s从题设中可知性能指标为数字指标，需将其转换为模从题设中可知性能指标为数字指标，需将其转换为模拟性能指标。拟性能指标。利用数字角频率和模拟角频率的映射关系利用数字角频率和模拟角频率的映射关系T T为采样周期。为采样周期。00.20.3T wpwp=0.2=0.2* *pi; pi; wsws=0.3=0.3* *pi;pi;rprp=1;rs=15;ts=0.01;Nn=128;=1;rs=15;ts=0.01;Nn=128;WpWp= =wp/tswp/t

8、s; Ws=; Ws=ws/tsws/ts; ; N,WnN,Wn=buttord(Wp,Ws,rp,rs,sbuttord(Wp,Ws,rp,rs,s););z,p,k=z,p,k=buttap(Nbuttap(N);); Bap,AapBap,Aap=zp2tf(z,p,k);=zp2tf(z,p,k);b,a=lp2lp(Bap,Aap,Wn);b,a=lp2lp(Bap,Aap,Wn); bz,azbz,az=impinvar(b,a,1/ts);=impinvar(b,a,1/ts);freqz(bz,az,Nn,1/ts)freqz(bz,az,Nn,1/ts)IIRIIR滤波器完

9、全设计函数滤波器完全设计函数在在MATLABMATLAB信号处理工具箱中提供了信号处理工具箱中提供了IIRIIR滤波器滤波器设计的完全工具函数，用户只要调用这些工具设计的完全工具函数，用户只要调用这些工具函数即可一次性完成设计，而不需像上面分步函数即可一次性完成设计，而不需像上面分步实现。实现。MATLABMATLAB提供的函数有：提供的函数有：ButterButter、cheby1cheby1、cheby2cheby2、ellipellip等。等。在使用这些函数设计数字滤波器时，数字频率在使用这些函数设计数字滤波器时，数字频率采用标准化频率（归一化频率）采用标准化频率（归一化频率）归一化频率

10、归一化频率频率的取值范围在频率的取值范围在0101之间，标准化频率之间，标准化频率1 1对应的数字对应的数字频率为频率为 ，对应的模拟频率为采样频率的一半。，对应的模拟频率为采样频率的一半。归一化的处理方法：归一化的处理方法： 归一化频率归一化频率 = =模拟频率模拟频率/ /采样频率一半采样频率一半巴特沃斯数字滤波器巴特沃斯数字滤波器格式：格式： b , a= b , a=butter(n,wbutter(n,wn n,ftype,ftype)其中，其中，n n为滤波器阶数，为滤波器阶数，w wn n为滤波器截止频率为滤波器截止频率ftypeftype为滤波器类型为滤波器类型 highhig

11、h为高通滤波器，截止频率为高通滤波器，截止频率w wn n stop stop为带阻滤波器，截止频率为带阻滤波器，截止频率w wn n=w1,w2=w1,w2 缺省时为低通或带通滤波器缺省时为低通或带通滤波器b,ab,a分别为滤波器传递函数分子和分母的系数向量。分别为滤波器传递函数分子和分母的系数向量。设计一个设计一个ButterworthButterworth高通数字滤波器，通带边界频率高通数字滤波器，通带边界频率为为300Hz300Hz，阻带边界频率为，阻带边界频率为200Hz200Hz，通带波纹小于，通带波纹小于1dB1dB，阻带衰减大于阻带衰减大于20dB20dB，采样频率为，采样频率

12、为1000Hz1000Hz。fs=1000;wp=300/(fs/2); ws=200/(fs/2);rp=1; rs=15; Nn=128;N,Wn=buttord(wp,ws,rp,rs);b,a=butter(N,Wn,high);freqz(b,a,Nn,fs)IIR滤波器直接设计法滤波器直接设计法经典设计法只限于几种标准的低通、高通、带通、带经典设计法只限于几种标准的低通、高通、带通、带阻滤波器，而对于具有任意形状或多频带滤波器的设阻滤波器，而对于具有任意形状或多频带滤波器的设计是无能为力的。计是无能为力的。直接设计法直接设计法 采用最小二乘法拟合给定的幅频响应，使设计的滤波采用最小

13、二乘法拟合给定的幅频响应，使设计的滤波器幅频特性逼近期望的频率特性。器幅频特性逼近期望的频率特性。MATLAB提供的工具函数提供的工具函数函数调用格式函数调用格式 b,a=yulewalk(n,f,m)说明说明 n：滤波器阶数；：滤波器阶数； f：给定的频率点向量（标准频率），第一个频点必须：给定的频率点向量（标准频率），第一个频点必须为为0，最后一个必须为，最后一个必须为1； m：和频率向量对应的理想幅值响应向量；：和频率向量对应的理想幅值响应向量； b,a：所设计的滤波器分子和分母多项式系数向量。：所设计的滤波器分子和分母多项式系数向量。例：用直接法设计一个多频带滤波器例：用直接法设计一个

14、多频带滤波器已知：已知：f=0,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1; m=0,0,1,1,0,0,1,1,0,0,0MATLAB程序程序N=10;f=0:0.1:1;m=0,0,1,1,0,0,1,1,0,0,0;b,a=yulewalk(N,f,m);h,w=freqz(b,a,128);plot(f,m,b-,w/pi,abs(h),m-)title(直接法设计多频带滤波器直接法设计多频带滤波器)三、实验内容三、实验内容用脉冲响应不变法设计一个切比雪夫型数字带通滤波用脉冲响应不变法设计一个切比雪夫型数字带通滤波器，要求通带频率器，要求通带频率 ，通带最大衰减，