IIR数字滤波器设计2.doc

1、吉林建筑工程学院电气与电子信息工程学院数字信号处理课程设计报告设计题目： IIR数字滤波器设计 专业班级： 学生姓 学 号： 指导教师： 设计时间： 2012.01.022012.01.06 教师评语：成绩 评阅教师 日期 IIR数字滤波器设计报告一、 设计的作用、目的课程设计是理论学习的延伸，是掌握所学知识的一种重要手段，对于贯彻理论联系实际、提高学习质量、塑造自身能力等于有特殊作用。本次课程设计一方面通过MATLAB仿真设计内容，使我们加深对理论知识的理解，同时增强其逻辑思维能力，另一方面对课堂所学理论知识作一个总结和补充。有以下几方面作用：1、有利于基础知识的理解2、有利于逻辑思维的锻炼

2、在许多常规学科的日常教学中，我们不难发现这样一个现象，不少学生的思维常常处于混乱的状态。写起作文来前言不搭后语，解起数学题来步骤混乱，这些都是缺乏思维训练的结果。程序设计是公认的、最能直接有效地训练学生的创新思维，培养分析问题、解决问题能力的学科之一。即使一个简单的程序，从任务分析、确定算法、界面布局、编写代码到调试运行，整个过程学生都需要有条理地构思，这中间有猜测设想、判断推理的抽象思维训练，也有分析问题、解决问题、预测目标等能力的培养。3、有利于与其他学科的整合在程序设计中，我们可以解决其它学科有关问题，也利用其它课程的有关知识来解决比较抽象很难理解的知识。4、有利于治学态度的培养。程序设

3、计中，语句的语法和常量变量的定义都有严格的要求，有时输了一个中文标点、打错了一个字母，编译就不通过，程序无法正常运行。因此，程序设计初学阶段，学生经常会犯这样的错误，可能要通过几次乃至十多次的反复修改、调试，才能成功，但这种现象会随着学习的深入而慢慢改观。这当中就有一个严谨治学、一丝不苟的科学精神的培养，又有一个不怕失败、百折不挠品格的锻炼。课程设计的目的就是要锻炼大家运用所学知识解决实际问题的能力，让我们在实践中更加深刻的了解滤波器。数字信号处理课程设计的目的是为了使学生加深对所学的知识的理解，培养学生的专业素质，提高其利用知识处问题的能力，为今后专业课程的学习、毕业设计和工作打下良好的基础

4、。使学生能比较扎实地掌握数字信号处理的基础知识和基本理论，掌握数字信号处理及有关设备的分析、开发等基本技能，受到必要工程训练、初步的科学研究方法训练和实践锻练，增强分析问题和解决问题的能力，了解数字信号处理的新发展。二、设计任务及要求1.掌握双线性变换法及脉冲响应不变法设计IIR数字滤波器以及窗函数法设计FIR数字滤波器的原理、具体方法及计算机编程2.观察双线性变换法、脉冲响应不变法及窗函数法设计的滤波器的频域特性，了解各种方法的特点3.用MATLAB画出三种方法设计数字滤波器的幅频特性曲线，记录带宽和衰减量，检查结果是否满足要求。三、设计内容IIR数字滤波器的设计内容：已知通带截止频率,通带

5、最大衰减,阻带截止频率，阻带最小衰减，T=1ms，按照以上技术要求，用脉冲响应不变法和双线性变换法设计巴特沃斯数字低通滤波器，并观察所设计数字滤波器的幅频特性曲线。四、设计原理4.1数字低通滤波器的设计原理1.滤波器的概念：滤波器是对输入信号起滤波的作用的装置。当输入、输出是离散信号，滤波器的冲激响应是单位抽样响应 h（n）时，这样的滤波器称作数字滤波器（ DF ）。DF是由差分方程描述的一类特殊的离散时间系统。2.数字滤波器的系统函数与差分方程:系统函数差分方程 对上式进行 Z反变换，即得滤波器的功能与实现滤波就是对输入序列 x(n)进行一定的运算操作，从而得到输出序列y(n)。实现滤波从运

