结合双线性变换法设计一个数字切比雪夫带通IIR滤波器

课程设计任务书学生姓名：蒋立豪专业班级：通信1303指导教师：魏勤工作单位：信息工程学院题目:利用MATLAB仿真软件系统结合双线性变换法设计一个数字切比雪夫带通IIR滤波器。初始条件：1MATLAB编程的基础知识2《信号与系统》滤波器的相关知识3《数字信号处理》的相关知识要求完成的主要任务:（包括课程设计工作量及技术要求，以及说明书撰写等具体要求）1、课程设计工作量：1周。2、技术要求：1）设计一个数字滤波器，要求运用切比雪夫原理和双线性变换法2）具有带通的功能，能够滤除小频率的信号和大频率的信号3）能够让频带内的信号通过4)确定设计方案，按功能模块的设计程序，写出总体程序，并阐述基本原理。5)查阅至少5篇参考文献。按《武汉理工大学课程设计工作规范》要求撰写设计报告书。全文用A4纸打印，图纸应符合绘图规范。时间安排：序号阶段内容所需时间1方案设计1天2软件设计2天3系统调试1天4答辩1天合计5天指导教师签名：年月日系主任（或责任教师）签名：年月日摘要面对庞杂繁多的原始信号,如何提取所需信号、抑制不需要的信号这就需要使用滤波器。滤波器的作用主要是选择所需频带的信号内容而抑制不需要的其他频带的信号内容。数字滤波器因其精度高、可靠性好、灵活性大等优点,在语音信号处理、信号频谱估计、信号去噪、无线通信中的数字变频以及图像处理等工程实际应用中都很广泛。MATLAB是“矩证实验室”（MATrixLABoratoy）的缩写，它是一种以钜阵运算为基础的交互式程序语言，专门针对科学﹑工程计算机绘图的需求。与其他计算机语言相比，其特点是简洁和智能化，适应科技专业人员的思维方式和书写习惯，使得编程和调试效率大大提高本次课程设计将完成一个数字切比雪夫带通IIR滤波器的设计，利用双线性变换和无限冲激响应IIR原理完成设计，并利用MATLAB进行仿真。关键字：滤波器、数字信号处理、噪声、IIR、带通、MATLABAbstractFacedwithahugevarietyoforiginalsignal,howtoextractthedesiredsignal,suppressthesignaldoesnotneedtousethefilter.Thefunctionofthefilteristoselectthesignalcontentofthedesiredfrequencybandandsuppressthesignalcontentoftheotherfrequencybandwhichisnotneeded.Digitalfilterbecauseoftheadvantageofitshighaccuracy,goodreliabilityandgreatflexibility,inspeechsignalprocessing,signalspectrumestimation,signaltonoise,wirelesscommunicationinthedigitalfrequencyconversionandimageprocessingetc.inactualengineeringapplicationisverywidespread.Matlabistheabbreviationof"momentoflaboratorymedicine"(matrixLABoratoy)itisakindoftohugeoperationalmatricesbasedinteractiveprogramminglanguage,specificallyfortheneedsofscienceandengineeringcomputergraphics.Comparedwithothercomputerlanguages,itischaracterizedbyitssimplicityandintelligence,thewayofthinkingandwritinghabits,whichmakestheprogramminganddebuggingefficiencygreatlyimproved.ThecurriculumdesignwillbecompletedadigitalChebyshevbandpassIIRfilterdesign,usingbilineartransformandinfiniteimpulseresponseIIRprincipleofthecompletionofdesign,andissimulatedwithMATLAB.Keywords:filter,digitalsignalprocessing,noise,IIR,band-pass,MATLAB目录摘要IAbstractII1前言12数字滤波器22.1数字滤波器介绍22.2IIR数字滤波器32.2.1IIR数字滤波器的特点32.2.2IIR数字滤波器与FIR数字滤波器的区别43切比雪夫滤波器43.1概述43.2切比雪夫滤波器的种类54双线性变换法85IIR数字滤波器设计原理与方法125.1IIR数字滤波器设计原理125.2IIR数字滤波器设计的基本方法146IIR带通滤波器的MATLAB设计166.1IIR带通滤波器的设计步骤166.2IIR带通滤波器设计166.4仿真结果20滤波器性能仿真20滤波器性能验证217.用FDATOOL设计数字滤波器268总结与体会309致谢32参考文献33附录：源程序341前言在现代通信系统中，由于信号中经常混有各种复杂成分，所以很多信号的处理和分析都是基于滤波器而进行的。数字滤波器可以理解为是一个计算程序或算法，将代表输入信号的数字时间序列转化为代表输出信号的数字时间序列，并在转化过程中，使信号按预定的形式变化。数字滤波器有多种分类，根据数字滤波器冲激响应的时域特征，可将数字滤波器分为两种，即无限长冲激响应(IIR)滤波器和有限长冲激响应(FIR)滤波器。但是，传统的数字滤波器的设计使用繁琐的公式计算，改变参数后需要重新计算，从而在设计滤波器尤其是高阶滤波器时工作量很大。利用MATLAB信号处理箱(SignalProcessingToolbox)可以快速有效地实现数字滤波器的设计与仿真。本文设计一个IIR数字带通滤波器。IIR数字滤波器具有无限宽的冲激响应，与模拟滤波器相匹配，所以IIR滤波器的设计可以采取在模拟滤波器设计的基础上进一步变换的方法。其设计方法主要有经典设计法、直接设计法和最大平滑滤波器设计法。在对滤波器实际设计时，整个过程的运算量是很大的。设计完成后对已设计的滤波器的频率响应要进行校核，要得到幅频、相频响应特性，运算量也是很大的。平时所要设计的数字滤波器，阶数和类型并不一定是完全给定的，很多时候要根据设计要求和滤波效果不断地调整，以达到设计的最优化。在这种情况下，滤波器设计就要进行大量复杂的运算，单纯的靠公式计算和编制简单的程序很难在短时间内完成。利用MATLAB强大的计算功能进行计算机辅助设计，可以快速有效地设计数字滤波器，大大地简化了计算量。滤波器的功能是让一定频率范围内的信号通过，而将此频率范围之外的信号加以抑制或使其急剧衰减。当干扰信号与有用信号不在同一频率范围之内，可使用滤波器有效的抑制干扰。用LC网络组成的无源滤波器在低频范围内有体积重量大，价格昂贵和衰减大等缺点，而用集成运放和RC网络组成的有源滤波器则比较适用于低频，此外，它还具有一定的增益，且因输入与输出之间有良好的隔离而便于级联。由于大多数反映生理信息的光电信号具有频率低、幅度小、易受干扰等特点，因而RC有源滤波器普遍应用于光电弱信号检测电路中数字滤波器的实用型式很多，大略可分为有限冲激响应型和无限冲激响应型两类，可用硬件和软件两种方式实现。在硬件实现方式中，它由加法器、乘法器等单元所组成，这与电阻器、电感器和电容器所构成的模拟滤波器完全不同。数字信号处理系统很容易用数字集成电路制成，显示出体积小、稳定性高、可程控等优点。数字滤波器也可以用软件实现。软件实现方法是借助于通用数字计算机按滤波器的设计算法编出程序进行数字滤波计算。带通滤波器是指能通过某一频率范围内的频率分量、但将其他范围的频率分量衰减到极低水平的滤波器，与带阻滤波器的概念相对。一个模拟带通滤波器的例子是电阻-电感-电容电路(RLCcircuit)。这些滤波器也可以用低通滤波器同高通滤波器组合来产生.2数字滤波器2.1数字滤波器介绍数字滤波器由数字乘法器、加法器和延时单元组成的一种算法或装置。数字滤波器的功能是对输入离散信号的数字代码进行运算处理，以达到改变信号频谱的目的。数字滤波器是一个离散时间系统（按预定的算法，将输入离散时间信号转换为所要求的输出离散时间信号的特定功能装置）。基本工作原理是利用离散系统特性对系统输入信号进行加工和变换,改变输入序列的频谱或信号波形,让有用频率的信号分量通过,抑制无用的信号分量输出。与模拟滤波器相比,数字滤波器除了具有数字信号处理的固有优点外,还有滤波精度高(与系统字长有关)、稳定性好(仅运行在0与l两个电平状态)、灵活性强等优点。应用数字滤波器处理模拟信号时，首先须对输入模拟信号进行限带、抽样和模数转换。数字滤波器输入信号的抽样率应大于被处理信号带宽的两倍，其频率响应具有以抽样频率为间隔的周期重复特性，且以折叠频率即1/2抽样频率点呈镜像对称。为得到模拟信号，数字滤波器处理的输出数字信号须经数模转换、平滑。数字滤波器具有高精度、高可靠性、可程控改变特性或复用、便于集成等优点。数字滤波器在语言信号处理、图像信号处理、医学生物信号处理以及其他应用领域都得到了广泛应用。数字滤波器有低通、高通、带通、带阻和全通等类型。它可以是时不变的或时变的、因果的或非因果的、线性的或非线性的。数字滤波器根据其冲激响应函数的时域特性，可分为两种，即无限长冲激响应(IIR)数字滤波器和有限长冲激响应(FIR)数字滤波器。2.2IIR数字滤波器IIR(InfiniteImpulseResponse)数字滤波器，又名“无限脉冲响应数字滤波器”，或“递归滤波器”。递归滤波器，也就是IIR数字滤波器，顾名思义，具有反馈，一般认为具有无限的脉冲响应。2.2.1IIR数字滤波器的特点IIR滤波器有以下几个特点：1、封闭函数IIR数字滤波器的系统函数可以写成封闭函数的形式。2、IIR数字滤波器采用递归型结构IIR数字滤波器采用递归型结构，即结构上带有反馈环路。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。由于运算中的舍入处理，使误差不断累积，有时会产生微弱的寄生振荡。3、借助成熟的模拟滤波器的成果IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，其设计工作量比较小，对计算工具的要求不高。在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。4、需加相位校准网络IIR数字滤波器的相位特性不好控制，对相位要求较高时，需加相位校准网络。2.2.2IIR数字滤波器与FIR数字滤波器的区别1、单位响应IIR数字滤波器单位响应为无限脉冲序列，而FIR数字滤波器单位响应为有限的；FIR滤波器，也就是“非递归滤波器”，没有引入反馈。这种滤波器的脉冲响应是有限的。2、幅频特性IIR数字滤波器幅频特性精度很高，不是线性相位的，可以应用于对相位信息不敏感的音频信号上；FIR数字滤波器的幅频特性精度较之于IIR数字滤波器低，但是线性相位，就是不同频率分量的信号经过FIR滤波器后他们的时间差不变，这是很好的性质。3、实时信号处理FIR数字滤波器是有限的单位响应也有利于对数字信号的处理，便于编程，用于计算的时延也小，这对实时的信号处理很重要。3切比雪夫滤波器3.1概述目的：构造一个模拟低通滤波器。为了从模拟滤波器出发设计IIR数字滤波器，必须先设计一个满足技术指标的模拟滤波器，亦即要把数字滤波器的指标转换成模拟滤波器的指标，因此必须先设计对应的模拟原型滤波器。模拟滤波器的理论和设计方法己发展得相当成熟，且有一些典型的模拟滤波器供我们选择，如巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等，这些典型的滤波器各有特点。这里介绍切比雪夫滤波器。切比雪夫滤波器是在通带或阻带上频率响应幅度等波纹波动的滤波器。在通带波动的为“I型切比雪夫滤波器”，在阻带波动的为“II型切比雪夫滤波器”。切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快，但频率响应的幅频特性不如后者平坦。切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动。这种滤波器来自切比雪夫多项式，因此得名，用以记念俄罗斯数学家巴夫尼提·列波维奇·切比雪夫。3.2切比雪夫滤波器的种类根据频率响应曲线波动位置的不同，切比雪夫滤波器可以分为以下两种：1、I型切比雪夫滤波器：在通带（或称“通频带”）上频率响应幅度等波纹波动的滤波器称为“I型切比雪夫滤波器”；2、 II型切比雪夫滤波器：在阻带（或称“阻频带”）上频率响应幅度等波纹波动的滤波器称为“II型切比雪夫滤波器”。3.3特点切比雪夫滤波器特点：误差值在规定的频段上等波纹变化。切比雪夫滤波器在过渡带比巴特沃斯滤波器的衰减快，但频率响应的幅频特性不如后者平坦。切比雪夫滤波器和理想滤波器的频率响应曲线之间的误差最小，但是在通频带内存在幅度波动。巴特沃兹滤波器在通带内幅度特性是单调下降的，如果阶次一定，则在靠近截止处，幅度下降很多，或者说，为了使通带内的衰减足够小，需要的阶次N很高，为了克服这一缺点，采用切比雪夫多项式来逼近所希望的。切比雪夫滤波器的在通带范围内是等幅起伏的，所以在同样的通常内衰减要求下，其阶数较巴特沃兹滤波器要小。切比雪夫滤波器的振幅平方函数为(3-1)式中Ωc为有效通带截止频率，表示与通带波纹有关的参量，值越大通带不动愈大。VN（x）是N阶切比雪夫多项式，定义为(3-2)切比雪夫滤波器的振幅平方特性如图所示：N为偶数，cos2()=1，得到min，，（3-3）N为奇数，cos2(，得到max，（3-4）图3-1切比雪夫滤波器的振幅平方特性有关参数的确定:a、通带截止频率：预先给定;b、与通带波纹有关的参数,通带波纹表示成（3-5）所以,,给定通带波纹值分贝数后，可求得。ｃ、阶数N:由阻带的边界条件确定。、A2为事先给定的边界条件,即在阻带中的频率点处,要求滤波器频响衰减到1/A2以上。(3-6)(3-7)(3-8)（3-9）因此，要求阻带边界频率处衰减越大，要求N也越大，参数N，给定后，查阅有关模拟滤波器手册，就可求得系统函数Ha(s)。4双线性变换法目的：将模拟带通滤波器转换成数字带通滤波器为了克服冲激响应法可能产生的频率响应的混叠失真，这是因为从S平面到Ｚ平面是多值的映射关系所造成的。为了克服这一缺点，可以采用非线性频率压缩方法，将整个频率轴上的频率范围压缩到-π/T～π/T之间，再用z=esT转换到Z平面上。也就是说，第一步先将整个S平面压缩映射到S1平面的-π/T～π/T一条横带里；第二步再通过标准变换关系z=es1T将此横带变换到整个Z平面上去。这样就使S平面与Z平面建立了一一对应的单值关系，消除了多值变换性，也就消除了频谱混叠现象，映射关系如图4-1。图4-1双线性变换的映射关系为了将S平面的整个虚轴jΩ压缩到S1平面jΩ1轴上的-π/T到π/T段上，可以通过以下的正切变换实现（4-1）式中,T仍是采样间隔。当Ω1由-π/T经过0变化到π/T时，Ω由-∞经过0变化到+∞，也即映射了整个jΩ轴。将式（4-1）写成(4-2)将此关系解析延拓到整个S平面和S1平面，令jΩ=s，jΩ1=s1，则得(4-3)再将S1平面通过以下标准变换关系映射到Z平面z=es1T从而得到S平面和Z平面的单值映射关系为：（4-4)（4-5)式（4-4）与式（4-5）是S平面与Z平面之间的单值映射关系，这种变换都是两个线性函数之比，因此称为双线性变换。双线性变换符合映射变换应满足的两点要求。首先,把z=ejω，可得（4-6）即S平面的虚轴映射到Z平面的单位圆。其次，将s=σ+jΩ代入式（4-5），得因此（4-7）由此看出，当σ<0时，|z|<1；当σ>0时，|z|>1。也就是说，S平面的左半平面映射到Z平面的单位圆内，S平面的右半平面映射到Z平面的单位圆外，S平面的虚轴映射到Z平面的单位圆上。因此，稳定的模拟滤波器经双线性变换后所得的数字滤波器也一定是稳定的。双线性变换法优缺点：双线性变换法与脉冲响应不变法相比，其主要的优点是避免了频率响应的混叠现象。这是因为S平面与Z平面是单值的一一对应关系。S平面整个jΩ轴单值地对应于Z平面单位圆一周，即频率轴是单值变换关系。这个关系重写如下：（4-8）上式表明，S平面上Ω与Z平面的ω成非线性的正切关系，如图4-2所示。由图4-2看出，在零频率附近，模拟角频率Ω与数字频率ω之间的变换关系接近于线性关系；但当Ω进一步增加时，ω增长得越来越慢，最后当Ω→∞时，ω终止在折叠频率ω=π处，因而双线性变换就不会出现由于高频部分超过折叠频率而混淆到低频部分去的现象，从而消除了频率混叠现象。图4-2双线性变换法的频率变换关系但是双线性变换的这个特点是靠频率的严重非线性关系而得到的，由于这种频率之间的非线性变换关系，就产生了新的问题。首先，一个线性相位的模拟滤波器经双线性变换后得到非线性相位的数字滤波器，不再保持原有的线性相位了；其次，这种非线性关系要求模拟滤波器的幅频响应必须是分段常数型的，即某一频率段的幅频响应近似等于某一常数（这正是一般典型的低通、高通、带通、带阻型滤波器的响应特性），不然变换所产生的数字滤波器幅频响应相对于原模拟滤波器的幅频响应会有畸变，如图4-3所示。图4-3双线性变换法幅度和相位特性的非线性映射对于分段常数的滤波器，双线性变换后，仍得到幅频特性为分段常数的滤波器，但是各个分段边缘的临界频率点产生了畸变，这种频率的畸变，可以通过频率的预畸来加以校正。也就是将临界模拟频率事先加以畸变，然后经变换后正好映射到所需要的数字频率上。5IIR数字滤波器设计原理与方法5.1IIR数字滤波器设计原理利用MATLAB信号处理工具箱中的滤波器设计和分析工具(FDATool)可以很方便地设计出符合应用要求的未经量化的IIR数字滤波器。需要将MATLAB设计出的IIR数字滤波器进一步分解和量化，从而获得可用FPGA实现的滤波器系数。IIR数字滤波器的最通用的方法是借助于模拟滤波器的设计方法。模拟滤波器设计已经有了相当成熟的技术和方法，有完整的设计公式，还有比较完整的图表可以查询，因此设计数字滤波器可以充分利用这些丰富的资源来进行。对于IIR数字滤波器的设计具体步骤如下：(1)按照一定的规则将给出的数字滤波器的技术指标转换为模拟低通滤波器的技术指标。(2)根据转换后的技术指标设计模拟低通滤波器G(s)(G(s)是低通滤波器的传递函数)。(3)再按照一定的规则将G(s)转换成H(z)(H(z)是数字滤波器的传递函数)。若设计的数字滤波器是低通的，上述的过程可以结束，若设计的是高通、带通或者是带阻滤波器，那么还需要将高通、带通或带阻数字滤波器的技术指标转换为低通模拟滤波器的技术指标，然后设计出低通G(s)，再将G(s)转换为H(z)。IIR滤波器运算结构通常由延时、乘以系数和相加等基本运算组成，可以组合成直接型、正准型、级联型、并联型四种结构形式，都具有反馈回路。由于运算中的舍入处理，使误差不断累积，有时会产生微弱的寄生振荡。IIR数字滤波器在设计上可以借助成熟的模拟滤波器的成果，如巴特沃斯、契比雪夫和椭圆滤波器等，有现成的设计数据或图表可查，其设计工作量比较小，对计算工具的要求不高。在设计一个IIR数字滤波器时，我们根据指标先写出模拟滤波器的公式，然后通过一定的变换，将模拟滤波器的公式转换成数字滤波器的公式。IIR数字滤波器的相位特性不好控制，对相位要求较高时，需加相位校准网络。在MATLAB下设计IIR滤波器可使用Butterworth函数设计出巴特沃斯滤波器，使用Cheby1函数设计出契比雪夫I型滤波器，使用Cheby2设计出契比雪夫II型滤波器。与FIR滤波器的设计不同，IIR滤波器设计时的阶数不是由设计者指定，而是根据设计者输入的各个滤波器参数（截止频率、通带滤纹、阻带衰减等），由软件设计出满足这些参数的最低滤波器阶数。在MATLAB下设计不同类型IIR滤波器均有与之对应的函数用于阶数的选择。IIR数字滤波器是一种离散时间系统，其系统函数为（5-1）假设M≤N，当M＞N时,系统函数可以看作一个IIR的子系统和一个(M-N)的FIR子系统的级联。IIR数字滤波器的设计实际上是求解滤波器的系数和，它是数学上的一种逼近问题，即在规定意义上（通常采用最小均方误差准则）去逼近系统的特性。如果在S平面上去逼近，就得到模拟滤波器；如果在z平面上去逼近，就得到数字滤波器。设计高通、带通、带阻等数字滤波器通常可以归纳为如图所示的两种常用方法。图5-1数字滤波器设计的两种方法方法1：首先设计一个模拟原型低通滤波器，然后通过频率变换成所需要的模拟高通、带通或带阻滤波器，最后再使用冲激不变法或双线性变换成相应的数字高通、带通或带阻滤波器。方法2：先设计一个模拟原型低通滤波器，然后采用冲激响应不变法或双线性变换法将它转换成数字原型低通滤波器，最后通过频率变换把数字原型低通滤波器变换成所需要的数字高通、带通或带阻滤波器。方法一的缺点是，由于产生混叠失真，因此不能用冲激不变法来变换成高通或阻带滤波器，故一般采用第二种方法进行设计。本课程设计先构造一个切比雪夫模拟低通滤波器，然后将模拟低通滤波器转换成模拟带通滤波器，最后利用双线性变换将模拟带通滤波器转换成数字带通滤波器。5.2IIR数字滤波器设计的基本方法IIR数字滤波器的设计方法有两类,一类是借助于模拟滤波器的设计方法设计出模拟滤波器,利用冲激响应不变法或双线性变换法转换成数字滤波器,再用硬件或软件实现；另一类是直接在频域或时域中进行设计,设计时需要计算机作辅助工具。随着MATLAB软件尤其是MATLAB的信号处理工作箱的不断完善,不仅数字滤波器的计算机辅助设计有了可能,而且还可以使设计达到最优化。IIR数字滤波器设计的基本步骤如下:(1)根据任务,确定性能指标。在设计一个滤波器之前,首先根据工程实际的需要确定滤波器的技术指标如：边界频率：ωp，ωs，ωc；阻带最小衰减As和通带最大衰减Rp；(2)将数字滤波器的技术指标转换成模拟滤波器指标。利用冲激响应不变法与双线性变换法进行频率间的转换，主要是边界频率Wp与Ws的转换。（3)用模拟滤波器设计方法得到模拟滤波器的传输函数Ha(s)；可借助巴特沃斯(Butterworth)滤波器、切比雪夫(Chebyshev)滤波器、椭圆(Cauer)滤波器、贝塞尔(Bessel)滤波器等，这些滤波器都有严格的设计公式、现成的曲线和图表供设计人员使用。(4)映射实现。利用双线性变换法将模拟滤波器Ha(s)转换成数字滤波器H(z)。(5)用有限精度算法实现这个系统函数H(z)(包括选择运算结构、选择合适的字长、有效数字处理方法)。(6)用适当的软、硬件技术实现。包括采用通用计算机软件、数字滤波器硬件或者软硬件结合，确定DF采用的结构将会影响其精度、稳定性、经济性及运算速度等很多重要性质。6IIR带通滤波器的MATLAB设计6.1IIR带通滤波器的设计步骤根据以上IIR数字滤波器设计方法，下面运用双线性变换法基于MATLAB设计一个IIR带通滤波器。IIR带通滤波器的设计框图如下：模拟低通滤波器的设计模拟频带变换：低通变高通、带通等滤波器变换：模拟变为数字IIR滤波器模拟低通滤波器的设计模拟频带变换：低通变高通、带通等滤波器变换：模拟变为数字IIR滤波器设计的IIR带通滤波器是从低通变换过来的，利用的是双线性变换以及切比雪夫II滤波器的原型，其具体的设计流程如图所示。首先根据题目要求确定带通滤波器的技术指标，先要进行频率的预畸变，并且归一化频率，再设计出切比雪夫II模拟低通滤波器，并求出其阶数等相关参数。其次利用双线性变换法设计数字带通滤波器，，再调用函数进行双线性变换，并求出分子、分母的系数向量。最后通过画图求出其幅频响应、相频响应、幅度特性曲线与零极点，并画出波形图。最后进行验证，看所设计的滤波器能否达到要求的指标，若能达到，则说明该滤波器设计符合要求。6.2IIR带通滤波器设计根据以上IIR数字滤波器设计方法，下面运用双线性变换法基于MATLAB设计一个IIR带通滤波器，其中带通的中心频率为wp0=0.55π，;通带截止频率wp1=0.45π,wp2=0.65π;通带最大衰减Ap=1dB;阻带最小衰减As=40dB;阻带截止频率ws2=0.75π(1)确定性能指标在设计带通滤波器之前,首先根据工程实际的需要确定滤波器的技术指标:通带截止频率wp1=0.45π,wp2=0.65π;阻带截止频率ws1=0.3π，ws2=0.75π；阻带最小衰减As=40dB和通带最大衰减Ap=1dB;中心频率wp0=0.55π。(2)频率预畸变用Ω=2/T*tan(w/2)对带通数字滤波器H(z)的数字边界频率预畸变,得到带通模拟滤波器H(s)的边界频率主要是通带截止频率Wp1,Wp2;阻带截止频率Ws1，Ws2的转换。双线性变换法一般T=2s。通带截止频率Wp1=(2/T)*tan(wp1/2)Wp2=(2/T)*tan(wp2/2)阻带截止频率Ws1=(2/T)*tan(ws1/2)Ws2=(2/T)*tan(ws2/2)(3)模拟带通性能指标转换成模拟低通性能指标BW=Wp2-Wp1；%带通滤波器的通带宽度W0=Wp1*Wp2；WP=1；%归一化处理WS=WP*(W0^2-Ws1^2)/(Ws1*BW)；(4)模拟低通滤波器的构造借助切比雪夫(Chebyshev)滤波器得到模拟低通滤波器的传输函数Ha(s)。(5)模拟低通滤波器转换成模拟带通滤波器调用lp2bp函数将模拟低通滤波器转化为模拟带通滤波器。(6)模拟带通滤波器转换成数字带通滤波器利用双线性变换法将模拟带通滤波器Ha(s)转换成数字带通滤波器H(z)。(7)输入信号检验滤波器性能输入不同频率的正弦波，观察输出波形，检验滤波器性能。6.3Malab程序（1）数字滤波器部分clear%数字滤波器的技术指标wp1=0.45*pi;wp2=0.65*pi;wp0=0.55*pi;ws1=0.3*pi;ws2=0.75*pi;Ap=1;As=40;T=2;%带通到低通的频率变换Wp1=(2/T)*tan(wp1/2);Wp2=(2/T)*tan(wp2/2);Ws1=(2/T)*tan(ws1/2);Ws2=(2/T)*tan(ws2/2);Wp=(2/T)*tan(wp0/2);BW=Wp2-Wp1;%带通滤波器的通带宽度W0=Wp1*Wp2;WP=1;%归一化处理WS=WP*(W0^2-Ws1^2)/(Ws1*BW);%切比雪夫模拟低通原型滤波器设计[N,Wn]=cheb1ord(WP,WS,Ap,As,'s');[B1,A1]=cheby1(N,Ap,Wn,'s');%模拟低通原型滤波器幅频特性曲线(dB)[h1,w1]=freqs(B1,A1);subplot(3,2,1);plot(w1/pi,20*log(abs(h1)));gridon;xlabel('w(rad)');ylabel('|H(jw)|.dB');title('模拟低通滤波器幅频特性曲线');%由模拟低通原型滤波器变换为模拟带通滤波器[B2,A2]=lp2bp(B1,A1,Wp,BW);[h2,w2]=freqs(B2,A2);%模拟带通滤波器幅频特性曲线(dB)subplot(3,2,2);plot(w2,20*log(abs(h2)));axis([0,3,-400,50]);gridon;xlabel('w(rad)');ylabel('|H(jw)|.dB');title('模拟带通滤波器幅频特性曲线');%双线性变换：由模拟滤波器向数字滤波器的变换[B3,A3]=bilinear(B2,A2,0.5);[h,w]=freqz(B3,A3,64);phz=unwrap(angle(h));%数字带通滤波器幅频响应曲线subplot(3,2,3);plot(w/pi,abs(h));xlabel('w(rad)');ylabel('|H(z)|');title('数字带通滤波器幅频特性曲线');%数字带通滤波器幅频响应曲线(dB)subplot(3,2,4);plot(w/pi,20*log10(abs(h)));axis([-1,2,-250,50]);xlabel('w(rad)');ylabel('|H(z)|.dB');title('数字带通滤波器幅频特性曲线');%数字带通滤波器相频特性曲线(dB)subplot(3,2,5);plot(w/pi,phz);gridon;xlabel('w(rad)');ylabel('H(z)');title('数字带通滤波器相频特性曲线');（2）输入正弦波检验性能部分%输入正弦波波验证滤波器特性n=0:600;t=n/11000;x1=2*sin(2*pi*2750*t);%正弦波信号figure;subplot(121);plot(x1);gridon;%500Hz正弦波波形axis([0,10*pi,-5,5]);xlabel('t(s)');ylabel('x1');title('正弦波信号');y1=filter(B3,A3,x1);%数字滤波函数输出subplot(122);plot(y1);gridon;%数字滤波器输出波形axis([0,10*pi,-3,3]);xlabel('f(hz)');ylabel('y');title('数字滤波器输出波形');注:应输入一系列不同频率的正弦波进行验证，只需将x1=2*sin(2*pi*2750*t)中的2750频率值改变即可，这里取500HZ、2750HZ和5000HZ进行验证。6.4仿真结果6.4.1滤波器性能仿真源程序设计了模拟低通滤波器、模拟带通滤波器与数字带通滤波器等滤波器，对各部分滤波器的性能仿真如下，下面五个图分别为模拟低通原型滤波器幅频特性曲线、模拟带通滤波器幅频特性曲线、数字带通滤波器幅频响应曲线（有两个）、数字带通滤波器相频特性曲线，可以看到各部分滤波器波形基本满足设计要求。图6-3滤波器性能仿真6.4.2滤波器性能验证为了验证滤波器性能，可以输入一些频率不同的正弦波，观察通阻状态，因为通带截止频率wp1=0.45π,wp2=0.65π;阻带截止频率ws1=0.3π，ws2=0.75π，中心频率wp0=0.55π，所以可以将测试频率分别定在通带内，阻带内，和截止区，通带内分别输入2750HZ,2650HZ，阻带内分别输入1700HZ,和3650HZ，截止区分别输入1000HZ和4000HZ进行验证仿真后的结果如下：1)通带内验证仿真（1）输入为2750HZ的正弦波信号图6-4输入2750HZ正弦波波形（2）输入为2650HZ的正弦波信号图6-5输入2650HZ正弦波波形从仿真结果可以看出，频率在通带内的信号在误差允许的范围内能够无损耗的通过滤波器。2）阻带内验证仿真（1）输入为1700HZ的正弦波信号图6-6输入1700Hz正弦波波形（2）输入为3650HZ的正弦波信号图6-7输入3650Hz正弦波波形从仿真可以看出，频率在阻带内的信号能够通过滤波器，但有明显的损耗。3）截止区内验证仿真（1）输入为1000HZ的正弦波信号图6-8输入1000Hz正弦波波形（2）输入为4000HZ的正弦波信号图6-9输入4000Hz正弦波波形由仿真图可以看出，频率在截止区的信号，在误差允许的范围内可近似的认为被滤波器滤除掉了4）仿真图总分析 频率在通带内的信号在误差允许的范围内能够无损耗的通过滤波器。频率在阻带内的信号能够通过滤波器，但有明显的损耗频率在截止区的信号几乎无输出，在误差允许的范围内可近似的认为被滤波器滤除掉了。由于软件的精度要求以及所调用函数的特性，在误差允许范围内，上述的滤波器特性的仿真结果基本符合题目的要求，并且在一定的程度上设计的指标远远好于题目的要求。因此，所设计的滤波器达到了题目的要求。7.用FDATOOL设计数字滤波器 还可以可以使用Matlab中的FDATOOL组件设计滤波器。具体操作如下。在CommandWindow中输入“fdatool”回车，图7-1fdatool窗口图7-2幅值响应图图7-3相频响应图图7-4阶跃响应图图7-5

