IIR数字滤波器的结构与设计.docx

1、IIR数字滤波器的结构与设计.3 IR数字滤波器的结构与设计IR滤波器系统函数的极点可以在单位圆内的任何位置，实现IIR滤波器的阶次较低，所用的存储单元较少,效率高，又由于IIR数字滤波器能够保留一些模拟滤波器的优良特性，因此应用很广。设计II数字滤波器的方法主要有基于冲激响应不变法的II数字滤波器设计，基于双线性Z变换法的数字滤波器设计,数字高通、带通及带阻IR滤波器设计,基于MATLAB函数直接设计II数字滤波器。.3.基于冲激响应不变法的IIR数字滤波器设计冲激响应不变法的设计原理是利用数字滤波器的单位抽样响应序列H（z）来逼近模拟滤波器的冲激响应(t）。按照冲激响应不变法的原理,通过模

2、拟滤波器的系统传递函数G（),可以直接求得数字滤波器的系统函数（z)，其转换步骤如下:1)利用=（可由关系式推导出)，将，转换成，而,不变;) 求解低通模拟滤波器的传递函数G();3)将模拟滤波器的传递函数G(s)转换为数字滤波器的传递函数H（z)。尽管通过冲激响应不变法求取数字滤波器的系统传递函数比较方便,并具有良好的时域逼近特性,但若G(s）不是带限的，或是抽样频率不高,那么在中将发生混叠失真,数字滤波器的频率响应不能重现模拟滤波器的频率响应。只有当模拟滤波器的频率响应在超过折叠频率后的衰减很大时,混叠失真才很小，此时采样脉冲响应不变法设计的数字滤波器才能满足设计的要求,这是冲激响应不变法

3、的一个严重的缺点。【例4-1】设模拟低通巴特沃斯滤波器,通带纹波为dB，通带上限角频率=0.2，阻带下限角频率03，阻带最小衰减=15B，根据该低通模拟滤波器,利用冲激响应不变法设计响应的数字低通滤波器，并绘出设计后的数字滤波器的特性曲线。实现例4-1的MTLA程序如例程41所示，程序的运行结果如图4所示。例程4-利用冲激响应不变法设计数字低通滤波器%利用模拟巴特沃斯滤波器设计数字低通滤波器冲激响应不变法p=0.*i;w=0.*pi；Rp=；s15;T=1;性能指标Rip10(-Rp/20)；Atn=1(As）;OmgP=w*T;OmgSwsT;N,OmgC=bord(OmP,OgS，Rp，A

4、s,)； %选取模拟滤波器的阶数s,dsbtter(,Omg,s); %设计出所需的模拟低通滤波器,a=ipinvr(c,s,T）； %应用脉冲响应不变法进行转换 求得相对、绝对频响及相位、群迟延响应b,mag,pha，gd,wfreqz_m(b，a）;%下面绘出各条曲线ublt(2，2，1);plo(w/p，ag);tite（幅频特性);label(（/i);label(H(w）|）;axis(,1,0,);set（gca,XickMod,manual,XTick，0 0.0.30.51);et(ca，YTckode,anal，YTik, Atn p 1);gridsubplot(2,2，)

5、;plot(w/i,db);itl（幅频特性(d)）；lab(w（/p)；labl(B);xis（0,1,40,5);set(gc,XkMode,manul,XTick,0 0.2 0. 0.5 )；set（gca,TickMoe,mnual,YTik,-40 -As Rp 0);gidsbp(2，2,)；pt（w/pi,ha/pi);titl(相频特性);xlabe（w(/i）);ybel(pha（/i）;ais（0,1,-1,)；set(ca，XTickMode,mn,XTik,0 0.2 0.3 0.5 1);gridsublo(2，2,4);pt(w/pi，grd);title(群延时

6、）;lae(w(/pi)；ylb（ple);axi(0，1,2);se（gc，XTckMode，aual，XTik, 0.2.3 0.5)；grd%所使用的M文件函数fncindb，mag，pa，grd，w=fre_m(b,a)滤波器的幅值响应(相对、绝对)、相位响应及群延迟%age:b，ag,pha,d，w=req_(b，a)%50点对应0,p%w采样频率; b系统函数H(z)的分子项(对R,b=)% 系统函数H(z)的分母项(对FR,1),wfrz(，a,500）; %500点的复频响应ag=abs（H); %绝对幅值响应d=20lg10(ag/max（mag)）; %相对幅值响应phaa

7、ngle(H)； %相位响应gd=grpdlay(b,a,w); %群延迟响应图 利用冲激响应不变法设计的数字低通滤波器的特性曲线.32基于双线性Z变换法的IR数字滤波器设计由于的频率映射关系是根据推导的,所以使j轴每隔2T便映射到单位圆上一周，利用冲激响应不变法设计数字滤波器时可能会导致上述的频域混叠现象。为了克服这一问题,需要找到由s平面到z平面的另外的映射关系,这种关系应保证：1）s平面的整个j轴仅映射为z平面单位圆上的一周;2) 若G（s)是稳定的，由G（)映射得到的(z)也应该是稳定的；3) 这种映射是可逆的,既能由G（s)得到H()，也能由H(z）得到G（s）;） 如果G（j0）,

8、那么。双线性Z变换满足以上4个条件的映射关系，其变换公式为 (4-0)双线性Z变换的基本思路是：首先将整个s平面压缩到一条从/s到Ts的带宽为2/T的横带里，然后通过标准的变换关系将横带变换成整个z平面上去,这样就得到s平面与z平面间的一一对应的单值关系。在ALAB中，双线性变换可以通过biliear函数实现,其调用格式为:Bz，Az=bilinear(B,s）；其中B，A为模拟滤波器传递函数G()的分子分母多项式的系数向量，而Bz,Az为数字滤波器的传递函数H（）的分子分母多项式的系数向量。【例4-】使用双线性Z变换法设计一个低通数字滤波器,给定的数字滤波器的技术指标为fp=00Hz，fs=

9、3H，p3，p=20dB,抽样频率为Fs00Hz。实现例4-2的MALAB程序如例程4-2所示，程序的运行结果如图4-所示。例程4- 基于双线性变换法设计数字滤波器cleaal;f0;fs=300;Fs=1000；rp=3;rs=20;%p=2pi*f/Fs;s2*p*fs/Fs;Fs=Fs/Fs;%Frstl o finish frqe prewrig；ap=tan(wp/2);was=tan(ws/2);,n=utod（ap,wa,r,r,)%oe:s!z,p,k=butap(n）;bp,azpt(z,p，),as=lp2lp(bp,a，wap)%Note:s(2/Ts)(z1)(z+1）

10、;Ts=,hatis fs=1,s=0.;bz,az=bilnear(s，as,F/）h,w=frqz(z,az,256，s*100);plot（w，as(h);grid on;【运行结果】n = 2w =033bp = 01p 1.000 1.412.000s= 0.10a= 0000 .5951056bz 0.06754 0.067az =1.0000-.43.418图-2基于双线性Z变换法设计的数字滤波器的特性曲线.3.3数字高通、带通及带阻IIR滤波器设计除了低通数字滤波器之外,实际中还常常需要高通、带通及带阻数字滤波器，这三种数字滤波器的设计步骤如下:步骤1:将数字滤波器H(z)的技

11、术指标p和s，通过=tn（2）转变为模拟滤波器()的技术指标p和,作归一化处理后,得到p=1，sp；步骤2：化解为模拟原型滤波器G(s）的技术指标;步骤3：设计模拟原型滤波器G(p）；步骤:将G(p）转换为模拟滤波器的转移函数G(s）；步骤5:将G（s)转换成数字滤波器的转移函数H（z）,=(1)(z+1)。所谓原型滤波器是指归一化的低通滤波器。本节主要讨论通过IIR数字滤波器的原型转换设计法和IIR数字滤波器的直接设计方法来设计数字高通、带通及带阻滤波器其转换方法主要有3种:一是直接由模拟低通滤波器转换成数字高通、带通或带阻滤波器;二是先由模拟低通滤波器转换成模拟高通、带通或带阻滤波器,然后

12、再把它转换成相应的数字滤波器;三是将模拟低通滤波器先转换成数字低通滤波器,再通过变量代换变换成高通、带通或带阻滤波器。1. 直接由模拟低通滤波器转换成数字高通、带通或带阻滤波器 (1)模拟低通滤波器转换成数字高通滤波器【例4-3】利用巴特沃斯模拟滤波器设计数字高通滤波器,要求通带截止频率为0.6,通带内衰减不大于1B,阻带起始频率为0.4，阻带内衰减不小于15dB,采样周期Ts。该例的MATLB实现如例程-3所示,该数字滤波器的幅频特性曲线如图3所示。例程4-3基于巴特沃斯模拟滤波器设计数字高通滤波器W=.6*i;0.4*pi；Ap;As=5；N,wnbuttrd(Wi，Wspi,A,A);

13、%计算巴特沃斯滤波器阶次和截止频率，a=tter（N,wn,igh）； %频率变换法设计巴特沃斯高通滤波器db,mag，pha,rd，wfrqz_m（b,a); %数字滤波器响应Plt(wpi，mg);itl（数字滤波器幅频响应|H(ejOmega)图43 基于巴特沃斯模拟滤波器设计的数字高通滤波器幅频响应曲线(2)模拟低通滤波器转换成数字带通滤波器【例-】利用巴特沃斯模拟滤波器设计数字带通滤波器，要求通带上下截止频率为.4、.3，通带内衰减不大于3d，阻带上下起始频率为0.5、0.2,阻带内衰减不小于18dB。该例的MTLA实现如例程4-所示，该数字滤波器的幅频特性曲线如图4 所示。例程4-

14、基于巴特沃斯模拟滤波器设计数字带通滤波器W=0.3*p,0.4*pi；Ws0.2*pi，0.5*pi;Ap=3；As=8；N,wn=ttord（Wpp,Ws/pi,Ap,As); %计算巴特沃斯滤波器阶次和截止频率b,a=buter(N,wn，bapa)； %频率变换法设计巴特沃斯带通滤波器b,m，pha,grd,w=feqz_m(b,a); %数字滤波器响应Plot(w/pi,g);l(数字滤波器幅频响应|H(ejOega)）图44 基于巴特沃斯模拟滤波器设计的数字带通滤波器的幅频响应曲线（3) 模拟低通滤波器转换成数字带阻滤波器【例4-5】利用巴特沃斯模拟滤波器设计数字带阻滤波器,要求通带

15、上下截止频率为0和.，通带内衰减不大于1d，阻带上下起始频率为7和0.4，阻带内衰减不小于30d。该例的MATLB实现如例程4-5所示，该数字滤波器的幅频特性曲线如图4-5所示。例程4-5基于巴特沃斯模拟滤波器设计数字带阻滤波器Wp=0.pi,07*pi;Ws=0*，0.*pi；A；A=0；N,n=buttod(W/pi,Ws/pi，Ap,As）; %计算巴特沃斯滤波器阶次和截止频率,a=uer(N,wn,stop）; 频率变换法设计巴特沃斯带通滤波器d,mg,pha，grd,wreqz_m(，a); 数字滤波器响应t(w/p,a）;it(数字滤波器幅频响|H(ejOmega)| ）图4-基于

16、巴特沃斯模拟滤波器设计数字带阻滤波器2. 先由模拟低通滤波器转换成模拟高通、带通或带阻滤波器，然后再把它转换成相应的数字滤波器首先由频率变换将模拟低通原型滤波器变换成模拟高通、带通或带阻滤波器，然后根据双线性Z变换将模拟高通滤波器变换成相应的高通、带通或带阻数字滤波器。其模拟低通原型滤波器变换成模拟高通、带通或带阻滤波器的关系式如表-1所示。表4-1 模拟低通原型滤波器变换成模拟高通、带通或带阻滤波器频率变换模拟滤波器类型频率变换公式高通带通带阻3. 将模拟低通滤波器先转换成数字低通滤波器,再通过变量代换变换成高通、带通或带阻滤波器该方法中的模拟低通滤波器先转换成数字低通滤波器,可以利用冲激响

17、应不变法和双线性变换法来实现，再通过变量代换变换成数字高通、带通或带阻滤波器。【例-】试用双线性变换法设计一个带通椭圆滤波器使其幅频特性逼近于一个具有如下技术指标的模拟带通滤波器：Wp1=0Hz,Wp2=2Hz，在通带内的最大衰减为0.5dB，在阻带内的最大衰减为5dB,抽样频率为100Hz。MATAB实现如例程46所示，该数字滤波器的幅频特性曲线如图4-6所示。wp1=10;wp2=20;Fs100;=5；s50;p1=2*p*w1;w=2*wp2;w=wp-wp;Wo=sqt(wp)；z,k=eipap(7,r,rs）;A,B，C,D=zss(z,p,k)；A,B，Ct,=pbp(A，B,C，,Wo,Bw)；At，1，Ct,Dt1=bilir(At，B,C,D,Fs);n,den=ss2tf(At，Bt1，t1,Dt);H，=feqz(nu，dn);pot（W*F/（i）,abs（H);gid;xlabel（频率Hz);ylabel（幅值)；图4-6 椭圆滤波器的频率响应