第八章 有限冲激响应数字滤波器及设计.ppt

（3）改变 的值只会影响 坐标的比例、窗谱的主瓣宽度及窗函数的绝对值大小，而不会改变肩峰的相对值。增加窗函数的长度 ，只能减小窗函数 的主瓣宽度和各旁瓣宽度，但不能改变主瓣和旁瓣的相对比值，从而使 的通带和阻带内波动起伏变密，但 的相对振荡幅度不减小，这种现象称为吉布斯（Gibbs）效应。例如，对于矩形窗函数，当增加窗宽度 时，过渡带宽度 将随之减小，振荡起伏变密，但最大肩峰却总是8.95％，阻带最小衰减为 ，这在工程上往往满足不了要求，改善阻带衰减特性只能是改变窗函数。 * 8.2.3常用窗函数 对窗函数一般有两个方面的要求： （1）主瓣尽可能窄，以使设计出的滤波器具有较陡的过渡带； （2）旁瓣尽可能少，即应使其能量尽可能集中在主瓣内，使设计出的滤波器肩峰和余振较小，阻带衰减较大。 对任一具体窗函数来说，这两项要求相互矛盾，无法同时满足，只能根据具体的设计指标选择较为合适的窗函数。 以下介绍的窗函数均为偶对称函数，都具有线性相位特性。设窗的宽度为N ， N可为奇数或偶数，且窗函数的对称中心点在(N-1)/2 处。因此，均为因果函数。 * 1.矩形窗 长度为N 的矩形窗定义为 * 2. 三角窗（或巴特利特(Bartlett)窗） * 3. 升余弦窗 式中A,B,C为常数。 升余弦窗的频率特性比矩形窗有很大改善。根据A,B,C的不同取值，将得到汉宁窗、汉明窗和布莱克曼窗三种不同的升余弦窗。 * （1）汉宁（Hanning）窗—升余弦窗 此种情况下A=0.5，B=0.5，C=0。汉宁窗的定义为 * （2）汉明（Hamming）窗—改进的升余弦窗 此种情况下A=0.54，B=0.46，C=0。汉明窗的定义为 * （3）布莱克曼（Blackman）窗—二阶升余弦窗 此种情况下A=0.42，B=0.5，C=0.08。长度为N 的布莱克曼窗定义为 * 4.凯塞（Kaiser）窗 上面几种窗函数都是以牺牲主瓣宽度来换取对旁瓣的抑制，主瓣宽度与旁瓣衰减无法协调兼顾。而凯塞窗则能全面反映窗函谱主瓣宽度与旁瓣衰减之间的互换关系，从而实现以同一种窗类型来满足不同性能需求的目的。凯塞窗能够在同等性能下，实现最陡峭的过渡带。凯塞窗是一组参数可调的由零阶贝塞尔函数构成的窗函数，定义如下: * * 由于贝塞尔函数的复杂性，为了方便设计，凯塞提出了经验公式。 * 表8-2 凯塞窗参数对滤波器的性能影响 * 表8-3 常用的六种窗函数的比较 窗函数 旁瓣峰值衰减 主瓣宽度 加窗后过渡带宽?? 阻带最小衰减 矩形窗 -13 4?/N 1.8?/N ?21 三角窗 -25 8?/N 4.2?/N -25 汉宁（升余弦窗） -31 8?/N 6.2?/N -44 汉明（改进升余弦窗） -41 8?/N 6.6?/N -53 布莱克曼（二阶升余弦窗） -57 12?/N 11?/N -74 凯塞窗 β=7.865 -57 10?/N -80 * * 图8-8 常用窗函数的幅度特性 * 图8-9 理想低通加窗后的幅度特性 8.2.4窗函数法设计步骤 （1）根据给出的频率响应函数 Hd(ejω)，经过傅里叶反变换得到hd(n) ，如果要求的滤波器的频率响应 Hd(ejω)存在过渡带，则设计中所使用的截止频率 ωc由通带频率ωp和阻带频率ωs按下式求出 ωc=(ωp+ωs)/2 （2）选择窗函数，根据所允许的过渡带宽 ，按表8-3估计序列 的长度 N。 （3）计算数字滤波器单位冲激响应 * （4）用选择的窗函数对 h(n)进行加窗 (5)计算滤波器的频率响应 检验其是否符合要求，如不符合要求修改有关参数，重复上述步骤直到满意为止。 * 1.数字低通滤波器的设计 例8-1设计一个FIR低通滤波器，所希望的频率响应为 如果取N=21 ，观察加窗后对滤波器幅频特性的影响。 * * 图8-11 用矩形窗和汉明窗设计的FIR低通滤波器 2.数字高通、带通和带阻滤波器的设计 例8-4 设计一个FIR高通滤波器，所希望的频率响应为 取 N=21，试用矩形窗和汉明窗设计FIR。 * * 图8-17 例8-4用矩形窗和汉明窗设计的FIR高通滤波器 8.3 利用频率采样法设计FIR滤波器 前面讨论的窗函数设计法是在时域内，用有限长冲激响应 去逼近所要求的理想冲激响应 ，然后用窗函数对加以修正，得到频率响应 逼近理想的频率响应 。 本节所讨论的频率采样法则是在频域内，以有限个频率响应采样，去近似所要求的理想频率响应 的方法。 * 8.3.1设计思路 设计指标通常是在频域给出理想的频率响应 Hd(ejω)和允许的误差。由一个长度为N 的序列可以用 N个频率采样值唯一确定，所以可以直接从频域出发，对理想频响采样。 * h(n)为偶对