dsp论文-二阶带通滤波器 22

PAGE1信号处理课程设计题目二阶带通滤波器学院通信与电子工程学院专业班级电子097班学生姓名方长鑫指导教师王发智2012年12月25日PAGEPAGE22摘要滤波是将信号中特定波段频率进行滤除的操作，是抑制和防止干扰的一项重要措施。其在信号处理中使用非常广泛。所谓滤波器，就是对特定频率的频点或该频点以外的频率进行有效滤除的电路。模拟滤波器的理论和设计方法相当成熟，设计时可以选用典型的模拟滤波器，如巴特沃斯滤波器、契比雪夫滤波器、椭圆滤波器、贝赛尔滤波器等，它们有严格的设计公式和现成的曲线和图表可供参考，巴特沃斯滤波器具有单调下降的幅频特性；切比雪夫滤波器的幅频特性在通带和阻带内有波动，可以提高选择性。MATLAB语言及其工具箱，其将一个优秀软件包的易用性与可靠性，通用性和专业性，以及一般目的的应用和高深的专业应用近乎完美地集成在一起。并凭借其功能强大，技术先进和应用之深广，使其逐渐成为国际性的计算标准，为世界各地超过20万名科学家和工程师所采用。利用MATLAB软件，可以使滤波器的设计比较简单易于实现，因此，本文借助MATLAB软件来进行应用研究。本文是基于模拟低通滤波器的原型，通过巴特沃斯、切比雪夫I、切比雪夫II型的具体算法来实现模拟带阻滤波器的设计，并通过和直接利用MATLAB软件中已经成熟的函数实现的方法相比较，以验证理论学习过程中的相关算法的正确性，同时也显示该软件在仿真中的强大功能，表明模拟滤波器的设计已经非常成熟和高效率，而且结果稳定、可靠。我们坚信：随着MATLAB软件在信息技术中的不断推广和使用，将会使科技工程人员从各类烦琐的计算工作中解脱出来，有时间斯考和研究更重要的问题，这必将为科技的进步起到巨大的推动作用。关键字：MATLAB软件，巴特沃斯滤波器，切比雪夫滤波器

AbstractFilteringistosignalaparticularbandfrequencyfilteringoperation,istopreventinterferencesuppressionandoneofthemostimportantmeasures.Thesignalprocessingiswidelyusedin.Theso-calledfilter,isonaparticularfrequencyfrequencypointorthefrequencypointoutsideofthefrequencyeffectivelyfiltercircuit.Analogfiltertheoryanddesignmethodisquitemature,designcanchoosetypicalanalogfilter,suchasbutterworthfilter,chebyshevfilters,ellipticfilter,beisaierfilterandsoon,theyhavestrictdesignformulaandready-madecurveandchartsforreference,butterworthfilterhasthemonotonedeclineamplitudefrequencycharacteristics;Chebyshevfiltertheamplitude-frequencycharacteristicsofthepassbandandstop-bandfluctuationsin,canimprovetheselective.MATLABlanguageanditstoolbox,whichwouldbeagoodpackageofusabilityandreliability,generalityandprofessional,andgeneralpurposeapplicationsandadvancedprofessionalapplicationalmostperfectlyintegratesintogether.Andwithitspowerful,advancedtechnologyandapplicationofthedepth,makeitsbecomeinternationalcalculationstandardforallovertheworldmorethan200000scientistsandengineersuse.UsingMATLABsoftware,canmakethedesignoffilterrelativelysimpleeasytorealize,therefore,thispapermakesuseoftheMATLABsoftwaretocarryouttheapplicationresearch.Thispaperisbasedonthesimulationoflowpassfilterprototype,throughthebutterworth,chebyshevI,chebyshevtypeIIconcretealgorithmtosimulatethebandstopfilterdesign,andthroughanddirectuseofMATLABsoftwarehasbeenmaturefunction．Keywords：MATLAB，Butterworthfilter，Chebyshevfilter

目录TOC\o"1-3"\u摘要 2Abstract……………………..31模拟带通滤波器设计 51.1概述 51.2滤波器的分析 51.3工作原理 52二阶低通滤波器的设计及制作 72.1二阶带通滤波器设计 72.1.1原理图设计分析（见附录） 72.1.2电路PCB（见附录） 72.1.3中心频率的计算 82.1.4芯片介绍 93基于MatlabGUI的数字滤波器设计思路及实现 113.1GUI界面设计概述 113.1.1、控件风格和外观 123.1.2

