FIR低通滤波器加窗效应分析

1、目录一.摘要 2二.引言 3三.FIR滤波器设计 4 3.1 线性相位FIR滤波器的条件与特点 3.2 用窗函数法设计FIR滤波器的基本原理 3.3 用窗函数法设计FIR滤波器的一般步骤 3.4 FIR滤波器加窗效应分析 3.5 几种常用窗函数简介 四 MATLAB仿真滤波实现 14 4.1 MATLAB软件简介 4.2 设计中主要用到的MATLAB函数 4.3 实验程序及结果分析五 心得体会与总结 21六 参考文献 22一.摘要 数字滤波器一词出现在60年代中期。由于电子计算机技术和大规模集成电路的发展，数字滤波器已可用计算机软件实现，也可用大规模集成数字硬件实时实现。数字滤波器是一个离散时

2、间系统（按预定的算法，将输入离散时间信号(对应数字频率)转换为所要求的输出离散时间信号的特定功能装置）。应用数字滤波器处理模拟信号（对应模拟频率）时，首先须对输入模拟信号进行限带、抽样和模数转换。数字滤波器输入信号的数字频率（2*f/fs,f为模拟信号的频率，fs为采样频率，注意区别于模拟频率),按照奈奎斯特抽样定理，要使抽样信号的频谱不产生重叠，应小于折叠频率（ws/2=)，其频率响应具有以2为间隔的周期重复特性，且以折叠频率即=点对称。为得到模拟信号，数字滤波器处理的输出数字信号须经数模转换、平滑。数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。数字滤波器在语言信号处

3、理、图像信号处理、医学生物信号处理以及其他应用领域都得到了广泛应用。数字滤波器有低通、高通、带通、带阻和全通等类型。它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。应用最广的是线性、时不变数字滤波器，以及f.i.r滤波器。22二.引言 随着信息技术的迅猛发展，数字信号处理已成为一个极其重要的学科和技术领域。在通信、语音、图像、自动控制和家用电器等众多领域得到了广泛的应用。数字滤波是数字信号处理的重要环节，它在数字信号处理中占有着重要的地位,它具有可靠性好、精度高、灵活性大、体积小、重量轻等优点。随着数字技术的发展，数字滤波器越来越受到人们的重视，广泛地应用于各个领域。数字滤波器的

4、输入输出信号都是数字信号，它是通过一定的运算过程改变输入信号所含频率成分的相对比例或者滤除某些频率成分来实现滤波的，这种运算过程是由乘法器、加法器和单位延迟器组成的。数字滤波器是数字信号处理技术的重要内容，其对数字信号进行的最常见处理是保留数字信号中的有用频率成分和去除信号中的无用频率成分。按照时间域的特性，数字滤波器可以分为无限冲激脉冲响应数字滤波器（IIR滤波器）和有限冲激脉冲响应数字滤波器（FIR滤波器）。 三 FIR滤波器设计3.1 线性相位FIR滤波器的条件与特点 FIR DF的系统函数无分母，为，系统频率响应可写成：，令，H(w)称为幅度函数，称为相位函数。这与模和幅角的

5、表示法有所不同，H(w)为可正可负的实数，这是为了表达上的方便。 线性相位的FIR滤波器是指其相位函数满足线性方程： （是常数）根据群时延的定义，式中表示系统群时延，表示附加相移。线性相位的FIR系统都具有恒群时延特性，因为为常数，但只有0的FIR系统采具有恒相时延特性。第一类FIR DF的特点：恒相时延，相位曲线是过原点的曲线；可通过h(n)灵活设计幅度函数的零点位置；幅度函数对频率轴零点偶对称，对点偶对称。第二类FIR DF的特点：恒相时延，相位曲线是过原点的直线；幅度函数对频率轴零点偶对称；幅度函数对频率轴点奇对称。由的连续性，点一定是幅度函数的零点。即时，在z=-1处有零点；因此这类滤

6、波器不适合高通或带阻滤波器。第三类FIR DF的特点：恒群时延，有附加相移，相位曲线是截距为、斜率为的直线；幅度函数对零频点奇对称，零频是的零点；对奇对称，也是的零点。第四类FIR DF的特点：恒群时延，有附加相移，相位曲线是截距为、斜率为的直线；幅度函数对零频点奇对称，零频是的零点；对偶对称。3.2 用窗函数法设计FIR滤波器的基本原理 设所希望得到的滤波器的理想频率响应为。那么FIR滤波器的设计就在于寻找一个传递函数去逼近。在这种逼近中最直接的一种方法是从单位取样响应序列着手，使逼近理想的单位取样响应。我们知道可以从理想频率响应通过傅里叶反变换来得到，即： 但是一般来说，这样得到的单位取样

