FIR數字滤波器的设计与实现.doc

FIR数字滤波器设计与实现 一.摘 要：数字滤波器是一种具有频率选择性的离散线性系统，在信号数字处理中有着广泛的应用。其中FIR滤波器是一种常用的滤波器，它在保证幅度特性满足技术要求的同时，很容易做到严格的线性相位特性，在语音分析、图像处理、雷达监测等对信号相位要求高的领域有着广泛的应用，能实现IIR滤波器不能实现的许多功能。 二.关键词：FIR 窗函数 系统函数 MATLAB 三.内容提要： 数字滤波器的功能就是把输入序列通过一定的运算变换成输出序列，因此数字滤波器的结构系统中就必须包括一定数量和性能的运算器件和运算单元，而运算器件和运算单元的配置必须由数字滤波器的结构特点和性能特点来决定，因此在进行FIR数字滤波器的设计之前，有必要介绍和总结FIR数字滤波器的基本结构和相关特性（包括频响曲线（幅度和相位），单位冲激响应等），在介绍完其基本结构和相关特性后，就进行FIR数字滤波器的设计和实现。 （一） FIR滤波器的基本结构 在讨论任何一种滤波器时，都要着重分析其系统函数，FIR滤波器的系统函数为： 。从该系统函数可看出，FIR滤波器有以下特点： 1)系统的单位冲激响应h(n)在有限个n值处不为零； 2)系统函数H(z)在|z|0处收敛，极点全部在z=0处(稳定系统)； 3)结构上主要是非递归结构，没有输出到输入的反馈，但有些结构中(例如频率抽样结构)也包含有反馈的递归部分。 1.FIR滤波器实现的基本结构有： 横截型（卷积型、直接型） a.一般FIR滤波器的横截型(直接型、卷积型)结构： 若给定差分方程为：? 。 则可以直接由差分方程得出FIR滤波器结构如下图所示： 这就是FIR滤波器的横截型结构，又称直接型或卷积型结构。 b.线性相位FIR滤波器的横截型结构 若h(n)呈现对称特性，即此FIR滤波器具有线性相位，则可以简化成横截型结构，下面分情况讨论： ①N为奇数时线性相位FIR滤波器实现结构如图所示： ?? ②N为偶数时线性相位FIR滤波器实现结构如图所示 我们知道IIR滤波器的优点是可利用模拟滤波器设计的结果，缺点是相位是非线性的，若需要线性相位，则要用全通网络进行校正，比较麻烦，而FIR滤波器的优点是可以方便地实现线性相位。 2)、级联型 将H(z)分解为若干个实系数一阶或二阶因子相乘： 实现结构如下图所示： 该结构图中有2L=M个延迟器，2L+1=M+1个乘法器，2L=M个加法器。 分析H(z)及结构图可以得出级联型的特点： ①每个基本节控制一对零点，便于控制滤波器的传输零点。 ②系数比直接型多，所需的乘法运算多。 3）频率取样型 ? 若FIR滤波器的冲激响应为有限长(N点)序列h(n)，则有如图所示的关系： ??? 因此，对h(n)可以利用DFT得到H(k)，然后利用内插公式： 来表示系统函数，这就为FIR滤波器提供了另外一种结构：频率抽样结构， ? 这种结构由两部分级联而成：分析系统函数 ? 其中级联的第一部分为：?????????????? ?? 这是一个梳状滤波器，它滤掉了频率 及其各次谐波。 ? 级联的第二部分为N个一阶网络并联而成，第k个一阶网络为：?? ??????????????????????? ????????它在单位圆上有一个极点： ??? 这是一个谐振频率 的无损耗谐振器。这个谐振器的极点正好与梳状滤波器的一个零点(i=k)相抵消，从而使这个频率上的频率响应等于H(k)。这样，N个谐振器的N个极点就和梳状滤波器的N个零点相抵消，从而在N个频率抽样点上的频率响应就分别等于N个H(k)值。 有上叙的理论分析基础可以得到FIR滤波器的频率抽样结构。 FIR滤波器的频率抽样结构如图所示： 频率抽样结构的特点是它的系数H(k)就是滤波器在 处的响应，因此控制滤波器的频率响应很方便。 频率抽样结构存在问题的问题是：在有限长情况下，系数量化后极点不能和零点抵消，使FIR系统不稳定。 解决方法： 在r圆上进行(r1但近似等于1)取样，即用r代替，使极点和相应的零点移到单位圆内。 （a）当N为偶数时的频率取样型结构如图所示。 ???????????? （b）当N为奇数时频率抽样型结构如图所示。 ??????? 4）快速卷积结构 若FIR滤波器的单位冲激响应h(n)是一个N1点有限长序列，输入x(n)是一个N2点有限长序列，那么输出y(n)是x(n)与h(n)的线性卷积，它是一个L＝N1+N2-1点的有限长序列。 我们知道，将x(n)补上L－N2个零值点，将h(n)补上L－N1个零值点，然后进行L点圆周卷积，就可以代替原x(n)与h(n)的线性卷积。 而圆周卷积可以用DFT和IDFT的方法来计算，这样我们得到FIR滤波器的快速卷积结构：