FIR数字滤波器设计与仿真

毕业论文(设计)题目FIR数字滤波器设计与仿真学生姓名学号院系专业指导教师二○一一年六月二日目录1.引言.......................................................11.1课题研究的目的和意义...................................11.2国内外发展现状.........................................21.3本文的研究工作.........................................22.数字滤波器..................................................32.1数字滤波器的分类.......................................32.2FIR滤波器特点..........................................32.3数字滤波器的设计.......................................43.FIR滤波器的设计............................................53.1FIR数字低通滤波器的设计................................53.2FIR数字高通滤波器设计..................................93.3FIR带通滤波器设计.....................................113.4FIR数字带阻滤波器设计.................................154.FIR数字滤波器的实例仿真...................................19总结.........................................................21参考文献......................................................22致谢.........................................................23Abstract......................................................241FIR数字滤波器设计与仿真宗铭南京信息工程大学电子信息工程系，南京210044摘要：数字滤波器可以利用软件和数字处理硬件两种方式实现，本文在matlab基础上采用了窗函数设计法、波纹逼近设计法、最小乘方设计法分别设计了低通、带阻、带通、高通FIR数字滤波器，比较了不同方法设计出滤波器的通带、阻带、过渡带以及阶数，通过分析大量的设计仿真指出了不同设计方法的特点：窗函数设计法设计方法简便，通带内稳定性高，过渡带宽度随滤波器的阶数变化而变化；波纹逼近法设计出的滤波器阶数较低，计算量较小；最小乘方法误差较低，通带内平稳度高并且最小阻带衰减高。关键词:FIR；Matlab；数字滤波器1.引言1.1课题研究的目的和意义数字信号处理（DigitalSignalProcessing，简称DSP）是一门涉及多门学科并广泛应用于很多科学和工程领域的新兴学科，它是指利用计算机或专用处理设备，以数字形式对信号进行采集、变换、滤波、估值、增强、压缩、识别等处理，以便提取有用的信息并进行有效的传输与应用，得到符合人们需要的信号形式。与模拟信号处理相比，数字信号处理具有精确、灵活、抗干扰能力强、可靠性高、体积小、易于大规模集成等优点。数字信号处理的实现在早期采用通用的处理器实现，但其速度低，价格昂贵。数字信号处理器（DigitalSignalProcessor）的出现则改变了这一状况，其内部带有专门用于数字信号处理的硬件，因此速度非常快，方便实时快速地实现各种数字信号处理算法。数字滤波器由数字乘法器、加法器和延时单元组成的一种算法或装置。数字滤波器的功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。与模拟处理相比，数字处理有灵活性强、精度高、处理成本低以及对环境没有特殊要求等优点。它不仅能实现模拟处理的大部分功能，而且还能完成模拟处理由于成本、可靠性等原因而无法具体实现的功能。数字滤波与模拟滤波相比有很多优点，它除了可避免模拟滤波器固有电压漂移、温度漂移和噪声等问题外，还能满足滤波器对幅度和相位的严格要求。数字滤波在通信、图像编码、语音编码、雷达等许多领域中也有着十分广泛的应用。数字滤波器是通过对抽样数据进行数学处理来达到频域滤波的目的。数字滤波器可以用软件或设计专用2的数字处理硬件两种方式来实现。用软件来实现数字滤波器优点是：随着滤波器参数的改变，很容易改变滤波器的性能。1.2国内外发展现状现代滤波理论研究的主要内容是从含有噪声的数据记录（又称为时间序列）中估计出信号的某些特征或信号本身。一旦信号被估计出，那么估计出的信号将比原信号会有高的信噪比。现代滤波器把信号和噪声都视为随机信号，利用它们的统计特征（如自相关函数、功率谱函数等等）导出一套最佳的估值算法，然后用硬件和软件实现。目前现代滤波器主要有：维纳滤波器、卡尔曼滤波器、线性预测器、自适应滤波器等，很多专家将基于特征分解的频率估计及奇异值分解算法都归入现代滤波器的范畴。GR雷达信号处理分析系统中的信号分析中的滑动平均谱和常规处理中的反褶积运算采用了现代处理的部分功能。MATLAB是一套用于工程计算的可视化高性能语言与软件环境。它集成了数值计算、符号计算、可视化建模、仿真和图形处理等功能。MATLAB语言运算符丰富，既具有结构化的控制语言，又能面向对象编程，语法限制不严格，具又强大的图形功能和功能强大的工具箱。它为数字滤波的研究和应用提供了一个直观、高效、便捷的利器。它以矩阵运算为基础，把计算、可视化、程序设计融合到了一个交互式的工作环境中。MATLAB推出的工具箱使各个领域的研究人员可以直观方便地进行科学研究、工程应用，其中的信号处理（signalprocessing）、图像处理（imageprocessing）、小波（wavelet）等工具箱为数字滤波研究的蓬勃发展提供了有力的工具。其中MATLAB的信号处理工具箱是专门应用于信号处理领域的专用工具箱，它的两个基本组成就是滤波器的设计与实现部分以及谱分析部分。工具箱提供了丰富而简便的设计，使原来繁琐的程序设计简化成函数的调用。只要以正确的指标参数调用相应的滤波器设计程序或工具箱函数，便可以得到正确的设计结果，所以使用非常方便。MATLAB也是大量的计算机仿真软件中的优秀代表，它在科学研究特别是电子信息科学中有着极为广泛的应用。MATLAB现已被广泛应用于数学、通信、信号处理、自动控制、神经网络、图形处理等许多不同学科的研究中。1.3本文的研究工作本文在MATLAB环境下完成FIR滤波器的设计仿真。采用窗函数设计法、最优等波动设计法、最小二乘法设计法分别设计低通、带阻、带通、高通FIR数字滤波器，通过大量仿真比较不同方法设计的滤波器幅频特性的通带平稳度、过渡带宽度、最小阻带衰减以及滤波器阶数，分析不同设计方法的特点。32.数字滤波器所谓数字滤波器，是指输入、输出均为数字信号，通过一定运算关系改变输入信号所含频率成分的相对比例或者滤除某些频率成分的器件。数字滤波器和模拟滤波器相比，因为信号的形式和实现滤波的方法不同，数字滤波器具有比模拟滤波器精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配等优点。2.1数字滤波器的分类数字滤波器按照不同的分类方法，有许多种类，但总起来可以分成两大类。一类称为经典滤波器，即一般的滤波器，特点是输入信号中有用的频率成分和希望滤除的频率成分各占有不同的频带，通过一个合适的选频率波器达到滤波的目的。例如，输入信号中含有干扰，如果信号和干扰的频带互不重叠，可滤除干扰得到纯信号。但对于一般滤波器如果信号和干扰的频带相互重叠，则不用完成对干扰的有效滤除，对时需要采用另一类所谓的现代滤波器，例如维纳滤波器、卡尔曼滤波器、自适应滤波器等最佳滤波器。这些滤波器可按照随机信号内部的一些统计分布规律，从干扰中最佳地提取信号。一般数字滤波器从功能上分类，和模拟滤波器一样，可以分成低通、高通、带通和带阻等滤波器。这种理想滤波器是不可能实现的，因为它们的单位脉冲响应均是非因果且是无限长的，我们只能按照某些准则设计滤波器，使之尽可能逼近它，着些理想滤波器可作为逼近的标准用。另外，需要注意的是数字滤波器的传输函数)(jeH都是以2为周期的，滤波器的低通频率带处于2的整数倍数，而高频频带处于的奇数倍附近，这一点和模拟滤波器是由区别的。数字滤波器按照滤波器的网络结构或者单位脉冲响应分类:可分为IIR滤波器(即无限长单位冲激响应滤波器)和FIR滤波器(即有限长单位冲激响应滤波器)。2.2FIR滤波器特点IIR数字滤波器的设计简单方便，特别是采用双线性变换法设计的数字滤波器没有频谱混叠问题，效果较好。在有些实际应用场合，例如数据传输、图象处理等对滤波器的线性相位特性要求颇为严格，所以在这场合中一般选用FIR数字滤波器。FIR数字滤波器可设计得具有严格的线性相位，而且其幅度特性可以随意设计。FIR滤波器的单位冲激响应h(n)是有限长序列，其系统函数的极点位于Z平面原点，因此FIR滤波器不存在稳定问题。另外，FIR滤波器还可以采用FFT方法实现其功能，从而大大提高了效率，因此，FIR数字滤波器日益引起人们的注意4对于线性相位的滤波器而言，通过滤波器的信号不存在相位失真，只存在固定单位的延迟，而这一点在数字通信、图象处理、语音合成等诸多应用中都极其重要。FIR数字滤波器具有线性相位特征的条件是它的单位冲激响应h(n)满足偶对称或者奇对称。若h(n)满足偶对称，则其相位特性)(为：21)(N（1）若h(n)满足奇对称，则其相位特性)(为：221)(N（2）从频域上看h(n)是长度为N的实序列，所以h(n)的N点DFT满足共轭对称性，即：H(k)=H*(N-k)而幅度函数则一定呈偶对称特性，即：)()(kNHkH（3）当h(n)为偶对称，信号通过该滤波器时，其通带就产生(N-1)/2点群延迟，滤波器的相频特性是相位过原点的一条直线，其斜率为[-(N-1)/2]，其幅频特性相当于对信号低通滤波；当h(n)为奇对称，信号通过该滤波器时，其通带不但产生(N-1)/2点群延迟，还对所有通带内各种频率的信号都有一个90°移相，这相当于信号先经过一个90°移相器，然后再作带通滤波。FIR滤波器灵活，易于实现某些特殊的应用，如构建微分或积分器。2.3数字滤波器的设计数字滤波器是一种用来过滤时间离散信号的数字系统，它是通过对抽样数据进行数学处理来达到频域滤波的目的。数字滤波器可以用软件或设计专用的数字处理硬件两种方式来实现。用软件来实现数字滤波器优点是：随着滤波器参数的改变，很容易改变滤波器的性能。IIR滤波器和FIR滤波器的设计方法是很不相同的。IIR滤波器设计方法有两类，经常用的一类设计方法是借助于模拟滤波器的设计方法进行的。其设计步骤是；先设计模拟滤波器得到传输函数)(sHA，然后将)(sHA按某种方法转换成数字滤波器的系统函数)(zH。这一类方法相对容易一些，这是因为模拟滤波器设计方法已经很成熟，它不仅有完整的设计公式，还有完善的图表供查阅；另外，还有一些典型的滤波器类型可供我们使用。另一类是直接在频域或者时域中进行设计的，由于要解连立问题方程，设计时需要计算机作辅助设计。FIR滤波器不能采用由模拟滤波器的设计进行转换得方法，经常用的是窗函数法和频率采样法。还有一种比较有效的方法是切毕雪夫等波纹逼近法，需通过计算机辅助设计完成。对于线性相位滤波器，通常采用FIR滤波器，其单位脉冲响应应满足一定条件，可以证明其相位特性在整个频带中是严格线性的，这是模拟滤波器无法达到的。当