《数字信号处理教程——MATLAB释义与实现》第七章.ppt

第七章  FIR滤波器设计  FIR滤波器设计 本章的主要内容为： 7.1 数字滤波器设计的基本要求 7.2 线性相位滤波器的特性 7.3 设计FIR滤波器的窗函数法 7.4 设计FIR滤波器频率样本法 7.5 FIR滤波器的最优设计法 7.6 FIR滤波器设计的一些深入问题 滤波器设计概述 从本章起将着重研究信号分析的逆问题，那就是给定了输入输出指标，设计系统。这种可以人为调整参数来满足给定要求的硬件或软件系统，称为滤波器。这种问题没有唯一的解，因为可能有多种类型的滤波器都能满足给定的要求。要人为地进行选择，就要考虑更多的因素,所以设计比之分析更加接近与工程实际。 有两种类型的数字滤波系统。本章将讨论FIR数字滤波器的几种基本设计算法；下一章讨论IIR数字滤波器的设计问题。 滤波器设计概述 本书中的设计大多数是针对选频型的滤波器。从幅频特性的要求而言，有低通、高通、带通和带阻滤波器；从相频特性来说，往往要求线性相位。在FIR滤波器设计中，还会涉及到像微分器或希尔伯特变换器之类的系统，尽管它们不是选频滤波器，但是也遵循选频滤波器的设计技术。 在计算技术高度发展的今天，几乎所有的电子设计软件中都把滤波器设计作为其模块之一。那种用查表、手工计算的方法在工程实践中已经濒于淘汰。 滤波器设计概述 数字滤波器的设计要经过三个步骤： 确定指标：在设计一个滤波器前，必须有一些指标。这些指标要根据应用确定。 模型逼近：一旦确定了技术指标，就可利用前面学习过的基本原理和关系式，提出一个滤波器模型来逼近给定的指标体系。这是滤波器设计所要研究的主要问题。 实现：上面的两步的结果得到的滤波器，通常是以差分方程、系统函数或脉冲响应来描述的。根据这个描述，可以如第六章所讨论的那样，用硬件或计算机软件来实现它。 滤波器设计概述 数字滤波器的主要指标： 在很多实际应用，象语音或音频信号处理中，数字滤波器用来实现选频功能。因此，指标的形式通常为频域中的幅度和相位响应。 幅度指标可以两种方式给出。第一种，叫做绝对指标，它提出了对幅度响应函数|H(jω)|的要求。这些指标一般可直接用于FIR滤波器。第二种方法叫做相对指标，它以分贝(dB)值的形式提出对幅频特性的要求。其值定义为： FIR滤波器设计概述 相对幅度指标： 由于定义中所包含的归一化，滤波器的相对幅频特性最高处的值为0dB，又由于定义中的负号，幅频特性小的地方，其dB值反而是正的。为了说明这些指标，以低通滤波器为例进行讨论。首先看下面两种指标图，它们分别用绝对指标和相对指标表示。 FIR滤波器设计概述 FIR滤波器设计概述 绝对指标：图a是典型的低通幅特性绝对指标： [0,ωp]段叫做通带，δ1是在理想通带中能接受的振幅波动（或容限）,0δ11。 [ωs,π]段叫做阻带，δ2是阻带中能接受的振幅波动（或容限）， 0δ21。 [ωp,ωs]叫做过渡带，在此段上对幅度响应通常没有限制，也可以给些弱限制。 相对指标（dB）：低通的典型相对技术指标如图b所示，三种频带的定义不变，只是δ1 ，δ2 换成了Rp，As。其中 Rp是通带波动的dB数; As是阻带衰减的dB数； FIR滤波器设计概述 绝对指标和相对指标间的变换关系： 可以用MATLAB语句fgp711准确地求出绝对指标和相对指标的对应数值。并画出如图中所示曲线。 FIR滤波器设计概述 FIR滤波器设计概述 相位指标主要是线性相位条件。 理想的滤波器相位θ应该和频率ω成正比，也就是θ(ω)= -τω。因此满足此式的系统对信号中所有的频率分量都具有相同的时间延迟τ (无量纲)，因此无波形失真。 略低一点的要求是相位θ和频率ω成线性关系，即θ(ω)= θ0-τω。 IIR滤波器无法做到线性相位，通常用群迟延g(ω)来评价。其定义为相位θ对频率ω的导数 ，它愈接近常数愈好。 FIR滤波器设计概述 上述指标是针对低通滤波器的，其他类型的选频滤波器，像高通和带通滤波器，它们的技术指标可作类似的规定，并且可以由低通滤波器演化而得。 本章讨论FIR数字滤波器的设计和逼近。FIR滤波器在设计和实现上具有如下的优越性： (a). 相位响应可为严格的线性； (b). 不存在稳定性问题； (c). 只包含实数算法，不需要递推运算。 7.2 线性相位滤波器的特性 线性相位条件要求滤波器分子系数bn满足对称性条件。设滤波器的系数长度为N，则这些系数应关于中心点τ=n0=(N-1)/2对称。偶对称时，h(n)=h(N-n-1)；而奇对称时，h(n)=-h(N-n-1)。再考虑到N可以为奇数或偶数，总共有四种类型的线性相位FIR滤波器。 在讨论线性相位滤波器频率响应时需要引进幅特性正负号的概念。以