实验三：用窗函数法设计fir数字滤波器1

《实验三：用窗函数法设计fir数字滤波器1》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

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

2、d(A2)69(/7)(2-2)//(«)就作为实际设计的FTR数字滤波器的单位脉冲响应序列，其频率响应函数为:HAej，=Yh(n)e，(2-3)z/=0式中，N为所选窗函数必(Z2)的长度。由第七章可知，用窗函数法设计的滤波器性能取决于窗函数仍(/7)的类型及窗口长度N的取值。设计过程屮，要根据对阻带最小衰减和过渡带宽度的要求选择合适的窗函数类型和窗口长度N。各种类型的窗函数可达到的阻带最小哀减和过渡带宽度见第七章。这样选定窗函数类型和长度N后，求出单位脉冲响应=•a)(n),并按式(2-3)求出是否满足要求，要进行验算。一般在尾部加零使长度满足于2的整数

3、次幂，以便用FET计算如果要观察细节，补零点数增多即可。如果不满足要求，则要重新选择窗函数类型和长度N,再次验算，直至满足要求。如果要求线性相位特性，则/2(…还必须满足h(n)=±h(N-1-/7)(2-4)根据上式中的正负号和长度N的奇偶性乂将线性相位滤波器分成四类。要根据设计的滤波特性正确选择其中一类。例如，要设计线性低通特征，可选择/z(n)=/z(7V_l-W—类，而不能选h(n)=-h(N-1-n)一类。四、实验内容1.复习用窗函数法设计FIK数字滤波器一节内容，阅读本实验原理掌握设计步骤。2.编写程序①编写能产生矩型窗、哈明窗、汉宁窗、莱克曼窗的

4、窗函数子程序。②编写主程序。主程序框图如阁2—1所示，仅供参考。其屮幅度特性要求用分贝dB表示。图2—1主程序框图如设：H(k)=DFT[h{n)](2-4)H(k)=HR(k)+jHf(K)(2-5)(2-6)Ik为使曲线包络更画图时，^lgH(k)打印幅度特性。第k点对应的频率N接近的幅度特性曲线，变换区间要选大些。例如窗口长度N=33时，可通过在/7(n)末尾补零的方法，使长度变为64,再进行4点DFT,则可得到更精确的幅度衰减特性曲线。3.上机实验内容川叫种窗函数设计线性相位低通FTR数字滤波器，截止频率6rad,N=256。绘制相应的幅频特性曲线

5、，观察3dB带宽和2CMB带宽以及附带最小衰减，比较四种窗函数对滤波器特性的影响。设计低通F1K数字滤波器时，一般以理想低通滤波特性为逼近函数，即(2-8)，cc-j^ejcaideo=s[nCt)c(n~a)7r(n-a)(2-9)五、思考题1.如果给定通带截止频率和阻带截止频率以及阻带最小哀减，如何用窗函数法设•线性相位低通滤波器？写出设计步骤。2.定性说明本实验中3dB截止频率的理论值在什么位置？是否等于理想低通的截止频率CD?3.如果要求用窗函数法设H•带通滤波器，且给定上、下边带截止频率为％和必2,试求理想带通的单位脉冲响应叫Z7)。六、实验报告要求

6、1.简述实验目的及原理。2.按照实验步骤及要求，比较各种情况下的滤波性能，说明窗口长度N和窗函数类型对滤波特性的影响。3.总结用窗函数法设讣FIR滤波器的主要特点。4.简要回答思考题。参考程序1:N=i叩ut(’窗宽度N=’)；k=input(’窗型：1.矩形窗，2.hanning(升余弦窗)，3.hamming(改进的升余弦窗),4.Blackman请选择：’)；subplot(2,2,1);w二pi/5;a:(N-l)/2;n:0:(N-l);m=n—a+eps;h=sin(wm)./(pi本m);ifk==lB:bartlett(N);elseifk==

7、2B=hanning(N);elseifk==3B=hamming(N);elseifk==4B=blackman(N);endendendendhd=h.(B’);stem(n,hd,’.’)；xlabel(’n’)；ylabcl(’h(n)’)；titleC在矩形窗下的N=33时h(n)函数’)；subplot(2,2,2);[H,m]=freqz(hd,[1],1024);mag=abs(H);db=20logl0((mag+cps)/max(mag));plot(m/pi,db);xlabel(’w/n’)；ylabel(’201og[H(ejw)]’

8、)；titleCh(n)的幅频特性’)