噪声信号功率谱分析图7-6零极点图可以看到利用双线性变换法得到的数字带通滤波器的幅频特性，符合切比雪夫I型滤波器的特点，通带内等波纹，阻带内单调的，且比较“陡峭”，特性较好。图4.6为其相频特性图，利用了函数unwrap函数使得形成的相位连续。图4.7为IIR带通滤波器系统函数对应的零极点分布图8总结与体会此次Matlab课程设计是在一定的理论基础之上进行的，在先修课程《信号与系统》中，大量有MATLAB设计方面的知识，已经接触了解了MATLAB的操作和功能，有了一定的基础。而以前做实验时，对它都是半知半解。只学习了MATLAB的一些简单基本的语句功能，根本没有认识到MATLAB在信号处理方面的强大功能。这次的课程设计认真的了解了这方面的知识，对MATLAB程序有了整体的认识。首先，在信号滤波系统中，有时因为模拟滤波器阶数太高，硬件占用空间太大为某些仪器的实现设置了障碍，而对于一些窄带情况下的低通滤波器用模拟手段往往很难实现。在这些情况下，数字滤波器将会是一个很好的解决办法。MATLAB信号处理工具箱提供了丰富而简便的设计、实现FIR和IIR的方法，使原来繁琐的程序设计简化成函数的调用，特别是滤波器的表达方式和滤波器之间的相互转换显得十分简便。其次，IIR数字滤波器的设计和模拟滤波器的设计有着紧密的关系。通常要先设计出适当的模拟滤波器，再通过一定的频带变换把它转换成为所需的数字IIR滤波器。此外，任何数字信号处理系统中也还不可避免地用到模拟滤波器，比如A/D变换器前的抗混叠滤波器和D/A变换器后的平滑滤波器，因此模拟滤波器设计也是很重要的。最后，在比较设计滤波器的方法上应该明确其技术指标以及某些参数的实际意义。比如本文用双线性变换法设计数字带通滤波器时，必须先将频率归一化，并且进行频率预畸变，然后设计模拟滤波器，再利用频率变换法将模拟低通变为模拟带通，最后经过双线性变换法将模拟带通变换为数字带通。如果不进行频率预畸变，那么设计出来的带通滤波器的幅频特性与相频特性将会产生很严重的畸变，使设计的结果不满足给定的要求，在实际中会造成很严重的危害。通过本次课程设计，我深深地感受到了自己所学到知识多么的匮乏，纸上学来终觉浅，觉知此事要躬行，通过这五天的学习，我明白了只学好课本上的知识是不够的，要通过图书馆和互联网等各种渠道来扩充自己的知识。通过自己想办法解决在实验过程中遇到的问题，调试修改程序来锻炼自己，扩充知识面和提高自己对知识的运用能力9致谢在这次课程设计的撰写过程中，我得到了许多人的帮助。

首先我要感谢我的老师在课程设计上给予我的指导、提供给我的支持和帮助，这是我能顺利完成这次报告的主要原因，更重要的是老师帮我解决了许多技术上的难题，让我能把系统做得更加完善。在此期间，我不仅学到了许多新的知识，而且也开阔了视野，提高了自己的设计能力。

其次，我要感谢帮助过我的同学，他们也为我解决了不少我不太明白的设计中的难题。同时也感谢学院为我提供良好的做课程设计的环境。

参考文献[1]VinayK.Ingle主编：《数字信号处理及其MATLAB实现》，电子工业出版社，1998年出版。[2]SanjitK.Miltra编著：《DigitalSignalProcessingLaboratoryUsingMatlab》，McGraw-Hill出版社，2000年出版。[3]刘泉,阙大顺主编《数字信号处理原理与实现》电子工业出版社.2005[4]景振毅，张泽兵，董霖.北京：中国铁道出版社，202