第7章有限脉冲响应数字滤波器设计.ppt

第7章 有限脉冲响应数字滤波器的设计 7.1 线性相位FIR数字滤波器的条件和特点 7.2 利用窗函数法设计FIR滤波器 7.5 IIR和FIR数字滤波器的比较 用 N 表示FIR滤波器单位脉冲响应h(n)的长度，其系统函数H(z)为 H(z)是z－1的N－1次多项式， 有N－1个零点，在原点z=0处有一个N－1重极点。 因此，H(z)永远稳定。 稳定和线性相位特性是FIR滤波器最突出的优点。 FIR滤波器设计任务： 选择有限长度的h(n)，使频率响应函数H(ejω)满足技术指标要求。 7.1 线性相位FIR数字滤波器的条件和特点 1. 线性相位FIR数字滤波器 对于长度为N的h(n) ，传输函数为 线性相位FIR滤波器是指θ(ω)是ω的线性函数，即 (2)式θ(ω)不具有线性相位，但以上二式都满足群时延（相位特性曲线的斜率）是一个常数，即 满足(1)式是第一类线性相位； 满足(2)式是第二类线性相位。 7.5 IIR和FIR数字滤波器的比较 IIR与FIR滤波器各有所长，在实际应用中应该考虑多方面因素加以选择，同时也要兼顾经济方面的要求以及计算工具的实际条件。 其幅度函数为 　　　　　　　　　　　　　　　　　　 其幅度函数由五部分组成，它们都是移位不同，且幅度也不同的WRg(ω)函数，使旁瓣再进一步抵消。旁瓣峰值幅度进一步增加，其幅度谱主瓣宽度是矩形窗的3倍。布莱克曼窗的四种波形如图7.2.8所示，参数为: ?n = ?57 dB  ΔB =12π/N ?s = ?74 dB 图7.2.8 布莱克曼窗的四种波形 6． 凯塞—贝塞尔窗（Kaiser-Basel Window） 　　以上五种窗函数都称为参数固定窗函数，每种窗函数的旁瓣幅度都是固定的。凯塞—贝塞尔窗是一种参数可调的窗函数，是一种最优窗函数。 　　　　　　　　　　　　　　　 （7.2.15） 式中  I0(β)是零阶第一类修正贝塞尔函数，可用下面级数计算：  一般I0(β)取15~25项，便可以满足精度要求。 ?参数可以控制窗的形状。一般?加大，主瓣加宽，旁瓣幅度减小，典型数据为4< ? <9。 当? =5.44时，窗函数接近哈明窗。 ? =7.865时，窗函数接近布莱克曼窗。 在设计指标给定时，可以调整?值，使滤波器阶数最低，所以其性能最优。凯塞（Kaiser）给出的估算β和滤波器阶数N的公式如下: （7.2.17） 式中，Bt=|ωs－ωp|, 是数字滤波器过渡带宽度。应当注意，因为式（7.2.17）为阶数估算，所以必须对设计结果进行检验。另外，凯塞窗函数没有独立控制通带波纹幅度，实际中通带波纹幅度近似等于阻带波纹幅度。凯塞窗的幅度函数为 （7.2.16） （7.2.18） 对?的8种典型值，将凯塞窗函数的性能列于表7.2.1中，供设计者参考。由表可见, 当? =5.568时, 各项指标都好于哈明窗。6种典型窗函数基本参数归纳在表7.2.2中，可供设计时参考。 表7.2.1 凯塞窗参数对滤波器的性能影响 表7.2.2 6种窗函数的基本参数 表中过渡带宽和阻带最小衰减是用对应的窗函数设计的FIR数字滤波器的频率响应指标。 图7.2.4 常用的窗函数 随着数字信号处理的不断发展，学者们提出的窗函数已多达几十种，除了上述6种窗函数外，比较有名的还有Chebyshev窗、Gaussian窗等。 MATLAB信号处理工具箱提供了14种窗函数的产生函数，下面列出上述6种窗函数的产生函数及其调用格式： wn=boxcar(N) 　 　％列向量wn中返回长度为N的矩形窗函数w(n) wn=bartle