基于MATLAB的数字滤波器的设计综述.docx

基于MATLAB的数字滤波器的设计综述.docx

《基于MATLAB的数字滤波器的设计综述.docx》由会员分享，可在线阅读，更多相关《基于MATLAB的数字滤波器的设计综述.docx（39页珍藏版）》请在冰点文库上搜索。

基于MATLAB的数字滤波器的设计综述

本科毕业论文（设计）

论文题目

：

基于MATLAB的数字滤波器的设计

姓名

：

学号

：

班级

：

02班

年级

：

2011级

专业

：

通信工程

系部

：

信息工程学院

指导教师

：

（副教授）

完成时间

：

2015年4月15日

作者声明

本毕业论文（设计）是在导师的指导下由本人独立撰写完成的，没有剽窃、抄袭、造假等违反道德、学术规范和其他侵权行为。

对本论文（设计）的研究做出重要贡献的个人和集体，均已在文中以明确方式标明。

因本毕业论文（设计）引起的法律结果完全由本人承担。

毕业论文（设计）成果归武昌工学院所有。

特此声明

作者专业

：

通信工程

作者学号

：

作者签名

：

年月日

基于MATLAB的数字滤波器的设计

TheDesignOfDigitalFilterbasedonMATLAB

2015年4月15日

摘要

本文研究基于MATLAB的数字滤波器的设计。

首先介绍了本课题研究的目的及意义，调研了国内外数字滤波器应用现状及发展趋势；接着对MATLAB软件的发展历史、语言特点、功能特性进行了介绍；其次对数字滤波器进行了简介，其中包括数字滤波器的基本概念、类型、设计步骤、基本结构等；然后本文重点描述了IIR和FIR数字滤波器的设计过程和设计方法。

对于IIR数字滤波器，详细的分析了双线性变换法、脉冲响应不变法：

对于FIR数字滤波器，详细分析了窗函数法、最小二乘法、最优化设计法。

通过MATLAB软件对五种不同方法编写的程序进行仿真，分析仿真结果，比较各种设计方法的适用范围和它的不可替代性。

基于MATLAB的数字滤波器的设计，将数字滤波器的设计集于一个软件包，十分方便、直观、准确和可靠。

借助于MATLAB设计滤波器的这种优越性，推动了滤波技术的高速发展。

基于MATLAB语言可扩充的特点，提高了编程的效率，节约了大量的编程时间。

本文所研究的内容，对通过MATLAB软件设计数字滤波器的发展起到一定的促进和推动作用。

关键字：

IIR数字滤波器；FIR数字滤波器；脉冲响应不变法；双线性变换；窗函数；最小二乘法；最优化设计

Abstract

ThispaperisbasedonthedesignofMATLABdigitalfilter.ItfirstIntroducesthepurposeandsignificanceofthetopicandresearchestheapplicationstatusanddevelopmenttrendsofthedigitalfilterathomeandabroad;SecondlyitpresentsthehistoricaldevelopmentoftheMATLABsoftware.Thenitgivesabriefintroductionofthedigitalfilter,includingthebasicconceptsofdigitalfilters,types,designstepsandthebasicstructure;BesidesitemphaticallydescribesthedesignprocessanddesignmethodsIIRandFIRdigitalfilter.ForIIRdigitalfilter,adetailedanalysisofthebilineartransformandimpulseresponseinvarianceisgiven.ForFIRdigitalfilter,adetailedanalysisofthewindowfunctionmethod,leastsquaresmethod,optimizeddesignisgiven.SimulatingtheprogramwritteninfivedifferentmethodsthroughMATLABsoftware,analysisthesimulationresultsandthencomparethescopeofapplicationofvariousdesignmethodsandtheirirreplaceableability.ThedigitalfilterdesignbasedonMATLAB,settingthedigitalfilterdesigninapackage,veryconvenient,intuitive,accurateandreliable.Thissuperioritypromotestherapiddevelopmentoffilteringtechnology.ThefeaturesthatMATLABlanguagefeaturescanbeexpandedimprovestheprogram'sefficiencyandsavesalotofprogrammingtime.ThecontentofthispaperplaysacertainroleinpromotingandfacilitatingthedevelopmentofthedigitalfilterbyMATLABsoftwaredesign.

