基于MatlabIIR数字滤波器设计

《基于MatlabIIR数字滤波器设计》由会员分享，可在线阅读，更多相关《基于MatlabIIR数字滤波器设计（24页珍藏版）》请在装配图网上搜索。

1、-摘 要在现代通信系统中，由于信号中经常混有各种复杂成分，所以很多信号分析都是基于滤波器而进展的，而数字滤波器是通过数值运算实现滤波，具有处理精度高、稳定、灵活、不存在阻抗匹配问题，可以实现模拟滤波器无法实现的特殊滤波功能。数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。实现IIR滤波器的阶次较低，所用的存储单元较少，效率高，精度高，而且能够保存一些模拟滤波器的优良特性，因此应用很广。Matlab软件以矩阵运算为根底，把计算、可视化及程序设计有机融合到交互式工作环境中，并且为数字滤波的研究和应用提供了一个直观、高效、

2、便捷的利器。尤其是Matlab中的信号处理工具箱使各个领域的研究人员可以直观方便地进展科学研究与工程应用。本文首先介绍了数字滤波器的概念，分类以及设计要求。接着利用MATLAB函数语言编程，用信号处理图形界面FDATool来设计滤波器以及Sptool界面设计的方法，并用FDATool模拟IIR数字滤波器处理信号。重点设计Chebyshev I型和Chebyshev II型数字低通滤波器,并介绍最优化设计。【关键字】IIR滤波器FDAToolSptoolSimulinkABSTRACTIn modern munication systems,Because often mi*ed with va

3、rious signal ple* ponents,So many signal analysis is based on filters, and the digital filter is realized through numerical putation, digital filters filter with high precision, stability and fle*ibility, dont e*ist, can realize the impedance matching simulating the special filter cannot achieve fil

4、ter function. Digital filter according to its impulse response function and characteristics of the time can be divided into two kinds, namely the infinite impulse response (IIR) digital filter and finite impulse response (FIR digital filters). The order of realizing IIR filter is used, low and high

5、efficiency less storage unit, high precision, and can keep some simulation characteristics of filter, so it is widely used. Matlab software based on matri* putation, the calculation, visualization and program design of organic integration to interactiveenvironment for digital filter, and the researc

6、h and application of provides an intuitive, efficient and convenient tool. Especially in the Matlab signal processing to all areas of research toolbo* personnel can easily for scientific research and engineering application. This paper introduces the concept of digital filter, classification and des

7、ign requirements. Then using MATLAB language programming, with functions of signal processing FDATool graphical interface design of interface design and Sptool filter, and FDATool analog signal processing IIR digital filter. Key design Chebyshev type I and II digital Chebyshev lowpass filter, and in

8、troduces optimization design.【Keywords】 IIR Filter FDATool Sptool Simulink目录前言第一章数字滤波器第一节数字滤波器的概念第二节数字滤波器的分类第三节数字滤波器的设计要求第二章 IIR数字滤波器设计方法第一节 IIR数字滤波器的设计步骤第二节用脉冲相应不变法设计IIR数字滤波器一、设计原理二、脉冲响应不变法优缺点第三节双线性变换法设计IIR数字滤波器一、设计原理二、双线性变换法优缺点第三章 IIR滤波器的MATLAB设计第一节 IIR数字滤波器的典型设计法第二节 IIR数字滤波器的直接设计法第三节 FDATool介绍和界面

9、设计第四节 FDATOOL设计IIR数字滤波器第五节 SIMULINK 仿真IIR滤波器总结致参考文献完毕语z.-前 言随着信息时代和数字世界的到来，数字信号处理已成为当今一门极其重要的学科和技术领域。目前数字信号处理在通信、语音、图像、自动控制、雷达、军事、航空航天、医疗和家用电器等众多领域得到了广泛的应用。在数字信号处理中起着重要的作用并已获得广泛应用的是数字滤波器DF，Digital Filter，根据其单位冲激响应函数的时域特性可分为两类：无限冲激响应IIRInfinite Impulse Response滤波器和有限冲激响应FIRFinite Impulse Response滤波器。

10、与FIR滤波器相比，IIR的实现采用的是递归构造，极点须在单位圆，在一样设计指标下，实现IIR滤波器的阶次较低，即所用的存储单元少，从而经济效率高。MATLAB是英文MATri* LABoratory(矩阵实验室)的缩写。它是美国的MathWorks公司推出的一套用于科学计算和图形处理可视化、高性能语言与软件环境。MATLAB的信号处理工具箱是专门应用于信号处理领域的专用工具箱，它的两个根本组成就是滤波器的设计与实现局部以及谱分析局部。工具箱提供了丰富而简便的设计，使原来繁琐的程序设计简化成函数的调用。只要以正确的指标参数调用相应的滤波器设计程序或工具箱函数，便可以得到正确的设计结果，使用非常