6、算上看,只需三种运算：加法、单位延迟、乘常数。实现的方法有两种： （1）利用通用计算机编程，即软件实现; （2）数字信号处理器（DSP）,即专用硬件实现。 3.数字滤波器结构的表示方法 方框图法、信号流图法 方框图表示法 流图表示法 三种基本运算的表示图 数字滤波器是离散时间系统，所处理的信号是离散时间信号。 一般时域离散系统或网络可以用差分方程、单位脉冲响应以及系统函数进行描述。滤波器的传递函数为: 4. 数字滤波器的分类 可分为: 低通(LP), 高通(HP),带通(BP), 带阻(BS) 滤波一般分为两类：模拟滤波器和数字滤波器。在时域的实现方法与方式上，它们是完全不同的两类系统，其差别

7、主要有两点：(1) 数字滤波器处理的是数字信号，而模拟滤波器处理的是连续信号。(2) 数字滤波器可以用数字硬件结构来实现，即用数字硬件构成专用数字计算机；也可以用编程序的方法利用通用计算机来实现，亦即软件实现。而模拟滤波器则是由分立元件组成的电路网络系统来实现的。数字滤波器是数字信号处理的一个重要组成部分。数字滤波实际上是一种运算过程，其功能是将一组输入的数字序列通过一定的运算后转变为另一组输出的数字序列，因此它本身就是一台数字式的处理设备。数字滤波器一般可以用两种方法实现：一种是根据描述数字滤波器的数学模型或信号流图，用数字硬件装配成一台专门的设备，构成专用的信号处理机；另一种方法就是直接利

8、用通用计算机，将所需要的运算编成程序让计算机来执行， 这也就是用软件来实现数字滤波器。 对于不同结构的滤波器，其所需要的存储单元及乘法次数是不同的，前者影响复杂性，后者影响运算速度。此外在有限精度（有限字长）的情况下，不同的运算结构的误差、稳定性是不同的。下面分别就IIR和FIR的结构进行讨论1.无限长单位冲激响（IIR）滤波器的基本结构:IIR滤波器的特点： 1、单位冲激响应h(n)是无限长的。 2、系统函数H(z)在有限Z平面上有极点存在。 3、结构上是递归型的，即存在着输出到输入的反馈。N阶IIR滤波器的输入输出关系可以用下式表示的N阶的差分方程来描述：系统函数为：IIR数字滤波器的结构

9、分类1.直接型（型）一个N阶的IIR滤波器的N阶的差分方程为：从这个差分方程表达式可以看出，系统的输出y(n)由两部分构成：第一部分是一个对输入x(n)的M阶延时链结构，每阶延时抽头后加权相加，构成一个横向结构网络。 第二部分 是一个对输出y(n)的N阶延时链的横向结构网络，是由输出到输入的反馈网络。 由这两部分相加构成输出，如图5-4。从图上可以看出，直接型结构需要M+N个延时器和M+N+1个乘法器。 直接型结构图 2.直接型:直接型结构又称为正准型结构。由图5-4，直接型结构的系统函数H(z)也可以看成是两个独立的系统函数的乘积。输入信号x(n)先通过系统H1(z)，得到中间输出变量y1(

10、n)，然后再把y1(n)通过系统H2(z)得到输出信号y(n)。 即 令:假设所讨论的IIR数字滤波器是线性非时变系统，显然交换H1(z)和H2(z)的级联次序不会影响系统的传输效果，即 其结构如图所示: 直接型的变形结构 直接型结构 3.级联型 若把N阶IIR滤波器的系统函数H(z)的分子和分母分别进行因式分解，得到多个因式连乘积的形式:式中：A为常数，ci和di分别表示H(z)的零点和极点。由于H(z)的分子和分母都是实系数多项式，而实系数多项式的根只有实根和共轭复根两种情况。将每一对共轭零点（极点）合并起来构成一个实系数的二阶因子，并把每两个的实根因子合成一个二阶因子，则可以把H(z)表

