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

数字信号处理实验报告 实验项目名称：用窗函数法设计FIR数字滤波器  实验题目名称：用窗函数法设计FIR数字滤波器 实验目的： 掌握用窗函数法设计FIR数字滤波器的原理和方法。 熟悉线性相位FIR数字滤波器特性。 了解各种窗函数对滤波特性的影响。 实验环境：计算机，MATLAB软件 实验原理： 数字滤波器的设计是数字信号处理中的一个重要内容。数字滤波器设计包括FIR（有限单位脉冲响应）滤波器与IIR（无限单位脉冲响应）滤波器两种。 与IIR滤波器相比，FIR滤波器在保证幅度特性满足技术要求的同时，很容易做到严格的线性相位特性。设FIR滤波器单位脉冲响应h(n)长度为N，其系统函数H（z）为： H（z）是z－1的N－1次多项式，它在z平面上有N－1个零点，原点z=0是N－1阶重极点，因此H（z）是永远稳定的。稳定和线性相位特性是FIR滤波器突出的优点。 FIR滤波器的设计任务是选择有限长度的h(n)。使传输函数H( )满足技术要求。FIR滤波器的设计方法有多种，如窗函数法、频率采样法及其它各种优化设计方法，本实验介绍窗函数法的FIR滤波器设计。 窗函数法是使用矩形窗、三角窗、巴特利特窗、汉明窗、汉宁窗和布莱克曼窗等设计出标准响应的高通、低通、带通和带阻FIR滤波器。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);if k==1??? B=bartlett(N);else if k==2??????? B=hanning(N);??? else if k==3??????????? B=hamming(N);??????? else if k==4??????????????? B=blackman(N);??????????? end??????? end??? endendhd=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/п);ylabel(20log[H(ejw)]);title(h(n)的幅频特性);grid;pha=angle(H);subplot(2,2,3);plot(m,pha);xlabel(n);ylabel(φ);title(h(n)的相频特性);subplot(2,2,4);plot(m,mag);xlabel(w);ylabel(H(ejw));title(h(n)的幅频特性);FIR滤波器可以得到严格的线性相位。 FIR滤波器的主要特点有：FIR滤波器主要采用非递归结构，不论在理论上还是在实际的有限精度运算中都不存在稳定性问题，运算误差也较小。 此外，FIR滤波器可以采用快速傅里叶变换算法，在相同阶数的条件下，运算速度可以快得多。FIR滤波器设计则一般没有封闭函数的设计公式。 窗口法虽然仅仅对窗口函数可以给出计算公式，但计算通阻带衰减等仍无显式表达式。 一般，FIR滤波器设计只有计算程序可循， 因此对计算工具要求较高。FIR滤波器则要灵活的多， 尤其是频率采样设计法更容易适应各种幅度特性和相位特性的要求， 可以设计出理想的正交变换、理想微分、线性调频等各种重要网络。 因而有更大适应性。