11、方便。第一章 数字滤波器第一节数字滤波器的概念滤波器是指用来对输入信号进展滤波的硬件和软件。数字滤波器是对数字信号实现滤波的线性时不变系统。数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。数字滤波实质上是一种运算过程，实现对信号的运算处理。数字滤波器和模拟滤波器相比，因为信号的形式和实现滤波的方法不同，数字滤波器具有比模拟滤波器精度高、稳定、体积小、重量轻、灵活、不要求阻抗匹配等优点。输入数字信号数字序列通过特定的运算转变为输出的数字序列，因此，数字滤波器本质上是一个完成特定运算的数字计算过程，也可

12、以理解为是一台计算机。描述离散系统输出与输入关系的卷积和差分方程只是给数字信号滤波器提供运算规则，使其按照这个规则完成对输入数据的处理。时域离散系统的频域特性:式1-1其中、分别是数字滤波器的输出序列和输入序列的频域特性或称为频谱特性,是数字滤波器的单位取样响应的频谱，又称为数字滤波器的频域响应。输入序列的频谱经过滤波后,因此，只要按照输入信号频谱的特点和处理信号的目的，适中选择，使得滤波后的满足设计的要求，这就是数字滤波器的滤波原理。第二节数字滤波器的分类按照不同的分类方法，数字滤波器有许多种类，但总起来可以分成两大类：经典滤波器和现代滤波器。经典滤波器的特点是其输入信号中有用的频率成分和希

13、望滤除的频率成分占有不同的频带，通过一个适宜的选频滤波器滤除干扰，得到纯洁信号，到达滤波的目的。但是，如果信号和干扰的频谱相互重叠，则经典滤波器不能有效地滤除干扰，最大限度地恢复信号，这时就需要现代滤波器，例如维纳滤波器、卡尔曼滤波器、自适应滤波器等最正确滤波器。现代滤波器是根据随机信号的一些统计特性，在*种最正确准则下，最大限度地抑制干扰，同时最大限度地回复信号，从而到达最正确滤波的目的。经典数字滤波器从滤波特性上分类，可以分为：低通滤波器、高通滤波器、带通滤波器、带阻滤波器。图1-1 各种理想滤波器的幅频特性数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)数字

14、滤波器和有限长冲激响应(FIR)数字滤波器。IIR 数字滤波器的特征是，具有无限持续时间冲激响应，需要用递归模型来实现，其差分方程为：式1-2系统函数为： 式1-3设计IIR滤波器的任务就是寻求一个物理上可实现的系统函数H(z)，使其频率响应H(z)满足所希望得到的频域指标，即符合给定的通带截止频率、阻带截止频率、通带衰减系数和阻带衰减系数。第三节数字滤波器的设计要求滤波器的指标常常在频域给出。数字滤波器的频响特性函数一般为复函数，所以通常表示为：式1-4其中，|称为幅频特性函数，(w)称为相频特性函数。幅频特性表示信号通过该滤波器后各频率成分的衰减情况，而相频特性反映各频率通过滤波器后在时间

15、上的延时情况。一般IIR数字滤波器，通常只用幅频响应函数|来描述设计指标，相频特性一般不作要求。IIR滤波器指标参数如下列图所示。图中，p和s分别为通带边界频率和阻带边界频率；1和2分别为通带波纹和阻带波纹；允许的衰减一般用dB数表示，通带所允许的最大衰减dB和阻带允许的最小衰减dB分别为p和s表示：式1-5 式1-6一般要求：当时，；当时，。图1-2低通滤波器的技术要求第二章 IIR数字滤波器设计方法IIR数字滤波器是一种离散时间系统，其系统函数为 式2-1假设MN，当MN时,系统函数可以看作一个IIR的子系统和一个(M-N)的FIR子系统的级联。IIR数字滤波器的设计实际上是求解滤波器的系

16、数和，它是数学上的一种逼近问题，即在规定意义上通常采用最小均方误差准则去逼近系统的特性。如果在S平面上去逼近，就得到模拟滤波器；如果在z平面上去逼近，就得到数字滤波器。第一节IIR数字滤波器的设计步骤IIR数字滤波器的设计一般有两种方法：一个是借助模拟滤波器的设计方法进展。其设计步骤是，先设计模拟滤波器，再按照*种方法转换成数字滤波器。这种方法比拟容易一些，因为模拟滤波器的设计方法已经非常成熟，不仅有完整的设计公式，还有完善的图表供查阅；另外一种直接在频率或者时域进展，由于需要解联立方程，设计时需要计算机做辅助设计。其设计步骤是：先设计过渡模拟滤波器得到系统函数，然后将按*种方法转换成数字滤波