7、响应往往都是无限长序列；而且是非因果的。我们以一个截止频率为的线性相应位理想低通为例来说明。设低通滤波器的时延为，即： 则这是一个以为中心的偶对称的无限长非因果序列。这样一个无限长的序列怎样用一个有限长序列去近似呢？最简单的办法就是直接截取它的一段来代替它。例如把到的一段截取来作为，但是为要保证所得到的是线性相位滤波器。必须满足的对称性，所以时延应该取长度的一半，即这种直接截取的办法可以形象地想象为：好比是通过一个“窗口”所看到的一段。中表达为和一个“窗口函数”的乘积。在这里，窗口函数就是矩形脉冲函数，即 但是一般来说，窗口函数并不一定是矩形函数，可以在矩形以内还对作一定的加权处理，因此，一般

8、可以表示为 这里就是窗口函数。按照复卷积公式，在时域中的乘积关系可表示成在频域中的周期性卷积关系，即可得所设计的FIR滤波器的频率响应:其中，为截断窗函数的频率特性。由此可见，实际的FIR数字滤波器的频率响应逼近理想滤波器频率响应的好坏，完全取决于窗函数的频率特性。如果w(n)具有下列形式：w(n)相当于一个矩形，我们称之为矩形窗。即我们可采用矩形窗函数w(n)将无限脉冲响应截取一段来近似为。3.3用窗函数法设计FIR滤波器的一般步骤（1）确定逼近理想滤波器的频率响应函数；（2）求出理想的单位抽样响应；（3）根据过渡带宽和阻带最小衰件等要求，选择窗函数，并确定窗口长度；（4）求所设计的FIR滤

9、波器的单位抽样响应；（5）计算频率响应，验算指标是否满足要求，若不满足，则要重新设计3.4 FIR滤波器加窗效应分析经过加矩形窗后所得的滤波器实际频率响应能否很好地逼近理想频率响应呢？下图给出了理想滤波器加矩形窗后的情况。理想低通滤波器的频率响应如图中左上角图，矩形窗的频率响应为左下角图。根据卷积定理，即得实际滤波器的频率响应图形为图中右图。 图（1）由图可看出，加矩形窗后使实际频率响应偏离理想频率响应，主要影响有两个方面：（1）在理想幅频特性陡直边缘处形成过渡带，过渡带宽取决于矩形窗函数频率响应的主瓣宽度。（2）过渡带两侧形成肩峰和波纹，这是矩形窗函数频率响应的旁瓣引起的，旁瓣相对值越大，旁

10、瓣越多，波纹越多。3.5 几种常用窗函数简介1.矩形窗矩形窗函数的时域形式可以表示为：它的频域特性为：2.汉宁窗函数汉宁窗函数的时域形式可以表示为： 它的频域特性为：其中，为矩形窗函数的幅度频率特性函数。汉宁窗函数的最大旁瓣值比主瓣值低31dB，但是主瓣宽度比矩形窗函数的主瓣宽度增加了1倍，为8/N。3.海明窗函数海明窗函数的时域形式可以表示为: 它的频域特性为:其中，为矩形窗函数的幅度频率特性函数。海明窗函数的最大旁瓣值比主瓣值低41dB，但它和汉宁窗函数的主瓣宽度是一样大的。4.布莱克曼窗 增加一个二次谐波余弦分量，可进一步降低旁瓣，但主瓣宽度进一步增加，增加N可减少过渡带。频谱的幅度函数

11、为： +0.045.三角窗函数三角窗是最简单的频谱函数为非负的一种窗函数。三角窗函数的时域形式可以表示为：当n为奇数时：当n为偶数时：它的频域特性为：三角窗函数的主瓣宽度为8/N，比矩形窗函数的主瓣宽度增加了一倍，但是它的旁瓣宽度却小得多。6.凯泽窗 以上几种窗函数，都是以增加主瓣宽度为代价来降低旁瓣。凯泽窗则可自由选择主瓣宽度和旁瓣衰减,如图（2）。 图（2）式中I0(x)是零阶贝塞尔函数，参数可自由选择，决定主瓣宽度与旁瓣衰减。越大，w(n)窗越窄，其频谱的主瓣变宽，旁瓣变小。四 MATLAB仿真滤波实现4.1 MATLAB软件简介MATLAB是由美国mathworks公司发布的主要面对科