Keywords:

IIRdigitalfilter;FIRdigitalfilter;Impulseresponseinvariance;doublelineartransformation;Windowfunction;leastsquaresmethod;optimaldesign

目录

1引言

1.1课题研究的目的和意义

目前，数字信号处理理论与应用已成为一门及其重要的高新技术学科。

数字信号处理技术日益成熟,广泛应用于各个领域。

伴随着社会和经济的伟大变革，如今人人享有通信自由。

从过去的普通电话到现在随时随地视频聊天，移动电话不再是特殊阶层使用，已成为大众日常生活的普通工具。

网络资源的共享，使得学生或是工作人士轻而易举的获取到自己想要的资料。

科技带给我们生活的极大便利，可能很多人会想到计算机发挥了很大作用，但是，数字信号处理技术发挥的作用也是惊人的。

如果没有这一技术的快速发展和应用，可能我们还没有那么快就可以足不出户而正常生活着。

数字信号处理（DSP:

DigitalSignalProcessing）经过几十年的不断发展，一套较为完整的的学科领域和理论体系现在逐渐形成了。

滤波含有过滤之意，将信号中无用部分滤除，得到想要的部分，而滤波却是信号处理技术的一个重要分支。

凡是信号经过的系统，比如硬件系统、算法程序等软件系统，都会起到滤波的作用。

人们往往熟知通过不同的电子元器件组成的模拟滤波器，但随着计算机技术、数字电路技术、大规模集成电路技术的发展，数字滤波技术越来越受人们亲赖。

因为数字滤波器较之模拟滤波器有如下优势：

（1）精度高。

（2）灵活性大。

（3）性能可靠。

（4）易于大规模集成，设计简单。

（5）并行处理。

（6）性能参数易变。

基于上述众多优点，数字滤波器在许多的应用领域逐步代替了传统的模拟滤波器，并且正在各种科学研究或工程技术领域里快速发展。

所以，对数字滤波器的研究和分析有着很重要的实际需要。

数字滤波已成功应用在语音处理、图像处理、医疗、电视、通信、雷达、航天、声纳、音乐、军事等领域。

而基于MATLAB软件实现数字滤波器的设计，可通过修改滤波器的参数十分方便的改变滤波器的特性，进而解决一些特定需求。

因此我们有必要对滤波器的设计方法进行探讨和研究，理解滤波器工作原理，优化设计方法，设计开发出稳定性好的滤波器系统。

数字滤波涉及大量的运算，离开了计算机及相应的软件，就不可能解决任何稍微复杂的实际应用问题。

MATLAB已成为国际公认的数字滤波的标准软件和开发平台。

借助于MATLAB软件设计数字滤波器设计不仅方便而且易于实现，现已成为人们研究和学习的热点。

因此基于MATLAB的数字滤波器的设计，不仅有着重要的理论意义，而且有着重要的现实意义。

1.2国内外发展现状和趋势

在近代电信装备和各类控制系统中，滤波器应用极为广泛，凡是有能力进行信号处理的装置都可以称为滤波器；在所有的电子产品部件中，使用最广，技术最复杂要算滤波器了。

对滤波器的研究和生产历来为各国所重视，因为，滤波器的优劣直接决定产品的优劣。

早在上世纪40年代末，开始了关于数字滤波器的可能性的讨论，并在20世纪50年代也有人讨论了数字滤波的问题。

在此期间，科学家们做了大量卓有成效的工作，这些工作对于后来数字滤波器的快速发展，发挥了非常重要的作用。

