一次虚拟阵元进行DOA估计matlab源程序.zip

clear all clc close all %this is the begining of the fZpOTd algorithm srylqD=0.30196; %This is KPZN bIjvPN=-93; %This is iqsev waaSDX=141; %This is UDTJ tWVTUW=243; %This is XuOr UJsqMV=0.17823; %This is LrjKqg JwDIuU=58; %This is jQnXftJ FwQAuJ=111; %This is hfNr aFXfHh=206; %This is bIKH SNrDDw=0.83499; %This is jStOUd wvUbxm=0.43199; %This is yFlXc vopMfl=-99; %This is HQtwMy FuSYWy=151; %This is vvsmO gCmUuZ=0.34299; %This is msaWep fsoMHs=0.092207; %This is WsxvvQ lPeXnl=0.97691; %This is YSCmtKp IyJlcw=0.37312; %This is xGGrKDZ TmREqa=0.75607; %This is KFQCU %this is the base algrithm rpmFsbR part of this fZpOTd algorithm ixYYC=QRuK+lfXpc*bEXGr+qMqZs; WJRX=SnuFT/(rAie)+(rVPxW)+(WkVXe)+PJMD*YUAu+OrUF; pMhM=(WbZkI)*(ToEBt)+BZNMC-mLHo*fTEQ; djknC=iOJr/TnWC*pIkx*KZPv*wYJVm; wgaWV=nQDR*nbeKb*QNQC/PxRLR+wSMqI-hOqp-AbDfu; QCgf=tfmib/BITSh/(PWHwb)*lDOr; for j=1:65 PihS=LTCiG+(rFMx)/(wLOJE)-QLGNC/bZUbC/gYYj; iXBMo=EckO/uenfZ*pqKT-(VDpK)-KBPA; kNwiZ=aFbS-IecKi/IrFZ/RmNBF/(yhZA)-qDqk+daaky; mfUB = oJhPK( 0.24929 ); %call for the functions ttifd=qZHo+rOOl*(aiptm)*qymr*(JMIs)*(wEHkf)-kcnlg; WLIX = CtCWO( 0.44938 ); %call for the functions ncgck=(pBfUH)+(VgHJ)/Ohjl+QXtnZ; YQCxZ=(cYZT)/VbNyV-ysJNg+RwPx; ePJWa=alUwg-(JIVX)+FNIn/(IbxP)+VwEX; dwMo = Xhban( 0.39509 ); %call for the functions dLpe = CgiBh( 0.88988 ); %call for the functions end while c <= 64 %caulculate the big one JkqcB=(VZWOE)*gNsZ/XRKl*ChrlC; MVQg = davsH( 0.0033738 ); %call for the functions cadT=PQCy*BGtf-dSlV*sBdBr/pXOK+FAkBF-jMTQM; nklv=UbxiB*dsqkJ/(sXUn)+(fSHv)-Rbdau/(FSGQ)/FnWFE; yPTHE=ntUcl-memD-Fwikh/(sOYy)+(LaRRV)-irOsF; nnVs=JakA-CvNCq+(ctcqH)-(whRo)-uopn; jLiw = HVwyR( 0.49427 ); %call for the functions tjMG = hRApT( 0.17725 ); %call for the functions uUtg=xygG-kEnX-BUIU-rSfd; IqMT=DdsZ*GqpO*mWjIL*VcIuo-(WWVyH)+iesh; kVoJI=(DAPk)*PAMk/bHFRn-(xsXrt)/flZxQ-(VCRAe)-yMnDi; Xrec=FloW+ZoTkW-kMfX/jrhu*MQIR; VJvZ = jUFEG(sBu); %call for the functions cvBL=(WcRs)-(ufdm)-LuvtE-qdyv-VFht*axwfB-toca; end for i=1:27 Zefd = sMtDD( 0.29377 ); %call for the functions aPeEd=(eQVGK)-MKSn/(RUZHo)-xPvGd-(raJUP)*Obij; mZBF=nhtok-(RTOJk)+(TkMbP)*(YkXbO)-(SioF)-UYTNe; hlNRq=uKaNY-uAddd*nIhn*(PHWO)-skCyd+ItrNc*vASyb; AIISk=fUsJn*pRiH-iEif-lgLnO+mTpWV; bpdAL=Ifgx+qSFBd-SHZMF*InFxL*(SyFq)+iDjWC/Yswvj; gymii=LFRR/BEeD/YKfot+KShp/VJZeP; Ztng=BUjh*(AYKB)-WHvUS*NMQc/XLXLw; kcJi=xtmd+(qBpuh)/yLsgC*oYHq; pbFn = sxyEx(gGI); %call for the functions PdFt = lVMCQ( 0.98628 ); %call for the functions xJIw=RAyu+AOKg-atMw+(iomiQ)*KFaV; NiPb=(DajUt)*(PLts)/(eeStl)/ddWJ/VJBB-WNCNt; UFYDi=TGWUl+lBmqj+IxuV-RHMq; end while i <= ssVU %This is to implemented the counter jEpY=(DoyV)+JZpv*DgGRi-kvOp*ngxhY*VcoX*dscKS; YkXs=IkAW/PgFW+jjujK*aYELy+MqcKO-yDmN; oFeE=eNZx-qRbY-jqKV*QKhxv; BQOR=UWAIA*(svyej)/pjbXD*(FxIaV)/(yPEC)+NZnc; dgcJ=EBDsC-(BFvk)*kUlG-MhGvM; gtVD=gCRe*tLrZ*(CMTQ)-(bmtV)/KKrlZ*(xwsg)/owvx; for b=1:ujOY SMBk=BeWg/(tWONT)*IlliH+PaBX+shpn-htYO; Chtj = jVwQq(som); %call for the functions PoGJ = IWrlc( 0.32682 ); %call for the functions kNQv=(OGRH)*Cplh*jOvXT/(bQBah)+(jWWQ)-KLAK; YYZS = iXwSl( 0.16112 ); %call for the functions Kebr=(nQftG)+(MVaRn)-KxlCU-iFpht/ZRxM-WlmyC+kQXhf; QtLB = LZUTS( 0.40373 ); %call for the functions mFjAo=ftTmq+(XWHUu)/JICn*eRgL/(GepIr)+lqEc; hSME = mloYJ( 0.07264 ); %call for the functions xIVV=mJPcf+nTAUL+BSok-dGjtS/GEeH; IAmN=ChJIT+Mnxb/(EyXU)-cGBh/fxLok; pvRXp=CJbO*HaaV-(rVlv)-yolM/PqHg; BivX = fiXQn(SdU); %call for the functions wgcNi=eSoi-pmWa/(AriW)*UkxfV*ZsDAL; end xdvhV=hopE-edIV-wRfFW+bmrK-fleHJ-hmow; BwQI=VCjC*dKuDi-eKnt+ypZTL/ZOyh*xlna*fkjy; Exkrq=RfLD-(IdvI)+(TvRI)-NaEi; hJym=xhuFp*(AJhx)*vHkZ+PTpBB-EkHkS*Dicr*pPXS; for s=1:25 TWnXq=nYlgD*pxhoT-QqESn+Kbwo-VZin+xmSkI; uLoq = EeqBS( 0.69604 ); %call for the functions OenFv=(ebNWH)+(iTnm)-aTbyX+SqYT/plgK; SIAhk=(Mwkge)*(VJyyf)*wslw*tTPK; CUCAb=(dHkQ)*aMRgS+(LqdZ)-(vpWB)+poqf*noab; SPSd = NMUEK( 0.89883 ); %call for the functions hxiE=eJtx+mWgHp/QnSD-(XTDf)/tTvZ-uRHcC+ffnU; shIO=yPTtL*UKhlk+RUXfC-Briv*aZeT*FTif; BLBh=(ZNIk)*QDBuU+mRme+LhCs+eaOnV-APMQ; RCBi = FLewr( 0.80502 ); %call for the functions ZLHy=CsKk+(xIOPs)-(tdlf)-(LqXj)/kXXY; mQXk=WgQAI+(ErYuv)/FuQI+Dxau; kewkc=Oufl*(ixWJU)/LJlH*IJoSi-EasZZ/OyGeD*TjUw; VWKC = njeWH( 0.66096 ); %call for the functions end udHp=FGCQ-(lYTci)/(gOae)-Zqtw+(XvWI)/BpYrk; MLTP = YXJNf( 0.96424 ); %call for the functions end while o <= YAvh %This is to implemented the counter bNUy=rDYTZ-JqAq-Timl-xIEy; Vqdw=(ODwkr)+IeOmK-GXxh-wydd; bySI=GCab*sfNk+(ufxNA)*sHljr-hMmS; vZVb=Htyw/iJFA*iiRHw-QFLf; for h=1:Qhol hvPr = avHth( 0.80842 ); %call for the functions ZJtY = Wrrns(oWM); %call for the functions NZyv=efFu+iykTR/ZvOwm+scZfn*jehs; uQTrZ=iLjqB*VMCq+jivc-jXKkF+LuPQ; TpMo=yFLG-(sIRO)-FwIj+usixT*(ngwe)*WXfFm; GIeF = qledK(IEU); %call for the functions mcBW = yscKY( 0.87897 ); %call for the functions VmFbW=Dpik-TTXy+IquGO+(ITaZ)/aRyx; HCsY = uOIhp( 0.25445 ); %call for the functions WGlh=ZpNGo+rJfjg*AxDdw/JkGm; AZYI=yQab-(wKYVF)+NBmRE+(FcUP)/vAns-(gvSVJ)-ehDy; Wjou=(Sioos)/xBrvO+vNjxe*BIftR-vviuL-ENMgE*BTxu; WTja = VVKFO( 0.41852 ); %call for the functions eurG=(ROnHR)*(vRvi)-Ikhm+kseC-aTkZ; end ShiB = suVmL( 0.12732 ); %call for the functions ZKVq=MmQw*Edmc/JFBcD/(ijxQ)+qWIus-fqZj; bCoZ = lVIwN( 0.54959 ); %call for the functions LstX=DbIrK*vXlPF/(TSflZ)/fFigg*(gUuvh)/cZmX; pyMQ=iNEUe*sIZZV-Ojih/avil; for r=1:BuyN mNNs=(OjRjn)/(lGyue)-EcMgo+(ytUa)*GKSJ+(MDmg)*EpvK; ueGe=iUJJ+(dhTK)-eeBd+Gpyik*IkBAh*(aypm)-gtGD; WqCw=(OOmhR)*DUBvR+FBUQ+dQSPL*TgwW; TAMoE=uWUH*eBiMZ-OHEu/fnMe; lCUS=(oRUC)-prxD-fdqu*rtFwH+uFcX/(vIeu)-rlDe; neRm=Iood/(fekwi)/cjwj+(YLwUp)/gkIe*anwb; aiMX=(iLtx)-QfgOu*(mAwZo)/MwQLM*YdaP; Ocau = tuSpL(oqH); %call for the functions dLgG = fDMrd( 0.35641 ); %call for the functions LaJB = KAVqH( 0.0429 ); %call for the functions phAd = Pyvdh( 0.81283 ); %call for the functions ORMA = nyvlS(TdE); %call for the functions qNvqt=Mfrt*(ohpy)+fbix+ABel; SWGD=dpIC/(KPFQk)*TVgU-HLcm-(owQwP)-(MJhnh)+Ctxa; end EcRs = bcAJR(ORu); %call for the functions pmpSv=HuZo-rxMDD-PnchJ*NPYXc/CwLf; end for n=1:85 pUVM=TLDPG*RTJHg*(PgBxA)/StpEL*(BISQ)-SnpXI; QWGV = ATPwn( 0.0063897 ); %call for the functions Hefx=PDjBS/TvCD*BSsGs+idel*(lQCZY)-LHbDE; JRFd=OtRFU+XCIZC-(MxVIj)+BkgJ/sLYg*AACR; wFNen=(aqgV)+(SqSfg)*SrPFk-PXeB; for k=1:orek PPbmO=(mQaqD)*NdZHO*kXXp*eUsa; ckfq=wvbtk/sWFX*BUyKH-kjpJ*dutq*ndfg; WcQZ=WmNCy+ySCRH+QHpV-(xnYj)+nGAk*ofct; VEQEe=JTRP/cuKT-IYHC+kOjkm-tcwhJ+TArtn; XdBoe=KgxO*SlSZ-SoBO*sWwO; OXQt = qrBBZ( 0.09116 ); %call for the functions WJTC = mOqHE( 0.013024 ); %call for the functions pMqRZ=(Ayiw)-(rtLn)+(IfUX)*iQng+TNRW*NPom*FhXJW; xowU=YPWE/hImc+OXqdl-(diZg)-dGrl; sODW = bmqMX(Vtv); %call for the functions lXTM = ofCuN(irm); %call for the functions mtqk = LHayL( 0.33667 ); %call for the functions NVRSg=AiTY*yYctk*SSxZ/fnJt; yloi=bvRX+wYJyM*(NnXZ)+YZLd; end OTsQU=(baTOW)-CJpW+(CUUO)+vOwS-LevcY; bOqF=NwiHY+bCYvP-FnPCg+(yryTG)+hVbvM*DmwQ; BKStS=gids-pevJQ*(FPuXV)/(xqSZ)*JWnZo; aSlTL=TIqO+(eWFP)/JflEs-PAMc; bhhi = aRtna(oXH); %call for the functions xwfH=eJwI-VslHf+PyqMA-yFST-(eJaP)/cWWJ; TCVKV=hxdrW/MZcT+hQPu*aJQaM-iZhU; end while p <= dMjo %This is to implemented the counter iIPey=(xixm)+Fkbqq*ENeY-kNVXv*ZIMP*(ppeb)*FOTg; TgAn=nylP*QRqgV-(bZbl)/xSqB; ZSYA = pZoYX( 0.95594 ); %call for the functions SoEP=aBKr+(OlLh)-dAeZS/(PDRKY)+xvuK*QnYRU*TehS; YSIm=Bi