11、示成多个实系数的二阶数字网络Hj(z)的连乘积形式。 且:若每一个实系数的二阶数字网络的系统函数Hj(z)的网络结构均采用前面介绍的直接型结构，则一个二阶子系统的结构为：则可以得到系统函数H(z)的级联型结构，如图所示4.并联型把传递函数H(z)展开成部分分式之和的形式，就可以得到滤波器的并联型结构。 当N=M 时，展开式为 和级联型结构的方法类似，将上式中的共轭复根部分两两合并得到实系数的二阶网络，则有： 由上式可知，滤波器可由E个一阶网络、F个二阶网络和一个常数支路并联构成， 其结构如图所示IIR数字滤波器的结构特点1.直接型：直接型比直接型结构延时单元少，用硬件实现可以节省寄存器，比直接

12、型经济；若用软件实现则可节省存储单元。但对于高阶系统直接型结构都存在调整零、 极点困难，对系数量化效应敏感度高等缺点。 2.级联型：在级联型结构中，每一个一阶网络只关系到滤波器的一个零点、一个极点；每个二阶网络只关系到滤波器的一对零点和一对极点。调整系数0j、1j和2j只会影响滤波器的第j对零点，对其他零点并无影响；同样, 调整分母多项式的系数1j和2j也只单独调整了第j对极点。因此，与直接型结构相比， 级联型结构便于准确地实现滤波器零、极点的调整。此外，因为在级联结构中，后面的网络的输出不会流到前面，所以其运算误差也比直接型小。 3.并联型：并联型结构使用的加法器，乘法器，延时单元基本与级联

13、结构相同。它的每个一阶节单独确定一个实数极点，每个二阶节确定一对共轭极点，各条支路互不影响；它只能独立的调整各极点的位置，不能单独调整零点的位置；此外，由于各基本节是并联的，故并联结构的误差比级联结构的运算误差小。并联型结构也可以单独调整极点位置，但对于零点的调整却不如级联型方便，而且当滤波器的阶数较高时，部分分式展开比较麻烦。在运算误差方面，由于各基本网络间的误差互不影响， 没有误差积累， 因此比直接型和级联型误差稍小一点。 IIR 数字滤波器几种结构的比较 直接I型和直接II型实现起来具有简单直观的特点。需要(M+N)个加法器和(M+N)个乘法器，直接II型比直接I型节省M个延时单元，在M

14、=N的情况下，需要N个延时单元。直接性的主要缺点在于差分方程的系数ak，bk对滤波器的性能控制不直接，同时由于其高度反馈性，容易出现不稳定或产生较大误差。 级联型结构的特点是每个二阶节是相互独立的，可分别通过调整各个 “ 零极点对 ” 来对滤波器性能进行较好的控制，且各二阶节的顺序可重排，能有效的减少有限字长效应。实现需要(M+N)个加法器、(M+N)个乘法器和N个延时单元。该结构应用最广泛。 并联型结构使用的加法器，乘法器，延时单元基本与级联结构相同。它的每个一阶节单独确定一个实数极点，每个二阶节确定一对共轭极点，各条支路互不影响；它只能独立的调整各极点的位置，不能单独调整零点的位置；此外，

15、由于各基本节是并联的，故并联结构的误差比级联结构的运算误差小。滤波器（filter），是一种用来消除干扰杂讯的器件，将输入或输出经过过滤而得到纯净的直流电。对特定频率的频点或该频点以外的频率进行有效滤除的电路，就是滤波器，其功能就是得到一个特定频率或消除一个特定频率。数字滤波器：与模拟滤波器相对应，在离散系统中广泛应用数字滤波器。它的作用是利用离散时间系统的特性对输入信号波形或频率进行加工处理。或者说，把输入信号变成一定的输出信号，从而达到改变信号频谱的目的。数字滤波器一般可以用两种方法来实现：一种方法是用数字硬件装配成一台专门的设备，这种设备称为数字信号处理机；另一种方法就是直接利用通用计算