到了上个世纪60年代中期，美国科学家库里，杜克总结了前人研究成果，经过长期研究，使得快速傅立叶变换的这一成果的到来。

在数字信号处理和数字滤波中，快速傅立叶变换的出现是学术界一直认可的一个科学里程碑。

自那时以来，各种快速傅立叶变换已经出现在数字滤波领域，在数字滤波器中开始形成了成套正规理论。

在这一时期，提出了各种各样的数字滤波结构，有的则以运算速度快见长，有的以运算误差是最小误差为特点，出现了数字滤波器的各种实现方法。

对递归和非递归两类滤波器作了全面的比较，统一了数字滤波器的基本概念和理论降数字滤波器的一个重要发展是对有限冲激响应和无限冲激响应关系的认识转化。

在初期，一般认为IIR滤波器比FIR滤波器具有更高的运算效率，因而明显地趋向于前者。

后来，当人们提出用快速傅利叶变换实现卷积运算的概念之后，发现高阶FIR滤波器也可以用很高的运算效率来实现，这就促使人们对高性能FIR滤波器的设计方法和数字滤波器的频域设计方法进行了大量的研究。

数字滤波器尽管在语音、声纳、地震的信号处理中曾经发挥过作用，但由于以前计算机主机的价格很昂贵，严重阻碍了数字滤波器的发展，随着计算机科学的飞速发展，数字滤波器越来越多的展示出它自身的优势，正在加速进入各种工业领域经过最近几十年的发展，线性滤波方法，如维纳滤波、卡尔曼滤波和自适应滤波得到了广泛的研究和应用。

一些非线性滤波方法，如小波滤波、同态滤波、中值滤波、形态滤波等都是现代信号处理的前沿课题，不但有重要的理论意义，而且有广阔的应用前景。

80年代致力于各类新型滤波器性能提高的研究并逐渐扩大应用范围。

90年代至今在主要致力于把各类滤波器应用于各类产品的开发和研制。

当然，对滤波器本身的研究仍在不断进行。

我国当时主要用于话路滤波和报路滤波，广泛使用滤波器是50年代后期的事。

我国数字滤波器在研制、生产和应用等方面正在快速发展，但由于缺少专门研制机构，经过半个世纪的发展，集成工艺和材料工业跟不上来，使得我国许多新型滤波器的研制应用与国际发展有一段距离。

模拟滤波器在国内的工程技术中应用中占了很大的比重，数字滤波器在生产应用中所占的比例还不多。

从应用比例来看，必须加大对这方面的研究投入，数字滤波器才能实现大规模应用。

在国内也很多学者在数字滤波领域展开长期的深入研究，随着我国科学技术的飞速发展，如山东大学的赖晓平教授、天津大学的王兆华教授。

无论是在理论上还是在工程技术领域，都有很多科研成果。

我国现有滤波器的种类和所覆盖的频率己基本上满足现有各种电信设备。

从整体而言，我国提出了一套有关数字滤波器设计的新算法和新设计方案设计空间方法，目前还处于研究阶段。

在理论研究方面，在研究用模糊数学来设计数字滤波器，在超大规模集成电路领域，对DF的研究也十分活跃。

对滤波器的性能要求越来越高，功能也越来越多。

对FIR滤波器而言，从理论上讲要尽可能的降低其阶数，因为随着阶数的增加，实现数字滤波器的运算量就越大，如果阶数不能降低，则对算法的运算效率就要提出更高的。

具体实现它们的硬件电子元器件必须向大规模集成化方向发展是一个重要的趋势要求，如果没有重大的理论突破，则只能在硬件寻找突破口，所以电子元件的集成化尤其重要。

我国数字滤波器研制和生产与上述要求相差甚远，电子工程和科技人员肩负着重大的历史责任，为缩短这个差距。

而且集成化在世界范围内也是一个重要的发展趋势。