对象的常规信息 123.1.3控件回调函数的执行 133.1.4、控件当前状态信息 133.1.5IIR滤波器简介 153.1.6IIR滤波器设计原理 153.1.7IIR滤波器设计思想 153.1.8利用完全设计法设计数字滤波器的步骤： 163.2运行和结果显示 164总结 20谢辞 21附录 22第1章模拟带通滤波器设计1.1概述带通滤波器是指能通过某一频率范围内的频率分量、但将其他范围的频率分量衰减到极低水平的滤波器，与带阻滤波器的概念相对。一个模拟带通滤波器的例子是电阻-电感-电容电路(RLCcircuit)。这些滤波器也可以用低通滤波器同高通滤波器组合来产生.无源滤波器通常由RLC元件组成，有源滤波器有电阻，电容和运算放大器组成。1.2滤波器的分析滤波器按元件分类分为无源滤波器和有源滤波器两种：（1）无源滤波器：由电感L、电容C及电阻R等无源元件组成。（2）有源滤波器：一般由集成运放与RC网络构成，它具有体积小、性能稳定等优点，同时，由于集成运放的增益和输入阻抗都很高，输出阻抗很低，故有源滤波器还兼有放大与缓冲作用，利用有源滤波器可以突出有用频率的信号，衰减无用频率的信号，抑制干扰和噪声，以达到提高信噪比或选频的目的，因而有源滤波器被广泛应用于通信、测量及控制技术中的小信号处理。滤波器的阶数越高，幅频特性越好，但RC网络节数越多，元件参数计算越繁琐，电路的调试越困难。任何高阶滤波器都可由一阶和二阶滤波器级联而成。对于N阶为偶数的高阶滤波器，可以由N/2节二阶滤波器级联而成；而N为奇数的高阶滤波器可以由(N-1)/2节二阶滤波器和一节一阶滤波器级联而成，因此一阶滤波器和二阶滤波器是高阶滤波器的基础。1.3工作原理一个理想的滤波器应该有一个完全平坦的通带，例如在通带内没有增益或者衰减，并且在通带之外所有频率都被完全衰减掉，另外，通带外的转换在极小的频率范围完成。实际上，并不存在理想的带通滤波器。滤波器并不能够将期望频率范围外的所有频率完全衰减掉，尤其是在所要的通带外还有一个被衰减但是没有被隔离的范围。这通常称为滤波器的滚降现象，并且使用每十倍频的衰减幅度dB来表示。通常，滤波器的设计尽量保证滚降范围越窄越好，这样滤波器的性能就与设计更加接近。然而，随着滚降范围越来越小，通带就变得不再平坦—开始出现“波纹”。这种现象在通带的边缘处尤其明显，这种效应称为吉布斯现象。除了电子学和信号处理领域之外，带通滤波器应用的一个例子是在大气科学领域，很常见的例子是使用带通滤波器过滤最近3到10天时间范围内的天气数据，这样在数据域中就只保留了作为扰动的气旋。在频带较低的剪切频率f1和较高的剪切频率f2之间是共振频率，这里滤波器的增益最大，滤波器的带宽就是f2和f1之间的差值。第2章二阶低通滤波器的设计及制作2.1二阶带通滤波器设计由上概论可知，任何高阶滤波器都可由一阶和二阶滤波器级联而成。对于N阶为偶数的高阶滤波器，可以由N/2节二阶滤波器级联而成；而N为奇数的高阶滤波器可以由(N-1)/2节二阶滤波器和一节一阶滤波器级联而成，因此一阶滤波器和二阶滤波器是高阶滤波器的基础。可是难度也会越加困难。所以，现在选择二阶电路来设计此带通滤波器。2.1.1原理图设计分析图2.1电路原理图2.1.2芯片介绍本次试验主要用到的元器件有LM358，其内部包括有两个独立的、高增益、内部频率补偿的双运算放大器，适合于电源电压范围很宽的单电源使用，也适用于双电源工作模式，在推荐的工作条件下，电源电流与电源电压无关。它的使用范围包括传感放大器、直流增益模块和其他所有可用单电源供电的使用运算放大器的场合。1.1.内部频率补偿直流电压增益高(约100dB)单位增益频带宽(约1MHz)2.电源电压范围宽：单电源(3—30V)；双电源(±1.5一±15V)低功耗电流，适合于电池供电3.低输入偏流低输入失调电压和失调电流共模输入电压范围宽，包括接地差模输入电压范围宽，等于电源电压范围4.输出电压摆幅大(0至Vcc-1.5V)图2.4LM358管脚图2.1.3中心频率的计算（a）电路图(b)幅频特性图2.2电路图以及幅频特性