16、机，将所需要的运算编成程序让通用计算机来完成，即利用计算机软件来实现。低通滤波器：(low-passfilter)低通滤波器是容许低于截至频率的信号通过， 但高于截止频率的信号不能通过的电子滤波装置。对于不同滤波器而言，每个频率的信号的减弱程度不同。当使用在音频应用时，它有时被称为高频剪切滤波器, 或高音消除滤波器。低通滤波器概念有许多不同的形式，其中包括电子线路（如音频设备中使用的hiss 滤波器、平滑数据的数字算法、音障（acoustic barriers）、图像模糊处理等等，这两个工具都通过剔除短期波动、保留长期发展趋势提供了信号的平滑形式。低通滤波器在信号处理中的作用等同于其它领域如金

17、融领域中移动平均数（moving average)所起的作用； 低通滤波器有很多种，其中，最通用的就是巴特沃斯滤波器和切比雪夫滤波器。IIR滤波器与FIR滤波器的分析比较IIR和FIR数字滤波器的设计方法，选择哪一种滤波器取决于每种类型滤波器的优点在设计中的重要性。为了能在实际工作中恰当地选用合适的滤波器，现将两种滤波器特点比较分析如下：(1) 选择数字滤波器是必须考虑经济问题，通常将硬件的复杂性、芯片的面积或计算速度等作为衡量经济问题的因素。在相同的技术指标要求下，由于IIR数字滤波器存在输出对输入的反馈，因此可以用较少的阶数来满足要求，所用的存储单元少，运算次数少，较为经济。例如，用频率抽

18、样法设计一个阻带衰减为20dB的FIR数字滤波器，要33阶才能达到要求，而用双线性变换法只需45阶的切比雪夫IIR滤波器就可达到同样的技术指标。这就是说FIR滤波器的阶数要高510倍左右。(2) 在很多情况下,FIR数字滤波器的线性相位与它的高阶数带来的额外成本相比是非常值得的。对于IIR滤波器，选择性越好，其相位的非线性越严重。如果要使IIR滤波器获得线性相位，又满足幅度滤波器的技术要求，必须加全通网络进行相位校正，这同样将大大增加滤波器的阶数。就这一点来看，FIR滤波器优于IIR滤波器。(3) FIR滤波器主要采用非递归结构，因而无论是理论上还是实际的有限精度运算中他都是稳定的，有限精度运

19、算误差也较小。IIR滤波器必须采用递归结构，极点必须在z平面单位圆内才能稳定。对于这种结构，运算中的舍入处理有时会引起寄生振荡。(4) 对于FIR滤波器，由于冲激响应是有限长的，因此可以用快速傅里叶变换算法，这样运算速度可以快得多。IIR滤波器不能进行这样的运算。(5) 从设计上看，IIR滤波器可以利用模拟滤波器设计的现成的闭合公式、数据和表格，可以用完整的设计公式来设计各种选频滤波器。一旦选定了已知的一种逼近方法(如巴特奥兹，切比雪夫等)，就可以直接把技术指标带入一组设计方程计算出滤波器的阶次和系统函数的系数(或极点和零点)。FIR滤波器则一般没有现成的设计公式。窗函数法只给出了窗函数的计算

20、公式，但计算通带和阻带衰减仍无显式表达式。一般FIR滤波器设计仅有计算机程序可资利用，因而要借助于计算机。(6) IIR滤波器主要是设计规格化、频率特性为分段常数的标准低通、高通、带通和带阻滤波器。FIR滤波器则灵活很多，例如频率抽样法可适应各种幅度特性和相位特性的要求。因此FIR滤波器可设计出理想正交变换器、理想微分器、线性调频器等各种网络，适应性很广。而且，目前已经有很多FIR滤波器的计算机程序可供使用。 表4-1、 两种滤波器特点比较分析FIR滤波器IIR滤波器设计方法一般无解析的设计公式，要借助计算机程序完成利用AF的成果，可简单、有效地完成设计设计结果可得到幅频特性(可以多带)和线性

