fir数字滤波器设计与实现.doc

FIR数字滤波器设计与实现 ———————————————————————————————— 作者： ———————————————————————————————— 日期： 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个加法器 2L=M个延迟器，2L+1=M+1个乘法器，2L 特点：可以分别控制每个子系统的零点 分析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点的有