12、学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学计算软件的先进水平。4.2设计中主要用到的MATLAB 函数4.2.1 fir1 函数 fir1 函数实现线性相位FIR 数字滤波器的窗函数法设计. 利用这一函数可以设计出标准的低通、高通、带通和带阻线性相位FIR滤波器。调用格式：(1) b=fir1(

13、n,wn) 用海明窗设计低通或带通的 FIR 滤波器. wn 是截止频率，以“ð 弧度”为单位。 fir1函数返回FIR滤波器的系数向量b。如果wn 是一个标量，则返回的是一个n 阶的低通FIR 数字滤波器。 (2) b=fir1(n,wn, 'type', 'window') type为滤波器的类型，当type=high 时，设计高通FIR滤波器；当type=stop 时，设计带阻FIR滤波器。向量window 用来指定窗函数类型，其长度为n+1，缺省时为海明窗。4.2.2 freqz 函数 为了检验所设计的数字滤波器是否正确，可以画出其幅频特性.

14、MATLAB 提供了求解数字滤波器频率响应的函数freqz ，其调用格式如下：h,w=freqz(b,a,n)返回数字滤波器的n点频率响应h和频率向量w。若n缺省，则自动取n为默认值512。4.2.3 plot函数plot(Y)如果Y是m×n的数组，以1:m为X横坐标，Y中的每一列元素为Y坐标，绘制n条曲线；如果Y是n×1或者1×n的向量，则以1:n为横坐标，Y为坐标表绘制1条曲线；如果Y是复数，则plot(Y)等效于plot(real(Y),imag(Y)；其它使用情况下，忽略坐标数据中的虚部。4.2.4 title函数 title就是给已经画出的图加一个标题.

15、title ('')里面的内容就是图片标题的名称。4.2.5 legend函数 legend(string1,string2,string3, .)分别将字符串1、字符串2、字符串3标注到图中，每个字符串对应的图标为画图时的图标。4.2.6 xlabel（） ylabel（）函数Xlabel（）函数给x坐标轴添加标注，ylabel（）函数给y坐标轴添加标注。4.2.7 hold on Hold on命令可以保护当前的坐标系，使以后再使用plot（）函数时将新的曲线叠印在原来的图上，用hold off则可以取消保护状态。4.2.8 axisaxis（【xmin xmax ymin

16、 ymax】）设置坐标轴的极限范围。4.3实验程序及结果分析（1）相同的窗口加不同的窗长(M=10,20,40)程序一、汉宁窗不同窗长wc=0.25;N=11 21 41;M=N-1;hl1=fir1(M(1),wc,hanning(N(1);hl2=fir1(M(2),wc,hanning(N(2);hl3=fir1(M(3),wc,hanning(N(3);h1,w11=freqz(hl1,1,512);h2,w12=freqz(hl2,1,512);h3,w13=freqz(hl3,1,512);plot(w11/pi,20*log(abs(h1)/max(abs(h1);hold on

17、plot(w12/pi,20*log(abs(h2)/max(abs(h2),':');hold onplot(w13/pi,20*log(abs(h3)/max(abs(h3),'-');grid on;axis(0 1 -350 10);title('汉宁窗不同窗长');xlabel('/'); ylabel('|H(ej)|');legend('M=10','M=20','M=40'); 图（3） 图（3）分别给出了M=10，M=20及M=40时的的幅频特性曲线

18、。由图可知当M的值逐渐增大，即窗口的长度逐渐变长时，过渡带的宽度逐渐变窄，但阻带内的波动幅度没有多大改善。说明调整窗口长度可以有效的控制过度带的宽度，但不改变旁瓣的相对宽度，不能加大阻带的衰减。 （2）不同窗口加相同的窗长程序二、同窗长的汉宁窗与矩形窗比较wc=0.25;N=21;M=N-1;hl1=fir1(M,wc,hanning(N);hl2=fir1(M,wc,boxcar(N);h1,w11=freqz(hl1,1,512);h2,w12=freqz(hl2,1,512);plot(w11/pi,20*log(abs(h1)/max(abs(h1);hold onplot(w12/pi,20*log(abs(h2)/max(abs(h2),':');grid on;axis(0 1 -350 10);title('同窗长的汉宁窗与矩形窗比较