21、相位(最大优点)只能得到幅频特性，相频特性未知，如需要线性相位，须用全通网络校准，但增加滤波器阶数和复杂性稳定性极点全部在原点(永远稳定)无稳定性问题有稳定性问题阶数高低结构非递归系统递归系统运算误差一般无反馈，运算误差小有反馈，由于运算中的四舍五入会产生极限环4.2变换方法的原理1.用脉冲相应不变法设计IIR数字滤波器 利用模拟滤波器来设计数字滤波器，也就是使数字滤波器能模仿模拟滤波器的特性，这种模仿可以从不同的角度出发。脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列h(n)模仿模拟滤波器的冲激响应ha(t)，即将ha(t)进行等间隔采样，使h(n)正好等于ha(t)

22、的采样值，满足h(n)=ha(nT)式中,T是采样周期。如果令Ha(s)是ha(t)的拉普拉斯变换，H(z)为h(n)的Z变换，利用采样序列的Z变换与模拟信号的拉普拉斯变换的关系得(4-1)则可看出，脉冲响应不变法将模拟滤波器的S平面变换成数字滤波器的Z平面，这个从s到z的变换z=esT是从S平面变换到Z平面的标准变换关系式。图1-1脉冲响应不变法的映射关系 由（4-1）式，数字滤波器的频率响应和模拟滤波器的频率响应间的关系为(4-2)这就是说，数字滤波器的频率响应是模拟滤波器频率响应的周期延拓。正如采样定理所讨论的，只有当模拟滤波器的频率响应是限带的，且带限于折叠频率以内时，即(4-3)才能

23、使数字滤波器的频率响应在折叠频率以内重现模拟滤波器的频率响应，而不产生混叠失真，即| (4-4)但是，任何一个实际的模拟滤波器频率响应都不是严格限带的，变换后就会产生周期延拓分量的频谱交叠，即产生频率响应的混叠失真，如图4-2所示。这时数字滤波器的频响就不同于原模拟滤波器的频响，而带有一定的失真。当模拟滤波器的频率响应在折叠频率以上处衰减越大、越快时，变换后频率响应混叠失真就越小。这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。对某一模拟滤波器的单位冲激响应ha(t)进行采样，采样频率为fs，若使fs增加，即令采样时间间隔（T=1/fs）减小，则系统频率响应各周期延拓分量之间相距更

24、远，因而可减小频率响应的混叠效应。图4-2脉冲响应不变法中的频响混叠现象 脉冲响应不变法优缺点：从以上讨论可以看出，脉冲响应不变法使得数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，也就是时域逼近良好，而且模拟频率和数字频率之间呈线性关系=T。因而，一个线性相位的模拟滤波器（例如贝塞尔滤波器）通过脉冲响应不变法得到的仍然是一个线性相位的数字滤波器。脉冲响应不变法的最大缺点是有频率响应的混叠效应。所以，脉冲响应不变法只适用于限带的模拟滤波器(例如，衰减特性很好的低通或带通滤波器)，而且高频衰减越快，混叠效应越小。至于高通和带阻滤波器，由于它们在高频部分不衰减，因此将完全混淆在低频响应中

25、。如果要对高通和带阻滤波器采用脉冲响应不变法，就必须先对高通和带阻滤波器加一保护滤波器，滤掉高于折叠频率以上的频率，然后再使用脉冲响应不变法转换为数字滤波器。当然这样会进一步增加设计复杂性和滤波器的阶数。2.用双线性变换法设计IIR数字滤波器脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-/T/T之间，再用z=esT转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的-/T/T一条横带里；第二步再通过标准变换关系z=es1T将此横带变换到整个Z平面上去。

26、这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，映射关系如图4-3所示图4-3双线性变换的映射关系为了将S平面的整个虚轴j压缩到S1平面j1轴上的-/T到/T段上，可以通过以下的正切变换实现（4-5）式中,T仍是采样间隔。当1由-/T经过0变化到/T时，由-经过0变化到+，也即映射了整个j轴。将式（4-5）写成将此关系解析延拓到整个S平面和S1平面，令j=s，j1=s1，则得再将S1平面通过以下标准变换关系映射到Z平面z=es1T从而得到S平面和Z平面的单值映射关系为：（4-6）（4-7）式（4-6）与式（4-7）是S平面与Z平面之间的单值映射关系，这

