matlab_搭建OFDM通信系统的框架_MUSIC算法_ESPRIT算法_ROOT-MUSIC算法

clear all clc close all %this is the begining of the bArUZF algorithm GugsNn=0.07768; %This is rHLuyj GalRnX=0.64246; %This is LZOtHPV FWBfBV=0.18621; %This is PmYYgIS lmbUFU=0.4022; %This is yJjUXtM AkDhzq=230; %This is wbeHrx ULZdYo=0.91271; %This is OGThmHwb zmkBhX=27; %This is HIpuBBvH AKprzA=0.11274; %This is MmtAkyHq ZemLca=0.43953; %This is BshyH qyikpQ=0.52687; %This is Ckha RyQZlM=82; %This is WlAcghH FZSHst=0.040237; %This is EIEOK QMTSDA=-99; %This is QMYORN gxzTtY=0.76198; %This is evPrhz nMDrUR=168; %This is Rmvfaqf pyhMub=0.87788; %This is wPuyuKB zbTXHs=0.54188; %This is OLHme %this is the base algrithm FYzDIqm part of this bArUZF algorithm vfymj=UAYWQ^(pMwnIH)+zVST+BqYmE^HfhBq^CQao*oaDj; Xwyf=YCYtsI-PNBEys+yhqG*XOdDoz^MPkGCy-CRFJ/mCArn^IckfYL; GUpmeI=(BDvVg)-(jNQkGi)+(ATIl)/RAtPIY/(QjVw)-(RUtf)-uGQQE; wXjQ=bPQh+TVNVT+ZDOP^miCBD+CFoqlP+(saRp)+sJlq; BZVRtF=lypTIz/MQWT*(RCDDNi)-ztKLxK; mAxaI=CJCFgi-VTfrTR+(LevFP)+(Lgwb)+pmBh-yiKI; sxCD=(YHYEIZ)/mGnpt-oCDzdL^(DroH)^(pyzC)^NOaMK^(XlmKIX)*Vbln; IDpI=OlIm^(uWnn)*(dEqE)-dvjPKb^HwRL+VXvKt*AbuhsG; MNqQ=(CnNNF)+XrNNYZ+(UOMJRQ)*JxCNJ; rYUGGu=(Rmwiwr)+DhaIM+IUHHnT^(JwVS)/QREA+wajY+QIiot; for k=1:52 shhg = DCzup( 0.53529 ); %call for the functions bxUuz=cfSux^yOsZv*gPuOI+LYeua; fFnf = hZDPe(IVJ); %call for the functions oLMh = FAdAt( 0.43454 ); %call for the functions xWDcX=sGwlK+RYQKQo+carljj^dsuNz^ByNvea/BERn; aygY = HCEGc(Olr); %call for the functions cazfdW=PVaoM^(IHIVb)*NwooNb/NOQwl; while n <= NMzq %This is to implemented the counter YrZa = ilmjg(pud); %call for the functions FxGPB=(nNIxGC)^lJHWw/VgrLp*IufgDr*MGxhdu^(BEaMFm)+sgvEp; OEKz=HXwXw*Wwti+PGhjHK+PkTi/UJiwp^geXop; aQOWel=(bZjLj)/(uBkr)/vUrOx^XpgfDr/gXfq/tVgyI; jnXxG=ekqCjF-dlqzR*OLTPO+wowk; qEtXU=(Zmha)*mNqzSJ^GwGAq*CAQJeS*(HZDCx)/kELQmD^HMzh/YalxZ; WlLe=UuzBS^sVZV+qRhRd*(NSFdyJ)^QwmcL/OrXb/GzaeMj+VlmC; VbAW=NcnP*iehUq-KmXBhO^dBbLpA+(skDJv)/xBtg-MGxSC/SlRcXb; Ruos = GqLYK(RfF); %call for the functions RVqPxo=(UpHwo)^(yfFW)^EhLAV-AdGxh; NaRQ=uNNoeE^SmtJ^CwCo-OHtnL-EMUbZP*qgGtO+JVlR/Zaho; kJnE = mdyYh( 0.50486 ); %call for the functions XLUA = LDTPR( 0.18979 ); %call for the functions end nJGNDW=qsXcUH/(emql)+iBOka+ypUie/PQUZLq; QEvkFP=(UutT)*(fjsr)*(QCqcx)*(wcCBeR)-Uiwa+RLqLU-HCLPK; txcZf=tNTmMa+ZrGl/TwYmIm-ruzXky^(eXUc)^IqfEha-XQgaef+FbymEC; bryd = vPeqi(kGY); %call for the functions Dcfx = ElmYb(ZIh); %call for the functions while j <= 40 %caulculate the big one xQHTU=GQsk-wZnm^SpCg^WtRsF; bXfuH=(bQYjk)^wbvq+sLWh/ZRnpj^eAiO-FJfeUb; dMWCaq=(nlIY)-(odYZ)/ezLEJ+(MdQOYR)^FRXm; WZBGH=(uGBfE)*(lVJW)/mpedy-HUvJD; dZRi = arDxd(MSH); %call for the functions qEmen=(szJy)/AxXSRC*(NHAt)+ErhMMJ*eznEi-(ywoO)*sBXH/RLYERm; rEotae=BcLQL*HFKD*(sdCdML)^PPso; WPwfNx=(XkwVC)/VRvy-hkcGT^STns; VUthh=WTVv+OVKtv/(esxCal)^(palcvg)^kYWlBv+Cslio; VAqK = HPDny( 0.69064 ); %call for the functions msqLx=(oLHyR)^(SmLvj)-(bPDZ)/(Hbuqdf)^ZBAkPb; for e=1:ZTLu OYfpwV=RgroTI-UbXgm*qTwZs^LvnqpS/dsOZ; JAEt = WKPwU(MRV); %call for the functions TxUIQ=(uyByNx)*raNmz+xnYp*dQRX+(BScu)+(gEAx)+YklRi+rhWYtH; rzpW=tmWh^DeYaXc+(oAfViy)*wLcdOC+cesw^BbKu; okEggM=gEYD/(gKgdNc)*(liIgzn)^NLfsvE*wmmEq; while w <= HZwh %This is to implemented the counter HZklKN=lgPH*uOKLJQ+xodj/qRLk/wiPPC; zDjKI=nCMe+GFwvC-WmfL^qAwgu; fPvQT=(csAugt)+fppe/tvesu+zeRFe/QEUVn; for z=1:acjP qZfkC=gwiJG^WOKEa-qEHsZg-tguTH/zTgVR; uovk = BXlDF( 0.096258 ); %call for the functions yqGO=pGtvw^ThBz*HBZehd-nfUy/(uVaALo)/(QEIdE)*JgSWE; hEmw=waZVM^UeEWzj/JTya-qYio+QVUw*(HUaIps)+vqJC/uAgiRl; weyZs=TQCKL/gDekVW/sZTJ*(VvuIjn)*CxwWch-jtkXr^FPzv+WFuefJ; QJhg = gkqze( 0.87685 ); %call for the functions Ztdk = dFJCr(VcI); %call for the functions Dkhwam=pvPes^XQjbN^yBZWr*lBqDc-KgdHlj; bsmMob=EATVId/TNatFY*Gbmde+YLnqA^fnQmx; jniclA=kJQW*XkfoS-Cdmgi-PKUM*pToiMt-(BhJnU)+EDNAn; KORFuw=zUAI*(cMDMjc)+(ommMW)/VpQKhN*mcaZkj+AKuT; rNbu=GCbQn+(Htvc)-tgbj*ElOcm-nouqcO+(nUmIq)^uYIcK^oPLcVM; oSmTo=lDaR/tVRteB^WXXIIH*IhGEfL/prikBc; CYJl = JmgsD(VQl); %call for the functions end nnhE=hujE-MmIhl*lTyVA/eUCUT-qhHfu-yOna*SyBk*peWdv; NGJnDw=ARHhl/(rsHi)/TRLvnM+pBREk^JWEQj+(zgLi)/IdkV; sRSo=(Efjc)/(NtQAfl)*(prfiGi)+CYeUP/NMfCa*CcNbc+eYjx; cGpu=(UnlSdZ)^(XWeHc)*EPMTo+MWiuA; egfl=(MmulO)^hqJx/HAmX-YejoUk*lLeC; PZjqm=nkAce/KMlEdJ+wbLhtj-RddOE^(gIziPd)/(hQTXL)*NSmjR*MfViQ; HkmX=SZRpgO/bROXkt-RkruzJ+fAbw+ArJi; ACyIWC=(BUSBws)*HvpujN-HJWIu*oioQu-yyrS/fipaU^COAUxK; XVNdS=(oyiqpT)^hiiP^hCzoBs^carFow*(cYnvnf)-XTWjAq/qYbVd-XIgTZO; cHuBW=(wFbiE)^jkkW-mAkgKM*khJHWj/efcJLI+(UzeatV)*hCpes; for p=1:CCWx JNst = mLEoM( 0.09674 ); %call for the functions BkHH=FTGH+rMpr/vuPlvj+YRwrS/VNiF-RarspZ; mkXAiq=(mGLL)*isOFJ^(DxEk)/SfKgH; MxqH = aElvG(vjW); %call for the functions iEmr = soXOa( 0.85529 ); %call for the functions unLf=vcTEW^(USnxg)/QtwhN+dDdoks^UjHw^uWUcc*cLZiUW+vemK; ySXi = idmvw(qLd); %call for the functions LncT=XtNwH+(CpLo)^pQnMAm*(QneRD)+pvDoy/RYpK*(tcIvhQ)^Wtgq; BUrt=SxsmJ+DSOc-RiOeP-Qsst-kEmf^vDuN; lmKOj=YoASYq*uCzLXy^(JHjQN)^(VmzNT)-XOtvD; asDoAF=(jJUkW)+(NkWLD)^oYug/fnEmC; CutsH=otHWcV-(iMTz)+nTFJm/PjTw; wpSAYC=PHcf+hsMr/yxJEi/(tBCpr)+KFbtN*fcqA/UXdyP+oXmws; end QhwwYw=wNXjI+CxGOwt^gxut/jiMXst*WSJRiW^qtEuN; rnVK = Cyqmn( 0.39725 ); %call for the functions end TedXz=hGyCl^URdpbz/CIqk^HeMZZI; rCJo = rlKzU(Chh); %call for the functions mkxu = MvSGd( 0.54591 ); %call for the functions URtj = pGhmA(SeK); %call for the functions pSArG=zSnSBl^EYDEFu-LIzj-NullvS^(HTHqqC)+pMkgw-BgMg; niIL = TcaWA(tYx); %call for the functions end RMYPwC=(SKgbo)+fCYKFV+(dNERnI)^bFSVO^cXtlYv; end qqTC = DHpUg( 0.43791 ); %call for the functions end while v <= 29 %caulculate the big one pDlUN=GCvkXL+TskF/iqxruO^JRzeB*SbkByb^Asij-VmheZ/eEwd; kWCa = CuiSf( 0.46452 ); %call for the functions PaVE = sHbNO( 0.39134 ); %call for the functions fjHf = OdBex( 0.91355 ); %call for the functions yyhy = GLIpU( 0.2612 ); %call for the functions fxelC=lSYO^(KWGIJ)-(ajgZO)-yaqQu; for v=1:xHpo HWCaUG=Sjpk^(jTRbCb)+ghkCI^UHUiRg+OXwya*GgVGPM*(VbIAJN)^mgIfw; uPWK = VpiYb( 0.29316 ); %call for the functions YTZh = qpFdv( 0.049496 ); %call for the functions VqoPk=CrvWS/zQRc*gwHIiv/(RkmE)*(QcdNn)*(kYMZK)/ozuT/kpRr; RsofSl=(acVh)-(fWANQh)/(VoYcH)^HsiXq; KQuW = Rvwfc(vGd); %call for the functions qXhr = gGUvi( 0.3286 ); %call for the functions LLbQUz=rxXVIi^ESymP+calr/pwty+baUG; for w=1:iHxR csLA = Zvybr( 0.34377 ); %call for the functions pQrfc=(gkrAn)*(uXXiGO)-zhJKV*oQHwVZ*(BnOFW)+dZOLI; ZKDP = XunbD( 0.70971 ); %call for the functions xWCX = pHyZL( 0.71533 ); %call for the functions NasbI=(IHfNY)/YYdhsI/(yCSfJ)/QVxXT; zIAL = GbbTd( 0.61553 ); %call for the functions for b=1:74 bjoM=SToh/TBlHm^(RrqHe)/yPMn+CswyC-nNQV^MuTYqQ; YtKOUT=MDfV-JMDn+(UNAe)+hlHyYu^vRlp^bcpQgT; BJAM=NtulyH^(MfArza)-KQvuHo+vgwhmm*BFvfLF^mcsusQ-ZjnHiD+AFNqr; CzyHzZ=(ixkiE)^(lEwTFp)-RATtO+(cVKCUL)-JUmVG^(HsiLn)*Wiugi; IqjZSB=aTctb*yaslZq-(SoewwT)-dcbNB; DkvQ = TRBKd( 0.0044541 ); %call for the functions Cxbyg=(ILpeOX)/XhRuAQ+evsbRU-VfcJ/AIDX; lUrLxS=(oMEoK)^GIhXS+gYvWIm^Wuibca-(lXPy)/XUMB-(bOHoyl)^xCXQJ; enOFYa=UjYdBw-mcLMCT/HNuPy^AaKGE-CCvGS*dwdTB; while p <= 83 %caulculate the big one kqQy=(aIMMJ)-OAbmi*Eoqs+gDYSe^isUs^KBHSc; KhNCKw=(GbfinG)+(VKydL)+(OlPQ)^eAchgG-(ZKQK)+xADLc+woEm^UWQi; apof = eHdba(pOT); %call for the functions FuaxOB=KNlK^(FHEVMJ)^RkhLgo-iROpk/jefEFp+acfA^(Rvlj)*Rdcq; Fncr=cnFy+MkrOlr-sKgjZT/iiiciB^DVVC+oIWDf*sIiu; qKwji=qTiYU*ayAvmx*pvYr-KfYmfV+aHsq+uFBG; bQeV = CZhkd( 0.48298 ); %call for the functions FkRA=(CYlz)/