FIR滤波器组的设计与Matlab仿真.docVIP

FIR滤波器组的设计与Matlab仿真 1 引言 在许多数字信号处理系统中，如语音或音频信号处理中，有限冲激响应（FIR）滤波器是最常用的组件之一，它完成信号预调、频带选择和滤波等功能。FIR滤波器虽然在截止频率的边沿陡峭性能上不及无限冲激响应（IIR）滤波器，但是却具有严格的线性相位特性，稳定性好，能设计成多通带（或多阻带）滤波器组，所以能够在数字信号处理领域得到广泛的应用。 Matlab是mathswork公司推出的一套高性能的数值计算和可视化软件，它集数值分析、信号运算、矩阵运算、信号处理和图形显示于一体，具有很强的绘图功能。利用它自带的丰富的函数和信号处理工具箱，可以很方便地设计和仿真FIR滤波器组。 2 用窗函数法设计FIR滤波器 一个离散时间系统H（z）=B（z）/A（z），若分母多项式A（z）的系数a1=…=aN=0,那么该系统即变成FIR系统，即 （1） 显然，系数b0，b1，…，bM即是该系统的单位抽样响应h(0) ，h(1) ，…，h(M) ，且当nM时，h(n)≡0。所对应的FIR滤波器的转移函数为 （2） H（z）的频率响应应近似为Hd(ejw)，而且是线性相位的。 FIR滤波器有许多行之有效的设计方法，如窗函数法、频率抽样法及最佳一致逼近法等，其中窗函数法很适合于带通滤波器组的设计。 一个截止频率为ωc的理想低通滤波器Hd(ejw)，其冲激响应hd(n)是无限长和非因果的。对hd(n)加窗（窗函数用w(n)表示，窗宽为有限值M）的结果，便得到一个FIR滤波器h(n)，它是对所要求的理想低通滤波器的逼近。如果窗函数w(n)和FIR滤波器h(n)的频率特性分别用W(ejw)和H(ejw)表示，则H(ejw)将等于W(ejw)和Hd(ejw)的卷积。由于窗函数频谱主瓣不是无限窄，因而FIR滤波器的频率特性在通带和阻带之间形成了一个过渡带。由于窗函数频谱中存在着旁瓣，因而FIR滤波器的频率特性在通带和阻带内产生了一些起伏振荡的波纹。 窗函数法设计的基本思想是对于给定的滤波器技术要求，选择滤波器长度M和具有最窄主瓣宽度和尽可能最小的旁瓣衰减的某个窗函数w(n)。常用的窗函数有矩形窗、巴特利特（Bartlett）窗、汉宁（Hanning）窗、哈明（Hamming）窗、布莱克曼（Blackman）窗和凯塞（Kaiser）窗，他们的主要区别在于参数主瓣3dB带宽B、最大边瓣峰值A（dB）和边瓣谱峰渐进衰减速度D（Db/oct）之间的差别。 3 FIR滤波器组的设计与Matlab仿真 3.1 FIR滤波器组的设计 用窗函数法设计FIR滤波器组的基本思路为: (1)根据所要求的分析频率范围和频率分辨率选择一组分析频率ωk及相对应的原型理想低通滤波器的截止频率ωck（k=0，1，…，L-1），L为滤波器组的通道数，得出理想带通滤波器组的复合频率响应。(2)用窗函数法近似实现各个理想带通滤波器。实际的复合频率响应等于理想复合频率响应与设计窗的频率特性的卷积。在用窗函数法设计FIR滤波器时，还需要知道通带内的纹波幅度δ，它可以根据阻带衰减量As的要求来确定，通常的典型值为 As=-20lgδ=40∽60 （3） 如果选择Kaiser窗，参数β的计算公式为 当As 50时，β=0.1102(As-8.7)； 当21As 50时，β=0.5845(As-21)0.5+0.07886(As-21)。 （4） 滤波器的阶数为 M=(As-7.95)/14.36△F+1 (5) 式中△F是归一化过渡带宽，即△F=△f?T，而 △f=△ω/2π （6） 过渡带宽△ω的选择根据对每个滤波器的频率分辨率要求来决定，必须满足限制条件 △ωk2ωck (7) 如果各通道采用相同的设计窗，保证中心频率及带宽的选择能覆盖-π到π频率范围（不管滤波器的中心频率和带宽如何分布），那么最终得到的实际复合响应一定是理想的，即具有平坦幅度和线性相位。 在数字语音处理中，语音信号取样频率f s=9.6kHz，频谱分析范围为200Hz到3200Hz，设计一FIR带通滤波器组，使阻