1、 毕业论文(设计)题 目 FIR 数字滤波器设计与仿真学生姓名 学 号 院 系 专 业 指导教师 二一一 年六月二日目 录1. 引言 11.1 课题研究的目的和意义 .11.2 国内外发展现状 21.3 本文的研究工作 22. 数字滤波器 32.1 数字滤波器的分类 32.2 FIR 滤波器特点 32.3 数字滤波器的设计 43. FIR 滤波器的设计 .53.1 FIR 数字低通滤波器的设计 53.2 FIR 数字高通滤波器设计 93.3 FIR 带通滤波器设计 .113.4 FIR 数字带阻滤波器设计 .154. FIR 数字滤波器的实例仿真 19总 结 .21参考文献 22致 谢 .23
2、Abstract 241FIR 数字滤波器设计与仿真宗铭南京信息工程大学电子信息工程系，南京 210044摘 要：数字滤波器可以利用软件和数字处理硬件两种方式实现，本文在 matlab基础上采用了窗函数设计法、波纹逼近设计法、最小乘方设计法分别设计了低通、带阻、带通、高通 FIR数字滤波器，比较了不同方法设计出滤波器的通带、阻带、过渡带以及阶数，通过分析大量的设计仿真指出了不同设计方法的特点：窗函数设计法设计方法简便，通带内稳定性高，过渡带宽度随滤波器的阶数变化而变化；波纹逼近法设计出的滤波器阶数较低，计算量较小；最小乘方法误差较低，通带内平稳度高并且最小阻带衰减高。关键词: FIR；Matl
3、ab；数字滤波器1. 引言1.1 课题研究的目的和意义数字信号处理（Digital Signal Processing，简称 DSP）是一门涉及多门学科并广泛应用于很多科学和工程领域的新兴学科，它是指利用计算机或专用处理设备，以数字形式对信号进行采集、变换、滤波、估值、增强、压缩、识别等处理，以便提取有用的信息并进行有效的传输与应用，得到符合人们需要的信号形式。与模拟信号处理相比，数字信号处理具有精确、灵活、抗干扰能力强、可靠性高、体积小、易于大规模集成等优点。数字信号处理的实现在早期采用通用的处理器实现，但其速度低，价格昂贵。数字信号处理器（Digital Signal Processor）
4、的出现则改变了这一状况，其内部带有专门用于数字信号处理的硬件，因此速度非常快，方便实时快速地实现各种数字信号处理算法。数字滤波器由数字乘法器、加法器和延时单元组成的一种算法或装置。数字滤波器的功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。与模拟处理相比，数字处理有灵活性强、精度高、处理成本低以及对环境没有特殊要求等优点。它不仅能实现模拟处理的大部分功能，而且还能完成模拟处理由于成本、可靠性等原因而无法具体实现的功能。数字滤波与模拟滤波相比有很多优点，它除了可避免模拟滤波器固有电压漂移、温度漂移和噪声等问题外，还能满足滤波器对幅度和相位的严格要求。数字滤波在通信、图像编码
5、、语音编码、雷达等许多领域中也有着十分广泛的应用。数字滤波器是通过对抽样数据进行数学处理来达到频域滤波的目的。数字滤波器可以用软2件或设计专用的数字处理硬件两种方式来实现。用软件来实现数字滤波器优点是：随着滤波器参数的改变，很容易改变滤波器的性能。1.2 国内外发展现状现代滤波理论研究的主要内容是从含有噪声的数据记录（又称为时间序列）中估计出信号的某些特征或信号本身。一旦信号被估计出，那么估计出的信号将比原信号会有高的信噪比。现代滤波器把信号和噪声都视为随机信号，利用它们的统计特征（如自相关函数、功率谱函数等等）导出一套最佳的估值算法，然后用硬件和软件实现。目前现代滤波器主要有：维纳滤波器、卡
6、尔曼滤波器、线性预测器、自适应滤波器等，很多专家将基于特征分解的频率估计及奇异值分解算法都归入现代滤波器的范畴。GR 雷达信号处理分析系统中的信号分析中的滑动平均谱和常规处理中的反褶积运算采用了现代处理的部分功能。MATLAB 是一套用于工程计算的可视化高性能语言与软件环境。它集成了数值计算、符号计算、可视化建模、仿真和图形处理等功能。MATLAB 语言运算符丰富，既具有结构化的控制语言，又能面向对象编程，语法限制不严格，具又强大的图形功能和功能强大的工具箱。它为数字滤波的研究和应用提供了一个直观、高效、便捷的利器。它以矩阵运算为基础，把计算、可视化、程序设计融合到了一个交互式的工作环境中。M
7、ATLAB 推出的工具箱使各个领域的研究人员可以直观方便地进行科学研究、工程应用，其中的信号处理（signal processing ） 、图像处理（image processing ） 、小波（wavelet ）等工具箱为数字滤波研究的蓬勃发展提供了有力的工具。其中 MATLAB 的信号处理工具箱是专门应用于信号处理领域的专用工具箱，它的两个基本组成就是滤波器的设计与实现部分以及谱分析部分。工具箱提供了丰富而简便的设计，使原来繁琐的程序设计简化成函数的调用。只要以正确的指标参数调用相应的滤波器设计程序或工具箱函数，便可以得到正确的设计结果，所以使用非常方便。MATLAB 也是大量的计算机仿真
8、软件中的优秀代表，它在科学研究特别是电子信息科学中有着极为广泛的应用。MATLAB 现已被广泛应用于数学、通信、信号处理、自动控制、神经网络、图形处理等许多不同学科的研究中。1.3 本文的研究工作本文在 MATLAB 环境下完成 FIR 滤波器的设计仿真。采用窗函数设计法、最优等波动设计法、最小二乘法设计法分别设计低通、带阻、带通、高通 FIR 数字滤波器，通过大量仿真比较不同方法设计的滤波器幅频特性的通带平稳度、过渡带宽度、最小阻带衰减以及滤波器阶数，分析不同设计方法的特点。32. 数字滤波器所谓数字滤波器，是指输入、输出均为数字信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤
9、除某些频率成分的器件。数字滤波器和模拟滤波器相比，因为信号的形式和实现滤波的方法不同，数字滤波器具有比模拟滤波器精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配等优点。2.1 数字滤波器的分类数字滤波器按照不同的分类方法，有许多种类，但总起来可以分成两大类。一类称为经典滤波器，即一般的滤波器，特点是输入信号中有用的频率成分和希望滤除的频率成分各占有不同的频带，通过一个合适的选频率波器达到滤波的目的。例如，输入信号中含有干扰，如果信号和干扰的频带互不重叠，可滤除干扰得到纯信号。但对于一般滤波器如果信号和干扰的频带相互重叠，则不用完成对干扰的有效滤除，对时需要采用另一类所谓的现代滤波器，例如维纳
10、滤波器、卡尔曼滤波器、自适应滤波器等最佳滤波器。这些滤波器可按照随机信号内部的一些统计分布规律，从干扰中最佳地提取信号。一般数字滤波器从功能上分类，和模拟滤波器一样，可以分成低通、高通、带通和带阻等滤波器。这种理想滤波器是不可能实现的，因为它们的单位脉冲响应均是非因果且是无限长的，我们只能按照某些准则设计滤波器，使之尽可能逼近它，着些理想滤波器可作为逼近的标准用。另外，需要注意的是数字滤波器的传输函数 都是以 为周期的，滤波器的低通频率带)(jeH2处于 的整数倍数，而高频频带处于 的奇数倍附近，这一点和模拟滤波器是由区别的。2数字滤波器按照滤波器的网络结构或者单位脉冲响应分类:可分为 IIR
11、 滤波器(即无限长单位冲激响应滤波器)和 FIR 滤波器(即有限长单位冲激响应滤波器)。2.2 FIR 滤波器特点IIR 数字滤波器的设计简单方便，特别是采用双线性变换法设计的数字滤波器没有频谱混叠问题，效果较好。在有些实际应用场合，例如数据传输、图象处理等对滤波器的线性相位特性要求颇为严格，所以在这场合中一般选用 FIR 数字滤波器。FIR 数字滤波器可设计得具有严格的线性相位，而且其幅度特性可以随意设计。FIR 滤波器的单位冲激响应 h(n)是有限长序列，其系统函数的极点位于 Z 平面原点，因此 FIR 滤波器不存在稳定问题。另外，FIR 滤波器还可以采用 FFT 方法实现其功能，从而大大
12、提高了效率，因此，FIR 数字滤波器日益引起人们的注意4对于线性相位的滤波器而言，通过滤波器的信号不存在相位失真，只存在固定单位的延迟，而这一点在数字通信、图象处理、语音合成等诸多应用中都极其重要。FIR 数字滤波器具有线性相位特征的条件是它的单位冲激响应 h(n)满足偶对称或者奇对称。若 h(n)满足偶对称，则其相位特性 为：)(（1）21N若 h(n)满足奇对称，则其相位特性 为：)(（2）从频域上看 h(n)是长度为 N 的实序列，所以 h(n)的 N 点 DFT 满足共轭对称性，即：H(k)=H*(N-k)而幅度函数则一定呈偶对称特性，即：（3）)()kHk当 h(n)为偶对称，信号通
13、过该滤波器时，其通带就产生(N-1)/2 点群延迟，滤波器的相频特性是相位过原点的一条直线，其斜率为-(N-1)/2，其幅频特性相当于对信号低通滤波；当h(n)为奇对称，信号通过该滤波器时，其通带不但产生(N-1)/2 点群延迟，还对所有通带内各种频率的信号都有一个 90移相，这相当于信号先经过一个 90移相器，然后再作带通滤波。FIR 滤波器灵活，易于实现某些特殊的应用，如构建微分或积分器。2.3 数字滤波器的设计数字滤波器是一种用来过滤时间离散信号的数字系统，它是通过对抽样数据进行数学处理来达到频域滤波的目的。数字滤波器可以用软件或设计专用的数字处理硬件两种方式来实现。用软件来实现数字滤波
14、器优点是：随着滤波器参数的改变，很容易改变滤波器的性能。IIR 滤波器和 FIR 滤波器的设计方法是很不相同的。IIR 滤波器设计方法有两类，经常用的一类设计方法是借助于模拟滤波器的设计方法进行的。其设计步骤是；先设计模拟滤波器得到传输函数 ，然后将 按某种方法转换成数字滤波器的系统函数 。这一类方法)(sHA)(sA )(zH相对容易一些，这是因为模拟滤波器设计方法已经很成熟，它不仅有完整的设计公式，还有完善的图表供查阅；另外，还有一些典型的滤波器类型可供我们使用。另一类是直接在频域或者时域中进行设计的，由于要解连立问题方程，设计时需要计算机作辅助设计。FIR 滤波器不能采用由模拟滤波器的设
15、计进行转换得方法，经常用的是窗函数法和频率采样法。还有一种比较有效的方法是切毕雪夫等波纹逼近法，需通过计算机辅助设计完成。对于线性相位滤波器，通常采用 FIR 滤波器，其单位脉冲响应应满足一定条件，可以证明其相位特性在整个频带中是严格线性的，这是模拟滤波器无法达到的。当然，也可以采用 IIR滤波器，但必须使用全通网络对其非线性相位特性进行相位校正，这样增加了设计与实现的复5杂性。一般情况下，数字滤波器是一个线性移不变离散系统，利用有限精度算法来实现。数字滤波器的设计一般包括：（1）按照任务的要求，确定滤波器的性能要求；（2）用一个因果稳定的离散线性移不变系统的系统函数去逼近这一性能要求。系统函
16、数有无限长单位冲激响应（IIR）系统函数及有限长单位冲激响应（FIR）系统函数两种；（3）利用有限精度算法来实现这个系统函数。这里包括选择运算结构、选择合适的字长（包括系数量化及输入变量、中间变量和输出变量的量化）以及有效数字的处理方法（舍入、截尾）等；（4）实际的技术实现，包括采用通用计算机软件或专用数字滤波器硬件来实现，或是用专用的或通用的数字信号处理器来实现。一般来说，滤波器的性能要求往往以频率响应的幅度特性的允许误差来表征。数字滤波器的频响特性函数 一般为复函数，所以通常表示为:)(jeH（4）)()(jjj eHe其中 称为幅频特性函数； 称为相频函数。幅频特性表示信号通过该滤波器后
17、)(j 各频率成分的衰减情况，而相频特性反映各频率通过滤波器后在时间上的延时情况。3. FIR 滤波器的设计FIR 滤波器的设计方法很多，文中选取通常使用的用窗函数法、波纹逼近法、最小乘方法设计 FIR 滤波器。采用不同的方法设计相同的滤波器（低通、高通、带阻、带通滤波器）并进行比较，分析不同设计方法的特点。3.1 FIR 数字低通滤波器的设计 分别采用窗函数法（凯泽窗、汉宁窗、汉明窗、布莱克窗窗） 、波纹逼近法、最小乘方法对 FIR 数字低通滤波器进行设计并比较。3.1.1 窗函数法凯泽窗设计如下：Fs=48000HzFpass=9500HzFstop=12000Hz6Apass=1dBAs
18、top=80dB图 1 凯泽窗函数法设计 101 阶 FIR 数字低通滤波器汉宁窗、汉明窗与布莱克曼窗设计如下：Fs=48000HzFc=10800Hz图 2 汉宁窗函数法设计 101 阶 FIR 数字低通滤波器7图 3 汉明窗函数法设计 101 阶 FIR 数字低通滤波器图 4 布莱克曼窗函数法设计 101 阶 FIR 数字低通滤波器比较三种不同窗函数设计的滤波器可以看出：在相同阶数下，凯泽窗最小阻带衰减最大，但过度带较宽。3.1.2 波纹逼近法设计如下：Fs=48000Hz;Fpass=9500Hz;Fstop=12000Hz;Apass=1dB;Astop=80dB;Filter ord
19、er 选择 Minimun order（matlab 自动选取滤波器最小阶数）8图 5 波纹逼近法设计 101 阶 FIR 数字低通滤波器3.1.3 最小乘方法设计如下：Fs=48000Hz;Fpass=9500Hz;Fstop=12000Hz;Apass=1dB;Astop=80dB;Filter order 选择 specify order：101图 6 最小乘方法设计 101 阶 FIR 数字低通滤波器通过比较窗函数法、波纹逼近法以及最小乘方法对相同低通滤波器的设计可以看出：窗函数法通带内平稳度高，过渡带较宽，阻带衰减随频率增加而减小，阶数较高；波纹逼近法在通带内有较小波纹，过渡带较宽，
20、阻带内衰减平稳，阶数较高；最小乘方法通带内平稳度高，过渡带较小，阻带衰减随频率减小。93.2 FIR 数字高通滤波器设计分别采用窗函数法（凯泽窗、汉宁窗、汉明窗、布莱克窗窗） 、波纹逼近法、最小乘方法对FIR 数字高通滤波器进行设计并比较。3.2.1 窗函数法设计如下：Fs=48000HzFpass=12000HzFstop=9600HzApass=80dBAstop=1dB图 7 凯泽窗函数法设计 102 阶 FIR 数字高通滤波器汉宁窗、汉明窗与布莱克曼窗设计如下：Fs=48000HzFc=10800Hz图 8 汉宁窗窗函数法设计 102 阶 FIR 数字高通滤波器10图 9 汉明窗函数法
21、设计 102 阶 FIR 数字高通滤波器图 10 布莱克曼窗窗函数法设计 102 阶 FIR 数字高通滤波器比较三种不同窗函数设计的滤波器可以看出：在相同阶数下，凯泽窗最小阻带衰减最大，但过度带较宽。3.2.2 波纹逼近法设计如下：Fs=48000Hz;Fpass=12000Hz;Fstop=9600Hz;Apass=80dB;Astop=1dB;11图 11 波纹逼近法设计 56 阶 FIR 数字高通滤波器3.2.3 最小乘方法设计如下：Fs=48000HzFpass=12000HzFstop=9600HzApass=80dBAstop=1dB图 12 最小乘方法设计 101 阶 FIR 数
22、字高通滤波器通过比较三种不同方法设计相同高通滤波器可以看出：窗函数法阻带衰减随频率增加而增加，过渡带较宽，通带平稳性高，阶数较高；波纹逼近法阻带衰减平稳，过渡带较宽，通带内有较小波纹，阶数较低；最小乘方法最小阻带衰减高，过渡带较宽，通带内平稳度高。3.3 FIR 带通滤波器设计分别采用窗函数法（凯泽窗、汉宁窗、汉明窗、布莱克窗窗） 、波纹逼近法、最小乘方法对12FIR 数字带通滤波器进行设计并比较。3.3.1 窗函数法设计如下：Fs=48000HzFstop1=8400HzFpass1=9600HzFpass2=12000HzFstop2=14400HzAstop1=60dBApass=1dB
23、Astop2=80dB图 13 凯泽窗函数法设计 101 阶 FIR 数字带通滤波器汉宁窗、汉明窗与布莱克曼窗设计如下：Fs=48000HzFc1=8400HzFc2=13200Hz13图 14 汉宁窗窗函数法设计 101 阶 FIR 数字带通滤波器图 15 汉明窗窗函数法设计 101 阶 FIR 数字带通滤波器图 16 布莱克曼窗函数法设计 101 阶 FIR 数字带通滤波器比较三种不同窗函数设计的滤波器可以看出：在相同阶数下，凯泽窗最小阻带衰减最大，但过度带较宽。143.3.2 波纹逼近法设计如下：Fs=48000HzFstop1=7200HzFpass1=9600HzFpass2=120
24、00HzFstop2=14400HzAstop1=60dBApass=1dBAstop2=80dB图 17 波纹逼近法设计 50 阶 FIR 数字带通滤波器3.3.3 最小乘方法设计如下：Fs=48000HzFstop1=8400HzFpass1=9600HzFpass2=12000HzFstop2=14400HzWstop1=60dBWpass=1dBWstop2=80dB15图 18 最小乘方法设计 101 阶 FIR 数字带通滤波器通过不同方法设计滤波器进行比较得出：窗函数法阻带衰减随频率变化，过渡带较宽，通带平稳度高；波纹逼近法阻带内平稳，过渡带较宽，通带内有较小波纹，阶数较低；最小乘
25、方法最小阻带衰减高，过渡带较宽，通带内平稳度高。3.4 FIR 数字带阻滤波器设计分别采用窗函数法（凯泽窗、汉宁窗、汉明窗、布莱克窗窗） 、波纹逼近法、最小乘方法对FIR 数字带阻滤波器进行设计并比较。3.4.1 窗函数法设计如下：Fs=48000HzFpass1=8400HzFstop1=9600HzFstop2=12000HzFpass2=14400HzApass1=0.5dBAstop=60dBApass2=1dB16图 19 凯泽窗函数法设计 74 阶 FIR 数字带阻滤波器汉宁窗、汉明窗与布莱克曼窗设计如下：Fs=48000HzFc1=8400Hz Fc2=13200Hz图 20 汉
26、宁窗函数法设计 74 阶 FIR 数字带阻滤波器图 21 汉明窗函数法设计 74 阶 FIR 数字带阻滤波器17图 22 布莱克曼窗函数法设计 74 阶 FIR 数字带阻滤波器比较三种不同窗函数设计的滤波器可以看出：在相同阶数下，凯泽窗最小阻带衰减最大，但过度带较宽。 3.4.2 波纹逼近法设计如下：Fs=48000HzFpass1=8400HzFstop1=9600HzFstop2=12000HzFpass2=14400HzApass1=0.5dBAstop=60dBApass2=1dB图 23 波纹逼近法设计 46 阶 FIR 数字带阻滤波器183.4.3 最小乘方法设计如下：Fs=480
27、00HzFpass1=8400HzFstop1=9600HzFstop2=12000HzFpass2=14400HzWpass1=0.5dBWstop=60dBWpass2=1dB图 24 最小乘方法设计 50 阶 FIR 数字带阻滤波器通过几种设计方法的比较可以得出：窗函数法通带内平稳度高，过渡带较宽，阻带衰减随频率变化，阶数较高；波纹逼近法通带内有较小波纹，过渡带较宽，最小阻带衰减高，阶数低；最小乘方法通带内平稳度高，过渡带较宽，最小阻带衰减高。 综合以上设计得出以下表格：表一几种不同设计法的特点比较方法 参数 通带内平稳度 过渡带宽度 最小阻带衰减 阶数凯泽窗 高 较宽 高 高汉宁窗 高
28、 一般 低 高汉明窗 高 一般 低 高布莱克曼窗 高 一般 低 高波纹逼近 有波纹 一般 不定 低最小乘方 高 一般 高 高194. FIR 数字滤波器的实例仿真假设原始信号是 x=sin(2*pi*70*t)+2*sin(2*pi*120*t)，抽样频率为 fs=1000Hz，信号被白噪声污染，实际获得信号是 xn=x+randn(size(t)，设计出一个 FIR 滤波器来恢复出原始信号。因为白噪声分布于整个频带，所以要在不衰减原信号的前提之下，对整个频带滤波，利用最小二乘法设计出一个多带滤波器。滤波器设计要求在65/500 75/500（即0.13 0.15）段、115/500 125/
29、500（即0.23 0.25）段频带内幅度为 l，在 0 60/500（即0 0.12）段、80/500 110/500（即0.16 0.22）段、130/500 1（即0.26 1）段频带内幅为 0。程序设计如下：t=0:1/fs:2; %取 2 秒的序列x=sin(2*pi*70*t)+2*sin(2*pi*120*t); %先生成相应信号xn=x+randn(size(t); %加入噪声f=0 0.12 0.13 0.15 0.16 0.22 0.23 0.25 0.26 1;m=0 0 1 1 0 0 1 1 0 0;b=firls(N,f,m);h,b=freqz(b,1);滤波器幅
30、频特性如图 25 所示：图 26 为滤波信号和原始信号比较：程序如下：t=0:1/fs:0.75; x=sin(2*pi*70*t)+2*sin(2*pi*120*t); xn=x+randn(size(t); subplot(3,1,1);plot(t,x); %原始信号subplot(3,1,2);plot(t,xn); %污染信号b=firls(N,f,m);xo=filter(b,1,xn); %滤波器对 xn 进行滤波subplot(3,1,3);plot(t,xo);20图 25 滤波器幅频特性图 26 滤波后信号比较21总 结本文在 matlab 环境下采用了窗函数设计法、波纹逼
31、近设计法、最小乘方设计法分别设计了低通、带阻、带通、高通 FIR 数字滤波器，采用控制变量法比较相同阶数下不同设计方法设计出的滤波器的通带平稳度、阻带最小衰减以及过渡带的宽度，通过对大量的设计仿真进行分析，指出了不同设计方法的特点：在同样滤波器阶数情况下，窗函数设计法设计方法相对简便，通带内稳定性高，过渡带宽度随窗函数种类变化，其中汉宁窗、汉明窗、布莱克曼窗函数法的过渡带相比凯泽窗较窄，但凯泽窗最小阻带衰减最高；波纹逼近法设计出的滤波器相对其他几种方法通带内不够平稳；最小乘方法相比于其他设计方法误差较低，通带内平稳度较高并且最小阻带衰减较高。通过分析推广到实际设计中去：对通带内平稳度要求较高时
32、优先选用窗函数法；对最小阻带衰减要求较高时可优先选择最小乘方法；要求阶数或计算量较低时可优先选择波纹逼近法。本文还存在可以深入之处，其一，选取的设计方法只是三种最常用的，实际设计中还有很多种方法，选取更多的方法可以得到更加细致的结论；其二，比较时选取的通带平稳度、过渡带宽度、最小阻带衰减与阶数，其实能够体现滤波器性能的还有主瓣宽度、旁瓣峰值等，增加一些基本参数进行比较可以使滤波器设计方法的比较更有说服力；其三，数字滤波器分为 FIR与 IIR，可以在这基础上用类似方法研究 IIR 滤波器，得出更加广泛的结论。22参考文献1郑君里.信号与系统M.北京：高等教育出版社，2000.2楼顺天基于 MA
33、TLAB的系统分析与设计-信号处理 M西安：西安电子科技大学出版社，19983程卫国MATLAB5.3 应用指南 M北京：人民邮电出版社，19994钟麟MATLAB 仿真技术与应用教程M北京：国防工业出版社，2004.5刘敏MATLAB 通信仿真与应用M北京：国防工业出版社，20016丁玉美数字信号处理M西安：西安电子科技大学出版社，20017何振亚数字信号处理的理论与应用M北京：人民邮电出版社，19838郭德才.基于 Matlab的 FIR数字低通滤波器的设计与仿真J.吉林大学，20099 陈明军. 改进窗函数在 FIR数字滤波器设计中的应用 J . 继电器, 2007.10 程明. 基于
34、MATLAB的数字滤波器分析J. 通信电源技术, 2006.11McClellan J H，Parks T WA Unified Approach to the Design of Optimum FIR Linear-Phase Digital FiltersIEEE TransSCircuit Theory，197312Steigilitz KComputer-Aided Design of Recursive Digital FiltersIEEE TransJAudio Electroacoust，June 1970 13Oppenheim A V，Schafer R W，Digita
35、l signal ProcessingMPrentice Hall，Inc，197514Stanley W DDigital Signal ProcessingMReston Publishing Company，Inc，1975 15Johan Baeten, Joris De Schutter. Appendix D: Image Processing Implementation on DSPJ.Journal of Mathematical Imageing and Vision. 2003.23致 谢转眼 2 个月过去了，我的毕业设计也告一段落，回顾这 2 个月，我收获了很多，这次设
36、计是整个大学期间学习知识综合运用的检验，也是对个人能力的一次锻炼。从开始的选题到最后的仿真都不容易，都是对我所积累知识的考验与挑战，与此同时，我巩固了大学所学理论知识，能力也得到了提升。这对于我今后踏上社会工作岗位都会起到积极的影响，我将终生受益。我的毕业设计是在行鸿彦教授的悉心指导下完成的。在毕业设计这几个月里，行教授给了我很大的帮助和指导。在此，我要向行教授表示我最衷心的感谢。还有各位同学的关爱与帮助，才使我的毕业设计得以顺利进行。在此，我向给与我大力帮助的各位同学表示我衷心的感谢！24The design and Simulation of FIR Digital FilterZong Ming(Department of Electronic the best method to deal with such fluctuations in the window method precision design difficulties; constrained minimum squares design method uses techniques which make the squared error minimum. Key words: FIR ; Matlab ; digital filter