由于缺少专门研制机构，集成工艺和材料工业跟不上来，虽然我国滤波器在研制、生产应用等方面已有一定进步，但使许多新型滤波器的研制应用与国际水平有一段距离。

1.3本文研究的主要内容

本次毕业设计中主要完成的内容包括：

（1）本课题研究的目的及意义，数字滤波器国内外的应用现状和发展趋势。

（2）MATLAB的简介，熟悉软件的功能和使用方法。

通过书籍阅读和文献搜索学习MATLAB和数字滤波器，并掌握MATLAB的使用和数字滤波器的设计方法。

（3）基于MATLAB的IIR数字滤波器的设计，对于IIR数字滤波器的设计，采用冲激响应不变法和双线性变换法，编写程序仿真。

（4）基于MATLAB的FIR数字滤波器的设计，对于IIR数字滤波器的设计，采用窗函数法，最小二乘法，最优化设计法，并进行仿真。

（5）通过仿真结果分析这几种设计数字滤波器的方法。

2MATLAB软件概述

2.1发展历程及语言特点

20世纪70年代中期，CleveMoler博士和其同事开发了调用EISPACK和LINPACK的FORTRAN子程序库。

其中，EIPACK用来求解特征值，LINPACK用来解线性方程。

同年后期，CleveMoler担任美国墨西哥大学计算机系的系主任，为了让学生方便的调用EISPACK和LINPACK，他设计了接口程序，并取名MATLAB（MAtrixLABoratory矩阵实验室）。

早期的MATLAB是用Fortran编写的，只能进行矩阵运算，绘图也只能用星号描点等简单形式，只提供了几十个内部函数。

虽然其功能如此简单，当作为免费软件推出以后，还是吸引了大批的使用者。

第一个商业化的MATLAB是于1984年由MATHWORKS公司推出的，该版本基于当时流行的DOS操作系统，本身也由C语言重新编写。

1992年具有划时代意义的MATLAB4.0版本发布，用户数剧增。

1994年的4.2版本扩充了4.0版本的功能，尤其在图形界面设计方面提供了新的方法。

1997年推出的5.0的版本允许了更多的数据结构，使其扩展为一种非常方便编程的语言工具。

1999年推出的MATLAB5.3版本在很多方面又进一步提高了MATLAB语言的功能。

2000年10月底推出了其全新的MATLAB6.0正式版，在核心数值算法、界面设计、外部接口、应用桌面等诸多方面有了极大的改进。

2004年推出7.0版本，新增12个新产品模块，升级了28个产品模块。

对MATLAB编程环境、代码效率、数据可视化、数学计算、文件I/O等方面进行了升级。

在这之后，基本每年都会稳定的发布两个更新版本。

现在的MATLAB支持各种操作系统，能够运行在十几个操作系统平台。

MATLAB再也不是一个简单的矩阵实验室了，逐渐演变成为一种基于算法开发、数据可视化、数据分析以及数值计算的计算机高级编程语言和交互式环境。

MATLAB语言是一种交互性的数学脚本语言，开发环境直接简洁，其语法与C/C++类似。

它支持包括逻辑、数值、文本、函数柄和异质数据容器在内的15种数据类型，每一种类型都定义为矩阵或阵列的形式。

在线性代数、数理统计、自动控制、数字信号处理、动态系统仿真等方面，MATLAB已表现了其不可替代的一面。

是科研工作人员和工程技术人员进行科学研究和生产实践的有利武器。

　　MATLAB语言具有下述显着特点：

（1）具有强大的矩阵运算能力

MatrixLaboratory（矩阵实验室），使得矩阵运算非常简单。

MATLAB软件是基于矩阵计算开发的，在其他编程语言中需要使用多个for语句才能完成的操作，在MATLAB中直接使用矩阵即可完成计算，因而MATLAB在数据计算分析，特别是对海量数据的处理方面表现出相比其他的编程语言更大的优势。

（2）多样化的操作途径

MATLAB语言为用户提供了多种操作方式选择。