27、种变换都是两个线性函数之比，因此称为双线性变换式（4-5）与式（4-6）的双线性变换符合映射变换应满足的两点要求。首先,把z=ej，可得（4-8）即S平面的虚轴映射到Z平面的单位圆。其次，将s=+j代入式（4-8），得因此由此看出，当0时，|z|0时，|z|1。也就是说，S平面的左半平面映射到Z平面的单位圆内，S平面的右半平面映射到Z平面的单位圆外，S平面的虚轴映射到Z平面的单位圆上。因此，稳定的模拟滤波器经双线性变换后所得的数字滤波器也一定是稳定的。双线性变换法优缺点双线性变换法与脉冲响应不变法相比，其主要的优点是避免了频率响应的混叠现象。这是因为S平面与Z平面是单值的一一对应关系。S平面整

28、个j轴单值地对应于Z平面单位圆一周，即频率轴是单值变换关系。这个关系如式（4-8）所示，重写如下：上式表明，S平面上与Z平面的成非线性的正切关系，如图4-4所示。由图7-7看出，在零频率附近，模拟角频率与数字频率之间的变换关系接近于线性关系；但当进一步增加时，增长得越来越慢，最后当时，终止在折叠频率=处，因而双线性变换就不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象，从而消除了频率混叠现象。 图4-4双线性变换法的频率变换关系但是双线性变换的这个特点是靠频率的严重非线性关系而得到的，如式（4-8）及图4-4所示。由于这种频率之间的非线性变换关系，就产生了新的问题。首先，一个线性相位的

29、模拟滤波器经双线性变换后得到非线性相位的数字滤波器，不再保持原有的线性相位了；其次，这种非线性关系要求模拟滤波器的幅频响应必须是分段常数型的，即某一频率段的幅频响应近似等于某一常数（这正是一般典型的低通、高通、带通、带阻型滤波器的响应特性），不然变换所产生的数字滤波器幅频响应相对于原模拟滤波器的幅频响应会有畸变，如图4-5所示。 图4-5双线性变换法幅度和相位特性的非线性映射对于分段常数的滤波器，双线性变换后，仍得到幅频特性为分段常数的滤波器，但是各个分段边缘的临界频率点产生了畸变，这种频率的畸变，可以通过频率的预畸来加以校正。也就是将临界模拟频率事先加以畸变，然后经变换后正好映射到所需要的数

30、字频率上。五、设计步骤5.1 Matlab设计步骤：目前，IIR数字滤波器设计最通用的方法是借助于模拟滤波器的设计方法。模拟滤波器设计已经有了一套相当成熟的方法，它不但有完整的设计公式，而且还有较为完整的图表供查询，因此，充分利用这些已有的资源将会给数字滤波器的设计带来很大方便，IIR数字滤波器的设计步骤是：(1)按一定规则将给出的数字滤波器的技术指标转换为模拟滤波器的技术指标；(2)根据转换后的技术指标设计模拟低通滤波器H(s)；(3)在按一定规则将H(s)转换为H(z)。若所设计的数字滤波器是低通的，那么上述设计工作可以结束，若所设计的是高通、带通或者带阻滤波器，那么还有步骤：(4)将高通

31、、带通或者带阻数字滤波器的技术指标先转化为低通滤波器的技术指标，然后按上述步骤(2)设计出模拟低通滤波器H(s)，再由冲击响应不变法或双线性变换将H(s)转换为所需的H(z)。s - z 映射的方法有：冲激响应不变法、阶跃响应不变法、双线性变换法等。下面讨论双线性变换法。双线性变换法8是指首先把s 平面压缩变换到某一中介平面s1 的一条横带(宽度为2T,即从- T到T) ,然后再利用的关系把s1平面上的这条横带变换到整个z 平面。这样s 平面与z 平面是一一对应关系, 消除了多值变换性, 也就消除了频谱混叠现象。s 平面到z 平面的变换可采用 (5-1) (5-2)令 ，有： (5-3)从s1