17、器的系统函数。这是因为模拟滤波器的设计方法已经很成熟，不仅有完整设计公式，还有完善的图表和曲线供查阅；另外，还有一些典型的优良滤波器类型可供我们使用。为了保证转换后的稳定且满足技术指标要求，对转换关系提出两点要求:(1) 因果稳定的模拟滤波器转换成数字滤波器，仍是因果稳定的。(2) 数字滤波器的频率相应模仿模拟滤波器的频响特性，s平面的虚轴映射为z平面的单位圆，相应的频率之间呈线性关系。利用模拟滤波器成熟的理论设计IIR数字滤波器的过程是：(1)确定数字低通滤波器的技术指标：通带边界频率、通带最大衰减、阻带截止频率、阻带最小衰减。(2)将数字低通滤波器的技术指标转换成相应的模拟低通滤波器的技术

18、指标。(3)按照模拟低通滤波器的技术指标设计过渡模拟低通滤波器。(4)用所选的转换方法，将模拟滤波器转换成数字低通滤波器系统函数。IIR数字滤波器的设计流程图如下：数字滤波器技术指标 模拟滤波器 技术指标 数字滤波器 模拟滤波器 模拟滤波器 设计方法变换变换图2-1 IIR数字滤波器的设计步骤流程图成熟的模拟滤波器设计方法主要有脉冲响应不变法和双线性变换法。第二节用脉冲相应不变法设计IIR数字滤波器一、设计原理利用模拟滤波器来设计数字滤波器，也就是使数字滤波器能模仿模拟滤波器的特性，这种模仿可以从不同的角度出发。脉冲响应不变法是从滤波器的脉冲响应出发，使数字滤波器的单位脉冲响应序列h(n)模仿

19、模拟滤波器的冲激响应ha(t)，即将ha(t)进展等间隔采样，使h(n)正好等于ha(t)的采样值，满足h(n)=ha(nT)式中,T是采样周期。如果令Ha(s)是ha(t)的拉普拉斯变换，H(z)为h(n)的Z变换，利用采样序列的Z变换与模拟信号的拉普拉斯变换的关系得(式2-2)则可看出，脉冲响应不变法将模拟滤波器的S平面变换成数字滤波器的Z平面，这个从s到z的变换z=esT是从S平面变换到Z平面的标准变换关系式。图2-2脉冲响应不变法的映射关系由2-2式，数字滤波器的频率响应和模拟滤波器的频率响应间的关系为(式2-3)这就是说，数字滤波器的频率响应是模拟滤波器频率响应的周期延拓。正如采样定

20、理所讨论的，只有当模拟滤波器的频率响应是限带的，且带限于折叠频率以时，即(式2-4)才能使数字滤波器的频率响应在折叠频率以重现模拟滤波器的频率响应，而不产生混叠失真，即(式2-5)但是，任何一个实际的模拟滤波器频率响应都不是严格限带的，变换后就会产生周期延拓分量的频谱交叠，即产生频率响应的混叠失真。这时数字滤波器的频响就不同于原模拟滤波器的频响，而带有一定的失真。当模拟滤波器的频率响应在折叠频率以上处衰减越大、越快时，变换后频率响应混叠失真就越小。这时，采用脉冲响应不变法设计的数字滤波器才能得到良好的效果。图2-3脉冲响应不变法中的频响混叠现象对*一模拟滤波器的单位冲激响应ha(t)进展采样，

21、采样频率为fs，假设使fs增加，即令采样时间间隔T=1/fs减小，则系统频率响应各周期延拓分量之间相距更远，因而可减小频率响应的混叠效应。二、脉冲响应不变法优缺点从以上讨论可以看出，脉冲响应不变法使得数字滤波器的单位脉冲响应完全模仿模拟滤波器的单位冲激响应，也就是时域逼近良好，而且模拟频率和数字频率之间呈线性关系=T。因而，一个线性相位的模拟滤波器例如贝塞尔滤波器通过脉冲响应不变法得到的仍然是一个线性相位的数字滤波器。脉冲响应不变法的最大缺点是有频率响应的混叠效应。所以，脉冲响应不变法只适用于限带的模拟滤波器(例如，衰减特性很好的低通或带通滤波器)，而且高频衰减越快，混叠效应越小。至于高通和带

22、阻滤波器，由于它们在高频局部不衰减，因此将完全混淆在低频响应中。如果要对高通和带阻滤波器采用脉冲响应不变法，就必须先对高通和带阻滤波器加一保护滤波器，滤掉高于折叠频率以上的频率，然后再使用脉冲响应不变法转换为数字滤波器。当然这样会进一步增加设计复杂性和滤波器的阶数。第三节双线性变换法设计IIR数字滤波器一、设计原理脉冲响应不变法的主要缺点是产生频率响应的混叠失真。这是因为从S平面到平面是多值的映射关系所造成的。为了克制这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率围压缩到-/T/T之间，再用z=esT转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的-/T/T一条横

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

