第课FIR数字滤波器的设计.ppt

第7章 FIR数字滤波器的设计 线性相位FIR滤波器的性质 窗函数法设计FIR滤波器 频率取样法设计线性相位FIR滤波器 线性相位FIR滤波器的优化设计 窗函数法设计FIR滤波器 频率取样法 积分平方误差定义为 由Parseval等式，e 2可表示为 可选择 h[k]= hd[k]， 0?k ?M 使积分平方误差最小。 吉伯斯(Gibbs)现象 0 0 0.5 1 0.2 p 0.4 p 0.6 p 0.8 p p W A (W) M=14 M=60 0.2 p 0 0 0.5 1 0.4 p 0.6 p 0.8 p p A (W) W 矩形窗对H(e jW)的影响 矩形窗的幅度函数为 W N p 2 N p 2 - N p 4 N p 4 - N p p - 主瓣 旁瓣 ) ( W W 矩形窗的幅度函数 将理想滤波器的频率响应表示为 则可得FIR滤波器的频率响应为 所以FIR滤波器的幅度函数为 ) ( W A W c W c W - p p - q p p - c W c W - ) ( q - W W ) ( q d A ) ( W A W c W c W - p p - q c W c W - p p - ) ( q d A ) ( q - W W q ) ( q d A ) ( q - W W c W - c W p p - ) ( W A W c W c W - p p - 由矩形窗截断产生的波峰大约是9%，所以阻带最小衰减为 20log10(9%)?-21dB。 用矩形窗设计的Wc=p/2 FIR滤波器的幅度响应 0 0.25 0.5 0.75 1 -40 -30 -21 -10 0 M=14 M=30 Gain db 常用窗函数 矩形窗 Ap? 0.82dB, As? 21dB Hann(汉纳)窗(w=hanning(M+1)) Ap? 0.056dB, As? 44dB 由Hanning窗设计的Wc=p/2 FIR滤波器的频响特性(M=38) 0 0.25 0.5 0.75 1 -80 -60 -44 -20 0 Square Hanning Gain dB 0 0.25 0.5 0.75 1 -80 -60 -52 -20 0 Square Hamming Hamming（哈明）窗( w=hamming(M+1) ) Gain db Blackman窗 ( w=blackman(M+1) ) 0 0.25 0.5 0.75 1 -100 -75 -60 -40 -20 0 Square Blackman Gain db Kaiser窗( w=kaiser(M+1,beta) ) b 是一可调参数?。I0(b ): the modified zeroth-order Bessel function. I0(b )可用幂级数表示为 一般求20项就能达到所需精度。 用Kaiser窗设计FIR滤波器的步骤： 1. 估计滤波器的阶数M A= -20log10(min{dp，ds }) 2. 估计b 1. 估计滤波器的阶数M 2. 估计b b = 0.1102(As-b) 3. 设定理想低通的截频 4. h[k]=hd[k]*w[k] 例: 用Kaiser窗设计一满足下列指标的I型线性相位FIR低通滤波器。 Wp=0.2p, Ws=0.4p，Ap=0.3dB, As=50dB Gain response of lowpass FIR filter 0 0.1 0.2 0.3 0.4 0.5 -80 -50 -30 0 Normalized frequency Gain, dB M=30 Ap=0.0105 dB As=50.7524 dB 例：用Hamming窗设计一个逼近截频为Wc的线性相位FIR 高通滤波器 选I型FIR(M为偶) k?0.5M k=0.5M * * 线性相位FIR滤波器的性质 线性相位系统的时域特性 线性相位系统的频域特性 线性相位系统H(z)的零点分布特性 FIR滤波器的定义 M阶(长度N=M+1) 的FIR数字滤波器 FIR滤波器的特点 1)h[k]在有限范围内非零，系统总是稳定的。 2)容易设计成线性相位 3)可利用FFT实现 4)运算量比IIR大 FIR滤波器设计指标 严格线性相位定义 例： 单频信号exp(jW0 k)通过线性相位(LTI)系统的响应