窗函数法设计线性相位FIR数字滤波器.pptVIP

* 三、常用窗函数 矩形窗 Ap =-20lg(1-dp) ? 0.82dB, As = -20lg(ds ) ? 21dB dp=ds=0.09 三、常用窗函数 Hann(汉纳/汉宁)窗(w=hanning(M+1)) Ap =-20lg(1-dp) ? 0.056dB, As=-20lg(ds ) ? 44dB dp=ds =0.0064 三、常用窗函数 Hamming(哈明/汉明)窗(w=hamming(M+1)) Ap? 0.019dB, As? 53dB dp=ds =0.0022 N 三、常用窗函数 Blackman (布莱克曼)窗(w=blackman(M+1)) Ap? 0.0017dB，As? 74dB dp=ds =0.0002 N 三、常用窗函数 常用窗函数性质 此处 N 为窗的长度。 例2：用非矩形窗设计截频为Wc =0.5p的线性相位理想低通滤波器，并将设计结果与利用矩形窗设计结果比较。 解： (1) 确定线性相位FIR滤波器类型：可选用I型或II型 fd(W)= -0.5MW (2) 确定理想滤波器的幅度函数Ad(W)和相位fd(W) (3) 计算IDTFT得hd[k] (4) 截断hd[k] 可以是汉宁窗、哈明窗或Blackman窗 矩形窗和汉宁窗设计的低通对比 0 0.25 0.5 0.75 1 -80 -60 -44 -20 0 Square Hanning 例2：用非矩形窗设计截频为Wc =0.5p的线性相位理想低通滤波器，并将设计结果与利用矩形窗设计结果比较。 矩形窗和哈明窗设计的低通对比 0 0.25 0.5 0.75 1 -80 -60 -52 -20 0 Square Hamming 例2：用非矩形窗设计截频为Wc =0.5p的线性相位理想低通滤波器，并将设计结果与利用矩形窗设计结果比较。 矩形窗和Blackman窗设计的低通对比 0 0.25 0.5 0.75 1 -100 -75 -60 -40 -20 0 Square Blackman 例2：用非矩形窗设计截频为Wc =0.5p的线性相位理想低通滤波器，并将设计结果与利用矩形窗设计结果比较。 四、Kaiser(凯泽)窗 b 是一可调参数，调节窗函数的形状。 I0(x ): 零阶第一类修正贝塞尔函数。 I0(x )可用幂级数表示为： 一般求20项就能达到所需精度。 定义 ( w=kaiser(N,beta) ) 四、Kaiser(凯泽)窗 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 ——b取不同值时窗的形状(M=20) 四、Kaiser(凯泽)窗 其中，A= -20lg (min{dp，ds }) N与b 的确定 四、Kaiser(凯泽)窗 Kaiser窗设计FIR滤波器的MATLAB函数 [M,Wc,beta,ftype] = kaiserord(f,a,dev) f: 表示需设计的FIR滤波器的频带。 a: B个元素的向量, 表示FIR滤波器在B个频带中的 幅度值。一般对通带取值为1，阻带取值为0。 dev: B个元素的向量, 表示FIR滤波器在B个频带中 的波动值。 返回参数M及beta，分别表示FIR滤波器阶数M及Kaiser窗的参数b。 返回参数Wc和ftype是函数fir1的调用参数。 h = fir1(M,Wc,ftype,window) M表示滤波器的阶数。 Wc表示理想FIR滤波器的B个频带。 ftype表示滤波器的类型，缺省值为空。 window是一长度为N+1的向量, 默认为Hamming窗。 四、Kaiser(凯泽)窗 获得FIR滤波器的单位脉冲响应 例3：用Kaiser窗设计满足下列指标的I型线性相位FIR低通滤波器。Wp=0.3p, Ws=0.5p，Ap=0.1dB, As=40dB。 解： (1) 由给定指标确定待逼近理想低通的截频Wc 由于理想低通滤波器的|H(ejW)|在截频Wc处收敛于0.5，因此常将截频Wc取在过渡带的中点 Wc =(Wp + Ws )/2=0.4p 解： (2) 由给定指标确定Kaiser窗的参数N和b A= -20lg(min{dp，ds })=As=40dB I型线性相位滤波器阶数必须是偶数，取M=24 例3：用Kaiser窗设计满足下列指标的I型线性相位FIR低通滤波器。Wp=0.3p, Ws=0.5p，Ap=0.1dB, As