北理工数字信号处理实验3 IIR数字滤波器设计

1、实验3 IIR数字滤波器设计1、 实验目的1、 掌握利用脉冲响应不变法和双线性变换法设计IIR数字滤波器的原理及具体方法2、 加深理解数字滤波器和模拟滤波器之间的技术指标转化3、 掌握脉冲响应不变法和双线性变换法设计IIR数字滤波器的优缺点及适用范围2、 实验设备与环境计算机、MATLAB软件环境3、 实验理论基础4、 实验内容解：（1）巴特沃斯Wn=0.2*pi;Rp=1;As=15;N=6;b,a=butter(N,Wn,'s');bz,az=impinvar(b,a);H,w=freqz(bz,az);subplot(221);plot(w/pi,abs(H);grid

2、on;xlabel('Omega(pi)');ylabel('|H(jOmega)|');subplot(222);plot(w/pi,20*log10(abs(H)/max(abs(H);grid on;xlabel('Omega(pi)');ylabel('|H(jOmega)|,dB');subplot(223);plot(w/pi,angle(H)/pi);grid on;xlabel('Omega(pi)');ylabel('Phase of H(jOmega)(pi)');t=0:0.0

3、1:30;h=impulse(b,a,t);subplot(224);plot(t,h);grid on;xlabel('t(s)');ylabel('Impulse Response');（2）切比雪夫 1Wn=0.2*pi;Rp=1;N=6;b,a=cheby1(N,Rp,Wn,'s');bz,az=impinvar(b,a);H,w=freqz(bz,az);subplot(221);plot(w/pi,abs(H);grid on;xlabel('Omega(pi)');ylabel('|H(jOmega)|

4、9;);subplot(222);plot(w/pi,20*log10(abs(H)/max(abs(H);grid on;xlabel('Omega(pi)');ylabel('|H(jOmega)|,dB');subplot(223);plot(w/pi,angle(H)/pi);grid on;xlabel('Omega(pi)');ylabel('Phase of H(jOmega)(pi)');t=0:0.01:30;h=impulse(b,a,t);subplot(224);plot(t,h);grid on;xlab

5、el('t(s)');ylabel('Impulse Response');（3）切比雪夫 2Wn=0.3*pi;As=15;N=6;b,a=cheby2(N,As,Wn,'s');bz,az=impinvar(b,a);H,w=freqz(bz,az);subplot(221);plot(w/pi,abs(H);grid on;xlabel('Omega(pi)');ylabel('|H(jOmega)|');subplot(222);plot(w/pi,20*log10(abs(H)/max(abs(H);gr

6、id on;xlabel('Omega(pi)');ylabel('|H(jOmega)|,dB');subplot(223);plot(w/pi,angle(H)/pi);grid on;xlabel('Omega(pi)');ylabel('Phase of H(jOmega)(pi)');t=0:0.01:30;h=impulse(b,a,t);subplot(224);plot(t,h);grid on;xlabel('t(s)');ylabel('Impulse Response');（4

7、） 椭圆Wn=0.2*pi;Rp=1;As=15;N=6;b,a=ellip(N,Rp,As,Wn,'s');bz,az=impinvar(b,a);H,w=freqz(bz,az);subplot(221);plot(w/pi,abs(H);grid on;xlabel('Omega(pi)');ylabel('|H(jOmega)|');subplot(222);plot(w/pi,20*log10(abs(H)/max(abs(H);grid on;xlabel('Omega(pi)');ylabel('|H(jOmega)|,dB');subplot(223);plot(w/pi,angle(H)/pi);grid on;xlabel('Omega(pi)');ylabel('Phase of H(jOmega)(pi)');t=0:0.01:30;h=impulse(b,a,t);subplot(224);plot(t,h);grid on;xlabel('t(s)');ylabel('Impulse Response');五、实验心得与