基于burg算法的谱估计研究及其matlab实现学士学位论文Word文档格式.docx
《基于burg算法的谱估计研究及其matlab实现学士学位论文Word文档格式.docx》由会员分享,可在线阅读,更多相关《基于burg算法的谱估计研究及其matlab实现学士学位论文Word文档格式.docx(20页珍藏版)》请在冰豆网上搜索。
开题报告
机电系电子信息工程专业
学号:
学生姓名:
指导教师:
(职称:
讲 师)
(职称:
)
XXXX年XX月X日
课题来源
功率谱估计在近30年中获得了飞速发展。
涉及到信号与系统、随机信号分析、概率统计、随机过程、矩阵代数等一系列学科,广泛应用于雷达、声纳、通信、地质、勘探、军事、天文、生物医学工程等众多领域。
科学依据(包括课题的科学意义;
国内外研究概况、水平和发展趋势;
应用前景等)
实际中,数字信号的功率谱只能用所得的有限次记录的有限长数据来予以估计,这就产生了功率谱估计这一研究领域。
功率谱的估计大致可分为经典功率谱估计和现代功率谱估计,针对经典谱估计的分辨率低和方差性能不好等问题提出了现代谱估计,AR模型谱估计就是现代谱估计常用的方法之一。
随着信号处理理论和电子技术的迅猛发展,谱估计分析在已经在国内外各个领域中得到广泛应用。
例如在现代军事领域中,人们通过分析雷达信号,根据回波信号的功率谱密度、谱峰的宽度、高度与位置,可以确定运动目标的位置、辐射强度和运动速度。
在引信系统中,通过对探测器接收到的信号进行谱分析,可以确定弹目之间的交会速度、方向,从而提高引信的打击精度和毁伤效率。
现代谱分析受算法的复杂度高和集成电路处理速度的影响,因此在引信系统中的应用会受到很大的限制。
但随着电子技术的大力发展,特别是DSP技术的不断发展,现代谱分析的方法在引信系统中的应用逐渐成为可能,并且极大地提高了无线电引信引战配合性能。
研究内容
①熟悉谱估计的发展历程;
②熟练掌握经典谱估计方法:
直接法和间接法、它们之间的关系、估计质量以及估计性能比较;
③熟练掌握现代谱估计方法:
信号建模、AR模型参数求解的Levinson-Durbin算法和BURG算法,阶数的确定方法和原则,稳定性以及对信号建模的讨论;
④熟悉熟练使用MATLAB仿真。
针对一个具体的随机信号,分别采用经典谱估计和现代谱估计方法估计出其功率谱,对经典谱估计和现代谱估计方法谱估计的分辨率和方差性能作一个综合评价;
⑤熟练使用MATLAB提供的图形用户界面(GUI)工具。
拟采取的研究方法、技术路线、实验方案及可行性分析
熟悉经典谱估计和现代谱估计相关的理论知识,并掌握MATLAB的使用方法,达到熟练使用并融会贯通。
应用有关理论知识在MATLAB上进行仿真实验,然后对实验结果进行对比和分析。
最后总结归纳实验结果,提出自己的观点。
研究计划及预期成果
研究计划:
2009年10月12日-2009年12月25日:
按照任务书要求查阅论文相关参考资料,并填写毕业设计开题报告书。
2010年1月11日-2010年3月5日:
进行毕业实习,并填写毕业实习报告。
2010年3月8日-2010年3月14日:
按照要求修改毕业设计开题报告。
2010年3月15日-2010年3月28日:
学习并翻译一篇与毕业设计相关的英文材料。
2010年3月29日-2010年4月11日:
MATLAB程序设计与MATLAB程序调试。
2010年4月12日-2010年4月18日:
GUI设计。
2010年4月19日-2009年5月21日:
毕业论文撰写和修改工作。
预期成果:
了解谱估计的发展进程,掌握各种算法,并熟练掌握MATLAB的基本操作。
训练自己综合运用所学知识的能力,培养自己独立分析和解决问题的能力,启发自己奋发向上的精神。
特色或创新之处
由于经典谱估计中将数据工作区外的未知数据假设为零,相当于数据加窗,导致其分辨率降低。
另外,功率谱密度原始定义中既无求均值又无求极限,所以使得经典谱估计方法的方差性能比较差,尤其是在数据纪录很短的情况下,这些问题更为突出。
而现代谱估计则不再简单地将观察区外的未知假设为而零,而是先将信号的观测数据估计模型参数,按照求模型输出功率的方法估计信号功率谱,回避了数据观测区以外的数据假设问题,因而其谱的分辨率得到提高。
可以看出现代谱估计性能优于经典谱估计。
已具备的条件和尚需解决的问题
已具备的条件:
可通过书籍和网络了解功率谱估计相关的理论知识,并已学习过MATLAB程序设计。
尚需解决的问题:
运用MATLAB进行仿真,并分析结果。
指导教师意见
该生查阅了大量的相关资料,设计方案可行,同意开题。
指导教师签名:
年月日
教研室(学科组、研究所)意见
教研室主任签名:
年月日
系意见
主管领导签名:
英文原文
ParametricspectralestimationonasingleFPGA
ABSTRACT
Parametric,modelbased,spectralestimationtechniquescanofferincreasedfrequencyresolutionoverconventionalshort-termfastFouriertransformmethods,overcominglimitationscausedbythewindowingofsampled,timedomain,inputdata.However,parametrictechniquesaresignificantlymorecomputationallydemandingthantheFourierbasedmethodsandrequireawiderrangeofarithmeticfunctionality;
forexample,operationssuchasdivisionandsquare-rootareoftennecessary.Thesearithmeticprocessesexhibitcommunicationbottleneckandtheirhardwareimplementationcanbeinefficientwhenusedinconjunctionwithmultipliers.Aprogrammable,bit-serial,multiplier/divider,whichovercomesthebottleneckproblemsbyusingadatainterleavingscheme,isintroducedinthispaper.ThisinterleavedprocessorisusedtoshowhowtheparametricModifiedCovariancespectralestimatorcanbeefficientlyroutedonafieldprogrammablegatearrayforreal-timeapplications.
1.INTRODUCTION
DuetoitseaseofhardwareandsoftwareimplementationtheshorttermfastFouriertransform(STFFT)iswidelyusedforspectralestimationandisknownastheconventionalmethod.However,thetechniquehasdrawbacksintermsofspectralresolutionandaccuracycausedbythefinitelengthoftheinputdatasequenceused.WindowingofinputdatacausesspectralbroadeningandGibb’sphenomenonofspectralleakagecanmasktheweakerfrequencycomponentsofthetruepowerspectraldensity(PSD)[1].Theseunwantedeffectscanbereducedbyusinglongerdatasequencelengths,sothatthetransformedsignalbecomesabetterrepresentationoftheinfinitedatasequence,butinreallifethisusuallyisnotfeasibleasthecharacteristicsoftheinputdatamaychangewithtime.Overshortperiodsoftimethedatasignalscanoftenbeassumedtoexhibitwide-sensestationarity,wherethesignalcharacteristicsareassumedapproximatelyconstantbutthespectralresolutionisthereforelimited.InattemptstoimprovethePSDestimation,windowingfunctions,BartlettorHanningforexample,canbeusedtoreduceside-lobelevelsbuttheselowerspectralresolutionbybroadeningthemainlobeofthePSD[2].
Modelbased,parametricspectralestimationtechniquescanalternativelybeused,wheretheunrealisticassumptionthatdataiszerooutsidethewindowofinterestisdropped[1].Eitherknowledgeoftheunderlyingprocessorreasonableassumptionsaboutthenatureoftheunobserveddataareusedtoimprovefrequencyresolutionovertheconventionalapproaches.ThecomputationalburdenofsuchprocessorsishowevermuchhigherthantheSTFFTandarithmeticfunctionssuchasdivisionandsquare-rootoftenbecomenecessary.Inthedivisionandsquare-rootnon-restoringalgorithmsthereisaninherentdependencythattheresultbitsmustbecomputedinamostsignificantbit(MSB)firstmanner,withthecomputationofabitdirectlydependentupontheresultofthepreviousone[3].Thisinterdependencymakesitdifficulttoefficientlyrealizesucharithmeticfunctionsinhardware,andimplementationsareusuallymuchslowerthanotherbasicfunctionssuchasmultiplication,additionandsubtraction.Communicationbottleneckscanthereforeeasilyoccurinsystolicarrayswheredifferenttypesofprocessorsareinterconnected.
ThedifficultieswithhardwareimplementationofparametricspectralestimatorshaveledtoapreferenceofsoftwareimplementationonhomogeneousDSPnetworks[4].However,highlevelsofprocessingcapacityhavenotbeenfullyreflectedinsystemthroughputsincetheincreasedcommunicationincurredasaresultofparallelismisconstrainedbycommunicationbusperformance.Thisrestrictstherangeofproblemsthatcanbecomputedinrealtimeandthesoftwareapproachmaysometimesbeinadequateforreal-timespectralestimation.
Inthispaper,hardwareimplementationofaparametricspectralestimatorisaddressed.Abit-serialprocessorcapableofdivisionandinnerproductstepcomputationisdevelopedbycombiningseparateprocessorsforthesefunctions.ThedesignusesahighlevelofpipeliningsothatdivisioncanbecomputedatahighrateandmultiplicationisperformedonaMSBfirstdatastream,eliminatingthebottleneckproblem.Thehighlevelofpipeliningallowsmanyindependentcomputationstobeperformedsimultaneouslyorinterleaved.TheuseoftheinterleavingschemeisdemonstratedbyimplementingthedesignofaModifiedCovariancetypeofparametricspectralestimator,toproduceafieldprogrammablegatearray(FPGA)basedsystemforthespectralanalysisofDopplersignalsfromultrasonicbloodflowdetectors.
2.MODIFIEDCOVARIANCESPECTRALESTIMATION
Themodelorderp=4ModifiedCovariance(MC)spectralestimator,proventobeoptimallycostefficientforthebloodflowapplicationwheremeanvelocityandflowdisturbanceareofinterest[5],involvessolvingthefollowinglinearsystemofcovariancematrixequations:
(1)
whereeachelement
isobtainedfrom:
(2)
forawindowoflengthNdatasamples.The
filterparameterestimatesareobtainedbysolutionofthelinearsystem
(1),usingtheCholesky,forwardeliminationandbacksubstitutionalgorithms.Thesignalwhitenoisevarianceestimate,
iscalculatedas:
(3)
andthepowerspectraldensity(PSD),
isobtainedfrom:
(4)
Hence,theMCspectralestimatormaybepartitionedontofourdifferentprogrammingmodules:
CMR-calculationoftheelementsofthecovariancematrixandright-handsidevector,5Nmultiplyaccumulatestakingintoaccountmatrixsymmetry.
Cholesky-solutionofthelinearsystemofequations,6divisionsand10innerstepproductsfornon-square-rootCholesky,4divisionsand12innerstepproductsforsolvingtriangularsystems.
WNV-calculationofthewhitenoisevariance,4multiplyaccumulates.
PSD-computationofthepowerspectraldensity,4NinnerstepproductsforazeropaddedDFT,NmultiplicationstofindabsolutevalueofDFTandN/2divisionsforthePSD.
Thenumberofsamples,overthefixedtimedurationwindowof10ms,isrequiredtobeeither64,128,256or512dependingonDopplersignalconditions.ImplementationofthealgorithminMatlabsoftwareprovedtobeinexcessofafactorof
timestooslowforreal-timeoperationandthataperformanceofupto13.5MFLOPS/sisrequired[4].ExecutiontimesofMCalgorithmimplementationusingvarioustopologiesofTexasInstrumentsTMS320C40DSPswithT8transputersasroutershavealsofallenshortofthereal-timerequirements,whereprocessingtimeisover150mstoolongintheworstcase[4][6].UseofasingleDSPinaPChostedsystemhasbeenshownsufficientforthesmallerNbutthespecificationofN=512couldnotbeachieved[4],thuspromptingconsiderationofthehardwareapproach.
3.BIT-SERIALINTERLEAVEDPROCESSOR
Studyofword-parallelsystolicimplementationsoftheMCmethodhasshownthemethodtoprovidemorethanadequatethroughputforthespecifiedreal-timebloodflowapplicationbutthecostofsuchasystemisveryhighintermsofarithmeticunits,communicationburdenandcontrolcomplexity[7].Forexample,asystolicarrayprocessorfornon-squarerootCholeskydecomposition[8]requires13processingelements(PEs),eachPEhaving2to6portsofeitherm(singleprecision)or2m(doubleprecision)lines,andcontrolisnecessarytoreversedatastreamsbeforebacksubstitution.Analternativewaytoapproachthehardwaredesigninvolvesconsiderationofbit-serialprocessingtechniques.
Thenatureofmultiplicationalgorithmsnormallyinvolvethecomputationofleastsignificantbits(LSBs)firstandbit-seria