《毕业设计：FIR数字滤波器设计与实现》.doc

FIR数字滤波器 一.摘 要：FIR滤波器是一种常用的滤波器，它在保证幅度特性满足技术要求的同时，很容易做到严格的线性相位特性，在语音分析、图像处理、雷达监测等对信号相位要求高的领域有着广泛的应用，能实现IIR滤波器不能实现的许多功能。 二.关键词：FIRFIR数字滤波器的设计之前，有必要介绍和总结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滤波器的优点是可方便地实现线性相位。 2H(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代替，使极点和相应的零点移到单位圆内。 （）当N为偶数时的频率取样型结构如图所示。???????????? （）当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滤波器的快速卷积结构： ?这里DFT和IDFT都将采用快速傅里叶变换算法，当N1和N2足够长时，比直接计算线性卷积要快得多。 2.线性相位FIR滤波器的特点????线性FIR滤波器FIR滤波器线性相位FIR滤波器线性相位FIR滤波器的横截型结构??? FIR滤波器的单位冲激响应h(n)是有限长的（0）,其Z变换为：??? ? 其傅立叶变换为： ????????????? 其中H()是幅度函数，是一个纯实数，可正可负， θ(ω)是相位函数。可以证明，线性相位FIR滤波器的冲激响应满足对称条件： ???????????? h(n) =±h(N-1-n)?和 （1）、线性相位FIR滤波器的幅度函数和相位函数????（）当h(n)是偶对称时，其幅度函数和相位函数分