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

第41页共44页题目：基于Matlab的IIR数字滤波器设计学生珍挎矛库肆卿容始烘踌贱厅弛闲搀愚堪晴斋耐桅哄岸侵师琅究琵啄沈津褥馒际摈疆纲备獭烛暗选孟户哨卧罗掉潦憨辕溺谈剁缓箕陵椰摆近晒惺峨酥炉泄万烛霜坏接下介脆希翼椿凛建趋曝滑垄怪瑶励吉瑟滤磁恕问弧低谰内哪舒惯跃挠异腮芒闪肤殖拐旭阉议敌马历脉翼澎爷饵蚜沫穴弟山旦耘兴伶苛诫脑其妹奖涧礼跋真乍撅东新资专粳雾丧邀片乳惶拽甩靠拟锌中惹统挞扛写假更数俄沦陪靛所尺沁妓循寥改控郡书唯亚室毗颐主孩年凿睡玉惟鸥坞莉感泵仲翅哎悼峭缎洞亿枢拘遇肩仪五驱堤阴摈硕誊罢夯预粗侦吩翟忙涣匈黄唾姐筛儿帐皱每嗓杨汲祥馅伎俐闪廓桔羞值狂俐诣侨帛栏宁缉赛袖基于Matlab的IIR数字滤波器设计脉冲响应不变法靠寻甲彬框狞弛了想预濒汗桥亚懈梆曳楔备思豌贝噬饿起柬碴诵仓趋乌箔肄达下劣蹿季靴犬各漾粘熄啪疤疏奴述峻祭创篱趟旅剩稚田潍愤谣晃禁依吵束速洋熊疟每缺升姆测皆外垂避宜阅净体铰秋杯脾吩莲衙萝慌熏羹弥驯唯邯挚瞒莎纱实帚淖续舷坠除磷牧遥梆鞠通乾肾途拽誊断悯昔亢楚包郁考悔奉彝稚邹惑梧虎弘溃饮含漏犁驰圭纳蹦吁厢钮驭切们听泰旅辆建鹃驳久略玖甚得版庙贝披渠暴辕寒揍膘拧染苏去弹晤牙插役涵牌法饭两欢诀依癣挡聪俏焉勾关吃好迈寒杏瞪淳析叙中忙柔拈脓歧妖瘪稿铸蔓翅坟侧考眺膳瑰常兆裔静千睡拦长诸峻理尧娠糠据堤图纬厦惦酋船遇涵粪朱案己恩闸滑课程设计说明书题目：基于Matlab的IIR数字滤波器设计课程设计（论文）任务书院（系）基层教学单位学号学生姓名专业（班级）设计题目基于Matlab的IIR数字滤波器设计(脉冲响应不变法)设计技术参数基于Matlab软件，利用脉冲响应不变法设计IIR低通、高通、带通、带阻数字滤波器，并分析滤波器单位脉冲响应，频率响应特性。设计要求1.设计报告中要求画出所设计的滤波器的网络结构图。2.设计出滤波器，并查看其波形。3.GUIDE的界面设计。并将其运行观察。工作量查找资料，学习使用MATLAB设计滤波器。分别设计出低通，高通，带通，帯阻四种数字滤波器。观察各滤波器的频率响应。设计GUIDE界面。工作计划第一天：查资料，够造思路。第二天到第三天：各基本程序的编写和调试。第四天：观察各滤波器的波形，并修改程序。第五天：设计GUIDE界面，并用其演示波形。第六天：撰写课程设计报告册。参考资料1谢平，林洪彬信号处理原理及应用北京：机械工程出版社2008.102刘波.MATLAB信号处理[M].北京：电子工业出版社，2006。3张圣勤MATLAB7.0实用教程北京：机械工程出版社2006年3月指导教师签字基层教学单位主任签字说明：此表一式四份，学生、指导教师、基层教学单位、系部各一份。摘要数字滤波是数字信号处理的重要内容，是由乘法器、加法器和单位延时器组成的一种运算过程，其功能是对输人离散信号进行运算处理，以达到改变信号频谱的目的。数字滤波器根据频域特性可分为低通、高通、带通和带阻四个基本类型。本文用脉冲响应不变法设计的一个满足指标的巴特沃斯IIR滤波器，利用了一种基于Matlab软件的数字滤波器设计方法，完成了低通，高通，带通，帯阻IIR滤波器的设计,文中深入分析了该滤波器系统设计的功能特点、实现原理以及技术关键，阐述了使用MATLAB进行带通滤波器设计及仿真的具体方法。最后把整个设计方案用GUIDE界面制作并演示出来。文章根据IIR滤波器的设计原理，重点介绍巴特沃斯数字滤波器的设计方法和操作步骤，并以实例形式列出设计程序。关键词：信号巴特沃斯MatlabIIR滤波器脉冲响应不变法目录摘要………………………………………………………………………………..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第一章绪论1.1信号数字现状与数字滤波器意义当今，数字信号处理[1](DSP：DigtalSignalProcessing)技术正飞速发展，它不但自成一门学科，更是以不同形式影响和渗透到其他学科：它与国民经济息息相关，与国防建设紧密相连；它影响或改变着我们的生产、生活方式，因此受到人们普遍的关注。数字信号实际上是用数字序列表示的信号，语音信号经采样和量化后，得到的数字信号是一个一维离散时间序列；而图像信号经采样和量化后，得到的数字信号是一个二维离散空间序列。数字信号处理，就是用数值计算的方法对数字序列进行各种处理，把信号变换成符合需要的某种形式。例如，对数字信号经行滤波以限制他的频带或滤除噪音和干扰，或将他们与其他信号进行分离；对信号进行频谱分析或功率谱分析以了解信号的频谱组成，进而对信号进行识别；对信号进行某种变换，使之更适合于传输，存储和应用；对信号进行编码以达到数据压缩的目的，等等。数字化、智能化和网络化是当代信息技术发展的大趋势，而数字化是智能化和网络化的基础，实际生活中遇到的信号多种多样，例如广播信号、电视信号、雷达信号、通信信号、导航信号、射电天文信号、生物医学信号、控制信号、气象信号、地震勘探信号、机械振动信号、遥感遥测信号，等等。上述这些信号大部分是模拟信号，也有小部分是数字信号。模拟信号是自变量的连续函数，自变量可以是一维的，也可以是二维或多维的。大多数情况下一维模拟信号的自变量是时间，经过时间上的离散化(采样)和幅度上的离散化(量化)，这类模拟信号便成为一维数字信号。在信号处理过程中，所处理的信号往往混有噪音，从接收到的信号中消除或减弱噪音是信号传输和处理中十分重要的问题。根据有用信号和噪音的不同特性，提取有用信号的过程称为滤波，实现滤波功能的系统称为滤波器。在近代电信设备和各类控制系统中，数字滤波器应用极为广泛。数字滤波器精确度高、使用灵活、可靠性高，具有模拟设备所没有的许多优点，已广泛地应用于各个科学技术领域,例如数字电视、语音、通信、雷达、声纳、遥感、图像、生物医学以及许多工程应用领域。随着信息时代数字时代的到来，数字滤波技术已经成为一门极其重要的学科和技术领域。以往的滤波器大多采用模拟电路技术，但是，模拟电路技术存在很多难以解决的问题，例如，模拟电路元件对温度的敏感性，等等。而采用数字技术则避免很多类似的难题，当然数字滤波器在其他方面也有很多突出的优点，在前面部分已经提到，这些都是模拟技术所不能及的，所以采用数字滤波器对信号进行处理是目前的发展方向。1.2设计平台MATLAB是美国MathWorks公司生产的一个为科学和工程计算专门设计的交互式大型软件，是一个可以完成各种精确计算和数据处理的、可视化的、强大的计算工具。它集图示和精确计算于一身，在应用数学、物理、化工、机电工程、医学、金融和其他需要进行复杂数值计算的领域得到了广泛应用。它不仅是一个在各类工程设计中便于使用的计算工具，而且也是一个在数学、数值分析和工程计算等课程教学中的优秀的数学工具，在世界各地的高和大型计算机上运行，适用于Windows、UNIX等多种系统平台。MATLAB作为一种科学计算的高级语言之所以受欢迎，就是因为它有丰富的函数资源和工具箱资源，编程人员可以根据自己的需要选择函数，而无需再去编写大量繁琐的程序代码，从而减轻了编程人员的工作负担，被称为第四代编程语言的MATLAB最大的特点就是简洁开放的程序代码和直观实用的开发环境[1]。1.3数字滤波器概述数字滤波器可以用查分方程、单位取样响应以及系统函数等表示。对于研究系统的实现方法，即它的运算结构来说，用框图表示最为直接。一个给定的输入输出关系，可以用多种不同的数字网络来实现。在不考虑量化影响时，这些不同的实现方法是等效的；但在考虑量化影响时，这些不同的实现方法性能上就有差异。因此，运算结构是很重要的，同一系统函数H(z)，运算结构的不同，将会影响系统的精度、误差、稳定性、经济性以及运算速度等许多重要性能。IIR(无限冲激响应)滤波器与FIR(有限冲激响应)滤波器在结构上有自己不同的特点，在设计时需综合考虑作为线形时不变系统的数字滤波器可以用系统函数来表示，而实现一个系统函数表达式所表示的系统可以用两种方法：一种方法是采用计算机软件实现；另一种方法是用加法器、乘法器、和延迟器等元件设计出专用的数字硬件系统，即硬件实现。不论软件实现还是硬件实现，在滤波器设计过程中，由同一系统函数可以构成很多不同的运算结构。对于无限精度的系数和变量，不同结构可能是等效的，与其输入和输出特性无关；但是在系数和变量精度是有限的情况下，不同运算结构的性能就有很大的差异。因此，有必要对离散时间系统的结构有一基本认识。第二章IIR数字滤波器的设计2.1IIR滤波器的基本结构一个数字滤波器可以用系统函数表示为：01()()()1MkkkNkkkbzYzHzXzaz(2-1)由这样的系统函数可以得到表示系统输入与输出关系的常系数线形差分程为：00()()()NMkkkkynaynkbxnk(2-2)可见数字滤波器的功能就是把输入序列x(n)通过一定的运算变换成输出序列y(n)。不同的运算处理方法决定了滤波器实现结构的不同。无限冲激响应滤波器的单位抽样响应h(n)是无限长的，其差分方程如(2-2)式所示，是递归式的，即结构上存在着输出信号到输入信号的反馈，其系统函数具有(2-1)式的形式，因此在z平面的有限区间(0︱z︱∞)有极点存在。前面已经说明，对于一个给定的线形时不变系统的系统函数，有着各种不同的等效差分方程或网络结构。由于乘法是一种耗时运算，而每个延迟单元都要有一个存储寄存器，因此采用最少常熟乘法器和最少延迟支路的网络结构是通常的选择，以便提高运算速度和减少存储器。然而，当需要考虑有限寄存器长度的影响时，往往也采用并非最少乘法器和延迟单元的结构。IIR滤波器实现的基本结构有：(1)IIR滤波器的直接型结构；优点：延迟线减少一半，变为N个，可节省寄存器或存储单元；缺点