工作原理：这种滤波器的作用是只允许在某一个通频带范围内的信号通过，而比通频带下限频率低和比上限频率高的信号均加以衰减或抑制。典型的带通滤波器可以从二阶低通滤波器中将其中一级改成高通而成。如图1（a）所示。

电路性能参数1、通带增益

2、中心频率

3、通带宽度

4、选择性

Q=f0/B5、本次滤波器的设计指标：此电路的优点是改变Rf和R4的比例就可改变频宽而不影响中心频率。由上所有公式，可以确定本次试验的设计指标。

（1）、通带中心频率

f0=11khz（2）、品质因数Q=1.1（3）、有此可得带宽B=f0/Q=10KHZ（4）、通带中心频率处的电压放大倍数：A0=2.42.2滤波器滤波效果为测试带通滤波器的性能，我用试验数据做出曲线图，与用MATLAB软件生成的曲线做比较。并且在实际设计的滤波器中，也有一些的电路干扰，所以曲线不会很平滑。表2.3带通滤波器测试数据电压峰峰值频率电压峰峰值频率电压峰峰值频率0.650.928.90.9612.80.635.30.9459.20.94513.10.6655.60.9659.50.93513.40.6855.90.989.80.9213.70.7056.20.9910.10.895140.7256.5110.40.87514.30.756.8110.70.8514.60.7757.11110.82514.90.8057.4111.30.80515.20.827.7111.60.7815.50.84580.99511.90.75515.80.8758.30.9912.20.7250.97512.50.69516.4图2.4带通滤波器实际测试幅频特性曲线图2.5带通滤波器理想幅频特性曲线第3章基于MatlabGUI的数字滤波器设计思路及实现3.1GUI界面设计概述GUI设计面板是GUI设计工具应用的平面，面板上部提供了菜单和常用工具按钮，左边提供了多种如命令按钮、单选按钮、可编辑文本框、静态文本框、弹出式菜单等。进行设计时,首先单击面板左边所需的控件,然后在右边的图形界面编辑区中再次单击某一恰当位置,这时将在该位上为图形界面添加相应的控件。一个图形界面的完成,除了设计其外观外,还有相当的一部分是通过属性的设来完成的。因此在设置这些属性时,要注意下面几个常用又很重要的属性设置：3.1.1、控件风格和外观（1）BackgroundColor：设置控件背景颜色，使用[RGB]或颜色定义。（2）CData：在控件上显示的真彩色图像，使用矩阵表示。（3）ForegroundColor：文本颜色（4）String属性：控件上的文本，以及列表框和弹出菜单的选项。（5）Visible：控件是否可见。3.1.2

对象的常规信息

（1）Enable属性：表示此控件的使能状态，设置为on”，表示可选，为“off”时则表示不可选。（2）Style：控件对象类型。（3）Tag：控件表示（用户定义）。（4）TooltipString属性：提示信息显示。当鼠标指针位于此控件上时，显示提示信息。（5）UserData：用户指定数据。（6）Position：控件对象的尺寸和位置。（7）Units：设置控件的位置及大小的单位（8）有关字体的属性，如FontAngle，FontName等。