用户可以编写代码实现各种功能，代码可重复利用，同时，不擅长编程的用户也可以通过MATLAB图形界面操作，完成MATLAB的相应功能。

（3）功能强大的工具箱

对于算法的开发，MATLAB提供了大量现成的函数，用户可以直接调用。

MATLAB软件对目前主流的算法都提供了现成的函数，并封装成一个个适用于不同领域的算法工具箱。

常用的工具箱包括数学类、经济类、信号处理类工具箱，同时随着MATLAB软件版本的更新，不断有功能更强大的工具箱添加到MATLAB中。

MATLAB语言的这一特点，将大大节省算法开发的时间，用户无须详细了解算法的细节，只要掌握MATLAB中算法相应函数的调用即可。

（4）良好的扩展能力

利用MATLAB语言编写的程序具有良好的扩展能力，可以方便地与各种编程语言链接。

用户可以方便地在MATLAB中调用其他语言已编写好的程序，同时在其他语言中也可以方便地调用MATLAB的程序。

MATLAB语言具有良好的接口编程技术。

（5）完善的帮助系统

完善的帮助系统是MATLAB的又一突出特点，MATLAB向用户提供了多种帮助途径，在1.4节中将详细介绍MATLAB强大的帮助系统。

通过MATLAB的帮助系统，用户可以获取MATLAB常用函数的使用方法及应用实例，而且这种帮助可以是实时的、在线的。

同时，为了便于用户更好地使用MATLAB软件，在MATLAB中的主要算法都是可以直接看到源代码的。

2.2MATLAB的应用与功能特性

MATLAB具有以下几个功能特性:

（1）功能强大的数值运算能力：

MATLAB有超过500种数学、统计、科学及工程方面的函数可供调用，这些函数使用方法简单自然，允许用数学形式的语言编写程序。

另外，用户也可以加入自己的函数使系统成为使用者所需要的环境。

MATLAB编程效率高，易学易懂，因此，该语言被通俗的称为演算纸式科学算法语言。

（2）强大的图形处理能力：

在MATLAB中数据的可视化非常方便，可以很容易的制作高品质的图形。

用MATLAB绘图十分方便，它的系列绘图函数均只需调用不同的绘图函数。

在调用绘图函数时调整自变量可绘出不变颜色的、线、复线或多重线。

这种为科学研究着想的设计是一般通用的编程语言所不及的。

（3）高级且简单的程序环境：

既有结构化的控制语句，又有面向对象的编程特性。

用MATLAB编程十分简单，所花的时间约为C++的几分之一，而且不需要编译及链接即可执行，并且它的语法限制不严格，可移植性好。

（4）丰富的工具箱及模块集：

这些功能强劲的工具箱提供了使用者在许多特别应用领域所需的函数。

同时系统还包含很多种特殊函数，可将基于MATLAB的算法与外部应用程序和语言进行集成。

（5）易于扩充：

除内部函数外，所有MATLAB的核心文件和工具箱文件都是既可读又可改的源文件，用户能根据自己的需要对这些源文件进行修改或加入自己编写的文件。

运用MATLAB，将大量的复杂的熟悉运算用MATLAB进行形象、直观的计算机模拟与仿真实现。

随着版本的不断升级，内容不断扩充，功能更加强大，从而被广泛应用于仿真技术、自动控制和数字信号处理等领域。

全球现有超过五十万的企业用户和上千万的个人用户，广泛分布在航空航天、金融财务、机械化工、电信、教育等各个行业。

3数字滤波器简介

3.1数字滤波器的概念

滤波，即过滤之意，并不是对所有的干扰信号进行滤除，而是进行一定的衰减，让某些信号成分通过，又阻止或衰减了另一些信号成分，滤波器可广义地理解为一个信号选择系统。

它按冲激响应h（n）的长度分类可分为有限冲激响应（FIR）数字滤波器和无限冲激响应（IIR）滤波器两种，冲激响应本来是用于模拟系统，指系统对冲激函数δ（t）的响应。

