《中南大学数字信号处理实验四》由会员分享，可在线阅读，更多相关《中南大学数字信号处理实验四（9页珍藏版）》请在金锄头文库上搜索。
1、 实验报告实验名称 用窗函数法设计FIR数字滤波器 课程名称 数字信号处理 姓名 成绩 班级 学号 日期 2014年5月24号 地点 综合实验楼机房 备注：1.实验目的（1）掌握用窗函数法设计FIR数字滤波器的原理和方法。（2）熟悉线性相位FIR数字滤波器的特性。（3）了解各种窗函数对滤波特性的影响。2.实验环境应用MATLAB 6.5软件操作系统：windows XP3.实验原理如果所希望的滤波器的理想频率响应函数为（），则其对应的单位脉冲响应为 （41）窗函数设计法的基本原理时用有限长单位脉冲响应序列h(n)逼近。用窗函数将截断，并进行加权处理，得到： （42）就作为实际设计的FIR数字滤
2、波器的单位脉冲响应序列，其频率响应函数为 （43）式中，N为所选窗函数的长度。由书本第七章可知，用窗函数法设计的滤波器性能取决于窗函数的类型及窗口长度N的取值。设计过程中，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。 这样选定窗函数类型和长度N后，求出单位脉冲响应，并按式（43）求出。是否满足要求，要进行验算。一般在h(n)尾部加零使长度满足2的整数次幂，以便用FFT计算。如果要观察细节，补零点数增多即可。如果补满足要求，则要重新选择窗函数类型和长度N，再次验算，直至满足要求。如果要求线性相位特性，则h(n)还必须满足： 根据上式中的正负号和长度N的奇偶性又将线性相
3、位FIR滤波器分为四类。要根据所设计的滤波器特性正确选择其中一类。例如，要设计线性相位低通特性，可选择一类，而不能选择一类。4.实验内容（1）用升余弦窗设计一线性相位低通FIR数字滤波器，截止频率。窗口长度N15，33。要求在两种窗口长度情况下，分别求出h(n),打印出相应的幅频特性和相频特性曲线，观察3dB带宽和20dB带宽。总结窗口长度N对滤波特性的影响。 设计低通FIR数字滤波器时，一般以理想低通滤波特性为逼近函数，即 其中 （2）n33，用四种窗函数设计线性相位低通滤波器。绘制相应的幅频特性曲线，观察3dB和20dB带宽以及阻带最小衰减，比较四种窗函数对滤波器特性的影响。5.实验结果及
4、分析 （1）N=15时的矩形窗： N=33时的矩形窗： (2) N=15时的汉明窗： N=33时的汉明窗： (3) N=15时的汉宁窗： N=33时的汉宁窗： (4) N=15时的布莱克曼窗： N=33时的布莱克曼窗：（1）矩形窗的频率响应主瓣宽度为4/N，第一副瓣比主瓣低13dB。（2）哈明窗是改进的升余弦窗，能量更加集中在主瓣中，主瓣的能量约占99.96%，第一旁瓣的峰值比主瓣小40dB,但主瓣宽度和汉宁窗相同，为8/N。（3）汉宁窗的幅度函数由三部分相加，使能量更集中在主瓣中，但代价是主瓣宽度加宽到8/N。（4）布莱克曼窗的幅度函数由五部分组成，它们都是移位不同，且幅度也不同的函数，使旁
5、瓣再进一步抵消，阻带衰减进一步增加，过渡带是矩形窗过渡带的3倍。（5）调整窗口长度N可以有效低控制过渡带的宽度，减少带内波动以及加大阻带的衰减只能从窗函数的形状上找解决方法，如果能找到的窗函数形状，使其谱函数的主瓣包含更多的能量，相应旁瓣幅度就减小了，旁瓣的减小可使通带、阻带波动减小，从而加大阻带衰减，但这样总是以加宽过渡带为代价的。6.思考题（1）如果给定通带截止频率和阻带截止频率以及阻带最小衰减，如何用窗函数法设计线性相位低通滤波器？写出设计步骤。答： 首先确定模拟带通滤波器的技术指标。 确定归一化低通技术要求 设计归一化低通 将低通转化为带通（2）如果要求用窗函数法设计带通滤波器，且给定
6、上、下边带截止频率为和，试求理想带通的单位脉冲响应。答：可以通过和求出这样便可以利用公式求出。7.实验结论用窗函数法设计FIR滤波器的主要特点：设FT为希望逼近的频响特性函数,H() =FTh(n)为用窗函数法设计的实际滤波器的频响函数。通常取H()相应的理想频响特性作为。知识要点如下：（1）希望逼近的理想滤波器频响函数的表达式。因为数字滤波器一般要求设计成线性相位特性，所以必须满足上述线性相位FIR滤波器的频域特点。（2）熟悉各种常用窗函数的技术指标和加窗后对滤波特性的影响，这样才能根据设计指标正确选择窗函数类型及其长度N。（3）检验设计结果：窗函数法的设计结果单位脉冲响应。而检验一般在频域
7、进行。所以要计算检验3dB截止频率和阻带最小衰减，其计算量相当大，必须用计算机进行。（4）熟悉窗函数设计法的特点：设计过程简单.方便实用。但边界频率不易精确控制。所以设计完以后，必须检验结果。附件-实验代码主程序：s = -1;while(s0) clc;N=input(请输入窗函数长度 N=);s=1;endclose all;i=0;wc=pi/4;while(s) n=0:N-1; hd=ideal(wc,N); k=input(请选择窗口类型:n1(boxcar)n2(hamming)n3(hanning)n4(blackman)n请选择：,s); k=str2num(k); if(k
8、=1) B=boxcar(N); string=Boxcar,N=,num2str(N); else if(k=2) B=hamming(N); string=Hamming,N=,num2str(N); else if(k=3) B=hanning(N); string=Hanning,N=,num2str(N); else if(k=4) B=blackman(N); string=Blackman,N=,num2str(N); end end end end h=hd.*(B); H,f=freqz(h,1,1024,whole); db=20*log10(abs(H)/max(abs(
9、H); pha=angle(H); i=i+1; figure(i) subplot(2,2,1); n=0:N-1; stem(n,h,.); xlabel(n); ylabel(h(n); title(实际低通滤波器的h(n); text(0.3*N),0.27,string); hold off; subplot(2,2,2); plot(f/pi,db); axis(0 1 -100 0); xlabel(w/pi); ylabel(20log|H(jw)|); title(衰减特性(dB); grid; subplot(2,2,3); plot(f,pha); hold on; n=
10、0:7; x=zeros(8); plot(n,x); title(相频特性); xlabel(频率w（rad); ylabel(相位（rad）); axis(0 3.15 -4 4); subplot(2,2,4); plot(f,abs(H); title(幅频特性); xlabel(频率W(rad); ylabel(幅值|H(jw)|); axis(0 3.15 0 1.5); text(0.9,1.3,string); s=input(继续实验请选择:n1(继续实验)n0(退出)n请选择:); if(s=1) N=input(请输入窗函数的长度 N=); endend子程序：function hd=ideal(w,N);alpha=(N-1)/2;n=0:N-1;m=n-alpha+eps;hd=sin(w*m)./(pi*m);