24、2-10式2-11式2-10与式2-11是S平面与Z平面之间的单值映射关系，这种变换都是两个线性函数之比，因此称为双线性变换式2-6与式2-10的双线性变换符合映射变换应满足的两点要求。首先,把z=ej，可得式2-12即S平面的虚轴映射到Z平面的单位圆。其次，将s=+j代入式2-12，得 式2-13因此 式2-14由此看出，当0时，|z|0时，|z|1。也就是说，S平面的左半平面映射到Z平面的单位圆，S平面的右半平面映射到Z平面的单位圆外，S平面的虚轴映射到Z平面的单位圆上。因此，稳定的模拟滤波器经双线性变换后所得的数字滤波器也一定是稳定的。二、双线性变换法优缺点双线性变换法与脉冲响应不变法相

25、比，其主要的优点是防止了频率响应的混叠现象。这是因为S平面与Z平面是单值的一一对应关系。S平面整个j轴单值地对应于Z平面单位圆一周，即频率轴是单值变换关系。这个关系如式2-10所示，重写如下： 式2-15上式说明，S平面上与Z平面的成非线性的正切关系，如图2-4所示。由图2-4看出，在零频率附近，模拟角频率与数字频率之间的变换关系接近于线性关系；但当进一步增加时，增长得越来越慢，最后当时，终止在折叠频率=处，因而双线性变换就不会出现由于高频局部超过折叠频率而混淆到低频局部去的现象，从而消除了频率混叠现象。图2-5双线性变换法的频率变换关系但是双线性变换的这个特点是靠频率的严重非线性关系而得到的

26、，如式2-12及图2-4所示。由于这种频率之间的非线性变换关系，就产生了新的问题。首先，一个线性相位的模拟滤波器经双线性变换后得到非线性相位的数字滤波器，不再保持原有的线性相位了；其次，这种非线性关系要求模拟滤波器的幅频响应必须是分段常数型的，即*一频率段的幅频响应近似等于*一常数这正是一般典型的低通、高通、带通、带阻型滤波器的响应特性，不然变换所产生的数字滤波器幅频响应相对于原模拟滤波器的幅频响应会有畸变，如图2-5所示。图2-6双线性变换法幅度和相位特性的非线性映射对于分段常数的滤波器，双线性变换后，仍得到幅频特性为分段常数的滤波器，但是各个分段边缘的临界频率点产生了畸变，这种频率的畸变，

27、可以通过频率的预畸变来加以校正。也就是将临界模拟频率事先加以畸变，然后经变换后正好映射到所需要的数字频率上。第三章 IIR滤波器的MATLAB设计MATLAB是矩阵实验室Matri* Laboratory之意。除具备卓越的数值计算能力外，它还提供了专业水平的符号计算，文字处理，可视化建模仿真和实时控制等功能。MATLAB的根本数据单位是矩阵，它的指令表达式与数学,工程中常用的形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完一样的事情简捷得多.在新的版本中也参加了对C,FORTRAN,c+ ,JAVA的支持.可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函

28、数库中方便自己以后调用，此外许多的MATLAB爱好者都编写了一些经典的程序，用户可以直接进展下载就可以用，非常的方便。模拟滤波器的理论和设计方法已开展得相当成熟，且有一些典型的模拟滤波器供我们选择，如巴特沃斯(Butterworth )滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer )滤波器、贝塞尔(Bessel )滤波器等，这些典型的滤波器各有特点。用MATLAB进展数字滤波器的设计的步骤如下:（1） 将设计指标归一化处理。如果采用双线性变换法，还需进展预畸变。（2） 根据归一化频率，确定最小阶数N和频率参数Wn。可供选用的阶数择函数有:buttord，cheblord，ch

29、eb2ord，ellipord等。（3） 运用最小阶数N设计模拟低通滤波器原型。模拟低通滤波器的创立函数有:buttap,cheblap, cheb2ap,ellipap和besselap，这些函数输出的是零极点式形式，还要用zp2tf函数转换成分子分母多项式形式。如果想根据最小阶数直接设计模拟低通滤波器原型，可用butter,chebyl,cheby2,ellip,bessel等函数，只是注意要将函数中的Wn设为1。（4） 根据第2步的频率参数Wn，模拟低通滤波原型转换模拟低通、高通、带通、带阻滤波器，可用函数分别是：lp21p,lp2hp,lp2bp,lp2bs。（5） 运用脉冲响应不变法

30、或双线性变法把模滤波器转数字滤波器，调用的函数是impinvar和bilinear。脉冲响应不变法适用于采样频率大于4倍截止频率的锐截止低通带通滤波器，而双线性变换法适合于相位特性要求不高的各型滤波器。（6） 根据输出的分子分母系数，调用函数buttord计算N和wc，有系数向量可以写出数字滤波器系统函数Z,再用freqz函数验证设计结果。第一节IIR数字滤波器的典型设计法设计巴特沃斯数字低通滤波器和椭圆数字低通滤波器，要求通带边界频fp=2.1kHZ，通带最大衰减Rp=0.5dB；阻带边界频率fs=8kHZ，阻带最小衰减Rs=30dB,采样频率为Fs=20kHZ。低通巴特沃斯滤波器设计步骤如

