fir滤波器设计外文翻译.docx

FIR滤波器设计技术 摘要这份报告列举了一些设计FIR滤波器所使用的技术。首先讨论了窗函数法和 频率取样法的优点和缺点。FIR数字滤波器也包含了许多优化设计的方法,这些优化技术减少了在频率采样时非采样频率点的误差频率.对于用于设计数字滤波器的技术,例如matlab，进行了简明扼要的探讨。介绍FIR滤波器的系统函数是一个的多项式，因FIR滤波器的频率响应是频率的实函数，也称其为零相位滤波器。N阶FIR滤波器的系统函数表示为H(z)= （1）FIR滤波器是十分重要的,可应用于精确线性相位相应。FIR滤波器的实现方式保证了它是一个稳定的滤波器。FIR滤波器的设计可分为两部分：近似问题实现问题解决近似问题,要通过四个步骤找出传递函数：在频域内找出期望的或最理想的反应选择滤波器的阶数（FIR滤波器的长度N)选择近似结果中较好的选择一种算法寻找最优的滤波器传递函数选择部分结构处理实现传递函数的形式可能是线路图或程序。本质上来说，有三种著名的FIR滤波器设计方法：窗函数法频率取样法滤波器的优化设计窗函数法在该方法中，Park87，Rab75，ProakisOO从理想的频率响应Hd (w)出发,一般来说,单位脉冲相应hd (n)的持续时间是无限的，所以在某种程度上说,它必须截断。n=M_l约束着FIR滤波器的长度M。以M-1截断的hd (n)乘以窗函数就得到了滤波器的单位脉冲响应。矩形窗口的定义为 （2）FIR滤波器的单位脉冲相应为 h(n)=hd(n) w(n) = （3）现在,多元化的窗函数w(n)与hd(n)相当于hd(w)与w（w）的卷积，其中，w（w）是窗函数的频域表示。因此Hd(w)与w(w)的卷积为FIR数字滤波器的截断后的频率响应 （4）频率响应也可以利用以下的关系式 (5) 由于非均匀收敛的傅里叶级数的不连续性，其自身的波纹前后有一种近似于不连续的频率响应，因此直接截断的hd(n)来获得h(n)将导致吉布斯现象。与此同时，利用（5）得到的频率响应在频域内有波纹的振荡。为了减少波纹，hd(n)不是乘以一个矩形窗口w(n)，而是乘以一个含有圆锥和逐渐衰减到零的窗口。作为主体的序列的hd(n)和w(n)在时域内的卷积相当于其在频域内的乘积，其效果是平滑的。 滤波器的窗函数的傅里叶系数对滤波结果的频率响应的影响如下：（i） 一个主要的结果就是过渡带的不连续的两边出现中断（ii） 过渡带的宽度取决于窗函数的频率响应的主瓣宽度（iii） 滤波器的频率响应是通过卷积关系得到的，可以肯定的是，由产生 的滤波器绝不是最佳的（iv）随着M的增加，其主瓣宽度降低从而降低了过渡带的宽度，但是这也过滤掉了更多的脉冲频率响应。（v） 窗函数消除边缘响应引起的效果，并以较低的旁瓣代价增加过渡带的宽度1.Bartlett三角窗：2.广义余弦窗3.Kaiser窗 表一 从Park87得到的系数 Window a b c Retangular 1 0 0 Hanning 0.5 0.5 0 Hamming 0.54 0.46 0 Blackman 0.42 0.5 0.08 Bartlett窗函数的设计减少了信息的误差，但其过渡带较宽。Hanning,Hamming和Blackman窗的使用会更好，它们可以用于复杂的余弦函数，并可以提供理想的光滑截断的脉冲响应和频率响应。研究的结果表明，最佳的窗函数可能是有一个参数的Kaiser窗，它可以实现衰减和过渡带宽度的妥协。 窗函数的主要优点是它们比起其它方法更加的简单，且易于使用。事实上，计算窗函数的明确的方程系数就可以成功的使用该方法。在使用窗函数来设计滤波器时，会遇到三个问题：（i） 该方法只适用于Hd(w)是绝对可积的情况，即只有（2）式可以评估。当Hd(w)是复杂的或不能轻易被评价的闭合形式，写出Hd(n)的数学表达式就变得困难了。 （ii） 使用窗函数的灵活性比较差，例如，在低通滤波器的设计中，通频带的边缘频率一般不能用窗口完全掠过不连续区域。因此理想低通滤波器的截止频率，是通带截止频率f1和阻带截止频率f2相关的一个频率响应。 （iii） 窗函数法在设计标准滤波器，例如低通、高通、带通，是很有用的。 但这也使其在语音、图像处理的程序上的应用是十分有限的。 频率取样技术 在该方法中，Rad75、Park87、Proakis00是按照前面的方法提供理想的频率响应。现在，是在给定的频率响应中取一系列等间隔的频率，用以得到N的取样。因此，采样频率的响应Hd(w)在其本质上是给了我们Hd(2pnk/N)。因此，利用该滤波器可以计算出下面的公式： (6) 现在使用上述的N阶滤波器响应、连续性频率响应作为计算差值采样频率响应的方法。其近似误差就会完全接近于零点采样频率的误差，并将它们之间的频率有限化。越平滑的频率响应将会越近似，存在于样本点间的差值误差会很小。 一种减小误差的方法是增加频率采样的样本数目Rad75。其它的改善其质量的方法是找出一组无约束变量的制定频率样品。这些无约束变量的值一般都是由计算机来优化的一些简单的函数逼近，其误差接近最小。例如，一个无约束变量可能会选择在低通滤波器的通带于阻带之间的过渡带上的频率响应。 有两套不同的频率，可以用于取样。一组样品是fk=k/N其中k= 0,1,.N-1其他的均匀间隔的频率样品可以采取fk=(k+1/2)/N其中k=0、1、.N-1。 选择第二种时会在指定目标的第二个可能的频率的响应中给我们额外的灵活性。因此一个给定的过渡带边缘频率可能接近表二的频率采样点。在此情况下，依托于表二的设计将采用优化设计程序。 在一张由RabinerandGoldRabi70提供的纸上，Rabiner提到了一种技术，即基于采样频率的理念来设计FIR滤波器。Rabiner对于这种方法的步骤的建议如下：(i) 随着样品数量的变化，其提供了相应量级的响应。给定N，设计师决定使用什么样的插值。(ii) 当它被 Rabiner发现时，他实验了N从15到256的设计，16N样品的H（w）导致可靠的运算，因此，16到1的插值方法可用。 (iii) 给定N值，滤波器的单位样品的响应会被确定，h(N)的计算公式是逆傅里叶变换。(iv) Rabiner建议用两种程序去获得频率响应值。它们是（a）h(n)是由N/2样品或者（N-1 ）/2样品去移除其锐利的边缘，然后再15N 的冲击响应样品周围的位置都是对称的脉冲响应。 （b）h(n)分布在n/2样品之中，15N的样品安置在两脉冲响应之间。（c）零增广序列使用FFT算法转换为插值频率的响应。 频率取样技术的优点 （i） 与窗函数法不同，该技术可用于任何的响应。（ii）这种方法在设计非标准滤波器时是非常有用的，它可以处理任何不规则形状的响应。 频率取样法也有一些缺点，即通过插值得到的频率响应只是理想频率的采样点的响应。在其他点，将会出现一些错误。FIRFilterDesignTechniques Abstract ThisreportdealswithsomeofthetechniquesusedtodesignFIRfilters.Inthebeginning,thewindowingmethodandthefrequencysamplingmethodsarediscussedindetailwiththeirmeritsanddemerits.DifferentoptimizationtechniquesinvolvedinFIRfilterdesignarealsocovered,includingRabinersmethodforFIRfilterdesign.Theseoptimizationtechniquesreducetheerrorcausedbyfrequencysamplingtechniqueatthenon-sampledfrequencypoints.AbriefdiscussionofsometechniquesusedbyfilterdesignpackageslikeMatlabarealsoincluded.Introduction FIRfiltersarefiltershavingatransferfunctionofapolynomialinzandisanall-zerofilterinthesensethatthezeroesinthez-planedeterminethefrequencyresponsemagnitudecharacteristic.TheztransformofaN-pointFIRfilterisgivenby H(z)= (1)FIRfiltersareparticularlyusefulforapplicationswhereexactlinearphaseresponseisrequired.TheFIRfilterisgenerallyimplementedinanon-recursivewaywhichguaranteesastablefilter. FIRfilterdesignessentiallyconsistsoftwoparts(i) approximationproblem(ii) realizationproblem Theapproximationstagetakesthespecificationandgivesatransferfunctionthroughfoursteps.Theyareasfollows: (i) Adesiredoridealresponseischosen,usuallyinthefrequencydomain.(ii) Anallowedclassoffiltersischosen(e.g.thelengthNforaFIRfilters).(iii) Ameasureofthequalityofapproximationischosen. (iv) Amethodoralgorithmisselectedtofindthebestfiltertransferfunction.Therealizationpartdealswithchoosingthestructuretoimplementthetransferfunctionwhichmaybeintheformofcircuitdiagramorintheformofaprogram. Thereareessentiallythreewell-knownmethodsforFIRfilterdesignnamely:(1)Thewindowmethod (2)Thefrequencysamplingtechnique(3)Optimalfilterdesignmethods TheWindowMethod Inthismethod,Park87,Rab75,Proakis00fromthedesiredfrequencyresponsespecificationHd(w),correspondingunitsampleresponsehd(n)isdeterminedusingthefollowingrelation Ingeneral,unitsampleresponsehd(n)obtainedfromtheaboverelationisinfiniteinduration,soitmustbetruncatedatsomepointsayn=M-1toyieldanFIRfilteroflengthM(i.e.0toM-1).Thistruncationofhd(n)tolengthM-1issameasmultiplyinghd(n)bytherectangularwindowdefinedas （2）ThustheunitsampleresponseoftheFIRfilterbecomes h(n)=hd(n) w(n)(3) Now,themultiplicationofthewindowfunctionw(n)withhd(n)isequivalenttoconvolutionofHd(w)withW(w),whereW(w)isthefrequencydomainrepresentationofthewindowfunction ThustheconvolutionofHd(w)withW(w)yieldsthefrequencyresponseofthetruncatedFIRfilter (4) Thefrequencyresponsecanalsobeobtainedusingthefollowingrelation (5) Butdirecttruncationofhd(n)toMtermstoobtainh(n)leadstotheGibbsphenomenoneffectwhichmanifestsitselfasafixedpercentageovershootandripple beforeandafteranapproximateddiscontinuityinthefrequencyresponseduetothe non-uniformconvergenceofthefourierseriesatadiscontinuity.Thusthefrequencyresponseobtainedbyusing(8)containsripplesinthefrequencydomain.Inordertoreducetheripples,insteadofmultiplyinghd(n)witharectangularwindoww(n),hd(n)ismultipliedwithawindowfunctionthatcontainsataperanddecaystowardzerogradually,insteadofabruptlyasitoccursinarectangularwindow.Asmultiplicationofsequenceshd(n)andw(n)intimedomainisequivalenttoconvolutionofHd(w)andW(w)inthefrequencydomain,ithastheeffectofsmoothingHd(w). TheseveraleffectsofwindowingtheFouriercoefficientsofthefilterontheresultofthefrequencyresponseofthefilterareasfollows: (i) AmajoreffectisthatdiscontinuitiesinH(w)becometransitionbandsbetweenvaluesoneithersideofthediscontinuity. (ii) Thewidthofthetransitionbandsdependsonthewidthofthemainlobeofthefrequencyresponseofthewindowfunction,w(n)i.e.W(w). (iii) Sincethefilterfrequencyresponseisobtainedviaaconvolutionrelation,itisclearthattheresultingfiltersareneveroptimalinanysense. (iv) AsM(thelengthofthewindowfunction)increases,themainlobewidthofW(w)isreducedwhichreducesthewidthofthetransitionband,butthisalsointroducesmorerippleinthefrequencyresponse. (v) Thewindowfunctioneliminatestheringingeffectsatthebandedgeanddoesresultinlowersidelobesattheexpenseofanincreaseinthewidthofthetransitionbandofthefilter. (vi)SomeofthewindowsPark87commonlyusedareasfollows:1. Bartletttriangularwindow:2. Generalizedcosinewindows3. Kaiserwindowwithparameter: Thegeneralcosinewindowhasfourspecialformsthatarecommonlyused.Thesearedeterminedbytheparametersa,b,c TABLEI Valueofcoefficientsfora,bandcfromPark87 Window A b c Retangular 1 0 0Hanning 0.5 0.5 0Hamming 0.54 0.46 0 Blackman 0.42 0.5 0.08 TheBartlettwindowreducestheovershootinthedesignedfilterbutspreadsthetransitionregionconsiderably.TheHanning,HammingandBlackmanwindowsuseprogressivelymorecomplicatedcosinefunctionstoprovideasmoothtruncationoftheidealimpulseresponseandafrequencyresponsethatlooksbetter.ThebestwindowresultsprobablycomefromusingtheKaiserwindow,whichhasaparameter.thatallowsadjustmentofthecompromisebetweentheovershootreductionandtransitionregionwidthspreading. Themajoradvantagesofusingwindowmethodistheirrelativesimplicityascomparedtoothermethodsandeaseofuse.Thefactthatwelldefinedequationsareoftenavailableforcalculatingthewindowcoefficientshasmadethismethodsuccessful. Therearefollowingproblemsinfilterdesignusingwindowmethod:(i) ThismethodisapplicableonlyifHd(w)isabsolutelyintegrablei.eonlyif(2)canbeevaluated.WhenHd (w)iscomplicatedorcannoteasilybeputintoaclosedformmathematicalexpression,evaluationofhd(n)becomesdifficult. (ii) Theuseofwindowsoffersverylittledesignflexibilitye.g.inlowpassfilterdesign,thepassbandedgefrequencygenerallycannotbespecifiedexactlysincethewindowsmearsthediscontinuityinfrequency.ThustheidealLPFwithcut-offfrequencyfc,issmearedbythewindowtogiveafrequencyresponsewithpassbandresponsewithpassbandcutofffrequencyf1andstopbandcut-offfrequencyf2.(iii) Windowmethodisbasicallyusefulfordesignofprototypefilterslikelowpass,highpass,bandpassetc.Thismakesitsuseinspeechandimageprocessingapplicationsverylimited. TheFrequencySamplingTechniqueInthismethod,Park87,Rab75,Proakis00thedesiredfrequencyresponseisprovidedasinthepreviousmethod.NowthegivenfrequencyresponseissampledatasetofequallyspacedfrequenciestoobtainNsamples.Thus,samplingthecontinuousfrequencyresponseHd(w)atNpointsessentiallygivesustheN-pointDFTofHd(2pnk/N).ThusbyusingtheIDFTformula,thefilterco-efficientscanbecalculatedusingthefollowingformula (6) NowusingtheaboveN-pointfilterresponse,thecontinuousfrequencyresponseiscalculatedasaninterpolationofthesampledfrequencyresponse.Theapproximationerrorwouldthenbeexactlyzeroatthesamplingfrequenciesandwouldbefiniteinfrequenciesbetweenthem.Thesmootherthefrequencyresponsebeingapproximated,thesmallerwillbetheerrorofinterpolationbetweenthesamplepoints. NowusingtheaboveN-pointfilterresponse,thecontinuousfrequencyresponseiscalculatedasaninterpolationofthesampledfrequencyresponse.Theapproximationerrorwouldthenbeexactlyzeroatthesamplingfrequenciesandwouldbefiniteinfrequenciesbetweenthem.Thesmootherthefrequencyresponsebeingapproximated,thesmallerwillbetheerrorofinterpolationbetweenthesamplepoints. OnewaytoreducetheerroristoincreasethenumberoffrequencysamplesRab75.Theotherwaytoimprovethequalityofapproximationistomakeanumberoffrequencysamplesspecifiedasunconstrainedvariables.Thebandbetweenthepassbandandthestopbandofalowpassfilter. Therearetwodifferentsetoffrequenciesthatcanbeusedfortakingthesamples.Onesetoffrequencysamplesareatfk=k/Nwherek=0,1,.N-1.Theothersetofuniformlyspacedfrequencysamplescanbetakenatfk=(k+1/2)/Nfork=0,1,.N-1.Thesecondsetgivesustheadditionalflexibilitytospecifythedesiredfrequencyresponseatasecondpossiblesetoffrequencies.Thusagivenbandedgefreque