用MATLAB来仿真FIR滤波器 1、矩形窗程序： clear all; N=64; w ....docVIP

下载本文档

37
0
约2.59千字
约 12页
2018-02-16 发布于江西
举报
版权申诉

用MATLAB来仿真FIR滤波器 1、矩形窗程序： clear all; N=64; w ....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文档。上传文档

用MATLAB来仿真FIR滤波器 1、矩形窗程序： clear all; N=64; w ...

用MATLAB来仿真FIR滤波器 1、矩形窗程序： clear all; N=64; w=rectwin(N); Wvtool(w); 仿真结果：从仿真上来看：主瓣宽度为4π/N，第一副瓣比主瓣低13dB，阻带衰减21dB，能量主要集中在旁瓣上，而主瓣的能量较低。而且过渡带比较窄。 2、三角形窗程序： clear all; N=64; w=triang(N); wvtool(w); 仿真结果：从上图上来看：主瓣宽度为8π/N，第一副瓣比主瓣低26dB，阻带衰减25dB。 3、汉宁窗程序： clear all; N=64; w=hann(N); wvtool(w); 仿真结果为：从上图上来看：主瓣宽度为8π/N，第一副瓣比主瓣低31dB，阻带衰减44dB。汉宁窗的幅度函数由三部分相加，使能量更集中在主瓣中，但代价是主办的宽度加宽了8π/N 4、汉明窗程序： clear all; N=64; w=hamming(N); wvtool(w); 仿真结果为：从上图上来看：主瓣宽度为8π/N，第一副瓣比主瓣低40dB，阻带衰减53dB。这种改进的升余弦窗，能量更加集中在主瓣中，主办的能量约占99.96%，第一旁瓣的峰值比主瓣校40 dB，但主瓣宽度和汉宁窗相同，仍为8π/N 5、布莱克曼窗程序： clear all; N=64; w=blackman(N); wvtool(w); 仿真结果为：从上图上来看：主瓣宽度为12π/N，第一副瓣比主瓣低57dB，阻带衰减74dB。使旁瓣在进一步抵消，阻带衰减进一步增加，过渡带是矩形窗过渡带的3倍。 6、凯塞窗程序： clear all; N=64; w=kaiser(64,1.5); wvtool(w); 仿真结果为：能量更加集中在主瓣上，过渡带增大。为依赖于N的参数，选择N可以产生各种过渡带宽和接近最优的阻带衰减。对于相同的N，此窗可以提供不同的过渡带宽，这是其他窗做不到的。 7、第一、二种情况理想低通滤波器的单位冲激响应hd(n)计算的程序： function hd=ideal_lp(wc,M); alpha=(M-1)/2; n=0:1:M-1; m=n-alpha+eps; hd=sin(wc*m)./(pi*m); 利用汉明窗设计第二种情况数字低通滤波器： Wp=0.2pi ws=0.4pi Ap=0.24dB As=50dB clear all; wp=0.2*pi; ws=0.4*pi; tr_width=ws-wp M=ceil(6.6*pi/tr_width)+1 n=0:1:M-1; wc=(ws+wp)/2; hd=ideal_lp(wc,M); w_ham=(hamming(M)); h=hd.*w_ham; [db,mag,pha,w]=freqz_m2(h,[1]); delta_w=2*pi/1000; Ap=-(min(db(1:1:wp/delta_w+1))) As=-round(max(db(ws/delta_w+1:1:501))) subplot(221) stem(n,hd) title(理想单位脉冲响应hd(n)); grid subplot(222) stem(n,w_ham) title(汉明窗w(n)); grid subplot(223) stem(n,h) title(实际单位脉冲响应hd(n));grid subplot(224) plot(w/pi,db) title(幅度响应(dB)); grid axis([0,1,-100,10]) function[db,mag,pha,w]=freqz_m2(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); 仿真结果为：用布莱克曼窗设计上面的要求：程序为： clear all; wp=0.2*pi; ws=0.4*pi; tr_width=ws-wp M=ceil(11*pi/tr_width)+1 n=0:1:M-1; wc=(ws+wp)/2; hd=ideal_lp(wc,M); w_bman=(blackman(M)); h=hd.*w_bman [db,mag,pha,w]=freqz_m2(h,[1]); delta_w=2*pi/1000; Ap=-(min(db(1:1:wp/delta_w+1))) As=-round(max(db(ws/delta_w+1:1:501))) subplot(221) stem(n,hd) title(理想单位