用窗函数设计FI数字滤波器.docVIP

下载本文档

2
0
约1.25千字
约 2页
2019-07-23 发布于浙江
举报
版权申诉

用窗函数设计FI数字滤波器.doc

1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

实验名：用窗函数设计FIR数字滤波器实验任务：选择合适的窗函数设计FIR数字低通滤波器，要求：通带ωp=0.2π,Rp=0.05dB;阻带ωs=0.3π，As=40dB。描绘实际滤波器的脉冲响应、窗函数及滤波器的幅频响应曲线和相频响应曲线。实验程序： 1.建立新文件,以ideal_lp命名，键入以下程序，保存文件 function hd=ideal_lp(wc,N); tao=(N-1)/2; n=[0:(N-1)]; m=n-tao+eps; hd=sin(wc*m)./(pi*m); 2.建立新文件，任意命名，键入以下程序，保存并运行程序 wp=0.2*pi;ws=0.3*pi; deltaw=ws-wp; N0=ceil(6.2*pi/deltaw);%汉宁窗 N=N0+mod(N0+1,2) windows=(hanning(N)); wc=(ws+wp)/2; hd=ideal_lp(wc,N); b=hd.*windows; [db,mag,pha,grd,w]=freqz_m(b,1); n=0:N-1;dw=2*pi/1000; Rp=-(min(db(1:wp/dw+1))) As=-round(max(db(ws/dw+1:501))) subplot(2,2,1),stem(n,b); axis([0,N,1.1*min(b),1.1*max(b)]);title(实际脉冲响应); xlabel(n);ylabel(h(n)); subplot(2,2,2),stem(n,windows); axis([0,N,0,1.1]);title(窗函数特性); xlabel(n);ylabel(wd(n)); subplot(2,2,3),plot(w/pi,db); axis([0,1,-150,10]);title(幅频响应); xlabel(频率(单位：\pi));ylabel(H(e^{j\omega})); set(gca,XTickMode,manual,XTick,[0,wp/pi,ws/pi,1]); set(gca,YTickMode,manual,YTick,[-100,-40,-3,0]);grid subplot(2,2,4),plot(w/pi,pha); axis([0,1,1.1*min(pha),1.1*max(pha)]);title(相频特性 ); xlabel(频率(单位：\pi));ylabel(\phi(\omega)); set(gca,XTickMode,manual,XTick,[0,wp/pi,ws/pi,1]); set(gca,YTickMode,manual,YTick,[-pi,0,pi]);grid 实验运行结果： N =63 Rp =0.0888 As =44 实验图象：