基于Matlab的IIR数字滤波器设计(脉冲响应不变法)

第0页共44页课程设计说明书题目：基于Matlab的IIR数字滤波器设计第1页共44页课程设计（论文）任务书院（系）基层教学单位学号学生姓名专业（班级）设计题目基于Matlab的IIR数字滤波器设计(脉冲响应不变法)设计技术参数基于Matlab软件，利用脉冲响应不变法设计IIR低通、高通、带通、带阻数字滤波器，并分析滤波器单位脉冲响应，频率响应特性。设计要求1.设计报告中要求画出所设计的滤波器的网络结构图。2.设计出滤波器，并查看其波形。3.GUIDE的界面设计。并将其运行观察。工作量查找资料，学习使用MATLAB设计滤波器。分别设计出低通，高通，带通，帯阻四种数字滤波器。观察各滤波器的频率响应。设计GUIDE界面。工作计划第一天：查资料，够造思路。第二天到第三天：各基本程序的编写和调试。第四天：观察各滤波器的波形，并修改程序。第五天：设计GUIDE界面，并用其演示波形。第六天：撰写课程设计报告册。参考资料1谢平，林洪彬信号处理原理及应用北京：机械工程出版社2008.102刘波.MATLAB信号处理[M].北京：电子工业出版社，2006。3张圣勤MATLAB7.0实用教程北京：机械工程出版社2006年3月指导教师签字基层教学单位主任签字说明：此表一式四份，学生、指导教师、基层教学单位、系部各一份。课程设计第-2-页共44页摘要数字滤波是数字信号处理的重要内容，是由乘法器、加法器和单位延时器组成的一种运算过程，其功能是对输人离散信号进行运算处理，以达到改变信号频谱的目的。数字滤波器根据频域特性可分为低通、高通、带通和带阻四个基本类型。本文用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，利用了一种基于Matlab软件的数字滤波器设计方法，完成了低通，高通，带通，帯阻IIR滤波器的设计,文中深入分析了该滤波器系统设计的功能特点、实现原理以及技术关键，阐述了使用MATLAB进行带通滤波器设计及仿真的具体方法。最后把整个设计方案用GUIDE界面制作并演示出来。文章根据IIR滤波器的设计原理，重点介绍巴特沃斯数字滤波器的设计方法和操作步骤，并以实例形式列出设计程序。关键词：信号巴特沃斯MatlabIIR滤波器脉冲响应不变法课程设计第3页共44页目录摘要………………………………………………………………………………..3目录……………………………………………………………………………….....4第一章绪论…………………………………………………………………………..51.1信号数字现状与数字滤波器意义……………………………………………….51.2设计平台…………………………………………………………………………61.3数字滤波器概述……………………………………………………………….....6第二章IIR数字滤波器的设计……………………………………………………….72.1IIR滤波器的基本结构………………………………………………………..…..72.2滤波器的性能指标……………………………………………………………….102.3IIR数字滤波器的设计方法………………………………………………..……112.4巴特沃斯滤波器。……………………………………………………….……….13第三章IIR频率响应滤波器的实例………………………………………….……….153.2用脉冲响应不变法设计IIR低通数字滤波器实例………………………………153.2用脉冲响应不变法设计IIR高通数字滤波器实例………………………………173.3用脉冲响应不变法设计IIR带通数字滤波器实例……………………………...193.4用脉冲响应不变法设计IIR帯阻数字滤波器实例……………….……………...213.5（附）滤波信号的输入...............................................243.6滤波的效果...............................................................................................................24第四章界面设计…………………………………………………………………….254.1主界面………………………………………………………………………………254.2软件功能及使用方法……………………………………………………………..26总结..................................................................27程序清单………………………………………………………………………………29课程设计第4页共44页第一章绪论1.1信号数字现状与数字滤波器意义当今，数字信号处理[1](DSP：DigtalSignalProcessing)技术正飞速发展，它不但自成一门学科，更是以不同形式影响和渗透到其他学科：它与国民经济息息相关，与国防建设紧密相连；它影响或改变着我们的生产、生活方式，因此受到人们普遍的关注。数字信号实际上是用数字序列表示的信号，语音信号经采样和量化后，得到的数字信号是一个一维离散时间序列；而图像信号经采样和量化后，得到的数字信号是一个二维离散空间序列。数字信号处理，就是用数值计算的方法对数字序列进行各种处理，把信号变换成符合需要的某种形式。例如，对数字信号经行滤波以限制他的频带或滤除噪音和干扰，或将他们与其他信号进行分离；对信号进行频谱分析或功率谱分析以了解信号的频谱组成，进而对信号进行识别；对信号进行某种变换，使之更适合于传输，存储和应用；对信号进行编码以达到数据压缩的目的，等等。数字化、智能化和网络化是当代信息技术发展的大趋势，而数字化是智能化和网络化的基础，实际生活中遇到的信号多种多样，例如广播信号、电视信号、雷达信号、通信信号、导航信号、射电天文信号、生物医学信号、控制信号、气象信号、地震勘探信号、机械振动信号、遥感遥测信号，等等。上述这些信号大部分是模拟信号，也有小部分是数字信号。模拟信号是自变量的连续函数，自变量可以是一维的，也可以是二维或多维的。大多数情况下一维模拟信号的自变量是时间，经过时间上的离散化(采样)和幅度上的离散化(量化)，这类模拟信号便成为一维数字信号。在信号处理过程中，所处理的信号往往混有噪音，从接收到的信号中消除或减弱噪音是信号传输和处理中十分重要的问题。根据有用信号和噪音的不同特性，提取有用信号的过程称为滤波，实现滤波功能的系统称为滤波器。在近代电信设备和各类控制系统中，数字滤波器应用极为广泛。数字滤波器精确度高、使用灵活、可靠性高，具有模拟设备所没有的许多优点，已广泛地应用于各个科学技术领域,例如数字电视、语音、通信、雷达、声纳、遥感、图像、生物医学以及许多工程应用领域。随着信息时代数字时代的到来，数字滤波技术已经成为一门极其重要的学科和技术领域。以往的滤波器大多采用模拟电路技术，但是，模拟电路技术存在很多难以解决的问题，例如，模拟电路元件对温度的敏感性，等等。课程设计第5页共44页而采用数字技术则避免很多类似的难题，当然数字滤波器在其他方面也有很多突出的优点，在前面部分已经提到，这些都是模拟技术所不能及的，所以采用数字滤波器对信号进行处理是目前的发展方向。1.2设计平台MATLAB是美国MathWorks公司生产的一个为科学和工程计算专门设计的交互式大型软件，是一个可以完成各种精确计算和数据处理的、可视化的、强大的计算工具。它集图示和精确计算于一身，在应用数学、物理、化工、机电工程、医学、金融和其他需要进行复杂数值计算的领域得到了广泛应用。它不仅是一个在各类工程设计中便于使用的计算工具，而且也是一个在数学、数值分析和工程计算等课程教学中的优秀的数学工具，在世界各地的高和大型计算机上运行，适用于Windows、UNIX等多种系统平台。MATLAB作为一种科学计算的高级语言之所以受欢迎，就是因为它有丰富的函数资源和工具箱资源，编程人员可以根据自己的需要选择函数，而无需再去编写大量繁琐的程序代码，从而减轻了编程人员的工作负担，被称为第四代编程语言的MATLAB最大的特点就是简洁开放的程序代码和直观实用的开发环境[1]。1.3数字滤波器概述数字滤波器可以用查分方程、单位取样响应以及系统函数等表示。对于研究系统的实现方法，即它的运算结构来说，用框图表示最为直接。一个给定的输入输出关系，可以用多种不同的数字网络来实现。在不考虑量化影响时，这些不同的实现方法是等效的；但在考虑量化影响时，这些不同的实现方法性能上就有差异。因此，运算结构是很重要的，同一系统函数H(z)，运算结构的不同，将会影响系统的精度、误差、稳定性、经济性以及运算速度等许多重要性能。IIR(无限冲激响应)滤波器与FIR(有限冲激响应)滤波器在结构上有自己不同的特点，在设计时需综合考虑作为线形时不变系统的数字滤波器可以用系统函数来表示，而实现一个系统函数表达式课程设计第6页共44页所表示的系统可以用两种方法：一种方法是采用计算机软件实现；另一种方法是用加法器、乘法器、和延迟器等元件设计出专用的数字硬件系统，即硬件实现。不论软件实现还是硬件实现，在滤波器设计过程中，由同一系统函数可以构成很多不同的运算结构。对于无限精度的系数和变量，不同结构可能是等效的，与其输入和输出特性无关；但是在系数和变量精度是有限的情况下，不同运算结构的性能就有很大的差异。因此，有必要对离散时间系统的结构有一基本认识。第二章IIR数字滤波器的设计2.1IIR滤波器的基本结构一个数字滤波器可以用系统函数表示为：01()()()1MkkkNkkkbzYzHzXzaz(2-1)由这样的系统函数可以得到表示系统输入与输出关系的常系数线形差分程为：00()()()NMkkkkynaynkbxnk(2-2)可见数字滤波器的功能就是把输入序列x(n)通过一定的运算变换成输出序列y(n)。不同的运算处理方法决定了滤波器实现结构的不同。无限冲激响应滤波器的单位抽样响应h(n)是无限长的，其差分方程如(2-2)式所示，是递归式的，即结构上存在着输出信号到输入信号的反馈，其系统函数具有(2-1)式的形式，因此在z平面的有限区间(0︱z︱∞)有极点存在。前面已经说明，对于一个给定的线形时不变系统的系统函数，有着各种不同的等效差分方程或网络结构。由于乘法是一种耗时运算，而每个延迟单元都要有一个存储寄存器，因此采用最少常熟乘法器和最少延迟支路的网络结构是通常的选择，以便提高运算速度和减少存储器。然而，当需要考虑有限寄存器长度的影响时，往往也采用并非最少乘法器和延课程设计第7页共44页迟单元的结构。IIR滤波器实现的基本结构有：(1)IIR滤波器的直接型结构；优点：延迟线减少一半，变为N个，可节省寄存器或存储单元；缺点：其它缺点同直接I型。通常在实际中很少采用上述两种结构实现高阶系统，而是把高阶变成一系列不同组合的低阶系统(一、二阶)来实现。(2)IIR滤波器的级联型结构；特点：系统实现简单，只需一个二阶节系统通过改变输入系数即可完成；极点位置可单独调整；运算速度快(可并行进行)；各二阶网络的误差互不影响，总的误差小，对字长要求低。缺点：不能直接调整零点，因多个二阶节的零点并不是整个系统函数的零点，当需要准确的传输零点时，级联型最合适。(3)IIR滤波器的并联型结构。优点：简化实现，用一个二阶节，通过变换系数就可实现整个系统；极、零点可单独控制、调整，调整α1i、α2i只单独调整了第i对零点，调整β1i、β2i则单独调整了第i对极点；各二阶节零、极点的搭配可互换位置，优化组合以减小运算误差；可流水线操作。缺点：二阶阶电平难控制，电平大易导致溢出，电平小则使信噪比减小。课程设计第8页共44页a、直接型b、并联型课程设计第9页共44页c、串联型图12.2滤波器的性能指标我们在进行滤波器设计时，需要确定其性能指标。一般来说，滤波器的性能要求往往以频率响应的幅度特性