实验四 IIR 滤波器设计 湖南大学.docVIP

实验四 IIR 滤波器设计 一、实验目的 认真复习滤波器幅度平方函数的特性，模拟低通滤波器的巴特沃思逼近、切比雪夫型逼近方法；复习从模拟低通到模拟高通、带通、带阻的频率变换法；从模拟滤波器到数字滤波器的脉冲响应不变法、双线性变换法的基本概念、基本理论和基本方法。 掌握巴特沃思、切比雪夫模拟低通滤波器的设计方法；利用模拟域频率变换设计模拟高通、带通、带阻滤波器的方法。 掌握利用脉冲响应不变法、双线性变换法设计数字滤波器的基本方法；能熟练设计巴特沃思、切比雪夫低通、带通、高通、带阻数字滤波器。 熟悉利用 MATLAB 直接进行各类数字滤波器的设计方法。 二、实验内容 a. 设计模拟低通滤波器，通带截止频率为10KHz，阻带截止频率为16KHz，通带最大衰减1dB，阻带最小衰减20dB。 (1) 分别用巴特沃思、切比雪夫I、切比雪夫II 型、椭圆型滤波器分别进行设计，并绘制所设计滤波器的幅频和相频特性图。 (2) 在通带截止频率不变的情况下，分别用n=3,4,5,6 阶贝塞尔滤波器设计所需的低通滤波器，并绘制其相应的幅频响应和相频响应图。 实验代码： 实验结果： b. 设计模拟高通滤波器，通带截止频率为2000Hz，阻带截止频率1500Hz，通带最大衰减为3dB，阻带最小衰减为15dB。 (1) 分别用巴特沃思、切比雪夫I 型滤波器首先设计模拟低通滤波器，再通过频率转换成高通滤波器，并分别绘制所设计滤波器的幅频和相频特性图。 (2) 直接用巴特沃思、切比雪夫I 型滤波器设计高通滤波器，并分别绘制所设计滤波器的幅频和相频特性图。 实验1）代码： 实验1）结果： 实验2）代码： 实验2）结果： c. 设计模拟带通滤波器，其通带带宽为B=2π×200rad/s，中心频率Ω0=2π×1000rad/s，通带内最大衰减为δ1=3dB，阻带Ωs1=2π×830rad/s，而Ωs2=2π×1200rad/s，阻带最小衰减δ2=15dB。 (1) 分别用巴特沃思、切比雪夫I 型滤波器首先设计模拟低通滤波器，再通过频率转换成带通滤波器，并分别绘制所设计滤波器的幅频和相频特性图。 (2) 直接用巴特沃思、切比雪夫I 型滤波器设计带通滤波器，并分别绘制所设计滤波器的幅频和相频特性图。 实验1）代码： clc; %通带带宽为B=2π×200rad/s，中心频率Ω0=2π×1000rad/s，通带内最大衰减为δ1=3dB? %阻带Ωs1=2π×830rad/s，而Ωs2=2π×1200rad/s，阻带最小衰减δ2=15dB。?%主要技术指标? wp=2*pi*1100; ws=2*pi*1200; wz=2*pi*1000; Rp=3;As=15;B=2*pi*200;%巴特沃斯带通滤波器? [N1,wc1]=buttord(wp,ws,Rp,As,s); [z1,p1,k1]=buttap(N1); b1=k1*real(poly(z1)); a1=real(poly(p1)) [numT1,denT1]=lp2bp(b1,a1,wz,B); [H1,w1]=freqs(numT1,denT1); subplot(211) plot(w1,abs(H1)); title(巴特沃斯带通滤波器幅频特性); subplot(212) plot(w1,angle(H1)); title(巴特沃斯带通滤波器相频特性);%切比雪夫带通滤波器? figure(2); [N2,wc2]=cheb1ord(wp,ws,Rp,As,s) [z2,p2,k2]=cheb1ap(N2,Rp) b2=k2*real(poly(z2)); a2=real(poly(p2)); [numT2,denT2]=lp2bp(b2,a2,wz,B); [H2,w2]=freqs(numT2,denT2); subplot(211) plot(w2,abs(H2)); title(切比雪夫I型带通滤波器幅频特性); subplot(212) plot(w2,angle(H2)); title(切比雪夫I型带通滤波器相频特性); 实验1）结果： 实验2）代码： clc; %通带带宽为B=2π×200rad/s，中心频率Ω0=2π×1000rad/s，通带内最大衰减为δ1=3dB? %阻带Ωs1=2π×830rad/s，而Ωs2=2π×1200rad/s，阻带最小衰减δ2=15dB。?%主要技术指标? wp=2*pi*1100; ws=2*pi*1200; wn=[2*pi*900,2*pi*1100]; Rp=3;As=15;B=2*pi*200;%巴特沃斯带通滤波器? [N1,wc1]=buttord(wp,ws,Rp,As,s); [b1,a1]=butter(N1,wn,s); [H1