基于频率采样法的FIR低通滤波器的设计.doc

实验总成绩： 装 订 线报告份数： 西安邮电大学 通信与信息工程学院 数字信号处理报告专业班级： 学生姓名： 学号(班内序号): 2013 年 06 月 20 日 基于频率采样法的FIR低通滤波器的设计一、实验目的（1）熟悉线性相位FIR滤波器的幅频特性。（2） 掌握设计FIR数字滤波器的原理和方法。（3） 掌握用频率采样法设计FIR滤波器原理和方法。 (4) 学会调用MATLAB函数设计与实现FIR滤波器。二、实验程序%T1=0.5%设计条件:wp=0.2pi;ws=0.3pi;Rp=0.25dB;Rp=50dB;M=40;alpha=(M-1)/2;l=0:M-1;w1=(2*pi/M)*l;Hrs=ones(1,5),0.5,zeros(1,29),0.5,ones(1,4);Hdr=1,1,0,0;wdl=0,0.25,0.25,1; k1=0:floor(M-1)/2);k2=floor(M-1)/2)+1:M-1;angH=-alpha*(2*pi)/M*k1,alpha*(2*pi)/M*(M-k2);H=Hrs.*exp(j*angH);h=real(ifft(H,M);db,mag,pha,grd,w=freqz_m(h,1);Hr,ww,a,L=Hr_type2(h);subplot(1,1,1);subplot(2,2,1);plot(w1(1:21)/pi,Hrs(1:21),o,wdl,Hdr);axis(0,1,-0.1,1.1);title(Frequency Samples:M=40,T1=0.5);set(gca,XTickMode,manual,XTick,0,0.2,0.3,1);set(gca,YTickMode,manual,YTick,0,0.5,1);gridxlabel(frequency in pi units);ylabel(Hr(k)subplot(2,2,2);stem(l,h);axis(-1,M,-0.1,0.3);title(Impulse Response);xlabel(n);ylabel(h(n)xa=0.*l;hold onplot(l,xa,k);hold offsubplot(2,2,3);plot(ww/pi,Hr,w1(1:21)/pi,Hrs(1:21),o);axis(0,1,-0.2,1.2);title(Amplitude Responde)xlabel(frequency in pi units);ylabel(Hr(w)set(gca,XTickMode,manual,XTick,0,0.2,0.3,1);set(gca,YTickMode,manual,YTick,0,0.5,1);gridsubplot(2,2,4);plot(w/pi,db);axis(0,1,-100,10);title(Magnitude Response);xlabel(frequency in units);ylabel(Decibels);set(gca,XTickMode,manual,XTick,0,0.2,0.3,1);set(gca,YTickMode,manual,YTick,-30,0);gridset(gca,YTickLabelMode,manual,YTickLabels,30;0)3、 运行结果四、实验心得 本次实验让我受益无穷.这次实验设计的题目是基于频域抽样法的FIR数字低通滤波器设计，通过学习，为实际设计已经打好了理论基础，在此基础上，我通过自己动手设计完成了课程设计要求。 通过这次实验，我体会到了MATLAB软件的强大功能，了解到它在各种工程计算中的重要作用,我们更进一步理解数字滤波器设计原理，学会了数字滤波器设计的方法和一般步骤，能够独立设计一个数字滤波器，实现了把理论知