数字FIR滤波器试卷.ppt

下载文档 降价啦

6
0
约7.77千字
约 79页
2017-07-05 发布于湖北
举报
版权申诉
保障服务

数字FIR滤波器试卷.ppt

1、本文档共79页，可阅读全部内容。
2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

clf whm=0.54-0.46*cos(2*pi*(0:N-1)/(N-1)); axis([0,N-1,-0.3,1.2]) hold on bar(0:N-1,whm,0) plot(0:N-1,whm,*) plot([0,N-1],[0,0]) xlabel(n) title(Hamming window) grid hold off pause clf hn_hm=hn.*whm; axis([0,N-1,-0.3,1.2]) hold on bar(0:N-1,hn_hm,0) plot(0:N-1,hn_hm,*) plot([0,N-1],[0,0]) xlabel(n) title(Hamming windowed H(n)) grid hold off pause clf i=1; Hf=0; for f=-fs/2:(fs/200):fs/2 Hf(i)=hn_hm*exp(-j*2*pi*f/fs*(0:N-1)); i=i+1; end Hf_abs=abs(Hf); f=-fs/2:(fs/200):fs/2; axis([-fs/2,fs/2,-0.3,1.2]) hold on plot(f,Hf_abs) plot([-fs/2,fs/2],[0,0],k) xlabel(f in kHz) title(H(f),Hamming windowed) f = [-fs/2,-fc,-fc, 0, fc,fc,fs/2]; Df=[0, 0, 1, 0.5,1, 0, 0]; plot(f,Df,k--) plot(fk,Dfk,k*) grid hold off pause clf L=201;P=(L-1)/2; delta_f=fs/L; fk=-fs/2*(1-1/L): delta_f : fs/2*(1-1/L); Dfk=((1+abs(fk)/fc)/2).*(sign(fk+fc)-sign(fk-fc))/2; wi=2*pi*(-P:P)/L; for k=-P:P d(k+P+1)=Dfk*exp(-j*wi*k)/L; end dk=real(d); k=-P:P; axis([-P,P,-0.3,0.8]) hold on bar(k,dk,0) plot([-P,P],[0,0]) xlabel(k) title(d(k),LN) grid hold off pause clf hn=dk(P-M+1:P+M+1); k=0:N-1; axis([0,N-1,-0.3,0.8]) hold on bar(k,hn,0) plot(k,hn,*) plot([0,N-1],[0,0]) xlabel(n) title(rectangular windowed h(n) from d(k),k=-P:P) grid hold off pause clf i=1; Hf=0; for f=-fs/2:(fs/200):fs/2 Hf(i)=hn*exp(-j*2*pi*f/fs*(0:N-1)); i=i+1; end Hf_abs=abs(Hf); f=-fs/2:(fs/200):fs/2; axis([-fs/2,fs/2,-0.3,1.2]) hold on plot(f,Hf_abs) plot([-fs/2,fs/2],[0,0],k) xlabel(f in kHz) title(H(f),rectangular windowed,LN) f = [-fs/2,-fc,-fc, 0, fc,fc,fs/2]; Df=[0, 0, 1, 0.5,1, 0, 0]; plot(f,Df,k--) grid hold off pause clf i=1; hn_hm=hn.*whm; Hf=0; for f=-fs/2:(fs/200):fs/2 Hf(i)=hn_hm*exp(-j*2*pi*f/fs*(0:N-1)); i=i+1; end Hf_abs=abs(Hf); f=-fs/2:(fs/200):fs/2; axis([-fs/2,fs/2,-0.3,1.2]) hold on plot(f,Hf_abs) plot([-fs/2,fs/2],[0,0],k) xlabel(f in kHz) title(H(f),Hamming windowed,LN) f = [-fs/2,-fc,-fc, 0, fc,fc,fs/2]; Df=[0, 0, 1, 0.5,1, 0, 0]; plot(f,Df,k--) grid h