用窗函数法设计FIR数字滤波器

1、用窗函数法设计FIR数字滤波器一、实验目的1 .掌握用窗函数法设计FIR数字滤波器的原理和方法。2 .熟悉线性相位FIR数字滤波器特征。3 .了解各种窗函数对滤波特性的影响。二、实验仪器微型计算机matlab软件三、实验原理和方法如果所希望的滤波器的理想频率响应函数为Hd(ej)，则其对应的单位脉冲响应为1ii，一.、hd(n)=-Hd(e)ed(2-1)窗函数设计法的基本原理是用有限长单位脉冲响应序列h(n)逼近hd(n)。由于hd(n)往往是无限长序列，且是非因果的，所以用窗函数(n)将hd(n)截断，并进行加权处理，得到：h(n)=hd(n)(n)(2-2)h(n)就作为实际设计的FIR

2、数字滤波器的单位脉冲响应序列，其频率响应函数Hd(ej)为：N1Hd(ej)=h(n)ej(2-3)n0式中，N为所选窗函数(n)的长度。由第七章可知，用窗函数法设计的滤波器性能取决于窗函数(n)的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度No各种类型的窗函数可达到的阻带最小衰减和过渡带宽度见第七章。这样选定窗函数类型和长度N后，求出单位脉冲响应h(n)=hd(n),(n),并按式(2-3)求出H(ej)。H(ej)是否满足要求，要进行验算。一般在h(n)尾部加零使长度满足于2的整数次哥，以便用FFT计算H(ej)。如果要观察细节，

3、补零点数增多即可。如果H(ej)不满足要求，则要重新选才i窗函数类型和长度N,再次验算，直至满足要求。如果要求线性相位特性，则h(n)还必须满足h(n)h(N1n)(2-4)根据上式中的正负号和长度N的奇偶性又将线性相位FIR滤波器分成四类。要根据设计的滤波特性正确选择其中一类。例如，要设计线性低通特征，可选择h(n)h(N1n)一类，而不能选h(n)h(N1n)一类。四、实验内容1 .复习用窗函数法设计FIR数字滤波器一节内容，阅读本实验原理掌握设计步骤。2 .编写程序 编写能产生矩型窗、哈明窗、汉宁窗、莱克曼窗的窗函数子程序。 编写主程序。主程序框图如图H(k)DFTh(n)(2-4)H(

4、k) HR(k) jHi(K)(2-5)H(k)|v'HR(k)H2(k)(2-6)一，一、2,2回图时，201gH(k)打印幅度特性。第k点对应的频率kNk。为使曲线包络更接近H(ej)的幅度特性曲线，DFT变换区间要选大些。例如窗口长度N=33时，可通过在h(n)末尾补零的方法，使长度变为64,再进行4点DFT则可得到更精确的幅度衰减特性曲线。3 .上机实验内容用四种窗函数设计线性相位低通FIR数字滤波器，截止频率c制相应的幅频特性曲线，观察3dB带宽和20dB带宽以及阻带最小衰减,"4 rad , N=256。绘比较四种窗函数对设计低通FIR数字滤波器时，一般以理想低通

5、滤波特性为逼近函数Hd(ej )，即hd(n)Hd(ej )Hd(ej其中0, c)ej ndsin(2-8)-c(n)(2-9)(n )滤波器特性的影响。五、思考题如何用窗函数法设计线性1 .如果给定通带截止频率和阻带截止频率以及阻带最小衰减,相位低通滤波器？写出设计步骤。2 .定性说明本实验中3dB截止频率的理论值在什么位置？是否等于理想低通的截止频下边带截止频率为 1和2，试求说明窗口长度 N和窗函数类型对3 .如果要求用窗函数法设计带通滤波器，且给定上、理想带通的单位脉冲响应h(n)。六、实验报告要求1 .简述实验目的及原理。2 .按照实验步骤及要求，比较各种情况下的滤波性能,滤波特性

6、的影响。3 .总结用窗函数法设计FIR滤波器的主要特点。4 .简要回答思考题。参考程序1:N=input('窗宽度N=');k=input('窗型：1.矩形窗，2.hanning(升余弦窗),3.hamming(改进的升余弦窗)，4.Blackman请选择：');subplot(2,2,1);w=pi/5;a=(N-1)/2;n=0:(N-1);m=n-a+eps;h=sin(w*m)./(pi*m);ifk=1B=bartlett(N);elseifk=2B=hanning(N);elseifk=3B=hamming(N);elseifk=4B=blackma

7、n(N);endendendendhd=h.*(B');stem(n,hd,'.');xlabel('n');ylabel('h(n)');title('在矩形窗下的N=33时h(n)函数')；subplot(2,2,2);H,m=freqz(hd,1,1024);mag=abs(H);db=20*log10(mag+eps)/max(mag);plot(m/pi,db);xlabel('w/n');ylabel('20logH(ejw)');title('h(n)的幅频特性'

8、;);grid;pha=angle(H);subplot(2,2,3);plot(m,pha);xlabel('n');ylabel('6');title('h(n)的相频特性');subplot(2,2,4);plot(m,mag);xlabel('w');ylabel('H(ejw)');title('h(n)的幅频特性');参考程序2：b=1;closeall;i=0;while(b);temp=menu('选择窗函数长度N','N=10','N=15&

9、#39;,'N=20','N=25','N=30','N=33','N=35','N=40','N=45','N=50','N=55','N=60','N=64');menu1=10,15,20,25,30,33,35,40,45,50,55,60,64;N=menu1(temp);temp=menu('选择逼近理想低通滤波器截止频率Wc','Wc=pi/4','Wc=pi/2&#

10、39;,'Wc=3*pi/4','Wc=pi','Wc=0.5','Wc=1.0','Wc=1.5','Wc=2.0','Wc=2.5','Wc=3.0');menu2=pi/4,pi/2,3*pi/4,pi,0.5,1,1.5,2,2.5,3;w=menu2(temp);n=0:(N-1);hd=ideal(w,N);%得到理想低通滤波器k=menu('请选择窗口类型:','boxcar','hamming','

11、;hanning','blackman');ifk=1B=boxcar(N);string='Boxcar','N=',num2str(N);elseifk=2B=hamming(N);string='Hamming','N=',num2str(N);elseifk=3B=hanning(N);string='Hanning','N=',num2str(N);elseifk=4B=blackman(N);string='Blackman','N=

12、9;,num2str(N);endendendendh=hd.*(B)'%得到FIR数字滤波器H,m=freqz(h,1,1024,'whole');%求其频率响应mag=abs(H);%得到幅值db=20*log10(mag+eps)/max(mag);pha=angle(H);%得到相位i=i+1;figure(i)subplot(2,2,1);n=0:N-1;stem(n,h,'.');axis(0,N-1,-0.1,0.3);holdon;n=0:N-1;x=zeros(N);plot(n,x,'-');xlabel('n

13、');ylabel('h(n)');title('实际低通滤波器的h(n)');text(0.3*N),0.27,string);holdoff;subplot(2,2,2);plot(m/pi,db);axis(0,1,-100,0);xlabel('w/pi');ylabel('dB');title('衰减特性(dB)');grid;subplot(2,2,3);plot(m,pha);holdon;n=0:7;x=zeros(8);plot(n,x,'-');title('相频特性');xlabel('频率(rad)');ylabel('相位(rad)');axis(0,3.15,-4,4);subplot(2,2,4);plot(m,mag);title('频率特性');xlabel('频率W(rad)');ylabel('