《用窗口法设计FIR数字滤波器.doc》由会员分享，可在线阅读，更多相关《用窗口法设计FIR数字滤波器.doc（4页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、【精品文档】如有侵权，请联系网站删除，仅供学习与交流用窗口法设计FIR数字滤波器.精品文档.一.实验目的了解一个实际滤波器设计过程，加深掌握用窗口发设计FIR数字滤波器的原理和窗函数对数字滤波器性能的影响。二.实验原理设FIR滤波器的单位冲激响应h (n)为一个N点序列，0 n N 1，则滤波器的系统函数为H(z)=h(n)*z-n。就是说，它有（N1）阶极点在z = 0处，有（N1）个零点位于有限z平面的任何位置。三.实验内容与要求1编写用改进的升余弦窗设计FIR线性相应低通数字滤波器的程序,已知wc=0.5,N=21.这里wc为低通滤波器的截止频率,N为滤波器的长度,滤波器的阶数为N-1.
2、2调试运行程序,要求在屏幕上显示出单位脉冲响应h(n)的数值.画出其幅度响应 |H |及 20lg |H |的曲线.3画出窗口函数w(n)及其频谱W的幅度|W |和20lg |W |的曲线.四.试验程序及图谱1.N=21;Window=hamming(N);wc=0.5*pi;b=fir1(N-1,wc/pi,window);freqz(b,1,512);2.N=21; wc=0.5*pi; r=(N-1)/2; n=0:N-1; hdu=sin(wc*(n-r)/pi./(n-r); if rem(N,2)=0; hdu(r+1)=wc/pi; end wn=hamming(N); h=hd
3、u.*wn; H=fft(h,512); w=2*0:511/512; subplot(3,1,1); plot(w,20*log10(abs(H); subplot(3,1,2); plot(hdu); subplot(3,1,3); plot(w,abs(H);3.N=21; wc=0.5*pi; r=(N-1)/2; n=0:N-1; hdu=sin(wc*(n-r)/pi./(n-r); if rem(N,2)=0; hdu(r+1)=wc/pi; end wn=hamming(N); h=hdu.*wn; H=fft(h,512); Wn=fft(wn,100); magWn=abs
4、(Wn); phaWn=angle(Wn); w=2*0:511/512; subplot(3,1,1); plot(wn); subplot(3,1,2); stem(magWn); subplot(3,1,3); stem(20*log10(magWn);一.实验目的1）掌握数字滤波器的计算机仿真方法。2）通过观察对实际心电图信号的滤波作用，获得数字滤波器的感性知识二.实验原理取h(n)为低通数字滤波器, N=21.截止频率=0.2pi采用海明窗设计出h(n)。若H(z)是IIR滤波器的系统函数，在通带内频率低于0.2pi是，最大衰减小于1dB；在阻带0.3pi,pi频率区间上，最小衰减大
5、于15dB。三.试验程序及图谱1.x=4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0;N=21;Window=hamming(N);wc=0.2*pi;b=fir1(N-1,wc/pi,window);y=filter(b,1,x);stem(y);2.x=4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0;n=0:55;A=0.2318;b1=A,A*2,A;a1=1,0.4164,-0.6006;h1=filter(b1,a1,x)b2=A,A*2,A;a2=1,0.3093,-0.1888;h2=filter(b2,a2,h1)b3=A,A*2,A;a3=1,0.2692,-0.0349;h3=filter(b3,a3,h2)stem(n,h3,.);FIR IIR