数字信号处理25-6.3FIR DF设计(二).pptVIP

1 频率抽样法基本原理 二者对比 H (e j? ) 频率抽样 内插公式 问题：理想滤波器的Hd (ej? )已知，设计一物理可实现的滤波器，使其频响特性逼近Hd (ej? )。 目标 Hd(e j? ) H(k) * 1 频率抽样法基本原理 抽样 恢复(内插) 内插函数 * 2 线性相位约束 抽样值H(k)的幅度和相位要满足线性相位约束条件。 * 2 线性相位约束 1、 第Ⅰ种类型——h(n)偶对称，N为奇数 如果抽样值也写成纯标量模与相角的形式，就可以得到抽样值 H(k)的约束为： * 2 线性相位约束 2、 第Ⅱ种类型——h(n)偶对称，N为偶数 抽样值 H(k)的约束为： * 2 线性相位约束 3、 第Ⅲ种类型——h(n)奇对称，N为奇数 抽样值 H(k)的约束为： * 2 线性相位约束 4、 第Ⅳ种类型——h(n)奇对称，N为偶数 抽样值 H(k)的约束为： * 3 频率抽样法设计线性相位FIR DF的步骤 (1) 根据所需设计的数字滤波器类型，确定线性相位FIR滤波器的类型。 (2) 获得Hd(ej?)在??[0,2?)区间的N个抽样点上的值Hd(k)，并使Hd(k) 满足线性相位条件，即 (3) 利用IDFT得到h(n) * 例：用频率抽样法设计一个逼近截频?c=0.5p的理想低通滤波器。 解： (1) 确定线性相位FIR滤波器的类型，I, II, III, IV型均可采用。 选用I型线性相位系统 (2) 获得Hd(ej?)在W?[0,2?)区间的N个抽样点上的值Hd(k) 取N=11 * * 解： (2) 获得Hd(ej?)在??[0,2?)区间的11个抽样点上的值Hd(k) (3) 利用IDFT得到h(n) 10 - 11 9 - 11 周期为11 例：用频率抽样法设计一个逼近截频?c=0.5p的理想低通滤波器。 * * 例：用频率抽样法设计一个逼近截频?c=0.5p的理想低通滤波器。 %频率抽样法设计I型线性相位低通FIR N=11;Wp=0.5*pi; m=0:floor(N/2); Wm=2*pi*m./N; Ad=[Wm=Wp]; Hd=Ad.*exp(-j*0.5*(N-1)*Wm); Hd=[Hd conj(fliplr( Hd(2:ceil(N/2)) ) )]; h=real(ifft(Hd)); w=linspace(0,2*pi,1000); H=freqz(h,[1],w); %plot(w/pi,20*log10(abs(H)));grid; plot(w/pi,abs(H));grid; * * 解： 设计结果 ?c=0.5p 0 0.4p 0.8p 1.2p 1.6p 2p 0 0.2 0.4 0.6 0.8 1 N=11 N=31 Hd(k)及频率抽样法设计的低通滤波器幅度函数 例：用频率抽样法设计一个逼近截频?c=0.5p的理想低通滤波器。 * * 解： 设计结果 ?c=0.5p 0 0.4p 0.8p 1.2p 1.6p 2p -25 -20 -15 -10 -5 0 N=11 N=31 频率抽样法设计的低通滤波器的增益响应 例：用频率抽样法设计一个逼近截频?c=0.5p的理想低通滤波器。 * * 4 过渡带抽样的优化设计 实际 H (e j?) 与目标 H d (e j?) 相比，抽样点处完全一致，非抽样点处 H (e j?) 值由内插值得到，这就与实际 H d (e j?) 可能存在误差，称为“逼近误差”。 逼近误差的大小取决于 H d (e j?) 曲线变化的快慢，H d (e j?) 变化越缓慢，内插值越接近于理想值，逼近误差就越小。反之，逼近误差就越大。同时产生尖峰和波纹。 * 4 过渡带抽样的优化设计 为此，可稍加一些过渡抽样点，这些点的最佳值由计算机算出，这样可以减小频带边缘突变，增大阻带最小衰减！ * 4 过渡带抽样的优化设计 * 例：利用频率抽样法设计一个II型 (N=64 )线性相位FIR低通滤波器，?p=0.5p，?s=0.6p。 解:没有加过渡带 * * 例：利用频率抽样法设计一个II型 (N=64 )线性相位FIR低通滤波器，?p=0.5p，?s=0.6p。 直接设计出的滤波器在阻带的衰减不会超过20dB * * 例：利用频率抽样法设计一个II型 (N=64 )线性相位FIR低通滤波器，?p=0.5p，?s=0.6p。 利用1个过渡点设计的FIR低通滤波器的幅度函数 一个过渡点 幅值0.38 * * 利用1个过渡点设计的FIR低通滤波器的增益响应 阻带衰减 43dB 上例的设计结果 ? p=0.5p， ?