用窗口法设计FIR数字滤波器.docVIP

一.实验目的 了解一个实际滤波器设计过程，加深掌握用窗口发设计FIR数字滤波器的原理和窗函数对数字滤波器性能的影响。 二.实验原理 　设FIR滤波器的单位冲激响应h (n)为一个N点序列，0 ≤ n ≤N —1，则滤波器的系统函数为H(z)=∑h(n)*z^-n。就是说，它有（N—1）阶极点在z = 0处，有（N—1）个零点位于有限z平面的任何位置。 三.实验内容与要求 1．编写用改进的升余弦窗设计FIR线性相应低通数字滤波器的程序,已知wc=0.5,N=21.这里wc为低通滤波器的截止频率,N为滤波器的长度,滤波器的阶数为N-1. 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=hdu.*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(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]频率区间上，最小衰减大于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