版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進行舉報或認領(lǐng)
文檔簡介
1Introduction
☆
FIRfilter:directdesignofDTfilterwiththeoftenaddedlinear-phaserequirement(1)WindowedFourierseriesapproach(§10.2)(2)Frequencysamplingapproach(Problem10.31,10.32)
(3)Computer-basedoptimizationmethod(§10.3)
Chap.10FIRDigitalFilterDesign1IntroductionChap.10FIRDi210.1PreliminaryConsiderationsForFIRsystem:realpolynomialapproximationifalinearphaseisdesired
10.1.1BasicApproachestoFIRDigitalFilterDesign210.1PreliminaryConsiderati310.1.2EstimationoftheFilterOrderKaiser’sFormulaForlowpassFIRfilterdesign:P397-398Bellanger’sFormulaHermann’sFormulaParametersseeP398.310.1.2EstimationoftheFil410.2DesignofFIRFiltersbyWindowing(P400)
10.2.1LeastIntegral-SquaredErrorDesignofFIRFilters410.2DesignofFIRFiltersb510.2.2ImpulseResponsesofIdealFiltersIdeallinearphaselowpassfilterIdeallinearphase
highpassfilter510.2.2ImpulseResponsesof6ImpulseResponsesofIdealFilters(II)Ideallinearphasebandpassfilter
Ideallinearphasebandstopfilter
6ImpulseResponsesofIdealFi7ImpulseResponsesofIdealFilters(III)IdealmultibandfilterIdealdiscrete-timeHilberttransformer
Idealdiscrete-timedifferentiator
7ImpulseResponsesofIdealFi8Gibbsphenomenon:
OscillatorybehaviorinthemagnituderesponseofcausalFIRfiltersdesignedutilizingtruncation10.2.3GibbsPhenomenon8Gibbsphenomenon:9mainlobesidelobeMainlobewidth--truncationperiodiccontinuousconvolution=N/210.2.3GibbsPhenomenon(II)9mainlobesidelobeMainlobewidt10-2/(N+1)10-2/(N+1)1110.2.3GibbsPhenomenon(III)1110.2.3GibbsPhenomenon(II12
Noscillatemorerapidly,buttheamplitudesofthelargestripples=constantFor,Nm
,sidelobe,
10.2.3GibbsPhenomenon(IV)(2)Fortheintegral
,oscillationwilloccurateachsidelobeofmovespastthediscontinuity(3)
ThemethodstoreduceGibbsphenomenon:taperingthewindowsmoothlytozeroateachend,butmasmoothtransitioninmagnitudespecifications12Noscillatemorerapi1310.2.4FixedWindowFunctions(1)Hanningwindow:
A=-B=1/2,C=0;Hammingwindow:A=0.54,B=-0.46,C=0Blackmanwindow:A=0.42,B=-0.5,C=0.08.Rectangularwindow:w[n]=u[n]–u[n–
N
–
1]
Hanning,Hamming,Blackman:
Bartlettwindow:triangular1310.2.4FixedWindowFunctio14P406Fig.10.6Commonlyusedfixedwindows10.2.4FixedWindowFunctions(II)N/2NRectangular
Hamming
HanningBartlett
Blackmann
w[n]
114P406Fig.10.6Commonlyu1510.2.4FixedWindowFunctions(III)P407Fig.10.750N1510.2.4FixedWindowFunctio1610.2.4FixedWindowFunctions(IV)Parameterspredictingtheperformanceofawindowmainlobewidth
relativesidelobelevel(dB)Sameripplesinpassbandandstopbandwidthoftransitionband1610.2.4FixedWindowFunctioTypeofwindowRelativeSidelobeLevel(dB)Main-lobewidthMinimumStopbandAttenuation(dB)TransitionBandwidthRect.13.34π/(N+1)20.91.84π/NBartlett26.58π/NHanning31.58π/N43.96.22π/NHamming42.78π/N54.56.64π/NBlackman58.112π/N75.311.12π/N1710.2.4FixedWindowFunctions(V)P408Table10.2’TypeofwindowRelativeSidelob1810.2.4FixedWindowFunctions(VI)ExampletoillustratetheeffectofwindowsN=50P4091810.2.4FixedWindowFunctio1910.2.4FixedWindowFunctions(VII)Computeimpulseresponseofthedesiredfilter(accordingtotheinverseFourierequation)(2)Determinethesuitablewindowbytheminimumstopbandattenuationand(3)DeterminethelengthofFIRbythetransitionwidth(4)ObtainthedesignedFIRfilter:StepsforFIRfilterdesign:1910.2.4FixedWindowFunctio20Example10.6Page410
DesignanFIRlowpassdigitalfilterwithspecifications:theattenuationofthestopbandshouldmorethan40dB;.2)AccordingtoTable10.2,wecouldselectHanning,hamming,Blackmanwindow,thenthebandwidthofthetransition
bandshouldsatisfy(for
Hanning)TypeI:N=32;TypeII:N=3310.2.4FixedWindowFunctions(VIII)1)i.e.Pleaseselectasuitablewindowfunctionanddeterminethesmallestlengthofthewindow.20Example10.6Page410Des2110.2.4FixedWindowFunctionsExampleShowthattheidealhighpasstransformerwithafrequencyresponsedefinedby(1)Determinetheimpulseresponseh[n],therelationofαandN?(2)Whattypeoflinear-phaseFIRfilter?(3)Writetheimpulseresponseh[n]
usingtheHannwindows-basemethod.Solution:2110.2.4FixedWindowFunctio2210.2.4FixedWindowFunctions2210.2.4FixedWindowFunctio2310.2.4FixedWindowFunctions(2)
IfNisevenwhen,thefilterhaslinearphaseisinteger,hd[n]isanti-symmetries,andh[n]=-h[N-n],thefilteristypeIII.IfNisoddisn’tinteger,hd[n]issymmetries,andh[n]=h[N-n],thefilteristypeII.2310.2.4FixedWindowFunctio2410.2.4FixedWindowFunctions(3)2410.2.4FixedWindowFunctio25with
=N/2.
β
controlstheside-lobeamplitudes(attenuation)
controlsthemainlobewidth
Predictionformula:–attenuation
s
=20log10δsβ–transitionregionwidthω=ωsωp
togetherwithattenuation
s
N(10.39’)10.2.5AdjustableWindowFunctions(P410)KaiserwindowN25with=N/2.(10.39’)10.2.526Amplitude0.305101520(10.41)(10.42)10.2.5AdjustableWindowFunctions(II)26Amplitude0.3051015227Kaiserwindowdesignexample(1)DeterminethewindowfunctionKaiserwindow:,Ni.e.,s=0.01,Assume:Question:IsitsuitableforNtobe23?27Kaiserwindowdesignex28Kaiserwindowdesignexample(II)(2)Thedesiredimpulseresponse28Kaiserwindowdesignex29Kaiserwindowdesignexample(III)(3)TheFIRfilterdesignedWhereN=24,=3.395TypeIlinearphaseFIR29Kaiserwindowdesignex3010.3CADofEquirippleLinear-PhaseFIRFiltersApproximationmethods:(2)LeastIntegral-SquaredapproximationWindowedFourierSeriesapproach(1)InterpolationFrequencysamplingapproach(3)ChebyshevapproximationEquirippleapproximationParks-McClellanAlgorithm3010.3CADofEquirippleLine3110.3CADofEquirippleLinear-PhaseFIRFilters(II)Weightederrorfunction:(10.47)or(10.62)(10.68)3110.3CADofEquirippleLine3210.3CADofEquirippleLinear-PhaseFIRFilters(III)ChebyshevorMinimaxcriterion:equirippleFIRfilterMinimizethepeakabsolutevalueofLinear-phaseFIRfiltersobtainedbythecriterionpolynomialapproximation3210.3CADofEquirippleLine3310.3CADofEquirippleLinear-PhaseFIRFilters(IV)AlternationTheorem:LetRbeaunionofdisjointclosedsubsetsofLetadesiredfunctionD(x)andweightingfunctionW(x)becontinuousonRDefinetheerrorfunction
E(x)=W(x)[PL(x)-D(x)]Maximumerror
3310.3CADofEquirippleLine3410.3CADofEquirippleLinear-PhaseFIRFilters(V)necessaryandsufficientconditionforPL(x)beingtheuniqueLthorderpolynomialundertheMinimax
criterioncanbeexpressedbythealternationtheorem:E(x)hasatleastL+2alterationsonF
,i.e.
xi,i=1,...,L’≥L+2suchthatxi
<xi+1,E(xi)=-E(xi+1),fori=1,...,L’-1andE(xi)
=±Emax,fori=1,...,L’3410.3CADofEquirippleLine3510.3CADofEquirippleLinear-PhaseFIRFilters(VI)Parks-McClellanAlgorithmIterativemethodtodeterminethealternationfrequenciesωi
andtheripple
1.initialize
ωi
pute—erpolateapolynomialbetweenthealternationpoints4.findthemaximum/minimumvaluesoftheerror5.if|E(ω)|≤
:stopelsecomputenewωi’asextremeofE(ω),andgoto2(elserecursive)3510.3CADofEquirippleLine3610.5FIRDigitalFilterDesignUsingMatlabOrderEstimation:kaiord()Kaiser’sFormulabellangord()Bellanger’sFormularemezord()Hermann’sFormulakaiserord()filterorderforKaiserwindow-baseddesign3610.5FIRDigitalFilterDes3710.5FIRDigitalFilterDesignUsingMatlab(II)Equiripplelinear-phaseFIRfilterdesign:remez()equirippleFIRfilterdesignusingParks-McClellanalgorithmExample10.15DesignanequirippleFIRfilterwithspecifications:3710.5FIRDigitalFilterDes3810.5FIRDigitalFilterDesignUsingMatlab(III)3810.5FIRDigitalFilterDes3910.5FIRDigitalFilterDesignUsingMatlab(IV)3910.5FIRDigitalFilterDes4010.5FIRDigitalFilterDesignUsingMatlab(V)4010.5FIRDigitalFilterDes4110.5FIRDigitalFilterDesignUsingMatlab(VI)4110.5FIRDigitalFilterDes42WindowingmethodforFIRfilterdesign:fir1()andfir2()Example10.15DesignaFIRlowpassfilterusingakaiserwindowwithspecifications:10.5FIRDigitalFilterDesignUsingMatlab(VI)42WindowingmethodforFIRfil4310.1,10.2,10.3estimationformula10.4multibandfilterimpulseresponse10.5truncationapproximation10.6,10.7idealdigitalHilberttransformation10.8idealdigitaldifferentiator10.9delay-complementarypair10.10,10.11,10.12,10.18inverseDTFT10.15,10.16,10.17windowingmethoddesign10.20fractionaldelayFIRfilter10.21idealcombfilter10.27,10.28differentfittingalgorithm10.29filtersharpening10.31~10.35frequencysamplingmethod10.40WDFT10.36~10.38Parks-McClellanalgorithmweightingfunctionExercises4310.1,10.2,10.3estimat44Introduction
☆
FIRfilter:directdesignofDTfilterwiththeoftenaddedlinear-phaserequirement(1)WindowedFourierseriesapproach(§10.2)(2)Frequencysamplingapproach(Problem10.31,10.32)
(3)Computer-basedoptimizationmethod(§10.3)
Chap.10FIRDigitalFilterDesign1IntroductionChap.10FIRDi4510.1PreliminaryConsiderationsForFIRsystem:realpolynomialapproximationifalinearphaseisdesired
10.1.1BasicApproachestoFIRDigitalFilterDesign210.1PreliminaryConsiderati4610.1.2EstimationoftheFilterOrderKaiser’sFormulaForlowpassFIRfilterdesign:P397-398Bellanger’sFormulaHermann’sFormulaParametersseeP398.310.1.2EstimationoftheFil4710.2DesignofFIRFiltersbyWindowing(P400)
10.2.1LeastIntegral-SquaredErrorDesignofFIRFilters410.2DesignofFIRFiltersb4810.2.2ImpulseResponsesofIdealFiltersIdeallinearphaselowpassfilterIdeallinearphase
highpassfilter510.2.2ImpulseResponsesof49ImpulseResponsesofIdealFilters(II)Ideallinearphasebandpassfilter
Ideallinearphasebandstopfilter
6ImpulseResponsesofIdealFi50ImpulseResponsesofIdealFilters(III)IdealmultibandfilterIdealdiscrete-timeHilberttransformer
Idealdiscrete-timedifferentiator
7ImpulseResponsesofIdealFi51Gibbsphenomenon:
OscillatorybehaviorinthemagnituderesponseofcausalFIRfiltersdesignedutilizingtruncation10.2.3GibbsPhenomenon8Gibbsphenomenon:52mainlobesidelobeMainlobewidth--truncationperiodiccontinuousconvolution=N/210.2.3GibbsPhenomenon(II)9mainlobesidelobeMainlobewidt53-2/(N+1)10-2/(N+1)5410.2.3GibbsPhenomenon(III)1110.2.3GibbsPhenomenon(II55
Noscillatemorerapidly,buttheamplitudesofthelargestripples=constantFor,Nm
,sidelobe,
10.2.3GibbsPhenomenon(IV)(2)Fortheintegral
,oscillationwilloccurateachsidelobeofmovespastthediscontinuity(3)
ThemethodstoreduceGibbsphenomenon:taperingthewindowsmoothlytozeroateachend,butmasmoothtransitioninmagnitudespecifications12Noscillatemorerapi5610.2.4FixedWindowFunctions(1)Hanningwindow:
A=-B=1/2,C=0;Hammingwindow:A=0.54,B=-0.46,C=0Blackmanwindow:A=0.42,B=-0.5,C=0.08.Rectangularwindow:w[n]=u[n]–u[n–
N
–
1]
Hanning,Hamming,Blackman:
Bartlettwindow:triangular1310.2.4FixedWindowFunctio57P406Fig.10.6Commonlyusedfixedwindows10.2.4FixedWindowFunctions(II)N/2NRectangular
Hamming
HanningBartlett
Blackmann
w[n]
114P406Fig.10.6Commonlyu5810.2.4FixedWindowFunctions(III)P407Fig.10.750N1510.2.4FixedWindowFunctio5910.2.4FixedWindowFunctions(IV)Parameterspredictingtheperformanceofawindowmainlobewidth
relativesidelobelevel(dB)Sameripplesinpassbandandstopbandwidthoftransitionband1610.2.4FixedWindowFunctioTypeofwindowRelativeSidelobeLevel(dB)Main-lobewidthMinimumStopbandAttenuation(dB)TransitionBandwidthRect.13.34π/(N+1)20.91.84π/NBartlett26.58π/NHanning31.58π/N43.96.22π/NHamming42.78π/N54.56.64π/NBlackman58.112π/N75.311.12π/N6010.2.4FixedWindowFunctions(V)P408Table10.2’TypeofwindowRelativeSidelob6110.2.4FixedWindowFunctions(VI)ExampletoillustratetheeffectofwindowsN=50P4091810.2.4FixedWindowFunctio6210.2.4FixedWindowFunctions(VII)Computeimpulseresponseofthedesiredfilter(accordingtotheinverseFourierequation)(2)Determinethesuitablewindowbytheminimumstopbandattenuationand(3)DeterminethelengthofFIRbythetransitionwidth(4)ObtainthedesignedFIRfilter:StepsforFIRfilterdesign:1910.2.4FixedWindowFunctio63Example10.6Page410
DesignanFIRlowpassdigitalfilterwithspecifications:theattenuationofthestopbandshouldmorethan40dB;.2)AccordingtoTable10.2,wecouldselectHanning,hamming,Blackmanwindow,thenthebandwidthofthetransition
bandshouldsatisfy(for
Hanning)TypeI:N=32;TypeII:N=3310.2.4FixedWindowFunctions(VIII)1)i.e.Pleaseselectasuitablewindowfunctionanddeterminethesmallestlengthofthewindow.20Example10.6Page410Des6410.2.4FixedWindowFunctionsExampleShowthattheidealhighpasstransformerwithafrequencyresponsedefinedby(1)Determinetheimpulseresponseh[n],therelationofαandN?(2)Whattypeoflinear-phaseFIRfilter?(3)Writetheimpulseresponseh[n]
usingtheHannwindows-basemethod.Solution:2110.2.4FixedWindowFunctio6510.2.4FixedWindowFunctions2210.2.4FixedWindowFunctio6610.2.4FixedWindowFunctions(2)
IfNisevenwhen,thefilterhaslinearphaseisinteger,hd[n]isanti-symmetries,andh[n]=-h[N-n],thefilteristypeIII.IfNisoddisn’tinteger,hd[n]issymmetries,andh[n]=h[N-n],thefilteristypeII.2310.2.4FixedWindowFunctio6710.2.4FixedWindowFunctions(3)2410.2.4FixedWindowFunctio68with
=N/2.
β
controlstheside-lobeamplitudes(attenuation)
controlsthemainlobewidth
Predictionformula:–attenuation
s
=20log10δsβ–transitionregionwidthω=ωsωp
togetherwithattenuation
s
N(10.39’)10.2.5AdjustableWindowFunctions(P410)KaiserwindowN25with=N/2.(10.39’)10.2.569Amplitude0.305101520(10.41)(10.42)10.2.5AdjustableWindowFunctions(II)26Amplitude0.3051015270Kaiserwindowdesignexample(1)DeterminethewindowfunctionKaiserwindow:,Ni.e.,s=0.01,Assume:Question:IsitsuitableforNtobe23?27Kaiserwindowdesignex71Kaiserwindowdesignexample(II)(2)Thedesiredimpulseresponse28Kaiserwindowdesignex72Kaiserwindowdesignexample(III)(3)TheFIRfilterdesignedWhereN=24,=3.395TypeIlinearphaseFIR29Kaiserwindowdesignex7310.3CADofEquirippleLinear-PhaseFIRFiltersApproximationmethods:(2)LeastIntegral-SquaredapproximationWindowedFourierSeriesapproach(1)InterpolationFrequencysamplingapproach(3)ChebyshevapproximationEquirippleapproximationParks-McClellanAlgorithm3010.3CADofEquirippleLine7410.3CADofEquirippleLinear-PhaseFIRFilters(II)Weightederrorfunction:(10.47)or(10.62)(10.68)3110.3CADofEquirippleLine7510.3CADofEquirippleLinear-PhaseFIRFilters(III)ChebyshevorMinimaxcriterion:equirippleFIRfilterMinimizethepeakabsolutevalueofLinear-phaseFIRfiltersobtainedbythecriterionpolynomialapproximation3210.3CADofEquirippleLine7610.3CADofEquirippleLinear-PhaseFIRFilters(IV)AlternationTheorem:LetRbeaunionofdisjointclosedsubsetsofLetadesiredfunctionD(x)andweightingfunctionW(x)becontinuousonRDefinetheerrorfunction
E(x)=W(x)[PL(x)-D(x)]Maximumerror
3310.3CADofEquirippleLine7710.3CADofEquirippleLinear-PhaseFIRFilters(V)necessaryandsufficientconditionforPL(x)beingtheuniqueLthorderpolynomialundertheMinimax
criterioncanbeexpressedbythealternationtheorem:E(x)hasatleastL+2alterationsonF
,i.e.
xi,i=1,...,L’≥L+2suchthatxi
<xi+1,E(xi)=-E(xi+1),fori=1,...,L’-1andE(xi)
=±Emax,fori=1,...,L’3410.3CADofEquirippleLine7810.3CADofEquirippleLinear-PhaseFIRFilters(VI)Parks-McClellanAlgorithmIterativemethodtodeterminethealternationfrequenciesωi
andtheripple
1.initialize
ωi
to
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
- 6. 下載文件中如有侵權(quán)或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 學校秋季讀書節(jié)活動方案計劃
- 倉庫服務質(zhì)量的提升路徑計劃
- 前臺文員的職業(yè)危機應對策略計劃
- 倉庫應急預案與事故處理計劃
- 人教版小學英語教案
- 冀教版七年級英語下冊單元測試題
- 備戰(zhàn)夏季無意中發(fā)現(xiàn)這個東西很好用
- 新版車間安全培訓試題及答案 完整版
- 項目管理人員安全培訓試題帶答案(完整版)
- 公司員工安全培訓試題附答案(突破訓練)
- 市政道路及設施零星養(yǎng)護服務技術(shù)方案(技術(shù)標)
- 2024中水北方勘測設計研究限責任公司校園招聘【重點基礎提升】模擬試題(共500題)附帶答案詳解
- 城市燃氣管道等老化更新改造項目初步設計說明
- 三年級數(shù)學上冊課件-5. 倍的認識 - 人教版(共16張PPT)
- 綜合實踐活動創(chuàng)意彩泥(課件)-2020-2021學年三年級下冊(14)ppt課件
- 中小學幼兒園安全技術(shù)防范系統(tǒng)要求.doc
- 礦產(chǎn)地球物理勘查項目野外驗收意見書
- 銀行系統(tǒng)論文:農(nóng)信社信貸違規(guī)的表現(xiàn)形式及防范對策
- 監(jiān)控弱電工程施工技術(shù)方案(完整版)
- 對外投資業(yè)務流程圖(doc 1頁).DOC
- 熱食類食品制售操作流程
評論
0/150
提交評論