IIR数字低通滤波器的设计

IIR数字低通滤波器的 设计 领导形象设计圆作业设计ao工艺污水处理厂设计附属工程施工组织设计清扫机器人结构设计 《数字信号处理课程设计 报告 软件系统测试报告下载sgs报告如何下载关于路面塌陷情况报告535n,sgs报告怎么下载竣工报告下载 》 题 目: IIR数字低通滤波器的设计 学 院: 专 业: 班 级: 姓 名: 指导教师: 2012年 6月24日 目 录 1.课程设计的目的及要求 ................................... 1 1.1课程设计的目的 .................................... 1 1.2课程设计的要求 .................................... 1 2.设计的理论基础 ........................................ 1 2.1数字滤波器的工作原理 .............................. 1 2.2 数字滤波器的基本特性 .............................. 2 2.3 数字滤波器的基本结构 .............................. 4 2.4 IIR和FIR滤波器的区别 ............................ 4 3.MATLAB概述 ............................................ 5 4.IIR数字低通滤波器的设计 ............................... 6 4.1 IIR数字滤波器的设计 方法 快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载 .......................... 6 4.2 IIR滤波器经典设计 ................................ 7 4.3 MATLAB仿真步骤 .................................. 11 4.4 程序清单 ........................................ 12 5. 总结 初级经济法重点总结下载党员个人总结TXt高中句型全总结.doc高中句型全总结.doc理论力学知识点总结pdf ................................................. 13 参考文献: ............................................. 14 1.课程设计的目的及要求 1.1课程设计的目的 (1)学会MATLAB的使用，掌握MATLAB的程序设计方法; (2)掌握数字信号处理的基本概念、基本理论和基本方法; (3)掌握MATLAB设计IIR滤波器; (4)学会用MATLAB对信号进行分析处理。 1.2课程设计的要求 滤波器的主要技术指标取决于具体的应用或相互间的相互关系。具体的有最 ,,p,大通带增益(即通带允许起伏);最大阻带增益;通带截止频率;阻带 ,s截止频率。 ,，2krad/sc本设计中要求的技术指标是:通带截止频率=2π，阻带截 R,,2,，4krad/sRpss止频率，通带衰减小于3dB，阻带衰减大于15dB,采 F,20000Hzs样频率。 2.设计的理论基础 2.1数字滤波器的工作原理 在数字滤波中，我们主要讨论离散时间序列。如图2.1所示。设输入序列为，，，，，，，，xn,nhn,n，离散或数字滤波器对单位抽样序列的响应为。因在时域离散信号和系统中所起的作用相当于单位冲激函数在时域连续信号和系统中所起的作用。 图2.1 数字滤波器原理 ，，yn数字滤波器的序列将是这两个序列的离散卷积，即 , ，，，，，，yn,hkxn,k ,k,, 同样，两个序列卷积的z变换等于个自z变换的乘积，即 1 ，，，，，， Yz,HzXz jT,用代入上式，其中T为抽样周期，则得到 z,e ,jTj,Tj,TY，，，，，，e,HeXe j,Tj,TX，，eY，，e式中和 分别为数字滤波器输入序列和输出序列的频谱，而 j,Tj,TH，，eX，，e为单位抽样序列响应，，的频谱。由此可见，输入序列的频谱经hn j,Tj,Tj,TH，，，，eXeX，，e过滤波后，变为 ，按照的特点和我们处理信号的目的， j,Tj,Tj,TH，，eH，，，，eXe选取适当的使的滤波后的符合我们的要求。 2.2 数字滤波器的基本特性 如果利用离散时间系统对数字信号(时间离散、幅度量化的信号)进行滤波则构成数字滤波器。 数字滤波器的差分方程表示为: NM y(n),by(n,i)，ax(n,k) ,,,,10kk ik系统函数表示: Mk,az,k Yz()k,0Hz(),,NXz()i,1,bz,ii,1 j,H(e)数字滤波器的特性通常用其频率响应函数来描述，包括幅度特性 j,j,H(e)arg(H(e))和相位特性。 按信号通过系统时的特性(主要是幅频特性)来分类:可以有低通、高通、带通和带阻四种基本类型。 (1)低通数字滤波器:图2.2所示 j,,,,H(e) ,,jc,H(e), ,0 ,,,,,,c, 2 ωj|H(e)| 2π-2π πω-π c -ω c π f/2 f ff cs-f s -f-f/2 css 图2.2 低通数字滤波器的频谱 (2)高通数字滤波器:图2.3所示 ,j,,,,H(e) ,,,jc,H(e),,0 ,,,,c, ωj|H(e)| -ππ-ω-2π ω 2πc c 图2.3 高通数字滤波器的频谱 (3)带通数字滤波器:图2.4所示 j,,,,,H(e) ,,,2j1,H(e),, , ,,,,0,,,,||,0,12, ωj|H(e)| -ω ωω2-π-ωπ-2π2 1 1 2π 图2.4 带通数字滤波器的频谱 (4)带阻数字滤波器:图2.5所示 j,,,,,,,H(e) , ,,,||,0,j12,H(e),,0 ,,,,,,12, 3 ωj|H(e)| 2π-π-ω ω2-2ππ 2 ω-ω1 1 图2.5 带阻数字滤波器的频谱 2.3 数字滤波器的基本结构 (1)直接?型 (2)直接?型 优点:延迟器减少一半，变为N个，可以节省存储器或存储单元。 缺点:其他缺点同直接?型 (3)IIR滤波器的级联型结构 优点:系统实现简单，只需一个二阶系统通过改变输入系数即可完成，极点位置可以单独调整，运算速度快，各二阶网络的误差互不影响，总的误差小，对字长要求低。 缺点:不能直接调整零点，因多个二阶的零点并不是整个系统函数的零点，当需要准确的传输零点时，级联型最合适。 (4)IIR滤波器的并联型结构 优点:简单实现，用一个二阶的，通过变换系数就可实现整个系统。零极点可单独控制，调整。各二阶零极点的搭配可互换位置，优化组合以减小运算误差，可流水线操作。 缺点:二阶电平难控制，电平大，容易导致溢出，电平小则是信噪比减小。 2.4 IIR和FIR滤波器的区别 按系统冲击响应(或差分方程)可以分成无限冲击响应 IIR和有限冲击响应FIR滤波器两类。这两种滤波器都可以现实各种频率特性要求，但它们在计算流程、具体特性逼近等方面是有差别的。 (1) FIR滤波器(非递归型): 4 N,1 y(n),h(m)x(n,m) ,m,0 N,1,nH(Z),h(n)Z ,n,0 (2) IIR滤波器(递归型) NM y(n),aky(n,k)，bkx(n,k) ,,,,10kk M,kbkz,Y(z)k,0H(z),, NX(z)k,1,akz,k,1 3.MATLAB概述 MATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、 工程 路基工程安全技术交底工程项目施工成本控制工程量增项单年度零星工程技术标正投影法基本原理 设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言(如C、Fortran)的编辑模式，代表了当今国际科学计算软件的先进水平。 MATLAB 产品族可以用来进行以下各种工作:数值分析，数值和符号计算，工程与科学绘图，控制系统的设计与仿真，通讯系统设计与仿真，财务与金融工程， (1)友好的工作平台和编程环境 MATLAB由一系列工具组成。这些工具方便用户使用MATLAB的函数和文件，其中许多工具采用的是图形用户界面。包括MATLAB桌面和命令窗口、历史命令窗口、编辑器和调试器、路径搜索和用于用户浏览帮助、工作空间、文件的浏览器。随着MATLAB的商业化以及软件本身的不断升级，MATLAB的用户界面也越来越精致，更加接近Windows的标准界面，人机交互性更强，操作更简单。 (2)简单易用的程序语言 Matlab一个高级的矩阵/阵列语言，它包含控制语句、函数、数据结构、输入和输出和面向对象编程特点。用户可以在命令窗口中将输入语句与执行命令同 5 步，也可以先编写好一个较大的复杂的应用程序(M文件)后再一起运行。 (3)强大的科学计算机数据处理能力 MATLAB是一个包含大量计算算法的集合。其拥有600多个工程中要用到的数学运算函数，可以方便的实现用户所需的各种计算功能。函数中所使用的算法都是科研和工程计算中的最新研究成果，而前经过了各种优化和容错处理。在通常情况下，可以用它来代替底层编程语言，如C和C++ 。在计算要求相同的情况下，使用MATLAB的编程工作量会大大减少。 (4)出色的图形处理功能 MATLAB自产生之日起就具有方便的数据可视化功能，以将向量和矩阵用图形表现出来，并且可以对图形进行标注和打印。高层次的作图包括二维和三维的可视化、图象处理、动画和表达式作图。可用于科学计算和工程绘图。 (5)应用广泛的模块集合工具箱 MATLAB对许多专门的领域都开发了功能强大的模块集和工具箱。一般来说，它们都是由特定领域的专家开发的，用户可以直接使用工具箱学习、应用和评估不同的方法而不需要自己编写代码。 (6)实用的程序接口和发布平台 新版本的MATLAB可以利用MATLAB编译器和C/C++数学库和图形库，将自己的MATLAB程序自动转换为独立于MATLAB运行的C和C++代码。允许用户编写可以和MATLAB进行交互的C或C++语言程序。 (7)应用软件开发(包括用户界面) 在开发环境中，使用户更方便地控制多个文件和图形窗口;在编程方面支持了函数嵌套，有条件中断等;在图形化方面，有了更强大的图形标注和处理功能，包括对性对起连接注释等。 4.IIR数字低通滤波器的设计 4.1 IIR数字滤波器的设计方法 IIR滤波器是一种数字滤波器，滤波器的系统函数如式 M,rbz,r,Y(z),n,0rH(z),,h(n)z, ,NX(z),k,0n1，az,k,1k 6 由于它的脉冲响应序列是无限长的，故称无限冲激响应滤波器。IIR滤h(n) 波器的设计就是根据滤波器某些性能指标要求，设计滤波器的分子和分母多项式。它和FIR滤波器相比优点是在满足相同性能指标要求条件下，IIR滤波器的阶数要明显低于FIR滤波器。但IIR滤波器的相位是非线性的。 IIR滤波器设计方法可分为三种:模拟滤波器变换(经典设计法)、直接设计法、参数模型设计法、最大平滑滤波器设计。本设计选用的是模拟滤波器变化法(经典设计法)。 4.2 IIR滤波器经典设计 基于模拟滤波器变换原理，首先是根据滤波器的技术指标设计出相应的模拟滤波器，然后再将设计好的模拟滤波器变换成满足给定技术指标的数字滤波器。这就是IIR滤波器设计的经典法。在MATLAB中，经典法设计IIR数字滤波器采用下面的主要步骤: IIR滤波器 模拟离散化 频率变换 模拟低通滤波原型 图4.1 IIR经典法设计步骤 第二步完成后，一个达到期望性能指标的模拟滤波器(低通、高通、带通和带阻)已经设计出来。第三步离散化主要任务就是把模拟滤波器变换成数字滤波器，即把模拟滤波器的系数映射成数字滤波器的系统函数。 H(S)H(z) 实现系统传递函数s域至z域映射有脉冲响应不变法和双线性映射两种方法。 (1)脉冲响应不变法。 设计过程是:按照技术要求设计一个模拟低通滤波器，得到模拟低通滤波器的传 ，，Hs输函数转换成数字低通滤波器的系统函数H(z)。 a ，，Hs 将传输函数从s平面传换到z平面的方法有多种，但工程上常用的是脉a 冲响应不变法和双线性变换法。我们先研究脉冲响应不变法。 ，，，，Hsht设模拟滤波器的传输函数为，相应的单位冲激响应是，aa ，，，，，，Hshtht=LT[]，LT[.]代表拉氏变换，对进行等间隔采样，采样间隔为T，aaa 7 ，，，，hnThnT得到，将h(n)= 作为数字滤波器的单位取样响应，那么数字滤波aa 器的系统函数H(z)便是h(n)的Z变换。因此脉冲响应不变法是一种时域上的转 ，，ht在采样点上等于。 换方法，它是h(n)a ，，Hs设模拟滤波器只有单阶极点，且分母多项式的阶次高于分子多项式的a ，，Hs阶次，将用部分分式表示: a NAis，，Hs，，Hs, ，式中为的单阶极点。 ,iaas,si,1i ，，，，Hsht将逆拉氏变换得到: aa Nsti，，ht,Aeu(t) ，式中u(t)是单位阶跃函数。 ,ai,1i ，，ht对进行等间隔采样，采样间隔为T,得到: a NsnTi，，hn,h(nT),Aeu(nT) ,ai,1i 对上式进行Z变换，得到数字滤波器的系统函数H(z): NAsTii，，Hz,，，sHsA ，的极点映射到z平面，其极点变成，系数不e,iaisT,1i1,ez,1i 变化。下面分析从模拟滤波器转换到数字滤波器，s平面和z平面之间的映射关 ，，ht系，从而找到这种转换方法的优缺点。这里我们以采样信号作为桥梁，推a导其映射关系。 ˆ，，ht，，设的采样信号用ht表示， aa ,ˆ，，ht,h(t),(t,nT) ,aan,,, ˆ，，对ht进行拉氏变换，得到: a ,,stˆˆHs,h(t)edt，，aa,,, ,,st,[h(t)(t,nT)]edt , a,,,,n ,snT,h(nT)ea,n 8 ，，h(nT)ht式中是在采样点t=n T时的幅度值，它与序列h(n)的幅度值相等，aa h(nT)它与序列h(n)的幅度值相等，即h(n)= ，因此得到: a ,snT,nˆ ，，Hs,h(n)e,h(n)z,H(z)sT,,az,esTnnz,e 上式表示采样信号的拉氏变换与相应的序列的Z变换之间的映射关系可用下式 sT表示: z,e ˆ，，htH(j,)我们知道模拟信号的傅立叶变换和其采样信号，，的傅立叶htaaa ˆˆˆ变换和其采样信号，，的傅立叶变换之间的关系满足式: htH(j,)H(j,)aaa ,1ˆH(j,),H(j,,jk,) ,aasTk,,, 将代入上式，得 s,j, 1ˆH(j,),H(s,jk,) ,aasTk 1H(z),H(s,jk,)由上式得到: sT,asz,eTk ，，ht,,2,/T上式表明将模拟信号的拉氏变换在s平面上沿虚轴按照周期延as迟后，再按照映射关系，映射到z平面上，就得到H(z)。 (2)双线性变换法 这种变换方法，采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到 sT?π/T之间，再用转换到z平面上。设Ha(s)，s=jΩ，经过非线性频率z,e s压缩后用Ha(s1)，=jΩ1表示，这里用正切变换实现频率压缩: 1 2 ，，,,tan0.5,T1T 式中T仍是采样间隔，当Ω1从-π/T经过0变化到π/T时，Ω则由-?经过0 s变化到+?，实现了s平面上整个虚轴完全压缩到平面上虚轴的?π/T之间的1 转换。这样便有 9 ,sT1221,zs,th，，0.5,T, 1,st1TT1，z sT再通过 转换到z平面上，得到: z,e ,121,zs, ,1T1，z 2，sT z,2,sT 上式称为双线性变换。从s平面映射到s平面，再从s平面映射到z平面，其映11 射情况如图4.2所示。由于从s平面到s平面具有非线性频率压缩的功能，因此1 s不可能产生频率混叠现象。另外，从平面转换到z平面仍然采用标准转换关系 1 sT1，s平面的 ?π/T之间水平带的左半部分映射z平面单位圆内部，虚轴z,e1 映射单位圆。这样，Ha(s)因果稳定，转换成的H(z)也是因果稳定的。下面分析模拟频率Ω和数字频率ω之间的关系。 j,s,j,,z,e令 ，有 ,,j2121,e,,tan ,j,,,j,T2T1，e 图4.2 双线性变化法映射 上式说明，s平面上Ω与平面上的ω成非线性正切关系，如图4.3所示。在ω=0附近接近线性关系;当ω增加是，Ω增加得愈来愈快;当ω趋近π时，Ω趋近于?。正是因为这种非线性关系，消除了频率混叠现象。 10 图4.3 双线性变换法的频率变换 ω与Ω之间的非线性关系是双线性变换法的缺点，直接影响数字滤波器频响逼真的模仿模拟滤波器的频响，幅度特性和相位特性失真的情况如图4.4所示。这种非线性影响的实质问题是:如果Ω的刻度是均匀的，则影射到z平面ω的刻度不是均匀的，而是随ω增加愈来愈密。 图4.4 双线性变换法幅度和相位特性的非线性映射 双线性变换法可由简单的代数公式将Ha(s)直接转换成H(z)，这是该变换法的优点。 4.3 MATLAB仿真步骤 (1)根据给定的性能指标和方法不同，首先对设计性能指标中的频率指标进行转换，转换后的频率指标作为模拟滤波器原型设计性能指标。 (2)估计模拟低通滤波器最小阶数和边界频率，利用MATLAB工具函数buttord、cheb1ord等。 (3)设计模拟低通滤波器原型，利用MATLAB工具函数buttap、cheb1ap等。 (4)由模拟低通原型经频率变换获得模拟滤波器(低通、高通、带通、带阻)，利用MATLAB工具函数lp2lp、lp2hp、lp2bp、lp2bs。 11 (5)将模拟滤波器离散化获得IIR数字滤波器，利用MATLAB工具函数 bilinear。 4.4 程序清单 ,数字性能指标给出通带截止频率，阻带截止频率，通带衰减，阻带,Rspp R,衰减等。数字频率和的取值范围为0,π，单位:弧度，而MATLAB工,ssp ,具函数常采用标准化频率，和的取值范围为0,1。 ,sp ,模拟性能指标给出通带截止频率，阻带起始频率，通带衰减，阻,Rspp R,,带衰减等。模拟频率和的单位均为弧度/秒。 scs ,，2krad/s技术指标:通带截止频率=2π，阻带截止频率c ,,2,，4krad/sR，通带衰减小于3dB，阻带衰减大于15dB,采样频率Rssp F,20000Hz。 s 程序如下: Wp= 2000*2*pi; Ws=3000*2*pi; Rp=3; Rs=15; Fs=10000; Ts=1/Fs; Nn=128; [N,Wn]=buttord(Wp,Ws,Rp,Rs,'s') [b,a]=butter(N,Wn,'s'); [z,p,k]=buttap(N); [Bap,Aap]=zp2tf(z,p,k); [b,a]=lp2lp(Bap,Aap,Wn); [bz,az]=bilinear(b,a,Fs) freqz(bz,az,Nn,Fs) 12 N =5 Wn =1.3387e+004 bz =0.0171 0.0854 0.1708 0.1708 0.0854 0.0171 az =1.0000 -1.2271 1.1622 -0.5176 0.1450 -0.0159 仿真结果如图4.5所示。 图4.5 数字低通滤波器仿真 5.总结 本次课程设计设计的是数字低通滤波器，在设计的过程中通过上网、看书来查阅与本次数字信号课程设计有关的资料，加深了对数字滤波器工作原理、设计方法以及相关的设计过程的了解，增强了认真思考和合作解决问题的能力，更重要的是动手操作MATLAB并实现了数字低通滤波器的模拟过程。再加上考试的时间将近，通过这这一方面的实践，将与之相关的知识点牢牢地记在心里。 总之，这次设计让我受益匪浅。 13 参考文献: [1] 程佩青，《数字信号处理教程》，北京，清华大学出版社，2011年 工业出版社，2005年 [2] 刘泉，阙大顺，《数字信号处理原理与实现》， 北京，电子 [3] 胡广书，《数字信号处理:理论算法与实现》，北京，清华大学出版社，2003年 [4] 姚天任，江太辉，《数字信号处理》，武汉，华中理工大学出版社，2000年 [5] 丁美玉，高西全，《数字信号处理》，西安，西安电子科技大学出版社，2009年 [6] Ingle.V.K，《数字信号处理(MATLAB版)》，西安，西安交通大学出版社，2008年 [7] 吴镇扬，《数字信号处理》，北京，高等教育出版社，2010年 [8] 张志涌，杨祖樱，《MATLAB教程》，北京，北京航空航天大学出版社，2010年 [9] Holly Moore，《MATLAB使用教程(第二版)》，北京，电子工业出版社，2010年 [10] 张琨，《MATLAB 2010从入门到精通》，北京，电子工业出版社，2011年 14 《数字信号处理》课程设计指导教师评语 评语: 指导教师: (签字) 年 月 日 15