发展到数字滤波器后，工程上仍沿用这个名称，与单位抽样响应和单位脉冲响应的说法通用。

在特性、结构、设计方法、运用场合等方面两种数字滤波器均不相同。

用系统函数表示为：

（式3.1）

由系统函数可以得到常系数线形差分程为：

（式3.2）

可见数字滤波器的功能就是把输入序列x（n）通过一定的运算变换成输出序列y（n）。

不同的运算处理方法决定了滤波器实现结构的不同。

无限冲激响应滤波器的单位抽样响应h（n）是无限长的，是递归式的，在结构上存在着输出信号到输入信号的反馈，因此在z平面的有限区间（0<︱z︱<∞）有极点存在。

3.2IIR数字滤波器的基本结构

IIR数字滤波器实现的基本结构有：

（1）IIR滤波器的直接型结构；

优点：

延迟线减少一半，变为N个，可节省寄存器或存储单元；

缺点：

其它缺点同直接I型。

通常在实际中很少采用上述两种结构实现高阶系统，而是把高阶变成一系列不同组合的低阶系统来实现。

（2）IIR滤波器的级联型结构；

优点：

系统实现简单，只需一个二阶节系统通过改变输入系数即可完成；极点位置单独调整；运算速度快（可并行进行）；各二阶网络的误差互不影响，总的误差小，字长要求低。

缺点：

不能直接调整零点，因多个二阶节的零点并不是整个系统函数的零点，当需要准确的传输零点时，级联型最合适。

（3）IIR滤波器的并联型结构。

优点：

简化实现，用一个二阶节，通过变换系数就可实现整个系统；极、零点可单独控制、调整，调整α1i、α2i只单独调整了第i对零点，调整β1i、β2i则单独调整了第i对极点；各二阶节零、极点的搭配可互换位置，优化组合以减小运算误差；可流水线操作。

缺点：

二阶阶电平难控制，电平大易导致溢出，电平小则使信噪比减小。

a、直接型b、并联型

c、串联型

图3.1IIR滤波器的基本结构

3.3FIR数字滤波器的基本结构

FIR滤波器的单位抽样响应为有限长度，一般采用非递归形式实现。

通常的FIR数字滤波器有横截性和级联型两种。

FIR滤波器实现的基本结构有：

（1）FIR滤波器的横截型结构

表示系统输入输出关系的差分方程为：

（式3.3）

   直接由差分方程得出的实现结构如图所示：

图3-2横截型（直接型﹑卷积型）

若h（n）呈现对称特性，即此FIR滤波器具有线性相位，则可以简化加横截型结构，下面分情况讨论：

图3-3、N为奇数时线形相位FIR滤波器实现结构图3-4、N为偶数时线性相位FIR滤波器实现结构

（2）FIR滤波器的级联型结构

将H（z）分解成实系数二阶因子的乘积形式：

（式3.4）

   这时FIR滤波器可用二阶节的级联结构来实现，每个二阶节用横截型结构实现。

如图所示：

图3.5FIR滤波器的级联结构

这种结构的每一节控制一对零点，因而在需要控制传输零点时可以采用这种结构。

3.4数字滤波器的性能指标

在进行滤波器设计时，需要确定其性能指标。

一般来说，滤波器的性能要求往往以频率响应的幅度特性的允许误差来表征。

以低通滤波器特性为例，频率响应有通带、过渡带及阻带三个范围。

在通带内：

 1-AP≤|H（ejω）|≤1|ω|≤ωc

在阻带中：

|H（ejω）|≤Astωst≤|ω|≤ωc

其中ωc为通带截止频率,ωst为阻带截止频率，Ap为通带误差,Ast为阻带误差。

图3.6低通滤波器频率响应幅度特性的容限图

与模拟滤波器类似，数字滤波器按频率特性划分为低通、