数字信号处理_Lecture 22.ppt

数字信号处理第二十二讲 中国地质大学(北京) 地球物理与信息技术学院 电子信息工程教研室 制作 第九章 FIR数字滤波器的设计Chapter 9 Designs of FIR Digital Filter §9.4 频率采样设计法9.4 Designing of Sampling Frequency 从上图可看出：最小阻带衰减大约为?30dB，还未达到?50dB。 在过渡带再增加一点并经反复试验后选取N=36，用优化法 求出k =6、7、31、32，|H(6)|=|H(32)|=0.5925，|H(7)|=|H(31)|=0.1099； 用同样的方式得出图 9-18的结果，阻带最小衰减达到-62dB。 满足设计要求。 *电子信息工程教研室 * §9.1 数字滤波器设计的基本概念 §9.2 用窗函数法设计线性相位FIR滤波器原理 §9.3 各种常见窗函数 §9.4 频率采样设计法 本讲的主要内容 频率采样设计法 频率采样设计法设计线性相位FIR滤波器例题 直接设计法 理想频率响应Hd(e j?)，在0到2?之间进行等间隔采样，即 滤波器频率特性的采样值H(k) (9-32) 得到一个N?1阶FIR滤波器，即 频域采样理论中X(z)的内插公式 (9-33) ?(?)是内插函数 (9-34) 将式(9-33)代入式(9-34)，即 (9-35) 得到FIR滤波器频率响应H(e j?) 从内插公式(9-33)可以看到： (2) 在各采样点之间的频率响应误差大小取决于理想频率响应曲线形状，理想频率响应特性变化越平缓，则内插值越接近理想值，逼近误差越小，如图 9-15(b)梯形理想频率特性所示。在理想频率特性的不连续点附近，就会产生肩峰和波纹，如图 9-15 (a)矩形理想频率特性所示。 (3) 在靠近通带边缘的逼近误差较大，而在通带内的逼近误差较小。 (1) 图 9-15 频率响应的采样 线性相位的FIR数字滤波器，则其采样值H(k)的幅度和相位一定要满足线性相位约束条件，即h(n)为奇对称或偶对称。 则根据式(9-36)可知，对于第I型及第II型来说，θk必须满足 (9-38) (9-37) (9-36) 即 对于第III型及第IV型来说，θk必须满足 即 (9-39) 对于第I型与第IV型Hk必须满足： Hk=HN?k (9-40) 对于第II型与第III型Hk必须满足： Hk＝?HN?k (9-41) 最优设计法 [例 8?3] 利用频率采样法来设计一个低通FIR数字滤波器，通带截止频率为?p=0.2?，阻带截止频率为?s=0.4?，阻带衰减相当于As=50dB，滤波器具有线性相位。 解：滤波器截止频率为?c＝0.3?，其理想频率特性是 频率采样后的H(k)序列，如图 9-16(a)所示。由于|H(k)|对称于?＝?，可将 [N/2]+1?k?N?1的图形略去。截止频率?c＝0.2?，采样间隔??=2?/N，取N=10，??=0.2?，使?p和?s正好落在采样点上。通带0? ? ?0.2?有2个样本，即k=0, 1；过渡带0.2? ? 0.4?没有样本；阻带0.4? ?? ??有4个样本，即 ? =(N?1)/2=4.5，相位相应为 将这些值代入式(9-37) 计算这些H(k)值及20log10|H(k)|，对H(k)进行傅氏反变换后得到h(n)，如图 9-16(b)、(c)、(d)所示。 图 9-16 N=10 ： (a) 理想频率响应 |Hd(e j?)|及其采样值|H(k)|； (b) 设计出的振幅频率特性 |H(e j?)|(dB)；(c) 设计出的冲激响应h(n)；(d) 设计出的幅度特性H(?) 由上可知：最小阻带衰减大约为?14dB，远未达到?50dB。下面用最优设计法在过渡带增加采样点，并增加N。当N=20时，用优化法求出k=3，|H(3)|=0.5，k=17，|H(17)|=0.5。其它H(k)为 ? =(N?1)/2=9.5，相位相应为 如图 9-17所示。 图 9-17 N=20 及在过渡带加一优化点： (a) 理想频率响应|Hd(e j?)|及其采样值|H(k)|；(b) 设计出的振幅频率特性|H(e j?)|(dB)；(c) 设计出的冲激响应h(n)；(d) 设计出的幅度特性H(?) 图 9-18 N=36 及在过渡带加一优化点： (a) 理想频率响应|Hd(e j?)|及其采样值|H(k)|；(b) 设计出的振幅频率特性|H(e j?)|(dB)；(c) 设