【精品】毕业设计（论文）—基于matlab的fir数字滤波器设计及其软件实现

基于MATLAB的FIR滤波器设计摘要在这个具有挑战性的时代，把计算机充分运用到教学及工程运算中，虽然具有重要的意义，而随着计算机技术的发展，计算机软件在工程设计领域应用越来越广。本设计研究是基于MATLAB的FIR数字滤波器的设计，是基于应用计算机软件编程。首先了解滤波器的基本工作原理，针对FIR滤波器，有各种不同的方法。FIR滤波器一般要求信号具有线性相位，同时由于脉冲响应无限长，因此具有永远稳定的特性。根据这个特点，。本文主要对FIR滤波器采用了窗函数法，最优设计法，最小二乘设计法，升余弦函数设计法和任意响应设计法来设计FIR滤波器，同时，还对滤波器进行阶数评估。滤波是信号处理中最基本有极为重要的技术，利用滤波器技术可以从复杂的信号中提取出所需要的信号，抑制不需要的信号。绝大多数传感器输出的信号，在使用过程中，都必须进行滤波，所以滤波器是具有一定传输选择特性的，对信号进行加工处理的装置，它允许输入信号中的一些成分通过，抑制或衰减另一些成分。其功能是将输入信号变换为人们所需要的输入信号。本设计主要给出了FIR数字滤波器的基本特性和设计方法。关键词MATLAB特性FIR滤波器设计ABSTRACTINTHISCHALLENGINGERA,THEFULLUSEOFCOMPUTERSTOTEACHINGANDENGINEERINGOPERATIONS,ALTHOUGHOFGREATSIGNIFICANCE,ANDWITHTHEDEVELOPMENTOFCOMPUTERTECHNOLOGY,COMPUTERSOFTWAREINENGINEERINGDESIGNAPPLICATIONSMOREWIDELYTHEDESIGNSTUDYISBASEDONMATLABFORFIRDIGITALFILTERDESIGN,ISBASEDONTHEAPPLICATIONOFCOMPUTERSOFTWAREPROGRAMMINGFIRSTUNDERSTANDTHEBASICWORKINGPRINCIPLEOFTHEFILTERFORFIRFILTERS,AVARIETYOFDIFFERENTWAYSGENERALREQUIREMENTSFORFIRFILTERSWITHLINEARPHASESIGNALS,ANDBECAUSEOFINFINITEIMPULSERESPONSE,SOITHASALWAYSSTABLECHARACTERISTICSACCORDINGTOTHISFEATUREINTHISPAPER,FIRFILTERSUSEDONTHEWINDOWFUNCTIONMETHOD,OPTIMALDESIGNMETHOD,LEASTSQUARESDESIGNMETHOD,DESIGNRAISEDCOSINEFUNCTIONANDDESIGNMETHODTODESIGNARBITRARYRESPONSEFIRFILTERS,ITISALSOTHEORDEROFTHEFILTEREVALUATIONSIGNALPROCESSINGFILTERISTHEMOSTBASICTHEREISAVERYIMPORTANTTECHNOLOGY,THEUSEOFTECHNOLOGYCANFILTERTHESIGNALFROMTHECOMPLEXBYTHENEEDTOEXTRACTTHESIGNAL,TOSUPPRESSUNWANTEDSIGNALSMOSTOFTHESENSOROUTPUTSIGNAL,INTHECOURSE,MUSTBEFILTERED,SOTHEFILTERISCHOSENWITHACERTAINTRANSMISSIONCHARACTERISTICS,THESIGNALPROCESSINGDEVICE,WHICHALLOWSSOMECOMPONENTSOFTHEINPUTSIGNALTHROUGHTHEINHIBITIONORATTENUATIONOTHERINGREDIENTSITSFUNCTIONISTOTRANSFORMTHEINPUTSIGNALISNEEDEDFORPEOPLETOINPUTSIGNAL目录摘要IABSTRACTII第一章绪论111数字滤波器定义112数字滤波器的分类113实际滤波器的设计指标214几种常见的特殊滤波器2第二章FIR滤波器特性321FIR滤波器简介322FIR滤波器的线性相位特性3221FIR滤波器的第一类线性相位4222FIR滤波器的第二类线性相位6223线性相位FIR滤波器的零点特性8第三章FIR滤波器的设计931用窗函数法设计FIR数字滤波器9311FIR数字滤波器窗函数设计法9312加窗的线性相位FIR数字滤波器设计函数FIR1和FIR21232FIR滤波器的最优设计法15321FIRS函数15322REMEZ函数1633FIR滤波器最小二乘设计法17331FIRCLS函数17332FIRCLS函数1934FIR滤波器升余弦函数设计法2035FIR滤波器的任意响应设计法21第四章FIR滤波器阶数估计2441KAISERORD函数2442REMEZORD函数25第五章小结27参考文献28致谢29第一章绪论当前，通信技术高速发展，业务范围不断扩大，人们对产品的需求迅速增长，滤波器在这些产品电路中就扮演着重要的角色，当然数字滤波器更影响着人们生活的方方面面。数字滤波器既可以是用数字硬件装配成德一台完成给定运算的专用数字计算机，也可将所需的运算编成程序，让通用计算机来执行。MATLAB是“矩阵实验室”（MATRIXLABORATOY）的缩写，它是由美国MATHWORKS公司于1984年正式推出的，是一种以矩阵运算为基础的交互式程序语言，专门针对科学、工程计算及绘图的需求。MATLAB是功能强大的科学及工程计算软件，它不但具有以矩阵计算为基础的强大数学计算和分析功能，而且还具有丰富的可视化图形表现功能和方便的程序计算能力。MATLAB的应用领域极为广泛，除数学计算和分析外，还被广泛地应用于自动控制、系统仿真、数字信号处理、图形图像分析、数理统计、人工智能、虚拟现实技术、通信工程、金融系统等领域，因此，MATLAB是面向21世纪的计算机程序设计及科学计算语言。11数字滤波器定义数字滤波器通常是指一个有限精度算法实现的离散线性是不变系统。通常用的数字滤波器一般属于选频滤波器，可以用下式表示|称为滤波器幅频响应，称为滤波器的相频响应。幅频响式中，|（）应表示信号通过该滤波器后各频率成分的衰减情况，而相频响应反映各频率成分通过滤波器后在时间上的延时情况。滤波器性能最容易通过它的幅频响应的形状来描述。滤波器在某个频率的幅度增益决定了滤波器对此频率输入的放大因子，增益可取任意值。增益高的频率范围，信号可以通过，称之为滤波器的通带，增益低的频率范围，滤波器对信号有衰减和阻塞作用，称这位滤波器的阻带。12数字滤波器的分类数字滤波器按照不同的分类方法，有许多种类，总体上可分为两大类。一类是经典滤波器，即一般的线性系统滤波器。另一类即所谓的现代滤波器，它的理论建立在随机信号处理的理论基础上，利用了随机信号内部的统计特性对信号进行滤波。从功能上可以分为低通、高通、带通和带阻滤波器，根据数字滤波器冲激响应的时域特性，可将数字滤波器分为两种，即无限长冲激响应（IIR）滤波器和有限长冲激响应（FIR）滤波器3。13实际滤波器的设计指标当滤波器形状为非理想时，要用一些参数指标来描述其关键特性。滤波器的通带定义了滤波器允许通过的频率范围。在阻带内，滤波器对信号严重衰减。和分别称为通带截止频率（或通带上限频率）和阻带截止频率（或阻带下限频率）。参数定义了通带波纹，及滤波器通带内偏移单位增益的最大值。参1数定义了阻带波纹，及滤波器阻带内偏离零增益的最大值。参数定义了过2渡带宽度，即阻带下限和通带上限之间的距离，。过渡带一般是单|调下降的，通带内和阻带内允许的衰减一般用单位DB表示，通带内允许的最大衰减用表示，阻带内允许的最小衰减用表示，它们分别定义为DB20201111DB2020112式中，是通带内的幅度最大值；是通带内的幅度最小值，是阻带内最大值。幅度下降到0707即时，W，此时3DB，称为3DB通带截止22频率。14几种常见的特殊滤波器全通滤波器、数字陷波器、梳状滤波器、最小相位滤波器第二章FIR滤波器特性21FIR滤波器简介FIR系统只有零点，因此这类系统不像IIR系统那样易取得比较好的通带和阻带衰减特性。但FIR系统有自己突出的优点，其一是系统总是稳定的，其二是易实现线性相位，其三是允许设计多通带（或多阻带）滤波器。FIR滤波器是指系统的单位冲击响应H（N）仅在有限的范围内有非零值的滤波器。N1阶滤波器的系统函数H（Z）可表示为（）10（）H（N）是的N1次多项式，它在Z平面上有N1个零点，原点Z0是N11阶重极点。因此，FIR滤波器永远稳定。FIR滤波器的频率响应为（）10（）由于H（）一般为复数，因此，可将其表示成（）|（）|（）22FIR滤波器的线性相位特性FIR滤波器的频率响应为（）10（）由于H（）一般为复数，因此，可将其表示成（）|（）|（）线性相位特性是指是W的线性函数，即（）（）（）式中，是常数。此时通过这一系统的各频率分量的时延为相同的常数，系统的群时延为（）即系统的时延是一个与W无关的常数，称系统H（Z）具有严格的线性相位。由于严格线性相位条件在数学上处理较为困难，因此在FIR滤波器设计中一般使用广义线性相位。若一个离散系统的频率响应H（）可以写为H（）（）（）（）（）其中，和是与W无关的常数，是W的是函数（可任意取值），则称（）为广义幅频响应，称为广义相频响应。（）上式中，为（）（）当0时，称是第一类线性相位，当（一般）时，称是第二类线性02相位。221FIR滤波器的第一类线性相位FIR滤波器满足第一类线性相位的条件是（1）由于H（N）的点数N可分为奇数和偶数两种情况，因而第一类线性相位又可分为两类。（1）H（N）为偶对称，N为奇数（）（）10（）320（）（12）（12）11232011212123202COS1212显然（）12满足第一类线性相位的条件，而（）3202（）COS12（12）由于N为奇数，因此N为整数。由于关于W0，2这些点12COS（）都是偶对称的，因此关于W0，2也是偶对称的，这种情况适合设（）计低通、高通、带通、带阻滤波器。（2）H（N）为偶对称，N为偶函数采用类似的方法，可以得到（）210（）210（1）（1）210（）（1）（12）2102（）COS12其相频响应与N为奇数时相同，为（）12而幅频响应为（）2102（）COS12由于N为偶数，因此为05的奇数倍。关于W奇对称，12COS12当W时，0，因此，0，即H（Z）在Z1处必然COS12（）有一个零点，而且，关于W奇对称这种情况不适合设计高通和带阻（）滤波器。222FIR滤波器的第二类线性相位FIR滤波器满足第二类相位的条件（）（1）（1）H（N）为奇对称，N为奇数H（N）的中间项H（）必须为零，因此12（）320（）112（）320（1）1223202SIN12此时的相位为212满足第二类线性相位的条件，而幅频响应为（）3202SIN12因为在W0，处都为0，因此在W0，处为0，SIN（），2（），2也即H（Z）在Z1处都有零点，并且关于W0，奇对称。因（），2此，该类型滤波器不适合设计低通、高通和带阻滤波器，只适合设计带通滤波器。（2）H（N）为奇对称，N为偶数（）1222102SIN12相位为212幅频响应为（）2102SIN12由于在W0,2处为0，因此在W0,2处为0，即SIN12（）H（Z）在Z1上有零点，且关于W0,2奇对称。H（N）为奇对称有两种情况，对于任何频率都有固定的相移，因此第二类线性2相位FIR滤波器一般用于正交移相或微分器。223线性相位FIR滤波器的零点特性线性相位FIR滤波器的单位脉冲响应具有对称特性，即（）（1）则（）10（）10（1）将MN1N代入上式，得（）10（）（1）（1）10（）因此，其系统函数具有以下特点（）（1）（1）可以看出，若Z及Z也必定是零点。1第三章FIR滤波器的设计IIR滤波器设计利用了模拟滤波器成熟理论以及设计图表，保持了典型模拟滤波器优良的幅频特性，但设计中只考虑了幅频特性，没有考虑相位特性，因此滤波器的相位特性通常是非线性的，要得到线性相位必须外加相位矫正网络，使滤波器变得复杂。而FIR滤波器在保证幅频满足技术要求的同时，更荣易做到严格的线性相位。FIR滤波器的系统函数为（）（）（）在Z平面上有N1个零点，并且Z0是N1阶重极点，位于R1的单位圆内，系统永远稳定。FIR滤波器的设计方法有窗函数法、频率采样法和契比雪夫等波纹逼近法31用窗函数法设计FIR数字滤波器311FIR数字滤波器窗函数设计法通常，理想数字滤波器的单位取样序列是无限长的，是非因果序列，在（）物理上是不可实现的或者难以实现的。为了可物理实现，需要对截取一（）段H（N）来近似代表，即H（N）（）（）（）窗函数主瓣宽度要尽可能窄，以获取最小的过渡带；旁瓣相对值尽可能小，以使得通带波纹小，并且阻带衰减大。（1）常用的窗函数为矩形窗、巴特利窗汉宁窗、海明窗、布莱克曼窗、凯泽、窗和契比雪夫窗。1）矩形窗（）（）矩形窗的频域响应为（）SIN2SIN2121其主瓣宽度为。调用WBOXCAR（N）函数42）巴特利窗（）21，0（1）2221，12（1）1其频率响应为（）2SIN4SIN22（12）其主瓣宽度为，调用WBARTLETT（N）函数83）汉宁窗（）051COS21（）其频域表达式为（）05（）025（2）（2）主瓣宽度为，调用WHANN（N）函数84）海明窗（）054046COS（21）（）其频域表达式为（）054（）023（21）023（21）主瓣宽度为，调用WHAMMING（N）函数85）布莱克曼窗（）04205COS21008COS41（）其频域表示式为（）042（）02521210044141其主瓣宽度为，调用WBLACKMAN（N）函数126）凯泽窗（）（）（）01式中，为第一类修正贝塞尔函数，窗函数的幅12112（）度函数为（）（0）2121（）COS其主瓣宽度为，调用WKAISER（N，BETA）函数107）契比雪夫窗调用WCHEBWIN（N，R）表31窗函数基本参数窗函数旁瓣峰值幅度/DB过渡带阻带最小衰减/DB矩形窗13421巴特利窗26825汉宁窗31844海明窗40853布莱克曼窗571274凯泽窗571080（2）FIR滤波器的设计步骤1）根据技术要求确定待求滤波器的单位取样响应或者频率响应。（）2）根据对过渡带和阻带衰减的要求，选择窗函数形式，并估计窗口长度N。3）计算滤波器的单位取样响应H（N），H（N）即为所求。（）（）4）检测所设计的滤波器性能，验算技术指标是否满足条件。312加窗的线性相位FIR数字滤波器设计函数FIR1和FIR21）FIR1。其调用格式为BFIR1（N，WN）BFIR1（N，WN，FTYPE）BFIR1（N，WN，WINDOW）BFIR1（N，WN，FTYPE，WINDOW）BFIR1（,NORMALIZATION）该函数实现加窗的线性相位FIR数字滤波器设计，可以设计标准低通、带通、高通和带阻滤波器。BFIR1（N，WN）得到N阶低通FIR滤波器，调用后返回维数为N1的行向量B，它是滤波器的系数。B与FIR滤波器的系统函数的关系为1（）（1）（2）1对于高通、带阻滤波器，N取偶数，为滤波器的截止频率，范围为（0,1）对于带通和带阻滤波器，且；对于多通带滤波器，W1，212，频段为01，2，3，41，12，23。用窗函数法设计线性相位FIR低通滤波器，通带截止频率05，阻带衰减不小于40DB，通带衰减不大于3DB。06下面采用汉宁窗设计滤波器DESIGNFIRFILTERWITHWINDOWSWP05PIWS06PIWDELWSWPNCEIL8PIWDELWN0506PI/2WINDOWHANNINGN1BFIRLN,WN/PI,WINDOWFREQZB,1,512图31采用汉宁窗设计低通滤波器用窗函数法设计多通带滤波器，归一化通带为。002，04，06，08,1由于高频段为通带，滤波器的阶数应为偶数，这里定位36，实现程序代码如下，滤波器幅频相频特性如图DESIGNMULTIBANDFIRFILTERWITHWINDOWSW02040608BFIR136,W,DC1FREQZB,1,512图32多通带FIR数字滤波器2）FIR2。其调用格式为BFIR2（N，F，M）BFIR2（N，F，M，WINDOW）BFIR2（N，F，M，NPT）BFIR2（N，F，M，NPF，WINDOW）BFIR2（N，F，M，NPF，LAP）BFIR2（N，F，M，NPF，LAP，WINDOW）该函数用于设计基于窗函数的任意相应的FIR滤波器，其频率响应由向量F和M决定，取值在之间；N为滤波器阶数；B向量为返回滤波器的系数；0,1WINDOW为窗类型，长度必须为N1，默认时为HANNING窗；NPT为对频率响应进行内插点数，默认时为512；LAP参数用于指定FIR2在重复频率点附近插入的区域大小。设计多通带滤波器，滤波器阶数为，比较理想滤波器和实际滤波器的频率响应。滤波器设计程序代码为DESIGNMULTIBANDFIRFILTERWITHFIR2M0011000110001100F001015020250304045050550607075080851N40BFIR2N,F,M,HAMMINGN1H,WFREQZB,1,128PLOTF,M,W/PI,ABSHXLABELFREQUENCYYLABELMAGNITUDEGRID图33用FIR2设计的多通带FIR数字滤波器32FIR滤波器的最优设计法FIRS是FIR1和FIR2函数的扩展，它采用最小二乘法，使指定频段内的理想分段线性函数与滤波器幅频响应之间的误差平方和最小。REMEZ函数则利用PARKSMCCLELLAN算法设计线性相位FIR滤波器。该算法能使期望频率响应和实际频率响应之间的最大误差最小。FIRS函数设计的滤波器在频率响应上出现等波纹滤波器。函数FIRS、REMEZ调用的语法规则相同。321FIRS函数其调用格式是BFIRS（N，F，A）BFIRS（N，F，A，W）BFIRS（N，F，A，FTYPE）BFIRS（N，F，A，W，FTYPE）该函数用于设计N阶FIR滤波器，其幅频特性由F和A向量确定，调用后返回长度为N1的滤波器系数向量B，且这些系数遵循以下偶对称关系（）（2），1,2，1F是频率点向量，范围为，频率点是逐渐增大，允许向量中有重复的频0,1率点；A是指定频率点的幅频响应，期望的频率响应由（F（K），A（K1）的连线组成，FIRLS（）则把F（K1）和F（K2）之间的频带视为过渡带。所以，所需要的频率响应是分段线性的，其总体平方误差最小。BFIRLS（N，F，A，W）则使用权系数W给误差加权。W的长度为F和A的一半。BFIRLS（N，F，A，W，FTYPE）调用形式中，参数FTYPE用于指定所设计的滤波器类型，FTYPEHILBERT，为奇对称的线性相位滤波器，返回的滤波器系数满足B（K）B（N2K），K1，2，N1；FTYPEDIFFERENTIATOR，则采用特殊加权技术，生成奇对称的线性相位滤波器，使低频段误差大大小于高频段误差。设计24阶反对陈滤波器，使其具有分段线性的通带，并画出理想滤波器和实际滤波器的频率响应。实现代码如下PERFORMANCEOFFILTERF034679M01000505BFIRLS24,F,M,HILBERTFORI126PLOTFIFI1,MI,MI1,HOLDONENDH,FFREQZB,1,512,2PLOTF,ABSHGRID图34理想和实际分段带通滤波器特性注虚线是理想滤波器322REMEZ函数该函数的调用格式与参数意义都与FIRLS函数一样，只是采用的算法不同，下面通过实力比较它们的区别分别使用REMEZ和FIRLS函数设计20阶带通滤波器，通常频率为WP104，WP206，实现代码如下DESIGNFIRFILTERWITHFIRLSANDREMEZF0030406071M001100N20BFIRLSN,F,MH,WFREQZBBBREMEZN,F,MHH,WFREQZBBAXESPOSITION,02020505PLOTW/PI,ABSH,B,W/PI,ABSHH,RXLABELFREQUENCYYLABELMAGNITUDEGRIDON图35两种最优FIR滤波器的比较33FIR滤波器最小二乘设计法约束最小二乘法就是在给定滤波器幅频响应最大允许波纹阀值约束条件下，使滤波器在整个幅频范围内误差平方最小化。FIRCLS和FIRCLS1两个函数，可以实现此滤波器的设计。331FIRCLS函数其调用格式为BFIRCLS（N，F，A，UP，LO）该函数返回长度N1的线性相位滤波器，期望逼近的频率分段恒定，由向量F和A确定，频率的上下限由参数UP和LO确定，长度与A相同。F中元素为临界频率，取值范围，且按递增顺序排列。0,1使用FIRCLS函数设计一带通滤波器，通带为，并用设计的滤波器对信02,04号滤波SIN21505SIN29002SIN2300（信号采样频率为600HZ）。实现代码如下DESIGNFIRFILTERWITHFIRCLST01/6001SIGSIN2PI15T05SIN2PI90T02SIN2PI300TPLOTT,SIGN50F002041A010UP002102001LO002098001BFIRCLSN,F,A,UP,LONEWSIGFFTFILTB,SIG；FIGUREFTT301350；NSNEWSIG301350ZNSINTERPNS,12ZNTINTERPFT,12PLOTZNT,ZNS时间/S时间/S图36滤波前后的波形332FIRCLS函数该函数采用约束最小二乘法设计基本的线性相位高通和低通滤波器，其调用格式为BFIRCLS（N，WO，DP，DS）BFIRCLS1（N,WO,DP,DS,HIGH）BFIRCLS1N,WO,DP,DS,WTBFIRCLS1N,WO,DP,DS,WT,HIGHBFIRCLS1N,WO,DP,DS,WP,WS,KBFIRCLS1N,WO,DP,DS,WP,WS,K,HIGHBFIRCLS1N,WO,DP,DS,DESIGNFLAGBFIRCLS（N，WO，DP，DS）函数返回长度为N1的线性相位低通FIR滤波器，截止频率为WO，取值01之间。通常幅度偏离10的最大值为DP，阻带偏离0的最大值为DS。BFIRCLS1（N,WO,DP,DS,HIGH）则返回高通滤波器，N必为偶数。BFIRCLS1N,WO,DP,DS,WP,WS,K采用平方误差加权，通常的权值比阻带的大K倍。WP为通带边缘频率，WS为阻带边缘频率，其中；如果要设计高通滤波器，则必须使。设计55阶低通滤波器，截止频率在03。设计代码如下DESIGNFIRFILTERWITHFIRCLS1N55WO03DP002DS0008BFIRCLS1N,WO,DP,DSH,WFREQZBPLOTW/PI,ABSHGRID图37FIRCLS1设计的低通滤波器34FIR滤波器升余弦函数设计法用FIRCOS函数，用于设计有光滑、升余弦过渡带的低通线性相位滤波器，其调用格式为BFIRCOS（N，FO，DF，FS）BFIRCOS（N，FO，DF，FS，BANDWIDTH）BFIRCOS（N，FO，DF）BFIRCOS（N，FO，R，FS，ROLLOFF）BFIRCOS（,TYPE）BFIRCOS（，TYPE，DELAY，WINDOW）FIRCOS（），其中，N为滤波器阶数；FO为低通滤波器的截止频率，DF为过渡带频宽，FS是采样频率，单位都是HZ；FO和FO都介于0之间；B为返回滤波器系22数向量。设计20阶升余弦滤波器，截止频率为250HZ，过渡带宽为100HZ，采样频率为1000HZ。其程序代码如下DESIGNFIRFILTERWITHFIRRCOSN20FO250DF100FS1000BFIRRCOSN,FO,DF,FSH,FFREQZB,1,512,FSAXESPOSITION,02020505PLOTF,ABSHXLABELFREQUENCYHZYLABELMAGNITUDEGRIDON；图38升余弦函数设计低通滤波器35FIR滤波器的任意响应设计法CREMEZ函数是设计FIR滤波器的另一函数，该函数可用来设计任意响应的复数和非线性相位等波纹FIR滤波器。它在指定的滤波器频率响应的方式上和其他滤波器设计函数有所不同，接收的是函数。CREMEZ函数采用扩展REMEZ交换算法，能使契比雪夫误差最小，其调用格式为BCREMEZ（N，F，FRESP）BCREMEZ（N，F，FRESP，W）BCREMEZ（N，F，FRESP，P1，P2，W）BCREMEZ（N，F，A，W）BCREMEZ（，SYM）BCREMEZ（，SKIPSTAGE2）BCREMEZ（，DEBUG）BCREMEZ（，LGRID）B，DELTA，OPTCRME（）其中，参数F是频率边界向量，范围是1到1。W是每个频带的加权系数，长度为F的一半。P1/P2是可选参数，FRESP则指定滤波器类型，可能取值为（1）LOWPASS、HIGHPASS、BANDPASS、BANDSTOP。比如以下引用形式BCRME（N，F，LOWPASS，）BCREMEZ（N，F，LOWPASS，D，）（2）MULTIBAND。该函数指定设计多通带滤波器，例如BCREMEZ（N，F，MULTIBAND，A，）BCREMEZ（N，F，MULTIBAND，A，D，）其中，A为在边界频带向量F各点上的期望幅值。（3）DIFFERENTIATOR。该函数指定设计线性相位微分器，其零频必须在过渡带，而频带加权设置与频率成反比，引用方式如下，其中FS用于确定微分器的响应斜率，默认值为1。BCREMEZ（N，F，DIFFERENTIATOR，FS，）BCREMEZ（N，F，DIFFERENTIATOR，FS，D，）（4）HILBFILT。该函数指定设计线性相位希尔伯特变换滤波器，其零频必须在过渡带，比如BCRME（N，F，HILBFILT，）BCREMEZ（N，F，HILBFILT，D，）BCREMEZ（，SYM）则使用所设计滤波器的冲击响应有对称性，SYM取值如下NONE，默认值，若F中有负值或FRESP未赋值则为默认值。EVEN，使单位冲击响应为实偶对称序列。ODD，使单位冲击响应为实奇对称序列。设计希尔伯特和差分滤波器时，取默认值。REAL，使频率响应为共轭对称。BCRME（，DEBUG）则将显示中间计算结果，DEBUG可为TRACE、PLOT、BOTH或者OFF之一，默认值为OFF。设计30阶任意响应多频带滤波器。滤波器的频率边界频率向量为F1，05，04，03，04，08；各频段的幅值响应为A3，1，2，2，2，1；各频段最优化权向量为W1，10，5设计代码如下COMPLEXANDNONLINEARPHASEFIRFILTERDESIGNN30F10504030408A312221W1105BCREMEZN,F,MULTIBAND,A,WH,WFFREQZB,1,512,WHOLEPLOTWF/PI1,FFTSHIFTABSHXLABELFREQUENCYYLABELMAGINITUDEGRIDON图39任意复响应滤波器频谱图第四章FIR滤波器阶数估计41KAISERORD函数其调用格式为N,WN,BETA,FTYPEKAISERORDF,A,DEVN,WN,BETA,FTYPEKAISERORDF,A,DEV,FSCKAISERORDF,A,DEV,FS,CELL该函数返回值在函数BFIR1（N，WN，KAISER（N1，BETA），FTYPE，NOSCALE）中需要使用的阶数N、归一化频带边缘WN、KAISER窗函数BETA、滤波器类型FTYPE的估计值。所设计的滤波器近似满足由输入参数F、A和DEV给定的性能。F为频率向量，A为F定义的各频带的理想频率响应幅值，F的长度是A的2倍减去2（是偶数），第一频带的起始频率为0，最后一个频带截止频率为FS。向量DEV的各元素分别是各频带允许的最大偏移量。而CKAISERORD（F，A，DEV，FS，CELL）返回的是向量组，用于FIR1函数的参数。设计低通滤波器，通带范围为01HZ，阻带范围为41500HZ，通带波纹为5，阻带衰减为40DB，系数采样频率FS8000HZ。设计程序如下DESIGNLOWPASSFILTERWITHKAISERORDFSAMP8000FCUTS10001500MAGS10DEVS005001N,WN,BETA,FTYPEKAISERORDFCUTS,MAGS,DEVS,FSAMPHHFIR1N,WN,FTYPE,KAISERN1,BETA,NOSCALEFREQZHH图41利用KAISERORD设计的低通滤波器幅频特性42REMEZORD函数REMEZORD函数为REMEZ函数选择滤波器阶数，在给定频域中的性能指标后，REMEZORD可产生近似的满足指标的最小阶数。函数调用格式如下N，FO，AO，WREMEZORD（F，A，DEV）N，FO，AO，WREMEZORD（F，A，DEV，FS）CREMEZORD（F，A，DEV，FS，CELL）该函数能找出近似阶数N、归一化频率边界FO、频带内幅值AO、以及权向量W，使由REMEZ函数构成的滤波器满足参数F、A、DEV指定的性能要求。F为频率向量，长度是参数A的2倍减去2（是偶数），而A为F各频带理想频率响应的幅值，向量DEV的元素分别是各频带允许的最大偏差。设计最小阶数的低通滤波器，通常截止频率为500HZ，阻带截止频率为600HZ采样频率为2000HZ，阻带最小衰减为40DB，通带最大衰减为3DB。程序代码如下DESIGNLOWPASSFILTERWITHKAISERORDRP3RS40FS2000F500600A10DEV10RP/201/10RP/20110RS/20N,FO,AO,WREMEZORDF,A,DEV,FSBREMEZN,FO,AO,WFREQZB,1,1024,FSTITLELOWPASSFILTERDESI