32、 平面到z 平面的变换,即 (5-4)代入上式，得到： (5-5) 一般来说，为使模拟滤波器的某一频率与数字滤波器的任一频率有对应关系，可引入代定常数c， (5-6)则 (5-7) 这种s 平面与z 平面间的单值映射关系就是双线性变换。有了双线性变换，模拟滤波器的数字化只须用进行置换。设滤波器的通带截止频率为wp，止带截止频率为ws，通带衰减为Rp，止带衰减为As，通带波纹为Rp。1、根据给定指标得出，低通模拟原型滤波器的阶数和截止频率。利用函数buttord, cheb1ord,cheb2ord。 n,wn=buttord(wp,ws,Rp,As,s) n,wn=cheb1ord(wp,ws

33、,Rp,As,s) 这里我们采用Buttord函数。2、求出低通模拟原型滤波器，利用函数Buttap,Cheb1ap,Cheb2ap。 z,p,k=Buttap(n) z,p,k=Cheb1ap(n,Rp) 由于所得的结果为零极点型，还必须转成b/a型，可用函数zp2tf。 bap,aap=zp2tf(z,p,k)3、将模拟低通原型滤波器经频率变换为所要求的模拟滤波器（低通，高通，带通，带阻）可用函数lp2lp，lp2hp，lp2bs， lp2bp。 b,a=lp2lp(bap,aap,wn) b,a=lp2bs(bap,aap,wn)4、求出滤波器的幅频，相频及冲激响应。 db,mag,ph

