1、1第一章 数字滤波器概述数字滤波器是指用有限精度算法实现的，完成信号滤波处理的离散时间线性非时变系统，其输入是一组由模拟信号取样和量化得到的数字量，其输出是经过变换的另一组数子量。数字滤波器具有高稳定性、高精度、灵活性大等突出特点。这章主要讲述了数字滤波器的基本概念、及数学建模，以便在后续章节中实现滤波算法做好准备 1。1.1 数字滤波器的基本概念滤波器是一种能使有用频率信号通过而同时抑制无用频率信号的电子装置，或者说是可以人为调整参数来满足给定要求的硬件或软件系统 1。 滤波器也就是一个线性移不变离散系统，其基本输入输出关系为：(1-1)h(n)*xmh()xy(n)Nm10其中 h(n)为
2、系统的单位冲激响应。在 Z 域上，有：(1-2)-n-zazabbH(z)10分析式（1.2）可知，只要分母多项式 ai（i=1,2,n）中有一个不为零，滤波器的的内部就存在反馈环。这种内部有反馈环的滤波器称为递归型滤波器。相反，所有系数 ai（i=1,2,n）都为零的滤波器称为非递归型滤波器由于非递归型滤波器的内部没有反馈环，所以这种滤波器总是稳定的。从滤波器的单位冲激响应来看，滤波器可以分为 IIR 滤波器和 FIR 滤波器。由于 IIR 滤波器的传递函数存在 0 之外的极点，所以单位冲激响应是无限持续的。因而 IIR 滤波器与递归型滤波器是一致的。而 FIR 滤波器的传递函数不存在 0
3、之外的极点，所以，本质上 FIR 滤波器与非递归型滤波器是一致的。但是，如果稳定的递归型滤波器与非递归型滤波器级联后，其极点与零点相互抵消，使滤波器在 0 之外不存在极点，这种滤波器也是 FIR 滤波器 2。这种情况下，由于滤波器内部存在反馈环，因而这种滤波器就成为递归型滤波器。频率采样2滤波器就是这种结构的滤波器。具有频率选择功能的数字滤波器按功能分类有低通滤波器（LPF） 、高通滤波器（HPF） 、带通滤波器（ BPF） 、带阻滤波器 3。另外，在全频带幅频特性固定的滤波器称为全通滤波器，可用于改善滤波器的相频特性。1.2 数字滤波器的频率特性1.2.1 幅频特性设采样周期为 T，令 z=
4、jwT，从 z 变换理论即可得到 H（z ）的频率特性，可表示为：(1-3)j(jj e|(|)(e|)(|)(e wTwTjwTjwTH其中，| H（e jwT）|表示数字滤波器的幅频特性；H（e jwT） 、(wT)表示相频特性 3。另外，H（e jwT）是数字滤波器的频率传递函数，它表示决定数字滤波器的幅频特性与相频特性的频率特性。各种数字滤波器的理想幅频特性示于图 1-1。1.2.2 相频特性设计数字滤波器时，不仅要考虑它的幅频特性，相频特性也是需要考虑的重要问题。特别是在线性相位滤波器中，相位的延迟不会造成波形失真，所以在波形要求很高时，滤波器具有线性相位是非常重要的 4。我们知道，
5、一般情况下要减小（增大）通带与阻带的波纹，会造成过渡带的截止特性变差（变好） 。根据滤波器的使用目的，其重点是不一样的，在滤波器的阶次确定后，合理的折中方案是必要的。但是当幅频特性要求很高，不能得到满意结果时，不仅要采用更高阶次 FIR 滤波器，还必须考虑采用 IIR 滤波器。1.3 数字滤波器的数学建模任何线性时不变集总参数离散系统都可用下列数学形式之一来表述其输入输出关系，在时域中，滤波器的功能可以用差分方程(1.4)或卷积公式(1.5)来描述：图 1-1 理想幅频特性 H（e jwT）3(1-4)()()( knyainxby(1-5)*hmh在变换域中可用系统传递函数（1.6）来表示：
6、(1-6)nzazbbZXHZY 10)()(这三种方法是等价的，从任何一个都能推导出其他两个 5。即使用同一类表示方法。也存在着许多等价的算法结构。不同的算法会影响系统的某些实际性能。滤波器的工程实现要用计算机的硬件或软件帮助完成，这需要考虑许多的问题：（1）计算的效率，即完成整个滤波所需要的乘法和加法次数；（2）需要的存储量；（3）滤波器系数的量化影响；（4）运算中的舍入和截断误差、饱和和溢出。不同的滤波器结构可以实现同样的传递函数，但不同的算法在满足上述的要求方面是有差别的，有的差别还很大，改变计算结构对提高计算效率和节省存储量能起很大的作用。1.4 本章小结本章主要讨论了数字滤波器的一
7、些基本概念。数字滤波器按单位冲激响应可以分为 FIR 和 IIR 滤波器，按频率选择功能的数字滤波器按功能分类有低通滤波器（LPF） 、高通滤波器（ HPF） 、带通滤波器（ BPF） 、带阻滤波器。数字滤波的作用是滤除信号中某一部分频率分量。信号经过滤波处理，就相当于信号频谱与滤波器的频率相应相乘。从时域来看，就好似输入信号与滤波器的冲激相应作卷积和。数字滤波器具有幅度特性和相位特性。4第二章 数字滤波器的设计准则及基本步骤上一章主要讨论的是数字滤波器的数学表示及一些特性和分类的问题，并没有解决如何去设计一个数字滤波器。众所周知，计算机或 DSP 芯片处理的是数字信号，那么对于一个数字滤波器
8、，要知道它设计准则。在这章主要讨论数字滤波器的设计准则和基本步骤。2.1 数字滤波器的设计准则滤波器设计误差是滤波器的期望响应与实际响应之间的差别。在 FIR 滤波器设计中通常采用 3 种准则 6。（1）最小平方（LS）近似，它是期望滤波器与设计滤波器频率响应之间平方误差的平均值。（2）切比雪夫（Chebjyshe）近似，它是频率响应指定区域内的最大误差。（3）巴特沃思（Butterworth）近似，它是基于泰勒级数近似的期望响应近似。大多数的有源滤波器设计算法都是以这 3 种近似法的一种和 Parks McCllan(最优 )算法。频率采用快速，简单，但其对整体频率响应的控制最弱。傅立叶设计
9、的频率响应中存在不合乎要求的震荡或过冲，可以使用窗口函数控制这些效应。加窗方法的重要优点在于其简单性，其主要缺点是缺少对通带和阻带纹波控制的灵活性。Parks-McCllan 算法生成一个线性相位 FIR 滤波器，并使切比雪夫误差最小，但是如果不利用计算软件，其设计非常慢 6。2.2 数字滤波器设计的基本步骤2.2.1 指标在设计一个滤波器之前，必须首先根据工程实际的需要确定滤波器的技术指标。在很多实际应用中，数字滤波器常常被用来实现选频操作。因此，指标5的形式一般在频域中给出幅度和相位响应。幅度指标主要以两种方式给出。第一种是绝对指标。它提供对幅度响应函数的要求，一般应用于 FIR 滤波器的
10、设计。第二种指标是相对指标。它以分贝值的形式给出要求。在工程实际中，这种指标最受欢迎。对于相位响应指标形式，通常希望系统在通频带中有线性相位。运用线性相位响应指标进行滤波器设计具有如下优点 7：（1）只包含实数算法，不涉及复数运算；（2）不存在延迟失真，只有固定数量的延迟；（3）长度为 N 的滤波器（阶数为 N-1），计算量为 N/2 数量级。因此，本文中滤波器的设计就以线性相位 FIR 滤波器的设计为例。 2.2.2 逼近确定了技术指标后，就可以建立一个目标的数字滤波器模型。通常采用理想的数字滤波器模型。之后，利用数字滤波器的设计方法，设计出一个实际滤波器模型来逼近给定的目标。2.2.3 性
11、能分析和计算机仿真 上两步的结果是得到以差分方程或系统函数或冲激响应描述的滤波器。根据这个描述就可以分析其频率特性和相位特性，以验证设计结果是否满足指标要求；或者利用计算机仿真实现设计的滤波器，再分析滤波结果来判断。一些 FIR 滤波器的设计方法可能比较简单，如傅立叶技术。但是，不存在对于所有滤波器特性都是最优的单一方法。滤波器系数可以通过一些合理的公式来生成。但是，采用最优方法来获得最优滤波器需要大量的计算。此外，在定点 DSP 处理器实现实时滤波器的主要问题就是有限精度和动态范围。分析有限长效应的理论过于复杂，不适与手工计算，所以通常依靠计算机软件来设计数字滤波器 7。2.3 本章小结本章
12、主要讨论了数字滤波器的设计准则，和设计数字滤波器的基本步骤。设计准则有最小平方（LS）近似、切比雪夫（Chebjyshe）近似、巴特沃思（Butterworth ）近似，这三种方法各有优缺点，在工程实现时由于现在有了高效的滤波器设计软件包（如 MATLAB）可供利用，所以 Parks-McCllan 算法被6广泛应用于大多数应用中。设计一个数字滤波器时有一些基本的步骤，首先要按照实际需要，确定滤波器的性能指标，寻找一个满足预定性能要求的离散线性非时变系统，然后通过模拟，验证所设计的系统是否符合给定的性能指标。第三章 滤波器的选择上一章我们解决了两个问题，那就是数字滤波器的设计准则及基本步骤。滤
13、波器包括模拟滤波器和数字滤波器。但在许多信号处理应用中，使用数字滤波器代替模拟滤波器是有利的。数字滤波器不仅具有高可靠性、高精度，而且它是用可编程的处理器实现的，因而容易改变滤波器的参数。在应用环境改变时，这个特性方便设计工程师升级或修改已有的滤波器特性。并且数字滤波器容易实现不同的幅度和相位频率特性指标，克服与模拟滤波器器件相关的电压漂移、温度漂移和噪声问题，因此这里选择使用数字滤波器。而数字滤波器的设计又分为 FIR（有限长单位冲激响应）滤波器和 IIR（无限长单位冲激响应）滤波器。在这一章中，就专门讨论滤波器的选择问题及 FIR 滤波器的定义和FIR 的几种基本结构。3.1 FIR 与
14、IIR 滤波器的比较从性能上说，IIR 滤波器可用较低的阶数获得高的选择性，所用的存储单元少，计算量小，效率高。但是这个高效率是以相位的非线性为代价的。选择性越好，则相位非线性越严重。FIR 滤波器要达到高的选择性，必须用较高的阶数。对于同样的滤波器设计指标，FIR 滤波器所要求的阶数可能比 IIR 滤波器高 510 倍，成本较高，信号延时也较大。如果按线性相位要求来说，则 IIR滤波器就必须加全通网络进行相位校正，同样要大大增加滤波器的阶数和复杂性。而 FIR 滤波器却可以得到严格的线性相位 1。从结构上看，IIR 滤波器必须采用递归结构来配置极点，并保证极点位置在单位圆内。由于有限字长效应
15、，运算过程中将对系数进行舍入处理，引起极点7的偏移。这种情况有时会造成稳定性问题，甚至产生寄生振荡。相反，FIR 滤波器只采用非递归结构，不论是在理论上还是在实际的有限精度运算中都不存在稳定性问题，因此造成的频率特性误差也较小。此外 FIR 滤波器可以采用快速傅立叶变换算法，在相同阶数的条件下，运算速度可以快的多 1。另外，也应看到，IIR 滤波器虽然设计简单，但主要是用于设计具有分段常数特性的滤波器，如低通、高通、带通及带阻等，往往脱离不了模拟滤波器的格局。而 FIR 滤波器则要灵活的多，尤其是他易于适应某些特殊应用，如构成数字微分器或西尔波特变换器等，因而有更大的适应性和广阔的应用领域。3
16、.2 FIR 滤波器的定义根据上面的比较，在这次设计中我们选择的是 FIR 数字滤波器，所以我们只简单介绍 FIR 滤波器的结构。 FIR 网络结构特点是没有反馈支路，即没有环路，其单位脉冲响应是有限长的 8。设单位脉冲响应 h(n)长度为 N，其系统函数 H（z）和差分方程分别为(3-1)10110)( NnzbzbzbZ(3-2)10110 )()()()( innxxxny其单位脉冲响应 h(n)是有限长的，按照式（3-1）的 Z 反变换，h( n)表示为滤波器的阶数 N-1。它的长度（等于系数的个数）为 N。FIR 滤波器总是稳定的，同 IIR 结构相比而言，相对简单一些。而且，FIR
17、 滤波器可设计成具有线性相位，这是某些应用所希望的。3.3 FIR 滤波器的结构3.3.1 直接型设 N=5（即四阶 FIR 滤波器） ，按照 H（Z）或者差分方程式画出的直接结(3-3)bn 0Nn-10 其他h(n)=8构图如图 3-1 所示，它可以用抽头延迟线实现 3。则可以写出：(3-4)4()3()2(1)-()(0 n-xb-n-xbnxby 根据此信号流图可以直接写出方程式（3-4） ，而且不必写出联立方程了。3.3.2 级联型将 H（Z）进行因式分解，并将共轭成对的零点放在一起，形成一个系数对称为实数的二阶形式。这样级联型网络结构是由一阶或二阶实系数因式级联构成的，其中每一个因
18、式都可以用直接型实现。设 FIR 网络系统函数 H（Z）如式（3-5） ，要求画出 H（Z）的直接型和级联型结构。(3-5)3431212)( aZa将 H（Z ）进行因式分解，得到：(3-6)() 211-edcbZ其级联型结构和直接型结构如图 3-2 所示。级联型结构每个一阶因式控制一个零点，每个二阶因式控制一对共轭零点，因此调整零点位置比直接型方便，但 H（Z）中的系数比直接型多，因而需要的乘法器多。在上例中直接型需要四个乘法器，而级联型则需要五个乘法器。分解的因式越多，需要的乘法器也越多。当 H（Z）的阶次高于三阶时，不易分解，需要用 Matlab 协助。3.3.3 线性相位形式ax(
19、n)z-1b z-1z-1decy(n) x(n)y(n)z-1 z-1 z-1a1 a2 a3 a4图 3-2 H（Z）的直接型和级联型结构图 3-1 直接型结构x(n)y(n)z-1 z-1 z-1 z-1 z-1b0 b1 b2 b3 b49一般的，理想的低通滤波器具有线性相位，如（3-7）式所示。 （3-eHjw)(7）其中 =0 或/2， 为时延的拍数。要得到线性相位，实际滤波器的脉冲响应序列应当具有对称性。利用这种对称性可以简化滤波器的结构。线性相位结构在本质上仍然是直接形式，只是缩减了乘法计算量。因此，传递函数多项式的组成形式上，线性相位结构仍等于直接形式。3.4 本章小结经过对
20、 FIR 和 IIR 的比较得知：FIR 滤波器可以得到严格的线性相位；FIR滤波器只采用非递归结构，有限精度的计算不会产生震荡；FIR 滤波器可以采用 FFT 算法,在相同阶数下，运算速度可以快的多；IIR 滤波器可以借助于模拟滤波器的成果；而设计 FIR 需要借助计算机辅助设计；FIR 易于适用某些特殊应用。这里，选择的是 FIR 滤波器。FIR 是没有反馈支路，即没有环路，其单位脉冲响应是有限长的。FIR 的结构有直接型、接连型、转置型和线形相位型等结果，在实际中要根据不同的要求来选择不同的结构。10第四章 FIR 滤波器的 Matlab 工具箱实现本章介绍了设计滤波器的 Matlab
21、工具箱 SPTool 和 FDATool，这两种数字信号处理工具箱都是一种图形用户界面，能够方便快捷的让我们设计自己的数字滤波器。这两种工具箱都能查看所设计滤波器的幅度响应、相位响应、群延迟、零点极点图、冲激响应和阶跃响应。4.1 Matlab SPTool 工具箱的 FIR 滤波器设计 4.1.1 SPTool 工具箱简介SPTool 是一个图形环境，提供了一种用于信号查看、滤波器设计和频谱分析的多图形环境。它为信号处理工具箱中的很多函数提供了易于使用的界面，只需要操纵鼠标就可以观察、分析和打印数字信号，分析、实现和设计数字滤波器，以及频谱分析等 6。4.1.2 用 SPTool 工具箱设计
22、 FIR 滤波器启动 Matlab 后在命令窗口输入 SPTool 回车就可以打开 SPTool 的主窗口。如图 41 所示，由 SPTool 的主窗口可以看出，SPTool 有 3 个列表框：Signals列表框，Filters 列表框和 Spectra 列表框，它们对应着 SPTool 工具中的 4 个功能模块 6：(1)信号浏览器；(2)滤波器设计器；(3)滤波器观察器；(4)谱观察器;11在用 SPTool 进行信号处理时，必须先载入数据或先设计好滤波器。按照不同的数据类型，数据载入可以分为信号数据载入、滤波器数据载入以及功率谱数据载入。首先在 Matlab 工作空间创建信号数据。这里
23、以低通为例，创建信号数据时的 Gen_Signal.m 文件如下 9 10：fs=1000; t=0:1/fs:1;x=sin(2*pi*5*t);xn=x+rand(size(t); B,A=butter(20,0.6);plot(t,xn)可以利用File菜单中的Import选项将其他信号、滤波器或频谱由Matlab工作区中导入到SPTool工作区内。如图4-2所示：图 4-1 SPTool 的主窗口12已经从工作区中选择数据 xn 导入到 Signals 列中，所导入的数据命名为sig1。这时单击 OK 按钮，将以此信号 sig1 更新原来的 SPTool 窗口。选择 Signals 列
24、中信号的名称，并单击列表下面的 View 按钮，就可以查看信号了。如图 4-3 给出了受干扰的正弦波 sig1。Signals Browser 窗口允许查看、测试和分析时域信号的特定区域。它包括有水平和垂直标记、带有轨迹和范围的标记，并显示信号的波峰和波谷，用于测试和对比。接下来的任务是设计一个 FIR 低通滤波器去滤除高频噪声分量。将滤波器导入到 Matlab 工作区中最简单的方法是单击 File|Import，并导入分子 B 和分母A 系数，如图 4-4 所示。并单击 View 来查看所导入滤波器的特性。Filter Viewer 窗口允许用户查看所设计滤波器的特性，包括滤波器的幅度响应、
25、相位响应、群延迟、零点极点图、冲激响应和阶跃响应。图 4-2 导入信号 xn图 4-3 信号 sig113SPTool 工具箱中还可以根据自己定义的规格设计 FIR 和 IIR 滤波器，可以通过单击 SPTool 中 Filters 列之下的 New 按钮来完成这一操作。这里选择SPTool/Filters 中的 Edit，对输入信号进行滤波设计，将出现如图 4-5 所示设计的 Filter Designer 窗口，这里也允许自己去选择一个滤波器设计算法、滤波器类型和频率规格。最后通过单击 SPTool 窗口中的 Apply 按钮，可以选择所设计的滤波器或所导入的滤波器 LSLp，并将其应用至
26、输入信号 sig3。如下图 4-6 所示：图 4-4 导入分子分母系数图 4-5 Filter Designer 窗口14图 4-6 信号 sig3将出现一个新的窗口，他允许我们选择滤波器的结构，并指定输出信号的名称为 sig3。由 SPTool 窗口中选择 sig1 和 sig3，并单击 Signal 列下面的 View 按钮，可以查看输入信号(见图 4-11)和滤波后输出信号的时域图。如图 4-7(a) 、 (b)所示(a)滤波前信号 (b)滤波后的信号 图 4-7 从滤波前后波形对比，我们可以看出，高频分量被滤除掉了。4.2 FIR 滤波器的 Matlab FDATool 工具箱的简介M
27、atlab 中包括一组用于设计、分析和实现 FIR 滤波器的函数。利用FDATool 来为实际应用设计、分析和量化 FIR 滤波器，FDATool 中整合了滤波器设计工具箱和信号处理工具箱。这一强大的工具将传统的命令行滤波器设计功能集成在一个图形窗口环境中。4.2.1 滤波器设计与 FDATool 工具箱简介在 Matlab 命令窗口中键入 FDATool 启动 FDATool。FDATool 打开默认设计模式，如图 4-8 所示。15图 4-8 FDATool 窗口FDATool 包含一个交互式 GUI，它允许使用者指定以下滤波器设计参数 6： （1） 滤波器类型可以为低通、高通、带通、带阻
28、、多频器、微分器、Hibert转换器、任意幅度、任意相位、升余弦。（2） 滤波器设计方法可以是用于 FIR 滤波器的等纹波、最小平方和加窗法。（3） 滤波器阶数可由使用者指定，也可以是由滤波器设计方法确定的最小阶数。（4） 窗口技术要求仅可用于采用傅立叶方法的 FIR 设计。其中共有 16 种窗口可供选择。（5） 频率设计要求可以是采样频率(fs) 、通带截止频率或阻带截止频率。 （6） 幅度技术要求可以是通带纹波和阻带纹波，单位可以为线性单位或 dB.4.2.2 有限长单位冲激响应滤波器的指标这里以一个具体例子为例利用 FDATool 设计 FIR 滤波器的过程。滤波器的技术要求给出如下：（
29、1）应使用等纹波设计方法以最低阶数设计一个带通 FIR 滤波器。（2）以 Hz 为单位的频率技术要求为 Fstop1=800Hz, Fstop2=1 500Hz，Fpass1=1 000Hz，Fpass2=1 300Hz,采样频率 Fs=8 000Hz。（3）低阻与高阻带纹波的幅度技术要求分别为 50dB（As1）和 60dB(As2)。16通带纹波限制为 1dB(Ap)。4.3 Matlab FDATool 工具箱的 FIR 滤波器设计有了带通滤波器的技术要求后，就可以进行 Matlab 软件设计了。首先打开MATLAB 进入命令窗口，在命令窗口中敲入“FDATool”命令，然后回车，就进入
30、了 FDATool 的默认设计模式如上图 4-7 所示。在 FDATool 的默认设计模式中有“Response Type”，可以从中选择“Lowpass”， “Bandpass”， “Highpass”等滤波器类型，这里以“Bandpass”为例。把技术指标填入相应的筐中，再点击“Design Filter”,设计的滤波器就如图 4-9 Bandpass 所示。完成滤波器的设计之后，可以分析滤波器的特性。如图 4-8 所示，可以通过选择 Analysis 菜单中所列的方法或单击工具栏按钮来访问分析方法。从Analysis 菜单中选择 Magnitude Response,在右上窗口中显示幅度
31、响应曲线，如图 4-10 所示，还可以选择 Phase Response 如图 4-11 所示：图 4-9 Bandpass 17FDATool 还提供了给出滤波器特性的附加显示。这些曲线包括相位响应、群延迟、冲激响应、阶跃响应、极零图和滤波器系数。FIR 滤波器的冲激响应显示所设计滤波器的系数。极零图显示极点和零点的位置。例如，相位响应和冲激响应分别显示于图 4-11 和图 4-12 中。相位响应显示的滤波器在通带内具有线性相位，这是 FIR 滤波器的一个重要特性。获得对称系数即可获得线性相位滤波器，如冲激响应中所示（图 4-12）已经利用 FDATool 根据给定技术要求设计了一个 FIR
32、 滤波器。滤波器系数以双精度浮点格式表示，需要 64 位字长。需要将这些系数量化为定点数（14位或 32 位） ，以便在定点 DSP 处理器上实现。如图 4-13 所示：在 Analysis 菜单中选择 Filter Coefficents 选项可以查看量化系数。在图 图 4-11 Phase Response Response 图 4-10 Magnitude Response 图 4-12 Impulse Response图 4-13 量化窗口184-14 Filter Coefficients 中可以显示以格式 Q.16 格式表示的量化系数。可以在列表的最左列显示符号+、- 和 0。符号
33、+ 表示量化系数向 +溢出，符号-表示量化系数向-溢出，符号 0 标志着系数已经被量化为零 6。系数量化后，在 Filter Design Add Breakpoint SSBX SXM ; SSBX FRCT STM #d_data_buffer,FIR_DATA_PSTL A,*FIR_DATA_PSTM #(d_data_buffer+K_FIR_BFFR-1),FIR_DATA_P fir_task: RPTZ A,#K_FIR_BFFR-1MACD *FIR_DATA_P-,COFF_FIR_START,ARET.end26连接文件：MEMORYPAGE 0:PROG: o= 100
34、h,l= 2000h PAGE 1:DATA1: o= 2600h, l= 1000hDATA2: o= 2100h, l= 100hDATA3: o= 2200h, l= 100hDATA4: o= 2300h, l= 100hDATA5: o= 2400h, l= 100hDATA6: o= 2500h, l= 100hSECTIONScoff_fir : PROG PAGE 0fir_prog : PROG PAGE 0fir_vars : DATA1 PAGE 1fir_coff : DATA2PAGE 1fir_bfr : DATA3 PAGE 127参考文献1 周霖.DSP 通信工
35、程技术应用. 长沙：国防工业出版社， 2004.1982 陈怀琛.数字信号处理教程.电子工业出版社，2004.783 程佩青.数字信号处理教程. 北京：清华大学出版社，1995.145149，1781824 胡广书.数字信号处理.北京：清华大学出版社，2003.995 郑君里.信号与系统 .高等教育出版社.2003.56576 郭森楙，颜允圣.数字信号处理器体系结构实现与应用.北京：清华大学出版社，2005.78957 Miroslav.信号处理滤波器 .电子工业出版社.2002.146 1498 EmmanueL.数字信号处理实践方法.电子工业出版社.2004.1992039 张威.MATL
36、AB 基础与编程入门 .西安电子科技大学出版社， 2004.15215810 王华.Matlab 在电信工程中的应用.中国水利水电出版社.2005.353711 苏涛，卢光跃等 .DSP 实用技术.西安：西安电子科技大学出版社，2002.15515612 郑红.TMS320C54x DSP 应用系统设计.北京航空航天大学出版社.2004.15515713 ICETEK-VC5416-A 评估板实验指导书14 郑红，吴冠 .TMS320C54X DSP 应用系统设计.北京：北京航空航天大学出版社，2002.13213515 赵红怡.DSP 技术与应用实例.电子工业出版社， 2003.454728
37、致 谢经过长时间的忙碌，本次毕业设计已经接近尾声，通过进行这次设计，我对 FIR 滤波器的设计及实现有了较为清晰的认识。 作为一个本科生的毕业设计，由于经验的匮乏，难免有许多不足的地方，如果没有指导老师的督促指导，以及一起学习的同学们的支持，想要完成这个设计是非常难的。 在这里首先要感谢我的指导老师曹晖老师。曹老师平日里工作繁忙，但在我做毕业设计的每个阶段，都给予了我悉心的指导。在设计过程中经常遇到一些问题，曹老师都给了我详细的讲解，使我顺利的完成了毕业设计。我要感谢所有帮助过我的同学，在做毕业设计的过程中他们给予我莫大的帮助。我还要感谢所有教我的老师及实验室的指导老师，是你们的谆谆教诲，使我拥有了丰富的专业知识。