3.1.3控件回调函数的执行（1）BusyAction：处理回调函数的中断。有两种选项：即Cancel：取消中断事件，queue：排队（默认设置）。（2）ButtonDownFcn属性：按钮按下时的处理函数。（3）CallBack属性：是连接程序界面整个程序系统的实质性功能的纽带。该属性值应该为一个可以直接求值的字符串，在该对象被选中和改变时，系统将自动地对字符串进行求值。（4）CreateFcn：在对象产生过程中执行的回调函数。（5）DeleteFcn：删除对象过程中执行的回调函数。（6）Interruptible属性：指定当前的回调函数在执行时是否允许中断，去执行其他的函数。3.1.4、控件当前状态信息（1）ListboxTop：在列表框中显示的最顶层的字符串的索引（2）Max：最大值。（3）Min：最小值。（4）Value：控件的当前值。应用MATLAB制作这样一个过程是非常方便的,我们可以通过GUI操作来看到。该环境下要设计一个界面友好的仿真软件，一般应完成以下两个步骤：（1）GUI界面设计。主要是通过不同的文本框、按钮等许多工具的使用，设计出一个图形用户界面。要清楚这个图形界面的功能是什么，即在图形界面上的操作会引发什么样的结果。（2）回调函数的设计。用户应根据设计好的图形界面的功能，针对各个不同的图形对象来编写出能够实现该功能的函数代码，确保这个图形界面能够完成所预定的功能。图3.1.1GUI设计图形界面其中，设计的分析菜单包括滤波器的幅频特性、相频特性、阶跃响应、脉冲相应和零相应和零计的分析菜单包括滤波器的幅频特性、相频特性、阶跃响应、脉冲相应和零极点子菜单，计算菜单包括阶数、零点、极点、增益、转移函数(分子)和转移函数(分母)子菜单。图3.1.2运行显示界面3.1.5IIR滤波器简介IIR滤波器即无限长脉冲响应（InfiniteImpulseResponse）滤波器，它具有反馈，一般认为具有无限的脉冲响应。IIR滤波器为非线性相位（双线性变换法），对于非线性相位会造成的影响。IIR数字滤波器方便简单，但它相位的线性，要采用全通网络进行相位校正。图象处理以及数据传输，都要求信道具有线性相位特性.有限冲击响应(FIR)滤波器具有很好的线性相位特性，因此越来越受到广泛的重视。3.1.6IIR滤波器设计原理滤波器的设计实质上是寻找一个既能物理实现，又能满足给定频率特性指标要求的系统传输函数。IIR滤波器一般采用递归型的结构，系统的输入输出服从N阶差分方程：（1-2）相应的传递函数：（1-3）设计IIR数字滤波器是要确定传点函数a、b或零极点c、d、A，使滤波器的频率特性满足给定的性能指标要求。设计原理主要包括两个方面：一是根据设计指标，先设计出相应的模拟滤波器再通过脉冲响应不变法或双线性变换法转换成对应的数字滤波器；二是选择一种优准则，如最小均方准则，再在先最误差此准则下求出滤波器传输函数的系数。根据设计理论，在MATLAB环境下设计IIR数字滤波器主要有四种方法：一是典型设计法,二是完全设计法；三是最优设计法；四是工具设计法。由于完全设计法程序简单，因此使用完全设计法设计滤波器。3.1.7IIR滤波器设计思想IIR滤波器设计思想是：利用已有的模拟滤波器设计理论，首先设计指标设计一个合适的模拟滤器，然后再通过脉冲响应不变法或双线性变换法，完成从模拟到数字的转换。常用的模拟滤波器有巴特沃斯滤波器、切比雪夫滤波器、椭圆滤波器、贝塞尔滤波器等，这些滤波器各有特点，供不同设计要求选用。滤波器的模拟数字变换，通常是复变函数的映射变换，也必须满足一定的要求。由于数字滤波器传输函数只与频域的相对值有关，故在设计时可先将滤波器设计归一化处理，设采样频率为Fs,归一化频率的计算公式是：（1-6）3.1.8利用完全设计法设计数字滤波器的步骤：将设计指标归一化处理。(1)根据归一化频率，确定最小阶数N和频率参数Wn。可供选用的阶数选择函数：buttord、cheb1ord、cheb2ord，ellipord.(2)运用最小阶数N设计模拟滤波器原型。根据最小阶数直接设计模拟低通滤波器原型，用到的函数有：buttord、cheb1ord、cheb2ord，ellip和bessel。如[B、A]=butter(N、Wn、“”)设计‘type’型巴特沃斯滤波器filter。N为滤波器阶数，Wc为截止频率，type决定滤波器类型，type=high,设计高通IIR滤波器，ftype=stop,设计带阻IIR滤波器。(3)再用freqz函数验证设计结果。IIR滤波器设计编程实现例如选择设计IIR的Butterworth低通滤波器，其中Fs=5000Hz,Fp1=500Hz,Fs=1200Hz.Rp=2dB；Rs=20Db.程序和效果图（图9）所示3.2运行和结果显示1运行和结果显示在图形界面下,按“运行”,出现如图15所示界面，选择要设计的滤波器选项，编辑框中输人要求设计的数字滤波器的性能指标,,选择“Run”命令按钮,出现如图所示的模拟低通原型滤波器幅频响应和相频响应曲线,最小阶次。当选择了“AutoRun”按钮，在“滤波器类型选择”旁的下拉菜单框中选择其他类型的滤波器会立刻出现此类型的滤波器的图形。例如选择设计IIR的Butterworth低通滤波器，其Fs=1000Hz，Fp1=100Hz，Fs1=300Hz,Rp=3dB,Rs=20dB,n=4,显示如图16图3.2.1Butterworth低通滤波器显示不同滤波器的显示结果1、选择设计IIR低通滤波器，其Fs=1000Hz，Fp1=100Hz，Fs1=300Hz,Rp=3dB,Rs=20dB,n=4,利用四种不同滤波器进行设计，如图17、18、19、20显示：图3.2.2Butterworth低通滤波器图3.2.3ChebyshevⅠ低通滤波器图3.2.4ChebyshevⅡ低通滤波器图3.2.5Ellipse低通滤波器2、选择设计FIR带通滤波器，其Fs=4000Hz，Fp1=900Hz，Fp2=1300Hz，Fs1=600Hz,Fs2=1500Hz,Rp=1dB,Rs=40dB,n=20,Rp=3dB,Rs=20dB,n=4,利用四种不同滤波器进行设计，如图21、22、23、24、25、26显示：图3.2.6Boxar窗带通滤波器图3.2.7Bartlett窗带通滤波器图3.2.8Blackman窗带通滤波器图3.2.9Hanning窗带通滤波器图3.2.10Hamming窗带通滤波器图3.2.11Kaiser窗带通滤波器第4章总结本次实训主要分为两部分，第一部分是模拟低通滤波器的硬件设计，主要的实训要求是设计中心频率为我们的学号后两位数乘以1K，通频带为固定值10000，这样可以确保我们的Q值为学号的后两位除以10.0。所以中心频率为11K(HZ).Q值为1.1。第二部分是软件的设计，通过MATLAB中的GUI界面设计IIR和FIR滤波器，可以设计虑波器为低通高通、带通和带阻等。第一周我们进行了硬件设计的实训，首先有设计方案，然后再确定电路图进行仿真做板，于是我就去借阅书籍，查找资料，从而了解模拟带通滤波器的工作原理以及实现方法，并结合以前所学的《模拟电子技术》中学到的有关截止频率的知识以及本学期学习的《数字信号处理》，最终能够确定自己的方案，在查资料时，使我知道滤波器的阶数越高，幅频特性越好,所以我选择了二阶的巴特沃斯的设计方案，设计方案确定后，就是进行电路的设计与仿真，有源滤波器主要由电容和电阻组成，所以在电路设计中我用了两个滑动变阻器，以便在硬件做好后可以对其进行调试。后来在同学的帮助和老师的指导下我成功的完成了硬件设计的实训。第二周就到了软件的学习与设计了，由于只是在上个学期的信号与系统中学到一点有关MATLAB的一些内容，所以对这个软件只是在这个学期的DSP实验中有所了解，但是GUI界面并没有接触过，这就需要在网上查找资料辅助自己学习GUI界面的设计所以对于我来说还是有一定的难度的，本身自己对软件方面就不是很在行，因此这个设计对我来说有一定的挑战性。不过就像事物都是有内在联系的，在设计编程的时候我们发现我们可以利用以前学的C语言、Java语言的一些基础语言来编写回调函数。也使我知道了GUI是实现人机交互的终结，具有强大的功能可以实现很多复杂的程序模块，使用它，需要具有一定的知识储备和必要的经验技巧，并且要充分利用好MATLAB的帮助文档，仔细研读HELP是最好的办法需要了解函数句柄等必要基础知识及其特有的功能，并会采用不同的使用手段实现相同功能的设计。我主要设计的是IIR的巴特沃斯和切比雪夫的高通滤波器。通过这次实训不仅通过硬件设计提高了我的动手能力，而且通过对软件的设计使我们巩固了数字信号处理中说滤波器的设计，也了解了利用MATLAB穿件图形用户界面GUI的知识，不管以后利用这个功能的机会多不多我们都提高学习新知识的能力，这也是我们在这次实训中的收获。只有学会了学习，获得了学习的能力，才能为我们今后的再学习打下基础，总之，这次实训使我受益匪浅。参考文献[1]程佩青.数字信号处理教程。北京:清华大学出版社,2000[2]刘树棠.数字信号处理—使用MATLAB[M].西安：西安交通大学出版社，2002[3]薛年喜.MATLAB在数字信号处理中的应用[M].北京清华大学出版社.007[4]【美】恒里Y-F拉姆著.模拟和数字滤波器设计与实现.人民邮电，1985[5]杨栓科.模拟电子技术基础.高等教育出版社.2010[6]GregWelch,GaryBishop.AnIntroductiontotheKalmanFilter[J].UNC-ChapelHill.May2003[7]沈晓晶,潘俊民.基于自适应Kalman预测器的算法[M].计算机仿真.Vol.21.2004,10[8]李天庆,张毅,刘志.基于时空图的轮廓线模型研究[M].计算机工程与应用.2004,4[9]王洪剑,孙志宏,彭思龙.基于GVF的骨架snake模型[M].计算机应用.Vol.24.2006,9[10]侯志强,韩崇昭.基于力场分析的主动轮廓模型[J].计算机学报.Vol.27.2007,[11]李丽勤,高焕文,周兴祥.Snake模型初始轮廓选取的研究[M].计算机工程与应用.2004附录1.滤波器软件设计程序functionAutoChoose(handles)Nn=128;DigitalFilter_value=get(handles.DigitalFilter,'Value');FilterDesign_value=get(handles.FilterDesign,'Value');%Windows_value=get(handles.Windows,'Value');FilterType_value=get(handles.FilterType,'Value');DisplayType_value=get(handles.DisplayType,'Value');Order_value=get(handles.Order,'Value');Rp_value=str2double(get(handles.Rp,'String'));Rs_value=str2double(get(handles.Rs,'String'));Fs_value=str2double(get(handles.Fs,'String'));Fp1_value=str2double(get(handles.Fp1,'String'));%Fp2_value=str2double(get(handles.Fp2,'String'));Fs1_value=str2double(get(handles.Fs1,'String'));%Fs2_value=str2double(get(handles.Fs2,'String'));wp1=2*Fp1_value/Fs_value;%wp2=2*Fp2_value/Fs_value;ws1=2*Fs1_value/Fs_value;%ws2=2*Fs2_value/Fs_value;%wp=[wp1,wp2];%ws=[ws1,ws2];if(DigitalFilter_value==1)if(FilterDesign_value==1)if((FilterType_value==1)||(FilterType_value==2))[n,Wn]=buttord(wp1,ws1,Rp_value,Rs_value)set(handles.MinOrderDisplay,'string',num2str(n))elseif((FilterType_value==3)||(FilterType_value==4))[n,Wn]=buttord(wp,ws,Rp_value,Rs_value)set(handles.MinOrderDisplay,'string',num2str(n))endendelseif(FilterDesign_value==2)if((FilterType_value==1)||(FilterType_value==2))[n,Wn]=cheb1ord(wp1,ws1,Rp_value,Rs_value)set(handles.MinOrderDisplay,'string',num2str(n))elseif((FilterType_value==3)||(FilterType_value==4))[n,Wn]=cheb1ord(wp,ws,Rp_value,Rs_value)set(handles.MinOrderDisplay,'string',num2str(n))endendelseif(FilterDesign_value==3)if((FilterType_value==1)||(FilterType_value==2))[n,Wn]=cheb2ord(wp1,ws1,Rp_value,Rs_value)set(handles.MinOrderDisplay,'string',num2str(n))elseif((FilterType_value==3)||(FilterType_value==4))[n,Wn]=cheb2ord(wp,ws,Rp_value,Rs_value)set(handles.MinOrderDisplay,'string',num2str(n))endendelseif(FilterDesign_value==4)if((FilterType_value==1)||(FilterType_value==2))[n,Wn]=ellipord(wp1,ws1,Rp_value,Rs_value)set(handles.MinOrderDisplay,'string',num2str(n))elseif((FilterType_value==3)||(FilterType_value==4))[n,Wn]=ellipord(wp,ws,Rp_value,Rs_value)set(handles.MinOrderDisplay,'string',num2str(n))endendendendendendelseif(DigitalFilter_value==2)if((FilterType_value==1)||(FilterType_value==2))[n,Wn]=buttord(wp1,ws1,Rp_value,Rs_value);set(handles.MinOrderDisplay,'string',num2str(n))elseif((FilterType_value==3)||(FilterType_value==4))[n,Wn]=buttord(wp,ws,Rp_value,Rs_value);set(handles.MinOrderDisplay,'string',num2str(n))endendendendMinOrder_value=get(handles.MinOrder,'Value');if(MinOrder_value==0)n=str2double(get(handles.Order,'String'))end数字滤波器IIR选择switchDigitalFilter_value%选择IIR滤波器case1%IIR中的Butterworth、Chebyshev1选择switchFilterDesign_value%选择设计的Butterworth滤波器case1%选择滤波器类型switchFilterType_valu