31、下：（1） 确定阶数N。 取N=4（2） 求极点， 归一化低通原型系统函数为由N=4直接查表得到：极点：归一化低通滤波器系统函数为式中，0.0000,0.0999,0.1914,0.0252（3） 将去归一化最终得到通过计算可以总结出过程太麻烦，而且容易出错，结果不直观。下面用M程序设计来实现用脉冲响应不变法设计的巴特沃斯数字低通滤波器的M程序如下：fp=2100;fs=8000;Fs=20000;Rp=0.5;Rs=30;T=1/Fs; %设计指标W1p=fp/Fs*2;W1s=fs/Fs*2;%求归一化频率N,Wn=buttord(W1p,W1s,Rp,Rs,s);%确定butterwor

32、th的最小介数N和频率参数Wnz,p,k=buttap(N); %设计模拟低通原型的零极点增益参数bp,ap=zp2tf(z,p,k); %将零极点增益转换成分子分母参数bs,as=lp2lp(bp,ap,Wn*pi*Fs);%将低通原型转换为模拟低通bz,az=impinvar(bs,as,Fs); %用脉冲响应不变法进展模数变换 sys=tf(bz,az,T); %给出传输函数H(Z)H,W=freqz(bz,az,512,Fs); %生成频率响应参数subplot(2,1,1);plot(W,20*log10(abs(H); %绘制幅频响应grid on; %加坐标网格*label(频率

33、/Hz); ylabel(振幅/dB);subplot(2,1,2); plot(W,abs(H); grid on; *label(频率/Hz); ylabel(振幅/H);运行后的波形如下：图3-1 典型滤波器在Matlab上运行波形运行结果：N =4bz = 0.0000 0.0999 0.1914 0.0252az=1.0000 -1.4336 1.0984 -0.4115 0.0627可以得出：只需编程，结果非常直观。双线性变换法设计步骤如下：1首先写出该滤波器的系统函数2利用双线性变换法转换，数字滤波器的系统函数为用双线性变换法设计椭圆数字低通滤波器的M程序如下： fs=20000

34、; wp=2*pi*2100/fs; ws=2*pi*8000/fs; Rp=0.5; Rs=30; Ts=1/fs; Wp=2/Ts*tan(wp/2);Ws=2/Ts*tan(ws/2); %按频率转换公式进展转换 N,Wn=ellipord(Wp,Ws,Rp,Rs,s); %计算模拟滤波器的最小阶数 z,p,k=ellipap(N,Rp,Rs);%设计模拟原型滤波器 Bap,Aap=zp2tf(z,p,k); %零点极点增益形式转换为传递函数形式 b,a=lp2lp(Bap,Aap,Wn); %低通转换为低通滤波器的频率转化 bz,az=bilinear(b,a,fs); %运用双线性变

35、换法得到数字滤波器传递函数 H,f=freqz(bz,az,512,fs); subplot(2,1,1); plot(f,20*log10(abs(H); title(N=2 频率响应); grid on; *label(频率/Hz); ylabel(振幅/dB);subplot(2,1,2); plot(f,abs(H); grid on; *label(频率/Hz); ylabel(振幅/H);图3-2 双线性变换法设计的椭圆数字低通滤波器运行波形运行结果：N=2bz= 0.1213 0.1662 0.1213az= 1.0000 -0.9889 0.4218第二节IIR数字滤波器的直接

36、设计法除了典型设计以外，MATLAB信号处理工具箱提供了几个直接设计IIR数字滤波器的函数，直接调用就可以设计滤波器，这为设计通用滤波器提供了方便。主要有以下几种方法：1.零极点累试法这种方法也称为零极点累试法。在确定零极点位置时要注意：(1)极点必须位于z平面单位圆，保证数字滤波器因果稳定；(2)复数零极点必须共轭成对，保证系统函数有理式的系数是实的。2.在频域利用幅度平方误差最小法直接设计IIR数字滤波器3.在时域直接设计IIR数字滤波器设计Butterworth滤波器用函数butter()，可以设计低通、高通、带通和带阻的数字和模拟滤波器，其特性是通带的幅度响应最大限度的平滑，但损失了截

37、止频率处的下降斜度。设计Chebyshev I型滤波器用函数chebyl()。可以设计低通、高通、带通和带阻的数字和模拟Chebyshev I型滤波器，其通带为等波纹，阻带为单调。Chebyshev I型滤波器的下降斜度比II型大，但其代价目是通带波纹较大。设计Chebyshev II型滤波器用函数cheby2()。可以设计低通、高通、带通和带阻的数字和模拟Chebyshev II型滤波器，其通带为单调，阻带等波纹。Chebyshev II型滤波器的下降斜度比I型小，但其阻带波纹较大。设计椭圆滤波器用函数ellip()，与chebyl, cheby2类似，可以设计低通、高通、带通和带阻的数字和

38、模拟滤波器。与Butterworth和chebyshev滤波器相比，ellip函数可以得到下降斜度更大的滤波器，得通带和阻带均为等波纹。一般情况下，椭圆滤波器能以最低的阶实现指定的性能指标。在使用各类滤波器函数时应当注意以下重点:(1)、阶数和固有频率的选择:N,Wn=buttord(Wp,Ws,Rp,Rs)可得到符合要求性质的滤波器的最小阶数N以及数字Butterworth滤波器的固有频率Wn(即3dB )。设计的要在通带的衰减不超过Rp,在阻带的衰减不小于Rs，通带和阻带有截止频率分别是Wp, Ws,它们是归一化的频率，围是0, 1，对应弧度。(2)、有关滤波器设计当中的频率归一化问题:信

39、号处理工具箱中经常使用的频率是Nyquist频率，它被定义为采样频率的一半，在滤波器的阶数选择和设计中的截止频率均使用Nyquist频率进展归一化处理。例如对于一个采样频率为1000 Hz的系统，400Hz的归一化即为400/500=0.8。归一化频率的围在0, 1之间。如果要将归一化频率转换为角频率，则将归一化频率乘以;如果要将归一化频率转换为Hz,则将归一化频率乘以采样频率的一半。(3)、设计一个N阶的低通Butterworth滤波器使用函数B,A=butter(N, Wn)，返回滤波器系数矩阵B,A。其中固有频率Wn必须是归一化频率。它的最大值是采样频率的一半。Fs缺省时默认为2Hz。如

40、果Wn=Wl,W2是一个两元素的向量，则函数将设计出一个2N阶的带通滤波器，通带为W1,W2。设计Chebyshev I型和Chebyshev II型数字低通滤波器，要求通带边界频率fp=2.1kHZ，通带最大衰减Rp=0.5dB；阻带边界频率fs=8kHZ，阻带最小衰减Rs=30dB,采样频率为Fs=20kHZ。Chebyshev I型的M程序如下：Fs=20000; %抽样频率20KHzFlp=2100; Fls=8000;Wp=2*Flp/Fs; %归一化的通带截止频率Ws=2*Fls/Fs; %归一化的阻带截止频率Rp=0.5; %通带最大衰减单位：dBRs=30; %阻带最小衰减单位

41、：dBN,Wn=cheb1ord(Wp,Ws,Rp,Rs); %返回最小阶数和截止频率b,a=cheby1(N,Rp,Wn); %返回H(z)的分子分母系数hw,w=freqz(b,a);subplot(2,1,1);plot(w/pi,20*log10(abs(hw);grid on;*label(/);ylabel(幅度dB)title(切比雪夫I型幅频响应);subplot(2,1,2);plot(w/pi,abs(hw);grid on;*label(/);ylabel(幅度H);运行后的波形如下：图3-3 Chebyshev I型运行波形运行结果：N=2b= 0.1007 0.201

42、4 0.1007a= 1.0000 -0.9872 0.4140Chebyshev II型M程序如下：Fs=20000; %抽样频率20KHzFlp=2100; Fls=8000;Wp=2*Flp/Fs; %归一化的通带截止频率Ws=2*Fls/Fs; %归一化的阻带截止频率Rp=0.5; %通带最大衰减单位：dBRs=30; %阻带最小衰减单位：dBN,Wn=cheb2ord(Wp,Ws,Rp,Rs); %返回最小阶数和截止频率b,a=cheby2(N,Rs,Wn); %返回H(z)的分子分母系数hw,w=freqz(b,a);subplot(2,1,1);plot(w/pi,20*log1

43、0(abs(hw);grid on;*label(/);ylabel(幅度dB)title(切比雪夫II型幅频响应);subplot(2,1,2);plot(w/pi,abs(hw);grid on;*label(/);ylabel(幅度H);运行后的波形：图3-4 Chebyshev II型运行波形运行结果：N=2b= 0.2357 0.4241 0.2357a= 1.0000 -0.2996 0.1950从频率响应图中可以看出:巴特沃斯滤波器具有单调下降的幅频特性，通带平滑;切比雪I型滤波器的幅频特性在通带有波动，阻带单调;chebyshev II型滤波器的幅频特性在阻带有波动，通带单调;

44、椭圆滤波器的选择性相对前三种是最好的，下降斜度比拟大，通带和阻带均为等波纹，同样的性能指标，椭圆滤波器可以最低的阶数来实现。这样根据不同的要求可以选用不同类波器。第三节 FDATool介绍和界面设计FDATool(Filter Design&Analysis Tool)是MATLAB信号处理工具箱里专用的滤波器设计分析工具，MATLAB 6.O以上的版本还专门增加了滤波器设计工具箱(Filter Design Toolbo*)。FDATool可以设计几乎所有的常规滤波器，包括FIR和IIR的各种设计方法。它操作简单，方便灵活。FDATool。界面总共分两大局部，一局部是Design Filte

45、r。在界面的下半部，用来设置滤波器的设计参数；另一局部则是特性区，在界面的上半局部，用来显示滤波器的各种特性。Design Filter局部主要分为：Filter Type(滤波器类型)选项，包括Lowpass(低通)、Highpass(高通)、Bandpass(带通)、Bandstop(带阻)和特殊的 FIR滤波器。Design Method(设计方法)选项，包括IIR滤波器的Butterwotth(巴特沃思)法、Chebyshev Type I(切比雪夫I型)法、Chebyshev Type II(切比雪夫II型)法、Elliptic(椭圆滤波器)法和FIR滤波器的Equiripple法、

46、Least-Squares(最小乘方)法、 Window(窗函数)法。Filter Order(滤波器阶数)选项，定义滤波器的阶数，包括Specify Order(指定阶数)和Minimum Order(最小阶数)。在Specify Order中填入所要设计的滤波器的阶数(N阶滤波器，Specify Order=N-1)。如果选择Minimum Order，则MATLAB根据所选择的滤波器类型自动使用最小阶数。Frequency Specifications选项，可以详细定义频带的各参数，包括采样频率和频带的截止频率。它的具体选项由Filter Type选项和Design Method选项决定

47、。例如Bandpass(带通)滤波器需要定义Fstop1(下阻带截止频率)、Fpass1(通带下限截止频率)、Fpass2 (通带上限截止频率)、Fstop2(上阻带截止频率)，而Lowpass(低通)滤波器只需要定义Fstop1、Fpass1。采用窗函数设计滤波器时，由于过渡带是由窗函数的类型和阶数所决定，所以只需定义通带截止频率，而不必定义阻带参数。Magnitude Specifications选项，可以定义幅值衰减的情况。例如设计带通滤波器时，可以定义Wstop1(频率Fstop1处的幅值衰减)、Wpass (通带围的幅值衰减)、Wstop2(频率Fstop2处的幅值衰减)。当采用窗函

48、数设计时，通带截止频率处的幅值衰减固定为6db，所以不必定义。 Window Specifications选项，中选取采用窗函数设计，该选项可定义，它包含了各种窗函数。基于MATLAB 函数的滤波器设计完成后，需要对已设计的滤波器的频率响应要进展校核。要得到幅频、相频响应特性，运算量也是很大的。而利用MATLAB 强大的信号处理界面工具进展计算机辅助设计，可以快速有效地设计数字滤波器，大简化了计算量。第四节FDATOOL设计IIR数字滤波器我们将以一个IIR 滤波器的设计实例来具体说明使用MATLAB 工具箱的方便。要求设计Chebyshev1 型低通滤波器满足指标，首先在Response T

49、ype 中选择Bandpass高通滤波器，然后在下面的Desigh Method 中选择IIR 类型，并且指定Filter Order 项中的阶数SpecifyOrder=10，由于是设计chebyshev 滤波器，其下面Option 就不必选择。然后在FrequencySpecifications 中选择Unit 为Hz，给出采样频率Fs=1000，；最后在Magnitude Specifications 中选择Unit 为dB，Apass=1。设置完成后点击Design Filter 即可得到所设计的IIR 滤波器。通过菜单项选择项Analysis 可以在特性区看到所设计的幅频响应、相频响

50、应、冲击响应和零极点配置等特性。设计完成后将结果保存，封装为filterl.fda 文件，下面是运用FDATool 设计滤波图3-5FDA设计主界面图3-6相频特性图3-7冲击响应特性从以上这些我们可以清晰明了的看到设计的该滤波器的幅频，相频和冲击响应特性：我们能够看到，图3-3中Chebshev1 型滤波器具有等纹波的通频带、阻带衰减的特性，总之，使用FDATOOL 工具包设计和分析滤波器，是非常方便易行的，而且交互性良好，不需要极其复杂的程序编制就可以实现。而且我们在最后的仿真模型设计中将载入这一模块进展信号滤波。第五节SIMULINK 仿真IIR滤波器SIMULINK是一个进展动态系统建

51、模、仿真和综合分析的集成软件包。它可以处理的系统包括：线性、非线性系统；离散、连续及混合系统；单任务、多任务离散事件系统。本章通过调用Simulink中的功能模块构成IIR数字滤波器的仿真框图，在仿真过程中，可以双击各功能模块，随时改变参数，获得不同状态下的仿真结果。首先建立信号传输仿真模块，先翻开 SIMULINK 工具箱，并且建立一个 Model，在这个空白 Model 中进展环节库及框图的建立，在 SIMULINK 菜单下找到 Source，双击 Source 图标，将正弦信号源和常量信号源拉到 Model 中，然后分别找到 DSP 离散正弦信号源,乘法 器和加法器以及示波器，然后把需要

52、用来设计 IIR 数字滤波器的模块 FDATOOL 都拉到 Model 中，把环节都布好后，把各环节的端口按框图连接起来。下面就来仿真 3.5节号处理的模型，以框图的方式直观地现实信号处理的过程，信号源为 f= sin(2*pi*t*40)+sin(3*pi*t*50)+sin(2*pi*t*100)，在 DSP离散信号模块中设定参数分别为 40HZ，75HZ，100HZ,在 simulink环境中，将 3.5中封装好的 filterl.fda滤波器文件载入 FDATOOL中，生成仿真框图和滤波效果如下图。图3-8仿真模型图3-9滤波效果结果显示，仿真模型中显示 40HZ,75HZ,100HZ

53、三个不同频率的正弦信相加后的波形，干扰后的混合信号是一个夹杂多频率的信号，我们所需要的有用信号40HZ在传输过程中不可防止的受到其他频率信号干扰，我们的目的就是得到所需要的 40HZ频段信号，尽可能滤除干扰信号，滤波效果的图中就是滤波后信号与原始抽样信号的的波形图，经滤波后，信号周期约为 0.025s,即频率为 40HZ，说明滤波效果很好。通过分析，说明 SIMULINK中各种非常有用的工具箱不仅对于设计 IIR数字滤波器非常有用，而且对于整个信号仿真处理具有相当可视化的效果，从仿真的效果来看，是到达了技术指标的要求。总结通过这几个实验，对设计数字滤波器的整个过程有了很好的掌握。其中对双线性变

54、换法，巴特沃斯设计模拟滤波器的运用，也更加清楚了。通过对数字带通滤波器的设计，熟悉了MATLAB的运行环境，初步掌握了MATLAB语言在数字信号处理中一些根本库函数的调用和编写根本程序等应用；熟悉了滤波器设计的一般原理，对滤波器有了一个感性的认识；学会了数字高通滤波器设计的一般步骤；加深了对滤波器设计中产生误差的原因以及双线性变换法优缺点的理解和认识。总之，使理论联系了实际，稳固并深化了对课本根本知识的认识和理解，使理论得以升华。致 本论文是在雷教师的亲切关心和悉心指导下完成的。他严肃的科学态度，严谨的治学精神，精益求精的工作作风，深深地感染和鼓励着我。从课题的选择到工程的最终完成，雷教师都始

55、终给予我细心的指导和不懈的支持。一直以来，雷教师不仅在学业上给我以精心指导，同时还在思想、生活上给我以无微不至的关心，在此谨向雷教师致以诚挚的意和崇高的敬意。在此，我还要感在一起愉快的度过毕业设计生活的同学，正是由于你们的帮助和支持，我才能克制一个一个的困难和疑惑，直至本文的顺利完成。特别感我的同学，她们对本课题做了不少工作，给予我不少的帮助。在论文即将完成之际，我的心情无法平静，从开场进入课题到论文的顺利完成，有多少可敬的师长、同学、朋友给了我无言的帮助，在这里请承受我诚挚的意!最后我还要感培，你们!参考文献1 倪养华数字信号处理与实现：交通大学，19982 飞思科技产品研发中心MATLAB

56、7辅助信号处理技术与应用：电子工业， 20053 楼顺天，博苗基于MATLAB的系统分析与设计一信号处理：电子科技大学，1998.4 祥，娜MATLAB仿真技术与应用：清华大学，20035 胡广书数字信号处理、理论、算法与实现M ：清华大学，19976 希林，肖明清一种LabWindows/CVI与MATLAB混合编程的实现方法J.微计算机信息，20057波MATLAB信号处理：电子工业，20068施阳等MATLAB语言工具箱：西北工业大学，19999丁玉美数字信号处理电子科技大学 第二版完毕语本文在对滤波器实际设计中用到了多种设计方法，与其他高级语言的程序设计相比，MATLAB 环境下可以更

57、方便、快捷地设计出具有IIR 数字滤波器，节省大量的编程时间，提高编程效率，且参数的修改也十分方便，还可以进一步进展优化设计。但是运用函数设计法，整个过程的运算量是很大的。设计阶数较高的IIR 滤波器时，计算量更大，设计过程中要改变参数或滤波器类型时都要重新计算。它需要反复的实验，而且需要设计者凭借经历设定参数，平时所要设计的数字滤波器，阶数和类型并不一定是完全给定的，很多时候要根据设计要求和滤波效果不断地调整，以到达设计的最优化。在这种情况下，滤波器设计就要进展大量复杂的运算，单纯的靠公式计算和编制简单的程序很难在短时间完成。因此，基于对象的信号处理工具FDATOOL 以及SPTOOL 界面设计滤波器，可以有效的的解决这一问题，它不仅减少了设计复杂度，而且还为用户提供了一个便于分析和观察的界面。本文利用Matlab的信号处理工具及其SIMULINK仿真工具包箱能够方便快捷地设计和实现各种滤波器，使信号波形更加直观，并且用FDATOOL工具箱以及专门用于滤波器设计来设计IIR数字滤波器，充分利用了MATLAB 的交互性好的特点，而且最后用到了SIMULINK 仿真工具箱来实现信号传输和滤波器的设计，将滤波器的设计置于一个新的平台，这对于研究信号的传输和处理有着极其重要的作用。z.