ch52FIR窗函数法.pptVIP

矩形窗和Blackman窗设计的低通对比 0 0.25 0.5 0.75 1 -100 -75 -60 -40 -20 0 Square Blackman 例：用非矩形窗设计截频为Wc =0.5p的线性相位理想低通，并将设计结果与利用矩形窗设计结果比较。 Kaiser(凯泽)窗( w=kaiser(N,beta) ) b 是一可调参数，调节窗函数的形状。 I0(x ): the modified zeroth-order Bessel function. I0(x )可用幂级数表示为 一般求20项就能达到所需精度。 定义 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 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滤波器的单位脉冲响应 例：用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 例：用Kaiser窗设计满足下列指标的I型线性相位FIR低通滤波器。Wp=0.3p, Ws=0.5p，Ap=0.1dB, As=40dB。 解： (3) 设计截频Wc =0.4p的I型线性相位FIR低通滤波器 fd(W)= -0.5MW 例：用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。 wp=0.3*pi; ws=0.5*pi; As=50; N=ceil((As-7.95)/(ws-wp)/2.285) N=N+mod(N,2) beta=0.1102*(As-8.7); w=kaiser(N+1,beta); wc=(wp+ws)/2; alpha=N/2; k=0:N; hd=(wc/pi)*sinc((wc/pi)*(k-alpha)); h=hd.*w; omega=linspace(0,pi,512); mag=freqz(h,[1],omega); magdb=20*log10(abs(mag)); plot(omega/pi,magdb); axis([0,1,-70,0]);grid; 利用MATLAB实现上例 Wp=0.3p, Ws=0.5p，Ap=0.1dB, As=40dB。 例：设计一满足下列指标的线性相位 FIR高通滤波器。 Wp=0.67p, Ws=0.53p，Ap=0.3dB, As=50dB。 解: 由As确定截断所用窗函数:Hamming窗或kaiser窗 1) 由过渡带宽度确定滤波器长度N 可选I型滤波器，取N=51；或选IV型滤波器，取N=50 2) 由给定指标确定待逼近理想高通的截频Wc =0.6p 一、采用Hamming窗截断