华侨大学自动化专业数字信号处理.pptVIP

窗函数的频谱 2 p N N 0 4 p N W p 3 p N 矩形窗的幅度函数W(W) 1. W(W)的主瓣宽度 4p/N 2. 旁瓣相对衰减为常数 =13dB 长度为N的矩形窗函数的频谱为 加窗对所设计滤波器的频率特性的影响 FIR滤波器的幅度函数为 结论 1. 在理想频率特性不连续点Wc附近形成一个过渡带，窗函数的主瓣宽度决定了过渡带的宽度。 如何提高阻带衰减? 选用旁瓣幅度较小的窗函数 2. 在过渡带的两边 ，A(W)出现最大的肩峰值，肩峰的两侧形成起伏振荡，其振荡幅度取决于旁瓣的相对幅度，而振荡的多少，取决于旁瓣的多少。 3. 窗函数长度N增大，过渡带减小，起伏振荡变密，但不会改变肩峰的相对值，例如在矩形窗情况下，通带和阻带最大波纹幅度大约为滤波器幅度的9%，则阻带最大衰减为20lg(9%)?-21dB 矩形窗 Ap =-20lg(1-dp) ? 0.82dB, As = -20lg(ds ) ? 21dB dp=ds=0.09 3.常用窗函数 Hann(汉纳)窗(也称为Hanning汉宁) Ap? 0.056dB, As? 44dB dp=ds =0.0064 Hamming(汉明)窗 Ap? 0.019dB, As? 53dB dp=ds =0.002 2 Blackman窗(布莱克曼) Ap? 0.0017dB，As? 74dB dp=ds =0.0002 常用窗函数性质 例：设计一满足下列指标的I型线性相位 FIR高通滤波器。 Wp=0.67p, Ws=0.53p，Ap=0.3dB, As=50dB。 解: 由As确定截断所用窗函数:Hamming窗或Blackman窗 1) 由过渡带宽度确定滤波器长度N 2) 由给定指标确定待逼近理想高通的截频Wc =0.6p 采用Hamming窗截断，设计过程如下 I型滤波器，取N=51 由于理想高通滤波器的截频Wc是过渡带的中点，因此通常取为 3) 设计截频Wc =0.6p的I型线性相位FIR高通滤波器 解: fd(W)= -0.5MW 截断，得I型线性相位FIR高通滤波器的单位脉冲响应 h[k]=hd[k]w51[k] 例：设计一满足下列指标的I型线性相位 FIR高通滤波器。 Wp=0.67p, Ws=0.53p，Ap=0.3dB, As=50dB。 用Hamming窗设计的Wc=0.6p ,N=51 的FIR HP滤波器的幅度响应 b 是一可调参数，调节窗函数的形状。 I0(x ): 零阶第一类修正贝塞尔函数，可用幂级数表示为 一般求20项就能达到所需精度。 定义 4. Kaiser(凯泽)窗 Kaiser(凯泽)窗——b取不同值时窗的形状 0 2 4 6 8 10 12 14 16 18 20 0 0.5 1 b=0 0 2 4 6 8 10 12 14 16 18 20 0 0.5 1 b=1 0 2 4 6 8 10 12 14 16 18 20 0 0.5 1 b=3 A= -20lg (min{dp，ds }) b 与M的确定 例：用Kaiser窗设计满足下列指标的I型线性相位FIR低通滤波器。Wp=0.3p, Ws=0.5p，Ap=0.1dB, As=40dB。 解： (1) 由给定指标确定待逼近理想低通的截频Wc Wc =(Wp + Ws )/2=0.4p (2) 设计截频Wc =0.4p的I型线性相位FIR低通滤波器 解： (3) 由给定指标确定Kaiser窗的参数b和M A= -20lg(min{dp，ds })=As=40dB I型线性相位滤波器阶数必须是偶数，取M=24 例：用Kaiser窗设计满足下列指标的I型线性相位FIR低通滤波器。Wp=0.3p, Ws=0.5p，Ap=0.1dB, As=40dB。 解： (4) 利用Kaiser窗将hd[k]截断 例：用Kaiser窗设计满足下列指标的I型线性相位FIR低通滤波器。Wp=0.3p, Ws=0.5p，Ap=0.1dB, As=40dB。 解： 单位脉冲响应 设计结果 增益响应 G(W)= 20lg |H(ejW)| 例：用Kaiser窗设计满足下列指标的I型线性相位FIR低通滤波器。Wp=0.3p, Ws=0.5p，Ap=0.1dB, As=40dB。 本章小结 一、线性相位FIR数字滤波器的特性 线性相位系统的定义 线性相位条件 线性相位系统的频域特性 线性相位系统H(z)的零点分布特性 二、窗函数法设计线性相位FIR数字滤波器 基本思想 Gibbs（吉伯斯）现象 常用窗函数 Kaiser(凯泽)窗 线性相位FIR数字滤波器的优化设计 FIR