数字信号处理实验三滤波器设计.docx

实验三滤波器设计一．实验内容1．认真复习滤波器幅度平方函数的特性，模拟低通滤波器的巴特沃思逼近、切比雪夫型逼近方法；复习从模拟低通到模拟高通、带通、带阻的频率变换法；从模拟滤波器到数字滤波器的脉冲响应不变法、双线性变换法的基本概念、基本理论和基本方法。2.掌握巴特沃思、切比雪夫模拟低通滤波器的设计方法；利用模拟域频率变换设计模拟高通、带通、带阻滤波器的方法。3.掌握利用脉冲响应不变法、双线性变换法设计数字滤波器的基本方法；能熟练设计巴特沃思、切比雪夫低通、带通、高通、带阻数字滤波器。4.熟悉利用MATLAB 直接进行各类数字滤波器的设计方法。二．实验内容a. 设计模拟低通滤波器，通带截止频率为10KHz，阻带截止频率为16KHz，通带最大衰减1dB，阻带最小衰减20dB。(1) 分别用巴特沃思、切比雪夫I、切比雪夫II 型、椭圆型滤波器分别进行设计，并绘制所设计滤波器的幅频和相频特性图。(2) 在通带截止频率不变的情况下，分别用n=3,4,5,6 阶贝塞尔滤波器设计所需的低通滤波器，并绘制其相应的幅频响应和相频响应图。（1）新建函数a1，代码如下:%滤波器主要技术指标OmegaP=2*pi*10000; OmegaS=2*pi*16000;Rp=1;As=20;%巴特沃斯低通滤波器%确定滤波器的阶次和边缘频率N1=ceil(log10((10.^(0.1*abs(Rp))-1)./(10.^(0.1*abs(As))-1))/(2*log10(OmegaP/OmegaS)))OmegaC1=OmegaP/((10^(.1*abs(Rp))-1)^(1/(2*N1)));[z0,p0,k0]=buttap(N1)b0=k0*real(poly(z0));a0=real(poly(p0))[H1,w1]=freqs(b0,a0);subplot(211)plot(w1,abs(H1));title(巴特沃斯滤波器幅频特性);subplot(212)plot(w1,angle(H1));title(巴特沃斯滤波器相频特性);%切比雪夫I型figure(2);[N2,OmegaC2]=cheb1ord(OmegaP,OmegaS,Rp,As,s)[z2,p2,k2]=cheb1ap(N2,Rp)b2=k2*real(poly(z2));a2=real(poly(p2))[H2,w2]=freqs(b2,a2);subplot(211)plot(w2,abs(H2));title(切比雪夫I型滤波器幅频特性);subplot(212)plot(w2,angle(H2));title(切比雪夫I型滤波器相频特性);%切比雪夫2型figure(3)[N3,OmegaC3]=cheb2ord(OmegaP,OmegaS,Rp,As,s)[z3,p3,k3]=cheb2ap(N3,As)b3=k3*real(poly(z3));a3=real(poly(p3))[H3,w3]=freqs(b3,a3);subplot(211)plot(w3,abs(H3));title(切比雪夫II型滤波器幅频特性);subplot(212)plot(w3,angle(H3));title(切比雪夫II型滤波器相频特性);%椭圆形滤波器figure(4)[N4,OmegaC4]=ellipord(OmegaP,OmegaS,Rp,As,s)[z4,p4,k4]=ellipap(N4,Rp,As)b4=k4*real(poly(z4));a4=real(poly(p4))[H4,w4]=freqs(b4,a4);subplot(211)plot(w4,abs(H4));title(椭圆形低通滤波器幅频特性);subplot(212)plot(w4,angle(H4));title(椭圆形低通滤波器相频特性);生成图像如下：(2)在matlab中新建函数a2，代码如下：wp=2*pi*10000; ws=2*pi*16000;Rp=1;As=20;n=3;[b,a]=besself(n,wp);[H,w]=freqs(b,a);figure(1);subplot(211)plot(w,abs(H));title(n=3阶贝塞尔低通滤波器幅频特性);subplot(212)plot(w,angle(H));title(n=3阶贝塞尔低通滤波器相频特性);n=3;[b,a]=besself(n,wp);[H,w]=freqs(b,a);figure(2);subplot(211)plot(w,abs(H));title(n=4阶贝塞尔低通滤波器幅频特性);subplot(212)plot(w,angle(H));title(n=4阶贝塞尔低通滤波器相频特性);n=3;