实验二 IIR滤波器的设计.docVIP

实验二 IIR滤波器的设计牟雪曼 一、实验仿真 1. 设计一个巴特沃色低通滤波器，其技术指标为：3dB带宽为，阻带截至频率为,阻带衰减大于30dB，采样间隔。在同一图中画出原型、冲激不变法和双线性变换法所设计的数字滤波器的幅频特性图。 Matlab程序 clear; %原型滤波器参数 [n,wn]=buttord(0.2,0.5,3,30); [b,a]=butter(n,wn); %T=10us时，模拟滤波器参数 fs=10^5; %Ts=10us [n,wn]=buttord(0.2*fs*2*pi,0.5*fs*2*pi,3,30,s); [b1,a1]=butter(n,wn,s); %冲击响应不变法和双线性不变 [bz,az]=impinvar(b1,a1,fs); [numd,dend]=bilinear(b1,a1,fs); [h0,w0]=freqz(b,a,fs); [h1,w1]=freqz(bz,az,fs); [h2,w2]=freqz(numd,dend,fs); plot(w0,20*log10(abs(h0)),b);hold on; plot(w1,20*log10(abs(h1)),r);hold on; plot(w2,20*log10(abs(h2)),g); legend(Protype,Impulse invariance,Bilinear); xlabel(Freq w（rad/s）); ylabel(|H(jw)|(dB)); title(Frequency Responses of Protype/Impulse invariance/Bilinear); hold off; grid on 2. 用双线性变换法分别设计如下低通、高通、带通和带阻四种数字滤波器。模拟滤波器的原型分别为切比雪夫滤波器和椭圆滤波器，采样率Fs=20kHz。 低通：通带０～４kHz，阻带5kHz，通带衰减Rp=0.5dB，阻带衰减As=10dB。 Matlab程序： function[]=lp() %lowpass filter fs=2*10^4; %fs=20k %chebyshev [n1,wn1]=cheb1ord(4000*2*pi,5000*2*pi,0.5,10,s); [b1,a1]=cheby1(n1,0.5,wn1,s); [bz1,az1]=bilinear(b1,a1,fs); %ellipse [n2,wn2]=ellipord(4000*2*pi,5000*2*pi,0.5,10,s); [b2,a2]=ellip(n2,0.5,10,wn2,s); [bz2,az2]=bilinear(b2,a2,fs); [h1,f1]=freqz(bz1,az1,512,fs); [h2,f2]=freqz(bz2,az2,512,fs); plot(f1,20*log10(abs(h1)),r);hold on; plot(f2,20*log10(abs(h2)),b); legend(chebyshev,ellipse); xlabel(Freq f(Hz)); ylabel(|H(jw)| (dB)); title(LowPass Filter); grid on; 高通：阻带０～４kHz，通带5kHz，通带衰减Rp=0.5dB，阻带衰减 As=10dB。 Matlab程序： %function[]=hp() clear; %highpass filter fs=2*10^4; %fs=20k %chebyshev [n1,wn1]=cheb1ord(5000*2*pi,4000*2*pi,0.5,10,s); [b1,a1]=cheby1(n1,0.5,wn1,high,s); [bz1,az1]=bilinear(b1,a1,fs); %ellipse [n2,wn2]=ellipord(5000*2*pi,4000*2*pi,0.5,10,s); [b2,a2]=ellip(n2,0.5,10,wn2,high,s); [bz2,az2]=bilinear(b2,a2,fs); [h1,f1]=freqz(bz1,az1,512,fs); [h2,f2]=freqz(bz2,az2,512,fs); plot(f1,20*log10(abs(h1)),b);hold on; plot(f2,20*log10(abs(h2)),r); legend(chebyshev protype,ellipse protype); xlabel(Freq f(Hz)); ylabel(|H(jw)| (dB)); title(HighPa