发布时间 : 星期六 文章基于Matlab的IIR数字滤波器设计(脉冲响应不变法)毕业论文更新完毕开始阅读
第一章 绪论
1.1信号数字现状与数字滤波器意义
当今，数字信号处理 (DSP：Digtal Signal Processing)技术正飞速发展，它不但自成一门学科，更是以不同形式影响和渗透到其他学科：它与国民经济息息相关，与国防建设紧密相连；它影响或改变着我们的生产、生活方式，因此受到人们普遍的关注。数字信号实际上是用数字序列表示的信号，语音信号经采样和量化后，得到的数字信号是一个一维离散时间序列；而图像信号经采样和量化后，得到的数字信号是一个二维离散空间序列。数字信号处理，就是用数值计算的方法对数字序列进行各种处理，把信号变换成符合需要的某种形式。例如，对数字信号经行滤波以限制他的频带或滤除噪音和干扰，或将他们与其他信号进行分离；对信号进行频谱分析或功率谱分析以了解信号的频谱组成，进而对信号进行识别；对信号进行某种变换，使之更适合于传输，存储和应用；对信号进行编码以达到数据压缩的目的，等等。数字化、智能化和网络化是当代信息技术发展的大趋势，而数字化是智能化和网络化的基础，实际生活中遇到的信号多种多样，例如广播信号、电视信号、雷达信号、通信信号、导航信号、射电天文信号、生物医学信号、控制信号、气象信号、地震勘探信号、机械振动信号、遥感遥测信号，等等。上述这些信号大部分是模拟信号，也有小部分是数字信号。模拟信号是自变量的连续函数，自变量可以是一维的，也可以是二维或多维的。大多数情况下一维模拟信号的自变量是时间，经过时间上的离散化(采样)和幅度上的离散化(量化)，这类模拟信号便成为一维数字信号。
在信号处理过程中，所处理的信号往往混有噪音，从接收到的信号中消除或减弱噪音是信号传输和处理中十分重要的问题。根据有用信号和噪音的不同特性，提取有用信号的过程称为滤波，实现滤波功能的系统称为滤波器。在近代电信设备和各类控制系统中，数字滤波器应用极为广泛。数字滤波器精确度高、使用灵活、可靠性高，具有模拟设备所没有的许多优点，已广泛地应用于各个科学技术领域, 例如数字电视、语音、通信、雷达、声纳、遥感、图像、生物医学以及许多工程应用领域。随着信息时代数字时代的到来，数字滤波技术已经成为一门极其重要的学科和技术领域。以往的滤波器大多采用模拟电路技术，但是，模拟电路技术存在很多难以解决的问题，例如，模拟电路元件对温度的敏感性，等等。而采用数字技术则避免很多类似的难题，当然数字滤波器在其他方面也有很多突出的优点，在前面部分已经提到，这些都是模拟技术所不能及的，所以采用数字滤波器对信号进行处理是目前的发展方向。
4 [1]
1.2 设计平台
MATLAB是美国MathWorks公司生产的一个为科学和工程计算专门设计的交互式大型软件，是一个可以完成各种精确计算和数据处理的、可视化的、强大的计算工具。它集图示和精确计算于一身，在应用数学、物理、化工、机电工程、医学、金融和其他需要进行复杂数值计算的领域得到了广泛应用。它不仅是一个在各类工程设计中便于使用的计算工具，而且也是一个在数学、数值分析和工程计算等课程教学中的优秀的数学工具，在世界各地的高和大型计算机上运行，适用于Windows、UNIX等多种系统平台。
MATLAB作为一种科学计算的高级语言之所以受欢迎，就是因为它有丰富的函数资源和工具箱资源，编程人员可以根据自己的需要选择函数，而无需再去编写大量繁琐的程序代码，从而减轻了编程人员的工作负担，被称为第四代编程语言的MATLAB最大的特点就是简洁开放的程序代码和直观实用的开发环境[1]。
1.3数字滤波器概述
数字滤波器可以用查分方程、单位取样响应以及系统函数等表示。对于研究系统的实现方法，即它的运算结构来说，用框图表示最为直接。
一个给定的输入输出关系，可以用多种不同的数字网络来实现。在不考虑量化影响时，这些不同的实现方法是等效的；但在考虑量化影响时，这些不同的实现方法性能上就有差异。因此，运算结构是很重要的，同一系统函数H(z)，运算结构的不同，将会影响系统的精度、误差、稳定性、经济性以及运算速度等许多重要性能。IIR(无限冲激响应)滤波器与FIR(有限冲激响应)滤波器在结构上有自己不同的特点，在设计时需综合考虑
作为线形时不变系统的数字滤波器可以用系统函数来表示，而实现一个系统函数表达式所表示的系统可以用两种方法：一种方法是采用计算机软件实现；另一种方法是用加法器、乘法器、和延迟器等元件设计出专用的数字硬件系统，即硬件实现。不论软件实现还是硬件实现，在滤波器设计过程中，由同一系统函数可以构成很多不同的运算结构。对于无限精度的系数和变量，不同结构可能是等效的，与其输入和输出特性无关；但是在系数和变量精度是有限的情况下，不同运算结构的性能就有很大的差异。因此，有必要对离散时间系统的结构有一基本认识。
5
第二章 IIR数字滤波器的设计
2.1 IIR滤波器的基本结构
一个数字滤波器可以用系统函数表示为：
(2-1)
由这样的系统函数可以得到表示系统输入与输出关系的常系数线形差分程为：
y(n)??aky(n?k)??bkx(n?k)
k?0k?0NM(2-2)
可见数字滤波器的功能就是把输入序列x(n)通过一定的运算变换成输出序列y(n)。不同的运算处理方法决定了滤波器实现结构的不同。无限冲激响应滤波器的单位抽样响应h(n)是无限长的，其差分方程如(2-2)式所示，是递归式的，即结构上存在着输出信号到输入信号的反馈，其系统函数具有(2-1)式的形式，因此在z平面的有限区间(0<︱z︱<∞)有极点存在。
前面已经说明，对于一个给定的线形时不变系统的系统函数，有着各种不同的等效差分方程或网络结构。由于乘法是一种耗时运算，而每个延迟单元都要有一个存储寄存器，因此采用最少常熟乘法器和最少延迟支路的网络结构是通常的选择，以便提高运算速度和减少存储器。然而，当需要考虑有限寄存器长度的影响时，往往也采用并非最少乘法器和延迟单元的结构。
IIR滤波器实现的基本结构有： (1)IIR滤波器的直接型结构；
优点：延迟线减少一半，变为N 个，可节省寄存器或存储单元； 缺点：其它缺点同直接I型。
通常在实际中很少采用上述两种结构实现高阶系统，而是把高阶变成一系列不同组合的低阶系统(一、二阶)来实现。 (2)IIR滤波器的级联型结构；
特点：
? 系统实现简单，只需一个二阶节系统通过改变输入系数即可完成； ? 极点位置可单独调整；
6
? 运算速度快(可并行进行)；
? 各二阶网络的误差互不影响，总的误差小，对字长要求低。
缺点：
不能直接调整零点，因多个二阶节的零点并不是整个系统函数的零点，当需要准确的传输零点时，级联型最合适。 (3)IIR滤波器的并联型结构。 优点：
? ?
简化实现，用一个二阶节，通过变换系数就可实现整个系统；
极、零点可单独控制、调整，调整α1i、α2i只单独调整了第i对零点，调整β1i、β2i则单独调整了第i对极点；
各二阶节零、极点的搭配可互换位置，优化组合以减小运算误差； 可流水线操作。
二阶阶电平难控制，电平大易导致溢出，电平小则使信噪比减小。
? ?
缺点：
?
a、直接型
7