数字信号处理第5章FIR数字滤波器的设计2.pptVIP

FIR数字滤波器 线性相位FIR数字滤波器的性质 窗函数法设计FIR数字滤波器 频率取样法设计线性相位FIR数字滤波器 线性相位FIR数字滤波器的优化设计 利用MATLAB设计FIR数字滤波器 窗函数法设计FIR数字滤波器 设计方法 吉伯斯现象 常用窗函数 设计方法 例：设计一个幅度响应能逼近理想带通滤波器的线性相位FIR滤波器。 例：设计一个幅度响应能逼近理想带通滤波器的线性相位FIR滤波器。 例：设计一个幅度响应能逼近理想带通滤波器的线性相位FIR滤波器。 例：设计一个幅度响应能逼近理想带通滤波器的线性相位FIR滤波器。 例：设计一个幅度响应能逼近理想带通滤波器的线性相位FIR滤波器。 吉伯斯现象 吉伯斯现象 吉伯斯现象 吉伯斯现象 吉伯斯现象 吉伯斯现象 吉伯斯现象 吉伯斯现象 常用窗函数 常用窗函数 常用窗函数 常用窗函数 Kaiser(凯泽)窗( w=kaiser(N,beta) ) Kaiser(凯泽)窗——b取不同值时窗的形状 * FIR窗函数法 问题：理想滤波器的Hd (ejW )已知，设计一物理可实现的滤波器，使其频响特性逼近Hd (ejW )。 hd [k]一般情况下是非因果无限长序列， 需对其进行截短和因果化处理。 设计方法 1．根据要求确定线性相位FIR滤波器的类型 2．确定理想滤波器的幅度函数Ad(W)和相位fd(W) fd(W)= -0.5MW +b 3．计算IDTFT得hd[k] 4．截短hd[k] h[k]= hd[k]， 0?k ?M 设计FIR的步骤 解： (1) 确定线性相位FIR滤波器类型：可选用I型或II型 (2) 确定理想滤波器的幅度函数Ad(W)和相位fd(W) fd(W)= -0.5MW 解： (3) 计算IDTFT得hd[k] (4) 截断hd[k] 长度为N=M+1的矩形窗 解： 单位脉冲响应 幅度函数 Wc1=0.3p，Wc2=0.7p，M=30 结论： 当Wc1=0，滤波器为理想低通 当Wc2=p，滤波器为理想高通 1. 理想低通和高通滤波器是带通滤波器的特殊情况 结论： 2. 理想带阻滤波器可通过带通滤波器得到 由于HBS(ejW)=1- HBP(ejW)，故 其中 0 0.25 0.5 0.75 1 -40 -30 -21 -10 0 M=14 M=30 矩形窗设计的FIR低通滤波器在不同阶数时的幅度函数(Wc=p/2) 所设计滤波器的幅度函数在通带和阻带都呈现出振荡现象，且最大波纹大约为幅度的9%，这个现象称为Gibbs 现象。 产生吉伯斯现象原因分析 由于所设计滤波器的单位脉冲响应 N=M+1 利用DTFT的性质可得所设计FIR滤波器的频率响应 H(ejW)逼近Hd(ejW)的好坏，取决于窗函数的频谱W(ejW) 窗函数的频谱 窗函数的频谱 2 p N N 0 4 p N W p 3 p N 矩形窗的幅度函数W(W) 1.W(W)的主瓣宽度 4p/N 2. 旁瓣相对衰减为常数 =13dB 矩形窗设计FIR滤波器的频率响应H(ejW) H(ejW)主要由主瓣的面积确定 H(ejW)的波动由旁瓣引起 矩形窗设计FIR滤波器的频率响应H(ejW) H(ejW)逐渐减小，形成了滤波器的过渡带 FIR滤波器过渡带的宽度和窗函数主瓣的宽度密切相关 矩形窗设计FIR滤波器的频率响应H(ejW) H(ejW)完全由旁瓣的面积确定 旁瓣的大小决定了FIR滤波器在阻带的衰减 矩形窗设计FIR滤波器的频率响应H(ejW) 矩形窗截断产生的波峰大约是9%，故用矩形窗设计出的滤波器阻带最大衰减 -20log10(9%)?21dB 结论 1. 窗函数的主瓣宽度决定了H(ejW )过渡带的宽度,窗函数长度N增大， 过渡带减小。 2. 旁瓣的大小决定了FIR滤波器在阻带的衰减 用矩形窗设计出的滤波器阻带最大衰减为 -20log10(9%)?21dB 如何提高阻带衰减? 选用旁瓣幅度较小的窗函数 矩形窗 Ap =-20lg(1-dp) ? 0.82dB, As = -20lg(ds ) ? 21dB dp=ds=0.09 Hann(汉纳)窗(w=hanning(M+1)) Ap? 0.056dB, As? 44dB dp=ds =0.0064 Hamming(哈明)窗(w=hamming(M+1)) Ap? 0.019dB, As? 53dB dp=ds =0.0022 Blackman窗(w=blackman(M+1)) Ap? 0.0017dB，As? 74dB dp=ds =0.0002 b 是一可调参数，调节窗函数的