数字信号处理FIR滤波器的设计课程介绍.ppt

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Kaiser(凯泽)窗 A= -20log10(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 例：用Kaiser窗设计满足下列指标的I型线性相位FIR低通滤波器。Wp=0.3p, Ws=0.5p，Ap=0.1dB, As=40dB 解： (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)| 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窗截断，设计过程如下 3) 确定线性相位FIR滤波器类型： 4) 设计截频Wc =0.6p的线性相位FIR高通滤波器 例：设计一满足下列指标的线性相位 FIR高通滤波器。 Wp=0.67p, Ws=0.53p，Ap=0.3dB, As=50dB 解: 由As确定截断所用窗函数:Hamming窗或kaiser窗 一、采用Hamming窗截断，设计过程如下 采用I型线性相位滤波器，N=51，M=50 fd(W)= -0.5MW 4) 设计截频Wc =0.6p的线性相位FIR高通滤波器 例：设计一满足下列指标的线性相位 FIR高通滤波器。 Wp=0.6