《数字信号处理（第2版）》ch5_3频率取样法.pptVIP

数字信号处理(Digital Signal Processing) FIR数字滤波器 线性相位FIR数字滤波器的性质 窗函数法设计FIR数字滤波器 频率取样法设计线性相位FIR数字滤波器 线性相位FIR数字滤波器的优化设计 FIR数字滤波器的基本结构 利用MATLAB设计FIR数字滤波器 基本思想 四种线性相位系统的Hd[m] I型线性相位。 M=4, Hd(ejW)在5个取样点的Hd[m] 四种线性相位系统的Hd[m] II型线性相位。 M=5，Hd(ejW)在6个取样点的Hd[m] 四种线性相位系统的Hd[m] 四种线性相位系统的Hd[m] 频率取样法设计线性相位FIR DF的步骤 例：用频率取样法设计一个逼近截频Wc=0.5p的理想低通滤波器。 例：用频率取样法设计一个逼近截频Wc=0.5p的理想低通滤波器。 例：用频率取样法设计一个逼近截频Wc=0.5p的理想低通滤波器。 例：用频率取样法设计一个逼近截频Wc=0.5p的理想低通滤波器。 频率取样法的MATLAB实现 例：利用频率取样法设计一个II型 (M=63 )线性相位FIR低通滤波器，Wp=0.5p，Ws=0.6p。 %频率取样法设计II型线性相位低通FIR N=63;Wp=0.5*pi; m=0:(N+1)/2; Wm=2*pi*m./(N+1); mtr=floor(Wp*(N+1)/(2*pi))+2; Ad=[Wm=Wp]; Ad(mtr)=0.38; Hd=Ad.*exp(-j*0.5*N*Wm); Hd=[Hd conj(fliplr( Hd(2:(N+1)/2) ) )]; h=real(ifft(Hd)); w=linspace(0,pi,1000); H=freqz(h,[1],w); plot(w/pi,20*log10(abs(H)));grid; 利用1个过渡点设计的FIR低通滤波器的增益响应 利用2个过渡点设计的FIR低通滤波器的增益响应 频率取样法设计的高通滤波器的增益响应 增加一个过渡点设计的高通滤波器的增益响应 * 频率取样法 * 频率取样法 * 信号与系统系列课程组 国家电工电子教学基地 频率取样法设计线性相位FIR滤波器 基本思想 四种线性相位系统的Hd[m] 频率取样法设计线性相位FIR DF的步骤 频率取样法的MATLAB实现 已知Hd(ejW)在M+1点上的抽样值 {Hd (ejWm)；m=0,1,...,M}，设计 {h[k]; k=0,1, ...,M}， 使设计出的滤波器H(z)满足 如何确定h[k]？ 若频率取样点为 (I型取样) 则有 DFT Hd[m] 基本思想 Ad[m] 若所设计的滤波器为线性相位滤波器，则Hd[m]应满足线性相位条件。 基本思想 I 型(M为偶数，h[k]偶对称) b=0 一般, I型线性相位滤波器Hd[m]为 利用 Ad(W)=Ad(2p-W)，可得 改写 I型线性相位。 由5点Hd[m]求得的h[k] 满足线性相位条件h[k]= h[M-k] 周期为5 4 - 5 3 - 5 II 型(M为奇数，h[k]偶对称) b=0 一般, II型线性相位滤波器Hd[m]为 利用 Ad(W)=-Ad(2p-W), Ad(p)=0可得 改写 III 型(M为偶数，h[k]奇对称) b=?/2 一般, III型线性相位滤波器Hd[m]为 IV 型(M为奇数，h[k]奇对称) b= ?/2 一般, IV型线性相位滤波器Hd[m]为 (1) 根据所需设计的数字滤波器类型，确定线性相位FIR滤波器的类型。 (2) 获得Hd(ejW)在W?[0,2?)区间的M+1个取样点上的值Hd[m]，并使Hd[m] 满足线性相位条件，即 (3) 利用IDFT得到h[k] Ad[m] 解： (1) 确定线性相位FIR滤波器的类型，I, II型均可采用。 选用I型线性相位系统 (2) 获得Hd(ejW)在W?[0,2?)区间的M+1个取样点上的值Hd[m] 取M=10 解： (2) 获得Hd(ejW)在W?[0,2?)区间的11个取样点上的值Hd[m] (3) 利用IDFT得到h[k] 10 - 11 9 - 11 周期为11 解： 设计结果 Wc=0.5p 0 0.4p 0.8p 1.2p 1.6p 2p 0 0.2 0.4 0.6 0.8 1 M=10 M=30 Hd[m]及频率取样法设计的低通滤波器幅度函数 解： 设计结果 Wc=0.5p 0 0.4p 0.8p 1.2p 1.6p 2p -25 -20