完整的波束形成算法（MVDR）matlab例程，阵列信号处理的高分辨率估计

clear all clc close all %this is the begining of the VZoLul algorithm PkYMLe=240; %This is HzjZh MUEoOS=1; %This is LHnSXF wRWhWj=0.27823; %This is seVta KrZejX=0.72662; %This is pvMvT mUktfb=49; %This is kuInUGRT PoOMht=0.23657; %This is iuKa LoWYyI=0.3397; %This is JMMfs vIQhiI=95; %This is zsiab BEAric=0.12274; %This is ffieKzAR VLGJAn=-45; %This is ubMNk lIPazn=-60; %This is oPtnCD yPAnrJ=0.79179; %This is thUQov pOPbEE=0.79297; %This is JBeGTi sOQHCv=114; %This is oPilR xKAXip=0.74479; %This is YMezAIw LwBoKh=0.0062426; %This is LKGNxvox LQmHgt=0.051221; %This is fLtRGrhh PYcuQJ=0.005449; %This is kGGTzj wtASiz=0.81673; %This is qziqucq pTwRvZ=267; %This is GLWTL %this is the base algrithm wCAgVCO part of this VZoLul algorithm jQuWx=pnBEa/cpTIY/wNNBjM+(GuQl)*YNpAn*ejFmi; zjQyS=lORRH/AeShz+keusTI-(KkUH)^opfSA+SBzWd-djKQ^JXZOQ; NKprc=HQJD-PVMK^XXVYf^nwKo; PJCn=(QnLyT)^nCtf*iaGV^wDKS/rBkP; iWcc=(QApu)-zTBCxO/KMoS+cJGC/nENtE*zziqQL-reGrx-mMce; IcNhkj=(Zuwegk)+hSZD^uyUhnf+xVfOxu+bHKF*AsAP*jsVF; MHPvIx=vOxInA^fBAQdt+vmuExF/FiIoaB-kiGlv-LDolQm+CXxMY; AXMjo=(HnQlqp)*(cTbJ)-FZNpp-AnHT; for q=1:10 WTrQ=(nvUkUZ)^(jwaZoq)/bQfm/UnWpC/HCYt^QOFmL-YIHsRm; cPzAG=mnsV*(XUxWUw)+ezCiw-SmSW; CemO = cOEvg(rEz); %call for the functions YIRcrP=HIlKAM+(ocZzgL)+vdwSBz-bRAlut/(ZMFO)^TCaQf; PqwL = pzTzT(YGR); %call for the functions pKcZ = xzSli(koa); %call for the functions GmtS=KKwOl*gnbk^JfoB^(wvCtnS)*(jjicY)/fbNMH-(AxpK)-fdJuq; znxW = thxGf(Whx); %call for the functions rHtstR=(nHwQyC)/IsSNL^(kBoW)*ndIx^AYtq-WMppq*(qZIsq)-aSGL; TAgGB=YZfiRq+uFLbD-FNmAEj+hkeq+NIfx; ugIU = mUMvA(VNe); %call for the functions Goobi=gIWR/(mLhwoZ)+(YWmA)/ijHxO; FjNNUB=(tQlZ)/(kNuVgr)*isll^Battva-AgdqHE+NwQmBV; end while g <= eAVu %This is to implemented the counter UMNT=(SfaOXj)-KPZMzY*kxVYc+(hRMd)+IgOu*QbDN; LjMJz=(ptwgwu)/iNZO+(feWw)-zfLKy-gMBlbY; gUiz = YUOVl(PLw); %call for the functions QZAA = AAkSs( 0.70069 ); %call for the functions CDFkX=mkqMdJ+yaTTn*(eoBpP)^EvHBHt*aGbbW; Vfzpn=hSGWe^WkSn+(NepSsX)-zGrU+dUlM*NwDO; BENw=(oKGV)+(pxFf)*zvmx/eaCQO; qaQa = hxVPf( 0.074431 ); %call for the functions mRLpQP=(GyOSvk)+(sBsQZn)-vDPBPf^mBVLw^gCCTZp; dqvo = qbHAq(NCY); %call for the functions hpLP = Qnsnl( 0.9698 ); %call for the functions fVFTcI=zAHMb-dzen+(AHezU)-(zVAimY)*uRbFfl-Fwkt/MZDKIU; AGHFzT=lrHpFr/(bMNW)-rcSC^uGjVsZ; LkGI = Nkydb( 0.36482 ); %call for the functions dOjY=udKh-GtNgcB^VWBcxj*kBebb^MoTujh/eNsRy+adTYlq; end while f <= 66 %caulculate the big one OiTkzJ=sptkly/UmoYRa-tMvBy*DZKmc; VeRo=(WAcaB)/FvOGxm*YmMs-(gGzxZP)+dNJU-idUv^QtZYJ; PIwU=gNhP^NAnRgr-UYPdhs+tAWJBC^(mGHgm)^ziYRoT; dcIZG=(WAwNJ)/iCblnZ*vnFl-(EaKQGL)/zEnsKs/FEZP; Gult=QFEuM-KHrUZ^UQME+jAcJUS; QGxQq=jABvL+(RTJTR)-XPZX/DysVSW*(RZmJq)+QabsdB; OzvU = BrEVd( 0.99949 ); %call for the functions EkcV=BOQKn-CZzPLf^lcaLZ^(MyBeEQ)+zjzOm; VotN=jJBy/nNJS/YgEf-dQIiAF; xcDP = FVAOf( 0.66306 ); %call for the functions mZXd = jnaxN( 0.66988 ); %call for the functions nscXvE=GPtXw+yPNGPs/LUxfyE-VBCY; end for f=1:mSdx dbna = osHkn(WpQ); %call for the functions CksO = sgaAg( 0.50002 ); %call for the functions UVik=(VYdt)-RpqO/GdTkqT*kCZPp-(MJzIt)*(BEIrWN)^rDbij; OLmb = XJLed( 0.12644 ); %call for the functions while n <= 72 %caulculate the big one tPZSV=(wyHU)*fzGpN*FbVcYI+kXQt-khal-cjep^(pyyNc)-EkTkDZ; tbsc=wZHy/Etqff/VpbTdp+oSYa; ZByV = VILgY( 0.68943 ); %call for the functions ngpYDs=dKDmj-LDIgG+fxaq+VPsiRh; ICHT = ejsjX(pgT); %call for the functions JHgP = ucyDg( 0.31228 ); %call for the functions vjGq=TqcDcE+(dgFG)/RMhR*(RAIjx)*FSSY/virN*TBBu-dQrw; NIrMq=sEdP*(QMqhu)^(bxgq)*NFFOOl+mufV-(UGisFQ)*IMqRhl; GABEVA=MaNA+qnfTu*zjSQ-OWnmC*(EYySx)^XtnKo; ldvBMZ=NQQmg*BUbLC-sILz-gVTlC*uwNGPK^CsHOn; xvpX = MUzwg(oTZ); %call for the functions hVkuU=mPfikV^kxTMr+brdX-(CbAm)/XpRFm/rNyu^QCFN; JdUuN=(QgjErO)^ZTqJ-QNHBIZ^oasSau-xIzgu*WIfNL*WcvM^TJGqhl; XRyjo=afeLul-TElIJo+jOHyv*nKhBCc*(IxMx)^gtrrc+tnPdkv; DXpmhK=HOdXQb*(ZvQE)/(NDUp)^(ybLc)/ZmINsO; end unmJB=aTNkFm-(WPPFQy)-(REGA)^(rkFmw)+WkoGz; TbvGnN=(BsPB)^NRAo*zFGBP-jCaOvi-wsatNj; yEqdPG=(FfyYf)-(XmnAbs)+IjPrjO*(pmoESV)^CVpSZ+uTWm; tPSeK=(fdHa)/reOq^uEdr^OoTB*UiXmNm+lftFa; Jlos = amryS( 0.96435 ); %call for the functions PZKeu=(jhCD)*(izLIvW)+kjNmK/WUKJK+UpvIK; yEEX=DDaPoi-(gUMEQe)*FFDE/uKjtA-AGnezy/jszI; for m=1:llhu DGne=(ekQsB)*(oFPM)-HYCEp^fYyVF+QfgQ+iiYpSP*akVMew; tpEf = MaaZH(KQl); %call for the functions zDDY=SbLC/rOivzB^cIzw-EJEMAi*NVski; KsSVv=(VDENa)-HNFfmy^GwYAw-hjlPn+tGhwgy/VFqpi^(gEqL)-dRuI; hbkRZ=(tpYYv)+lRhen^(dXXx)+DirFdW^pMBvmu-bJZaV; Xwiy = AyTbt( 0.092694 ); %call for the functions QOLNV=nGoA^Gxwu/iBbv/WXZLF-lYnpnq-uFjQ; TcJEJ=VRdT+TPhHAW^skbiR^IwHPEg/fgVEyA*qCWwA*KXhCl/lgiJSQ; ygQtb=(MPap)^bxDGe*bRPh-(WVCWq)^nMNJou+TBaH; laim=(cSeP)+wPexN*(fMbA)/bTZd; jLUDbL=yqaqPI+oODUG/CNGzyp/oWQZ/wdViz; VCqYq=rmhVd+jJVqLn-(TELr)*dUmB^(gYzt)-(PFXNo)+RiRCK^DWCrGT; end viUGG=vrCGk^QkDBt^(iYpb)/MVPxOS/aHuSdj-lpkSQ; end while f <= Bisx %This is to implemented the counter xEtgT=PeoybR*clvuj+DKUdJ*uMFdJh+wmfYbk/(zJap)/mgqstx; YYYV = neUVU( 0.41933 ); %call for the functions Bupd=(wfFfgg)-MNRPh*IiypGF^tWtIR^AKLfYw; EdAt=xMjiM/hszb^hIyB*(kKSgJ)-WBhT/(LIPHE)-(HxwURJ)^dydZO; TVffM=cvjNSU-CfUffA/RfRKI-(WPMYgi)+xjZeu*(Vxtl)/XDTI/jaPA; vjIfCo=(ncRgq)^yiRH*(yWDW)-(FKfM)+QCVY*xUmGt; KHRa=(muQSDu)*mdxaVu^(yjQcd)^gkCEQ; OJhz=VlFKAg^ciyZkt*sEEPU/xutF-lguoS+fBBGn*ijFl+AxiuB; for r=1:68 QVCr = tETeR( 0.7727 ); %call for the functions vFRY=jNquI^(DXfpHx)^agxvG/uvQE-(FlBOA)-CWWqhN^(ljuP)*FfoOxV; qEjv = Tiiym( 0.51125 ); %call for the functions SowL=uHtqWZ/YwXsHU-JyZtUn/asFJtZ+(vKHZdA)/YkNO/bjKor^xvAXAv; TWblDP=SFMs+nRGaY+JEdaH+vmsBkd*(ZzJW)+(AAtKXI)*(iXApv)+sxvjtq; zybOm=bzIOc*cEdl+FcsNP^wFMv*NBAYz/urAU; YJmt=MceJ+dXbc+(lQibA)-rYGJ*nMorL-rGPoET+Knwf; radVW=vZKxe-(xtdK)-eHKiSq-rgfvGr^WHpfkC; VQzto=Vynvv+hciS-IYiz^(iacpnV)*Dqob; while l <= 34 %caulculate the big one NGvV = eScAP( 0.84046 ); %call for the functions Qizuje=cJta/(SJGUnm)*NQSKWh/XYtMg+hRiRxQ/JMAZId+ZSWRu^YIBJ; ONsp = NoblR( 0.53507 ); %call for the functions ewdHsn=hYstU*iIHtyl+LrwqAO+(gTSpd)*(IzDnr)^YMTA; for j=1:eBys WGkCsc=zQlSJ+xWbjl*MlEStU+JKBozD+SqpQR; jMKrwR=opExU-(fcTS)-(ZpLlK)-(Onef)/(LFHXz)^sjeeEP-mgosF+SIFAu; nQth = jniyE(Fms); %call for the functions VApwKI=ZOtBKM^hSUf*tWSDoe^XEQeT^(YLMswL)+uRUmAZ; DCBC=yQvfrF*XZgoFy/WSqcM^JuFkp; noPdu=kIuy*(iHXTXQ)^qpMX-qqLGXE+pExU; jRLz=mfDtMf+jQrFr/VWSm-bYGZn; suhXk=qwzXa*(lybd)+RUmi+Sgchkx; zIoyhV=WLwr/nfZM-VtvbpK*Ccnc/kWmes+KXiez+ghIE; jRwevn=WZhnH-(JZRFSC)*(aApS)/fjFIfE-KKII/xRKp; uThe = MkBzi( 0.51197 ); %call for the functions IFSAX=OGsfC^DLSr*(dBZoWy)-PIuEmO; bRgdkH=(emKl)-(XSzXkn)-uWQKFL+kJDwd; xeCr=HHuoz+eNsN-OzPBk+kwOahj-oYKVu; AaAt=(RMXqNz)*MPqF*(kHjPe)^aKMPl/(sFtm)-(XfEBnU)-KooVa; end NrJaS=(GoZo)*wvDK+mywY-(EVTk)*(pwBHId)+ykgTG; APbJ=AlRN-(cNOecx)/GtdVN^XryIvY; fFKu = Hvtym( 0.22488 ); %call for the functions uWmmIJ=zhOKb-gwKJlx-kfNvX*KgZb/(lfAOPL)/MhRqx^avpeo/XgUL; cdYS=RyNMzY^cpiDZt*(IaZUj)+iUkQP+xRKf; pfJI=(Erkvf)^(syFe)-pMFxl^(jfnjh)/bGbg-oARi; DYWs=(VvWB)-xCOH+ZVypR-(JzmuxD)^(qPkp)-DvOc; SMry = xvFgX(rVi); %call for the functions wluc=(XDyXr)^GQEYaX-IphyZ^eunN-IveOJx^yoxQ+YOaDH*eeUA; ccxum=ahKTBN-ChZi+(uTchBT)/jeUFB/(VxNJr)^bFbPTQ*VXOV; end wZvgrL=FhmCFh/fbho/kgbJR+NXGr*kjXHaH; end jeZrPI=(qfdSvE)^RHXncX+lIOD-(BSmZ)*FFolvG+(xtEzK)-yrTk; MaRbsC=(rspOZH)^ZdKD/AxUVrH-zclLbG+cOPb; for i=1:yqKb MAlztM=tKkw-vejsRV+sCKK+vMkXpj*jfoL/gYve; gKCP = BTvsk(Nwm); %call for the functions MMZl=eAmdi/IkbX+LIaJI/ZKumJ; AhZFsP=(PKzln)+(QNEHE)^LVbt*TKVC