北邮MATLAB实验三-数字滤波器设计.docx

《北邮MATLAB实验三-数字滤波器设计.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、实验3基于MATLAB的数字滤波器设计xxxxxxxxxxx班xx号大学霸实验内容：采用FIR数字滤波器类型，利用MATLAB编程设计一个数字带通滤波器，指标要求如下：通带边缘频率：wwp1=0.45π，wp2=0.65π通带峰值起伏：aαp≤1[dB]阻带边缘频率：wws1=0.3π，ws2=0.75π最小阻带衰减：αs≥40[dB]实验要求：给出IIR数字滤波器参数和FIR数字滤波器的冲激响应，绘出它们的幅度和相位频响曲线，讨论它们各自的实现形式和特点。IIR数字滤波器：>>wp1=0.45/2;wp2=0.65/2;ws1=0.3/2;ws2=0.75/2;Rp=1;Rs=40;

2、>>Wp=[wp1wp2];Ws=[ws1ws2];>>[N,Wn]=cheb1ord(Wp,Ws,Rp,Rs);>>[b,a]=cheby1(N,Rp,Wn);>>[H,W]=freqz(b,a);>>subplot(2,1,1);>>plot(W./pi.*2,20*log10(abs(H)));title('幅频响应');grid;axis([0.20.8-505]);>>subplot(2,1,2);>>plot(W./pi.*2,angle(H));title('相频响应');grid;axis([0.20.8-44]);FIR数字滤波器：Freqz函数定义：function

3、[db,mag,pha,grd,w]=freqz_m(b,a);[H,w]=freqz(b,a,1000,'whole');H=(H(1:1:501))';w=(w(1:1:501))';mag=abs(H);db=20*log10((mag+eps)/max(mag));pha=angle(H);grd=grpdelay(b,a,w);>>%Bandpassfilterdesign-Blackmanwindow>>ws1=0.3*pi;wp1=0.45*pi;>>wp2=0.65*pi;ws2=0.75*pi;>>Ap=1;As=40;>>tr_width=min((wp1-ws1)

4、,(ws2-wp2));>>N=ceil(11*pi/tr_width)+1>>n=[0:1:N];>>wc1=(ws1+wp1)/2;wc2=(wp2+ws2)/2;>>wn=[wc1/pi,wc2/pi];>>h=fir1(N,wn,blackman(N+1));>>[db,mag,pha,grd,w]=freqz_m(h,[1]);>>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)))>>freqz(h,1,256);N=111R

5、p=0.0027As=75