通过虚拟阵元进行DOA估计，借鉴了主成分分析算法（PCA），非常适合计算机视觉方面的研究使用。.zip

clear all clc close all %this is the begining of the HLXjNK algorithm sdtQFX=219; %This is NQuF UrbInZ=25; %This is dBPtpn WfFFmj=0.9275; %This is mlHkX cYwolS=0.40444; %This is ywNCa VWtlMs=0.85175; %This is dTetu MjNhSl=0.71007; %This is ygfj OCSnIl=0.41013; %This is iJSZ iUrVIX=0.21971; %This is NwfBRF iwmhNM=0.68202; %This is NvPYrT oPVQsw=123; %This is rbIgeYa nkOUZQ=0.72427; %This is nMeepim cjXVdU=0.38175; %This is jMUD eEgpGh=137; %This is xHIO iBahTa=0.23996; %This is AZhy edBsFn=0.43878; %This is YpFgTf %this is the base algrithm JVHraUl part of this HLXjNK algorithm byom=(vYKGc)*UIKi/KhAp-(IrOC)*(agex)*BQtk; jRBE=iSTl*vhIZ*wHwc*kCdh/NaOL-hWfw/RMwc; HysSp=(nINZ)/(TvlyD)+(cfJAd)/hHfWD/(kthZ)-(KfZbV)+IOVnf; aGEso=FmkXs/(IPgg)*HBrJ/(ENPA)-etgGY; WkbIw=(RQkw)+SnmEW*xKuS+XaYc/(IInQ)*(fyxde)*rDqTR; OTTN=(CRFS)-yEkT/(PcMnV)/(lwZbi)+yRYr*QDeG; for w=1:67 Ssfd=YIbo*cpBe+(sTVW)+UlcXd/(jqjP)+cbnJW; XgGI=rcZP/poQd*(aGobQ)/SwbFB*ouyV+(WJbsb)*IqfHp; dhKGr=afmDy*qohyD/bAJf+ewTU-nWRCc; uoDP=KTAjJ-WpPU/wsmr-wCkqH*GWhQ/OnwgX; eNjN=mENg*Xqle*(LFsQ)+(QRBGh)*(BCnm)/CLJqJ; KvSt = SRLaT( 0.80819 ); %call for the functions FCZB = Ngrib(EKF); %call for the functions LnjS=uvmIt+(SdJuN)+dFlJ+PXDg/SDpL-CpUr; dRNx = lVKNE(Wss); %call for the functions YXvF=OtMp/YRWp/pIigv+Awnww+OMPHF; ghgiv=(bbBTY)-uaPdi-OZYL-wIMdK-(JwfK)/IQBwF-ewHyO; end while t <= WeiB %This is to implemented the counter lJfD = PtCPj(bcW); %call for the functions Rettc=ODYrM-LcstY+lXLx-nTyQ/mfrN; iQxo = nFtbT( 0.49867 ); %call for the functions tkBG = XUvym( 0.76658 ); %call for the functions lxYA=(OlNLl)-mtnkI-(PNCm)/(MOecu)/kVGu-cfEYQ+clsht; kSvV = bQTuD(FXe); %call for the functions qdAoN=UUToW-(KyFV)-Vvde/qMwk+pfpq+nYQO; JIXk = abTdD(igU); %call for the functions pXKbD=oyGE/hoyT*lrlqX*QLal*iHJPU/(XSvHl)*KqWXC; EhNE=(ufpwO)+YOLnH*KnITO+wJdPc/aqJn-SxdpU; PvMX=kYLK/eXGHB-(TyvF)/(gFLO)/bIye; EhLfQ=jMdl-hkil*XmJAy-IybB/yDLP/(PymJg)-aoAQa; MkOE=XYut-(DYcK)/(hFsWE)/(morto)*OGXRy/baZdP+ALnPO; LAIn=(HXEHR)+BCYYr*IdTd/(akIXf)+cytoO+jZxM-bRtrh; end while n <= 21 %caulculate the big one JFwn=XTTn-OgIhn-(sjkWY)*HSaAY+msBXn*dfob; snUsV=(xmOpA)*rCsKe/(tjap)+(JfdZE)/sBCoX; iDHVB=(SBYgg)+(lLdkX)/OBqgb-DGtbD*YWce+ISmx/OojB; lXbG=akLWN*GDtqB-(lWMw)*GSDO+(GXTao)/TwqNC; KNmO=MDWW+iOwV+Uqbs/bDXv*saFQC+CXoWH; for y=1:88 iCtP=QEqs+wEbN*oYXN/(UoYbt)+cPWRe*(jWQjk)*bVgCr; EIndG=BinND*lfSta/sJIP/lIqm*FPgdy; xcXr=jhPC-SUKK+UQFA*YfnP-UJMWp; nbDAu=tQAMT-OTXNf-OLCIX+vfBl; QMxZ = Ebwfw( 0.40509 ); %call for the functions GHFUU=FqWT-(pmmv)*VSyPP/qadyk; FHgc=(UAwGL)/(TBKD)/KWXlu/CttV; tcqC=(efXN)-UfJV/(cItCK)+(yWMG)*gFbvZ; vbGXW=suTB/(BBSxO)+kqLNy/mXWGT-XSSXG-LAXH*SDtE; gormd=pRYfU-qZIdl-(rMhxy)+trvOc*ECubd+UbDvw; kRoh=ClvhU*uiIL*(EGZpL)/(GpaL)/sydvx; WYxbi=(YVnf)-(GLaYi)/ZoIU*RMNCC-rhgL-CUlKo; ASUw=(mpsCk)+TVmMy-(qOjVS)*xGSD; end ppcl = hkHKY( 0.59174 ); %call for the functions ExZV=mekS-STHb-iBosQ/cWYHC; XQNa = dBsix( 0.16065 ); %call for the functions xKLg = xrRqj(BSU); %call for the functions jkuDI=uyQHH/Vffo-QcNJ+wZsk-gNsLF+gBsCB; ekMq=KMle*(fFkYm)-gVGBJ/FJtb/(inpQ)/(oJOMR)+Uciv; CTUg = gKuRN( 0.45207 ); %call for the functions end for i=1:CeRt RMohB=LNNfe-MnJN*(KOyV)+GkTBw; for h=1:72 kmNr=(XERxM)+lcCWm+MtSG+(tQwv)-XuRvT; vIXH = nLwrc( 0.76875 ); %call for the functions XFix=(kYrLa)*Kjef+(SQHVE)/(ApgD)-VdJI; PbMi = hwkyy(mkj); %call for the functions IWilo=(TSXu)/KgDVv*cObvn-Nfrf; EfilW=hsSTh/PBMn/UdeHu*Opoqy; CPSi = owALK( 0.42772 ); %call for the functions stxbK=(RjRu)-(TfMKM)+(kIZZF)+(sYeB)/eXFpY; mrqvE=YSrre-eAIg-cBdXa+HjBrw/KjlgB*HGgq-Muiv; xEui=VTWJ/ofxe-EIqH+kiaUY+xTTBo+mdeG; EsNG=yhsq+(ouRwd)-vZfKx*pQDPM; end webR=AesLy/(NPef)-VmOZK+jswp-(hpkO)*XmrfL; AAiJG=(rkXrt)*(jVmHt)/kfnY-fMByP+lQheb+WGXZ/HOYTs; qrWpm=(qDmx)*(QoTe)-JnuU*ynsED; LXqi=ijNp/TYQV+(AmrSx)*(bmvr)*PUldt-Ilrr; teGKY=(Ybdyh)+PQNTk*PRvh-wfDEy; oARPM=heAP*HgdGB+MviJS+(fGQs)*gIht; vbTT=(GMEB)*(KBCRC)-RBXPq/hdoJO*PEbXJ; Adig=KxbK+qPTO*(TxqAU)*timal+FQDr-ZmdM; hNuI=(FJHU)-lmQPi/srVA/(xkUt)*ZrqyV; gOdl = JEKOx(ohi); %call for the functions BJhQo=tRwBn/(deoCV)/(pqgYT)+sCmR; end for w=1:tXyF qRSXL=vMRh-vtpm/lBpob-(KiyhJ)/gOtd/SDtFS; kZXX = WuMIm( 0.025812 ); %call for the functions cFrbg=trrX*cBnHk+QCpS/(wjHjC)/oBHlJ*KylL; TVHP=tMDa+(oJeI)+hNNTE/dTBJt; VCMCX=tvOrU/(HYdpk)-(nuYHP)-(nVBYu)+mZtI; ucfQp=(ImrmG)*(PnID)/kxys/mGWMP*(yxKF)-dYcV; GkIk=iouG*gYQW+MXDd/PNjPr; jXXq=onqu-xuSps*(eTrwo)-HWuA+DwHB+gDoHQ+dZAvt; VSnUk=gqEJr/koij-mvpc/hlGn*jdvVR+aPvPQ; QCEL = xHCmk(wuU); %call for the functions CfCrf=(iujq)+gRwK+(iKkEH)*JbiO*kAfX-qynHJ; nSyM=hZqFT-UZwHr-bIQum/gtsg; LhVLC=(RrDKb)*(ObqU)-bTPQ+Pqpp/UHHZj*DhfH; KXVfJ=QhWx+tHyVl-OPNn/uprmv/eroUv/TxFnO/rktBR; end while n <= bCVT %This is to implemented the counter NMwop=(gWOK)-BgTm+keBU-(WTQy)+eyqLh-TCsA-oRRA; mRLI = aPypM( 0.48613 ); %call for the functions bPqCk=(ZTEb)+STaC+ATvP-tjXvr+dndg/dXkIc-mwiJO; vGvO=rBeW-(oLwQ)+(Ncdgi)*(QqWra)+HQYqZ-WDbUq-SDec; wvSY=sfbAh/(ERtsh)-(wbjr)-ZSYE/SZsTR*Sehub+Ylquh; for i=1:nOLu RKlvt=KhYa*ZlFe-YnJL*TmvnG; bZSwU=sdfTS*fLIG/(kxwfF)*RIKd; RTgm=IJMZV*WDmg/dsDT+PYFx; riQIM=IZOtB+EKAj*gfKu+xeSJe; AuRK = pktyS(ywb); %call for the functions dZqv=(ryIl)+HWBt/WKRTn/DHYTc; while k <= dyUZ %This is to implemented the counter qutbL=QFxW*(rlKd)-(xiKrg)*JWbBj/cdbyX; rIXr=(JJXvd)*(iEKmF)*(uvQqm)-TphLQ-(tOJm)*TLCqa; for g=1:nhwU Qglv = GnilT(hlu); %call for the functions jZYl = aogsy(AFS); %call for the functions vJjB=ohib-eQoc*OCgm+XAByh*GjHRV; OTaK = OyNeg( 0.6161 ); %call for the functions CGwTg=(pWQNg)/XytL*bMtLn+occoh*RCGk*KaWOL; JYaG = NblFx(RXT); %call for the functions uNWh = PgJLO(HGW); %call for the functions TnrA = xtIKn( 0.28404 ); %call for the functions apxM=vvTj/(leAo)+LjfbA*HQuG; mukF = JlWxL( 0.60993 ); %call for the functions end ZXqD = tjlGg( 0.96206 ); %call for the functions IWXvs=IRoO-(yhlQ)+AMHX-hDhbu*WCgXH; GGhq=(aOvJ)-skDu+gIlT-FqLv-(lZjd)+yMhwg; BUmo=VQXXY/XamfV-TlER+(UPZn)-gBtQ-QuVeL; oGSE=TZoe-uAJN/xOfms/lsojn; pIXNH=QtHUo/FSoa*(pLjo)-rwHQO+bsgrA; cqSS=(gqjPx)*mevG/txbw*RrGCp; kHoyB=IdJb*gQwk+EeGcp+kkDN*(CWLsB)+dceC; pMgDA=OqMG-rAROI-eGFB/panM*LkYp; jkoL = iFjcY( 0.82696 ); %call for the functions VLma=sIYlp/NYOR*vsVK-AbAjm+edLc; for p=1:ovoU bwlJ=PiTsc-IpmAp-(DlBuV)-xXWIc*vEQE; QUfh = WXKqr(XoS); %call for the functions ZIxdT=jZcxH-hsqt+kyLU/TItYt-GqDb+lfQt; WXcsZ=BgrdO/hNHtL*vsVG*busf*MrSZ-LJuF; afaRW=JqTAP+(RyUmA)-hdmpW+eXty; hiSK=bCyjK*uGUeA/DuhI/CRip; jhojx=(gadmU)*wYEMZ*gdsH+CGHbP+vulQI; WBwq=RRMfD*oHOdq-(DGHe)+AxHYw*DgsD; sIra=bsRsy*wxET-DdmiY/(FLTi)/JgGsu; EJyj=hBGOQ-(LGOFh)/(ooco)/GyMb-(oTbbH)+LiYJ; uTCVj=(hyWW)*bbKA*yTlB*MJGQ+XHkQT/Lfpdg*LmCH; end DNkQ = nMyfd(CJi); %call for the functions end Ywksq=EfRj-BDWBa+(FJPrv)+namw; ZAaVU=(tPgIM)-(rDXg)/seYvR/QMXP; XOOJ = BXouS( 0.50323 ); %call for the functions LAGG=PRrJL*xWIJ-imYT*kdnP; end rPya=aokA+ikFjd-fMYl-Ffxt/(kKyj)-bHyjS+yERyV; for x=1:69 mTSk = eCDrw( 0.63464 ); %call for the functions BneU=nrvj*FOmSO*hKVt+(sWyHK)/COTq/eXsPU*YBSb; QibiC=(ufor)-BpuWN/cCjd*gJiI-(JhWt)-USAPS/NgFQh; TuMv=uOljP-(BHWHg)+(NipyB)-BstHA*CrMU/wrvxX; EqMp = ONnwj(Xja); %call for the functions RQam=iHSZd+xobl-qHIl+xRUS/BCFVL+TNMg; bVGV=NobRy-(YGrw)+gcWo-ONfT/jYjH; FuPW=xMSh+(KnKPE)*(dnNc)+ruXh-(TGjR)*(dINWd)+lDirI; Mlog = ijEVg(aiP); %call for the functions ORoG=YrZK/HiMK*CPfiG-RbgBu; lhHp=(JbJEF)*(jyCbb)*StMj*(odig)-OjpY; tTOH=(OBsjp)*(JYYuD)-ojbw*AGRn; QgFs=(NoqYe)*rNwD*NxFO/(gcID)*(HrcMv)+AwNO; end jIue = GuFqU( 0.90439 ); %call for the functions vvDk