数字信号课件-8第七篇 有限脉冲响应数字滤波器设计.pptVIP

7.1 引言 7.2 线性相位FIR滤波器的特点 7.3 利用窗函数法设计FIR数字滤波器 7.4 利用频率采样法设计FIR数字滤波器 7.5 IIR与FIR数字滤波器的比较 改变N 相对形状几乎不发生变化，只能改变主辨宽度。（矩形窗只改变过度带，不能改变波动） 主峰1，肩峰0.0895 8.95% 窗函数的选择原则： 1.主辨带宽要尽可能小，以获得较陡的过渡带； 2.与主辨的幅度相比，旁辨应尽可能地小，以减少通带内、阻带内波动的最大振幅。 几种常用窗口函数！ 1. 矩形窗(Rectangle)： 主瓣宽度为4?/N 第一副瓣比主瓣低13dB 2.三角形窗: 主瓣宽度为8?/N 第一副瓣比主瓣低26dB 3.汉宁窗（Hanning) ： 三部分之和使旁辨相互抵消，较前面的窗口来说能量更集中在主辨，使主辨宽度增加到 。 主瓣宽度为8?/N 第一副瓣比主瓣低31dB 4.哈明(Hamming)窗: 主瓣宽度为8?/N 第一副瓣比主瓣低40dB 5.布莱克曼(Blackman)窗： 主瓣宽度为12?/N 第一副瓣比主瓣低57dB 几种时域窗的比较 常用窗函数的幅度特性 (a).矩形窗 (b).三角形窗 (c).汉宁窗 (d).哈明窗 (e).布莱克曼窗 理想低通滤波器加窗后的幅度特性 (a).矩形窗 (b).三角形窗 (c).汉宁窗 (d).哈明窗 (e).布莱克曼窗 五种窗函数的基本参数 窗函数 旁瓣峰值幅度/dB 近似过渡带宽 精确过渡带宽 矩形窗 -13 4π/N 1.8π/N 三角形窗 -25 8π/N 汉宁窗 -31 8π/N 哈明窗 -41 8π/N 布莱克曼窗 -57 12π/N 阻带最小衰减/dB 21 25 44 53 74 6.1π/N 6.2π/N 6.6π/N 11π/N 三、窗函数设计FIR滤波器的步骤 （1）根据性能要求确定待求滤波器的冲激响应hd(n)。 ? 如果根据通带阻带衰减和边界频率要求，可选用理想滤波器作为逼近函数，从而用理想滤波器的频率响应Hd(ejω)作傅立叶逆变换求出hd(n)。 ? 如果Hd(ejω)较复杂或不能用封闭公式表示，则不能用上式求出hd(n)。此时可对Hd(ejω)从ω=0到ω=2?（一个周期内）均匀采样M点，采样值为Hd(ej2?k/M)，k=0,1,…,M-1，则上式可近似写成 如果M选得较大，可以保证窗口内hM(n) 有效逼近hd(n)。实际上，一个周期内的 hM(n)可用M点IDFT(IFFT)求得。 （2）根据对过渡带及阻带衰减的要求，选择窗口函数的形式，并估计窗口长度N。 设待求滤波器的过渡带用??表示,它近似于窗口函数的主瓣宽度。因过渡带近似与窗口函数长度N成反比，N?C/ ??，A取决于窗口函数的形式。如矩形窗C=4?；哈明窗C=8?等。按过渡带及阻带衰减情况选择窗口函数的形式。原则是满足阻带衰减的前提下尽量选择主瓣窄的窗口函数。 （3）计算滤波器的冲激响应h(n)。 （4）验证技术指标是否满足要求。 设计出的滤波器频率响应用下式表示： 计算上式时可用FFT算法。 如果不满足要求，重复（2）——（4）。 例：设计第一类线性相位FIR高通数字滤波器，3dB截止频率 rad，阻带最小衰减 ，过渡带宽度 。用窗函数法设计。 解：根据设计要求，N必须取奇数。 (1) 确定逼近理想高通频响函数 : (2) 求hd(n) : 其中: (3) 选择窗函数类型，估算窗函数长度N： 根据阻带最小衰减 ，查P209表7.2.2 ，选择哈明窗。 但N必须为奇数，故取N=129。 * 第七章 有限脉冲响应数字滤波器的设计 * * 第七章 有限脉冲响应数字滤波器的设计 7.1 引言 FIR滤波器的冲激响应h(n)的非零值为有限长； FIR滤波器一定是稳定的； 极点 原点处有N－1阶重极点。 许多信号的传输都要求线性相位，FIR滤波器具有非常优良的线性相位特点，同时具有任意的幅度特点； FIR滤波器的设计也要求满足频率特性指标。 7.2 线性相位FIR滤波器的特点 一、线性相位条件 对于长度为N的h(n)，传输函数为 式中： ：幅度特性，是?的实函数，可能取负值。 ：相位特性。 1. ；? 为常数，第一类线性相位 2. ；?0为起始相位，第二类线性相位 严格说，2中的?(?)不具有线性相位。但1、2都满