34、a,w=freqs_m(b,a,Wmax)(P208） ha,x,t=impulse(b,a)(matlab自带） 直接画出幅频特性和相频特性：freqs(b,a) matlab自带5.程序流程框图：开始读入数字滤波器技术指标将指标转换成归一化模拟低通滤波器的指标设计归一化的模拟低通滤波器阶数N和Wn截止频率模拟域频率变换，将G(P)变换成模拟带通滤波器H(s)用双线性变换法将H(s)转换成数字带通滤波器H(z)输入信号后显示相关结果结束5.2 计算步骤：(1)根据已知条件，确定技术指标,根据技术指标用buttord函数求出最小滤波器阶数N和截止频率Wn：N,Wn=buttord(Wp,Ws,

35、Rp,Rs)T=1ms；Fs=1000通带截止频率:wp=2*0.2；阻带截止频率:ws=2*0.3；通带最大衰减p=1dB和阻带最小衰减:s=25dB。(2)用=2/T*tan(w/2)对低通数字滤波器H(z)的数字边界频率预畸变，得到低通模拟滤波器H(s)的边界频率：通带截止频率p、阻带截止频率s。 (3)设计模拟低通原型滤波器。用模拟低通滤波器设计方法得到模拟低通滤波器的传输函数Ha(s);借助巴特沃斯(Butterworth)滤波器。(4)利用双线性变换法将模拟低通滤波器Ha(s)转换成数字低通滤波器H(z)。用butter函数求出滤波器系数：b,a=butter(N,Wn,ftype

36、)。用Filter函数平滑滤波。六、 数字低通滤波器MATLAB编程及幅频特性曲线6.1 MATLAB语言编程1.双线性变换法fp=200;fs=300;Rp=1;Rs=25;T=0.001;Fs=1000;Wp=2*pi*fp;Ws=2*pi*fs;N,Wn=buttord(Wp,Ws,Rp,Rs,s); % 阶数N、截止频率Wnb,a=butter(N,Wn,s); % 模拟低通滤波器极零点numd,dend=bilinear(b,a,Fs); % 双线性变换法w=linspace(0,pi,1000);h=freqz(numd,dend,w);plot(w,20*log10(abs(h)

37、;axis(0 pi -55 5);grid;xlabel(omega/rad);ylabel(幅度/dB);title(双线性变换法IIR滤波器幅频特性曲线);2.脉冲响应不变法fp=200;fs=300;Rp=1;Rs=25;T=0.001;Fs=1000;Wp=2*pi*fp;Ws=2*pi*fs; N,Wn=buttord(Wp,Ws,Rp,Rs,s); % 阶数N、截止频率Wnb,a=butter(N,Wn,s); % 模拟低通滤波器极零点numd,dend=impinvar(b,a,Fs); % 脉冲响应不变法w=linspace(0,pi,1000);h=freqz(numd,d

38、end,w);plot(w,20*log10(abs(h);axis(0 pi -55 5);grid;xlabel(omega/rad);ylabel(幅度/dB);title(脉冲不变响应法IIR滤波器幅频特性曲线);6.2 幅频特性曲线七、体会及建议 我将把我本次课程设计的心得和体会简述如下：(1)类型寻优： 确定好所需要设计的滤波器用途后，首先对其进行性能需求分析，明确该系统应该达到的各种性能指标，其次，拟定多种滤波器类型，对这些方案采用Matlab进行仿真，进行综合分析和比较，选择出最佳的滤波器类型作为本设计方案，然后，依据其性能指标编写matlab程序，确定二阶节系数。(2)系统思

39、维： 必须有系统的设计思维，把每一个细节都放到整个系统中考虑，考虑整个系统设计的可行性、完整性、稳定性和功能的实现，这样才不会局限在细节上，才能快速的完成性能优越的硬件设计。(3)把握细节： 系统也是由细节构成的，在把握整个系统思维的基础上把握每一个细节，因为每一个细节都有可能决定整个系统的性能。在写MATLAB代码时，把每一个变量与过程考虑清楚才可能完成整个程序。(4)勇于尝试： 系统即使经过了非常严格的论证也仍然可能存在问题，或许面对问题一时没有很好的解决方案，但是有一个或许可行的大胆的想法，不要犹豫，试一试吧。勇于尝试往往能够找到更好的解决方法。(5)保持自信： 无论遇到什么困难，我都相

40、信自己一定能够找到解决的方法，有的时候只要稍微再用一点力、使一点劲，结果就会不一样。运用MATLAB软件进行滤波器设计，免去了大量计算，通过绘制图像，更加直观的理解了两种变换方法间的不同以及各自的优缺点。对滤波器有了一个更加感性的认识。利用MATLAB的强大运算功能，基于MATLAB的信号处理工具箱(Signal Processing Toolbox)的数字滤波器设计法可以快速有效地设计由软件组成的常规数字滤波器，设计方便、快捷，大大减轻了工作量。在设计过程中可以对比滤波器特性。随时更改参数，以达到滤波器设计的最优化。利用MATLAB设计数字滤波器在数字通信系统和计算机领域信号处理中，有着广泛

41、的应用前景,可关注。基于MATLAB的信号处理工具箱为数字滤波器设计带来了全新的实现手段,设计快捷方便,仿真波形直观。八、参考文献1. 高西全 数字信号处理 西安电子科技大学出版社 20082. 丁玉美 数字信号处理 西安电子科技大学出版社 20083. 张明照 应用MATLAB实现信号分析和处理 科学出版社20054. 余成波 数字信号处理及MATLAB实现 清华大学出版社20025. 罗军辉 MATLAB7.0在数字信号处理中的应用 机械工业出版社20066. 董长虹 Matlab信号处理与应用 国防工业出版社 20047. 刘令普 数字信号处理 哈尔滨工业大学出版社 2002 8. 邓重一 滤波技术的发展现状 中国仪器仪表出版社 2004 9. 王世一 数字信号处理 北京理工大学出版社 2005 10. 李行一数字信号处理 重庆大学出版社 2002 11. 刘正士，王勇，陈恩伟，葛运建 一种数字滤波器的设计方法及其应用 中国机械工程 200612. 张亚妮 基于MATLAB的数字滤波器设计 辽宁工程技术大学出版社 200513. 刘波 MATLAB信号处理 电子工业出版社 2006 14. 刘卫国，陈昭民，张颖 Matlab程序设计与应用 高等教育出版社 2002