实验11用matlab设计fir数字滤波器

《实验11用matlab设计fir数字滤波器》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、41彡树&实验报告课程名称：数字信号处理实验成绩：指导教师：院系：信息工程学院班级：电信二班学号：姓名：日期：、实验目的:1、加深对窗蛾数法设计FIR数字滤波器的基本原理的理解。2、学习用MATLAB语言的窗函数法编写设计FIR数字滤波器的程序。3、了解MATLAB语言有关窗函数法设计HR数字滤波器的常用函数用法。二、实验内容及步骤2、选择合适的愈函数设计FIR数字低通滤波器，要求：wp=0.271,Rp=0.05dB；ws=0.3兀，As=40dBo描绘该滤波器的脉冲响应、窗函数及滤波器的幅频响放曲线和相频响应曲线。分析：根裾设计指标要求，并査表

2、11-1,选择汉宁窗。程序清单如下：functionhd=ideal_lp(wc，N)wp=0.2*pi;ws=0<3*pi;deltaw=ws-wp;tao=(N-l)/2;n=[0:(N-l)];m=n-tao+eps;hd=sin(wc*m)

3、=angle(H);grd=grpdelay(b,a,w);wp=0.2*pi;ws=0.3*pi;dcltaw=ws-wp;wc=(ws+wp)/2;N0=ceil(6.6*pi/deltaw)；z=zo+mcd(NO+l,2)windows=(hanning(N)y;hd=ideal_lp(wcz)；b=hd.*windows;【db,mag,phaOQrd>VJ=freqzlm(bJ)；n=0:N';dw=2*pi/l000;Rp=-(min(db(l:wp/dw+l)))％郎錄is球落&As=-round(max(db(ws/dw+l:50

4、1)))％^ie*、sKeMMsubplot(2,2J)；stem(n,b)；axis(【o,N,1•1*min(b)，1•1*max(b)】)；dtlef將蜀葵^sgt-)；xlabelfn-)；ylabelch(ny)；subplot(2,2,2)；stem(n,windows)；axis(【o,N,o,1•1】)；title(M®燁忝w)；xlabelf5;ylabelfwd(ny)；subplot(2,2,3)；plot(w/pi,db)；axis(【0,l,-80,10】)；titlefgi12&^-Msw)；set(gca,-XTick

5、Mode7manual7XTick-,o,wp/pi,ws/pi,lJ)；set(gca,-YTickMode7manual7YTick-,【-50,-20,-3,0】)；gridsubplot(2,2,4)；plot(w/pi,pha)；axis(【0,l,-4,41)；titlef益商雞樹sw)；set(gca,-XTickMode7manual7XTick-,o,wp/pi,ws/pi,lJ)；z=67Rp=0.0706As=44相位频率响应频率(单位：JC)频率(单位：兀)3、用凯塞窗设id~个FIR数字高通滤波器，要求：wp=0.371,

6、Rp=O.ldB；ws=0.2ti，As=50dB。描绘该滤波器的脉冲响应、窗函数及滤波器的幅频响应曲线和相频响应曲线。程序清单如下：functionhd=ideal—lp(wc，N)tao=(N-l)/2;n=[0:(N-l)];m=n-tao+eps;hd=sin('vc*m)./(pi*m);function[db，mag,pha,grd,w]=freqz_m(b,a);[H,w]=freqz(b,a,1000,’whole’)；H=(H(l:501))’；w=(w(l:501))';raag=abs(H);db=20*logl0((mag+

7、eps)/max(mag));pha=angle(H);grd=grpdelay(b,a,w);wp=0.3*pi;ws=0.2*pi;deltaw=wp-ws;N0=ceil(6.6*pi/deltaw);N=N0+mod(N0+l，2)windows:(kaiser(N,7.865))';wc=(ws+wp)/2;hd=ideal_lp(pi,N)-ideal_lp(wc,N);b=hd.*windows;[db,mag,pha,grd,w]=freqz_m(b,1);n=0:N-l;dw=2*pi/1000;Rp=-(min(db(wp/dw

8、+l:501)))%检验通带波动As=-round(max(db(1:ws/dw+1)))%检验最小阻带衰减subplot