实验五 IIR滤波器设计.docVIP

实验五 IIR滤波器设计 一、实验目的： 1、各种模拟原型滤波器的滤波特性； 2、模数滤波器变换时的脉冲响应不变法和双线性变换法； 3、低通滤波器变换成其他类型滤波器的方法； 4、能够根据指标选择合适的原型滤波器和合适的方法设计IIR滤波器。 二、实验内容： 1、自定设计指标（通带截止频率、通带波纹、阻带截止频率、阻带衰减），选择合适的模拟原型低通滤波器和合适的设计方法（脉冲响应不变法或双线性变换法），设计符合要求的IIR滤波器；并验证设计好的滤波器是否满足要求。 2、将此低通滤波器映射为高通、带通或带阻滤波器，并验证变换后的滤波器是否满足指标。 3、求输入x(n)=[cos(w1n)+cos(w2n)]u(n)经过系统后的输出y(n)。其中w1位于通带内，w2位于阻带内，要求做一个两行两列的子图，第一个做x的时域波形，第二个做x的幅度谱，第三个做y的时域波形，第四个做y的幅度谱。 三、实验平台： MATLAB集成系统 、、 参考程序： 窗口设计法： ws1=0.2*pi;wp1=0.35*pi;wp2=0.65*pi;ws2=0.8*pi;As=60; tr_width=min((wp1-ws1),(ws2-wp2));M=ceil(11*pi/tr_width)+1; n=[0:1:M-1];wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2; alpha=(M-1)/2;m1=n-alpha;fc1=wc1/pi;hd1=wc1*sinc(fc1*m1); alpha=(M-1)/2;m2=n-alpha;fc2=wc2/pi;hd2=wc2*sinc(fc2*m2); hd=hd1-hd2; w_bla=(blackman(M));h=hd.*w_bla;b=h;a=[1]; [H,w]=freqz(b,a,1000,whole); mag=abs(H);db=20*log10((mag+eps)/max(mag)); pha=angle(H);grd=grpdelay(b,a,w); delta_w=2*pi/1000; Rp=-min(db(wp1/delta_w+1:1:wp2/delta_w)) As=-round(max(db(ws2/delta_w+1:1:501))) subplot(2,2,1);stem(n,-hd);xlabel(n);ylabel(hd(n)); axis([0,M-1,-1.2,2]); subplot(2,2,2);stem(n,w_bla);xlabel(n);ylabel(w(n)); axis([0,M-1,0,1.1]); subplot(2,2,3);stem(n,-h);xlabel(n);ylabel(h(n)); axis([0,M-1, -1.2,2]); subplot(2,2,4);plot(w/pi,db);xlabel(w);ylabel(db); axis([0,1, -139,10]); 接上面的程序： ws1=0.2*pi;wp1=0.35*pi;wp2=0.65*pi;ws2=0.8*pi;As=60; tr_width=min((wp1-ws1),(ws2-wp2));M=ceil(11*pi/tr_width)+1; n=[0:1:M-1];wc1=(ws1+wp1)/2;wc2=(ws2+wp2)/2; alpha=(M-1)/2;m1=n-alpha;fc1=wc1/pi;hd1=wc1*sinc(fc1*m1); alpha=(M-1)/2;m2=n-alpha;fc2=wc2/pi;hd2=wc2*sinc(fc2*m2); hd=hd1-hd2; w_bla=(blackman(M));h=hd.*w_bla; w1=0.5*pi;w2=0.1*pi; x=cos(w1*n)+cos(w2*n); nyb=n(1)+n(1);nye=n(length(x))+n(length(h)); ny=[nyb:nye];y=conv(x,h); s=length(n); w=0:pi/s:2*pi-pi/s; k=-s:s-1; X=x*(exp(-j*pi/s)).^(n*k); N=length(ny); g=0:pi/N:2*pi-pi/N; p=-N:N-1; Y=y*(exp(-j*pi/N)).^(ny*p); magX=abs(X);magY=abs(Y); subplot(2,2,1);stem(n,x);xlabel(n);ylabel(x(n));axis([0,M-1,-2.5,2.5]); subplot(2,2,2);plot(w/(2*pi),magX);xlabel(w);y