有限长单位脉冲响应（FIR）滤波器的设计方法.ppt

有限长单位脉冲响应（FIR）滤波器的设计方法   §4.1 线性相位FIR数字滤波器的特性    3. h(n)奇对称，N为奇数，h(n)=-h(N-1-n)    四种线性相位FIR滤波器 小结：    §4.2 窗口设计法（时域） 一.矩形窗口法       例题：用矩形窗设计一个线性相位低通滤波器 解： （1） 因此 （2）因为 偶对称，所以若N为奇数，则属于第一种线性相位滤波器，若N为偶数，则属于第二种线性相位滤波器。 （3）若为汉宁窗，则 §4.3 频率采样法 二.设计方法 三、?约束条件  因此，Hk 也必须满足奇对称性： 相位关系同上， 其它两种线性相位FIR数字滤波器的设计，同样也要满足幅度与相位的约束条件。? 四、逼近误差        §4.5 IIR与FIR数字滤器的比较 补充例题： 4.8 用频率采样法设计一个线性相位低通滤波器， ，边沿上设一点过渡带 。试求各点采样值 。 解： 由于N为偶数，故为第二种类型滤波器。 （注：第四种类型滤波器不能实现低通！） 有： 频率间隔为： 又因为 ，则 故低通通带内包含的点数为 ，所以取 为过渡点。因为必须满足 单位圆上的频响为： 这是一个内插公式。 式中 为内插函数 令 则 内插公式表明： 在每个采样点上， 逼近误差为零，频响 严格地与理想频响的采样值 H(k)相等； 在采样点之间，频响由各采样点的内插函数延伸迭加而形成，因而有一定的逼近误差，误差大小与理想频率响应的曲线形状有关，理想特性平滑，则误差小；反之，误差大。在理想频率响应的不连续点附近， 会产生肩峰和波纹。 N增大，则采样点变密，逼近误差减小。 图 频率采样的响应 例：设计一个FIR数字 LP 滤波器，其理想特性为 采样点数 N=33，要求线性相位。 解：根据P.142的表4.1，能设计低通线性相位数字滤波器的只有1、2两种，因N为奇数，所以只能选择第一种。 即 h(n)=h(N-1-n)， 幅频特性关于π偶对称，也即 HK 偶对称。 利用 HK 的对称性,求π～2π区间的频响采样值。 根据指标要求，在0～2π内有33个取样点，所以第k点对 应频率为 而截止频率 0.5π位于 之间 ，所以，k=0～8时，取样值为1；根据对称性， 故 k=25～32时，取样值也为1，因 k=33 为下一周期，所以0～π区间有9个值为 1的采样点，π～2π区间有8个值为 1 的采样点，因此： 将 代入内插公式，求H(ejω)： 考虑到8k25时 Hk=0，而其它k时，Hk=1,令 k=33-n，则 从图上可以看出，其过渡带宽为一个频率采样间隔 2π/33，而最小阻带衰减略小于20dB。 对大多数应用场合，阻带衰减如此小的滤波器是不能令人满意的。 增大阻带衰减三种方法： 1）加宽过渡带宽，以牺牲过渡带换取阻带衰减的增加。 例如在本例中可在k=9和k=24处各增加一个过渡带采样点H9=H24=0.5，使过渡带宽增加到二个频率采样间隔4π/33，重新计算的H(ejω)见图4.12(c)，其阻带衰减增加到约 -40dB。 2）过渡带的优化设计 根据H(ejω)的表达式，H(ejω)是Hk的线性函数，因此还可以利用线性最优化的方法确定过渡带采样点的值，得到要求的滤波器的最佳逼近（而不是盲目地设定一个过渡带值）。 例如，本例中可以用简单的梯度搜索法来选择H9、H24,使通带或阻带内的最大绝对误差最小化。