资源描述：
《用窗口法设计fir数字滤波器》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、一.实验目的了解一个实际滤波器设计过程，加深掌握用窗口发设计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．调试运行程序
2、,要求在屏幕上显示出单位脉冲响应h(n)的数值.画出其幅度响应
3、H
4、及20lg
5、H
6、的曲线.3．画出窗口函数w(n)及其频谱W的幅度
7、W
8、和20lg
9、W
10、的曲线.四.试验程序及图谱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);ifrem(N,2)~=0;hdu(r+1)=wc/pi;endwn=ha
11、mming(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);ifrem(N,2)~=0;hdu(r+1)=wc/pi;endwn=hamming(N);h=hdu.*wn';H
12、=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)。
13、若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);
14、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.6
15、006];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,'.');FIRIIR
