第六章 有限长单位冲激响应（FIR）数字滤波器的设计方法.pptVIP

第六章 有限长单位冲激响应（FIR）数字滤波器的设计方法(2) 尚勇 6.3 窗函数设计法 设计方法 设计方法 设计方法 设计方法 设计方法 设计方法 设计方法 设计方法 设计方法 设计方法 设计方法 各种窗函数 各种窗函数 各种窗函数——三角形(Bartlett)窗 各种窗函数——三角形(Bartlett)窗 各种窗函数——汉宁(Hanning)窗(又称升余弦窗) 各种窗函数——汉宁(Hanning)窗(又称升余弦窗) 各种窗函数海明(Hamming)窗(改进升余弦窗) 各种窗函数海明(Hamming)窗(改进升余弦窗) 各种窗函数——布拉克曼(blackman)窗 各种窗函数 窗函数的设计步骤 窗函数的设计步骤-设计举例 6.4 频率抽样设计法 6.4 频率抽样设计法 6.4 频率抽样设计法 减小设计误差的方法 6.5 IIR与FIR数字滤波器的比较 6.5 IIR与FIR数字滤波器的比较 6.5 IIR与FIR数字滤波器的比较 6.5 IIR与FIR数字滤波器的比较 当理想Hd(ej?)在某些频点（如通带截止频率处）变化剧烈时，会引起较大的设计误差。一种减小这种误差的方法是对Hk在这些频点附近增加过渡点 通过增加过渡点可有效减小由H(k)内插得到的H(ej?)在?c附近的波动，从而减小设计误差。但过渡点的增加会使设计的滤波器的过渡带增大，同时过渡点的增加会使参数N增大，这样会增加设备的成本。故这种方法较适用于窄带滤波器的设计[通带短，点数少N]，该方法便于在频率域观察设计结果，比较直观 IIR滤波器 系统函数的极点可位于单位圆内任意位置，故可用较低的阶数获得高的选择性，所用的存贮单元少，经济高效 但它的相位无法保证是线性相位 而且一般选择性越好，相位的非线性越严重 性能： FIR滤波器 可以获得严格的线性相位 但由于极点固定在原点，故只能用较高的阶数才能达到好的选择性 对于同样的滤波器指标，FIR滤波器所需的阶数一般是IIR的5-10倍，因此成本较高，且信号的时延较大，故FIR滤波器经常用在要求具有线性相位的场合 结构： IIR滤波器必须采用递归结构，极点必须位于单位圆内，否则系统将不稳定，这时有限字长的处理会引起寄生振荡 FIR滤波器主要采用非递归结构，不存在稳定性的问题，运算误差也较小。另外FIR滤波器便于采用FFT技术进行相关运算，故运算速度可以快很多 * * 顾名思义，窗函数设计法就是利用加窗技术来设计线性相位FIR滤波器的。 加窗处理是在时域进行的，是一种时域设计方法。 其基本思想是，在给定理想线性相位FIR滤波器的频率响应函数Hd(ej?)的前提下，用有限长的一个h(n)去逼近理想的hd(n)。这种逼近可以看成是对hd(n)进行加窗截断处理而得到的。 已知理想滤波器的频率响应 Hd(ej?)（一般为具有矩阵形状的频率响应特性）。要求设计一个FIR滤波器响应 逼近Hd(ej?) 是一个无限长序列，若用有限长的h(n)去逼近它就要对hd(n)作加窗截断处理。设所用的窗函数为w(n) 则有 由于 由上式可以看到，在衡量h(n)与hd(n)的相似程度时，w(n)函数起到了关键作用。我们最直接所能够想到的窗函数为矩形窗函数。 下面就以矩形窗函数为加窗函数，以低通滤波器为例来分析窗函数的设计方法 设低通滤波器的通带截止频率为?c，低通滤波器的群时延为α，此时的Hd(ej?)如下 在通带内|Hd(ej?)|=1，相位为线性相位-?α，Hd(ej?)对应的hd(n)为 hd(n)是一个中心点在α的偶对称无限长的非因果序列 令 w(n)=RN(n) 则 h(n)=hd(n)RN(n) 思考： hd(n)为线性相位，而要求设计的h(n)也应满足线性相位的要求。回顾上节中关于线性相位FIR Filter的单位样值响应h(n)的约束条件，可知h(n)必须是偶对称的，对称中心应为长度的一半，即(N-1)/2，因而必须要求α=(N-1)/2，这样的h(n)才是具有线性相位频率特性的 故有 将（6-27）式的hd(n)代入上式可得 这里，由于α=(N-1)/2，故h(n)是满足线性相位约束条件的 得到了h(n)，就可以通过傅氏分析，找出所求FIR滤波器的频率特性 下面就来求解 h(n) 对应的频率响应函数 H(ej?)，并比较 H(ej?) 和 Hd(ej?) 的异同 由（6-25）式可得 由式（6-30）式可知，H(ej?) 逼近 Hd(ej?) 的程度关键取决于 W(ej?) 即取决于窗函数的频谱特性 若w(n)=RN(n)，有 其中 将Hd(ej?)也写成(6-31)的形式，可得 其中 将 (6-31)、(6-32) 式代入 (6-30) 式可得 (6-33)式说明：实际FIR滤波器的幅度函数 H(?) 是理想