汉宁窗设计带阻FIR滤波器

《汉宁窗设计带阻FIR滤波器》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、4E0频率特性ogp、罟竖iiaiaaiii•ii■{•"：■"iiI•ii、iii了iii•*•}iIi■10005101520253035404550频率/Hz零极图i!rq••»iii•••••P!11••••••I••111111111111111T1i•i•i•ii#ii(,1Y111111111111J••F7•i•iii•i•••••IIJiL』7••••'•>ii」x107-1.5-105Re[z]°0.51.5x108Figure1FileEditViewInsertToolsDesktopWindowHelp口曰Pd図氈貯⑥娱謠▼屈□匡］■口频率特性1000-

2、10005101520253035404550频率/Hz4单位脉冲响应IIIIIIIInunuo(u)q40团Reg0.4单位脉神响应0.8

3、:::f:0.6:；；■：•0102030405060708042JO-2-4x107零极图1rid••••・•&11•••rr1i••・••1••.80'：iI1■■・p•1•1・•Fii••4I」••・•i4d••i>I」•ILAReM2x108»wp仁03pi;ws仁O・0pi;ws2=O®pi;wp2=O/Tpi;»N=ceil(8wpi/wdelta)80»wdelta=ws1-wp1;»wp1=0.3*pi;ws1=0.4*pi;

4、ws2=0.6tpi;wp2=0.7*pi;»subplot(3,1,2)9impz(b,1)»wp1=0.3*pi;ws1=0.4*pi;ws2=0.6*pi;wp2=0.7*pi;»wdelta=ws1-wp1;»N=ceil(8*pi/wdelta)80»wc=[(ws1+wp1)/2,(ws2+wp2)/2];»n=0:N-1;»b=fir1(N,wc/pi,'stop1,hanning(N+1));»[H,f]=freqz(b,1,512,100);

5、»figure(1)»subplot(3,1,1),plot(f,20*log10(abs(H)))»xlabelC频率/

6、Hza);ylabel(a振幅/dB();gridon;»titleC频率特性)；»subplot(3,1,2),impz(b,1)»xlabel(,n');ylabel(h(n)');gridon;»titlef单位脉冲响应•);»subplot(3,1,3),zplane(b,1)»xlabel('Re[z]);ylabel('lm[z]');g「idon;»title©零极图J;»f1=10；f2=25；f3=40；

7、»dt=0.02；t=0:dt:3；»x=sin(2*pi*frt)+cos(2*pi*f2*t)+cos(2*pi*f3*t);»y=fftfilt(b,x)

8、;»figure⑵»subplot(2,1,1),plot(t,x),title('输入信号)»xl3bel('时间/s')»subplot(2,1,2),plot(t,y),titleC输出信号')»xlabelC时间/s')wp1=0.3*pi;wsl=0.4*pi;ws2=0.6*pi;wp2=0.7*pi;»N=ceil(8*pi/wdelta)80»wdelta=ws1-wp1;»wp1=0.3*pi;ws1=0.4*pi;ws2=0.6*pi;wp2=0.7*pi;»subplot(3,1,2),impz(b,1)»wp1=0.3*pi;ws1=0.4*pi;ws2=0

9、.6*pi;wp2=0.7*pi;»wdelta=ws1・wp1;»N=ceil(8*pi/wdelta)N=80»wc=[(ws1+wp1)/2,(ws2+wp2)/2];»n=0:N-l;»b=fir1(N,wc/pi/stophanning(N+1));»[H,f]=freqz(b丄512,100);»figure(1)»subplot(3,l,l),plot(f,20*log10(abs(H)))»xlabel('频率/Hz');ylabel('振幅/dB');gridon;»title。频率特性J;»subplot(3,l,2),impz(b,1)»xlabel('n'

10、);ylabel('h(n)');gridon;»titled单位脉冲响应»subplot(3J,3),zplane(b,1)»xlabel(Re[z『);ylabelCIm[z『);gridon;»titleC零极图J;»fl=10;f2=25;f3=40;»dt=0.02;t=0:dt:3;»x=sin(2*pi*fl*t)4-cos(2*pi*f2*t)+cos(2*pi*f3*t);»y=fftfilt(b,x);»figure(2)»subplot(2,l,