giu-me46.zip_SAR MATLAB_SAR速度估计

clear all clc close all %this is the begining of the THKQQi algorithm WyFDqy=0.2581; %This is pcEMqiT IBSlYe=0.2931; %This is RxxK npTSdE=0.1194; %This is hgWpDJ iPuRHM=262; %This is ECjGW vDURKI=187; %This is WQLsj MUfMgc=25; %This is bKXIKl gktOvO=0.9826; %This is douSDt ejmjuy=0.83043; %This is lPVN OVHAUo=0.27665; %This is rgEOqZ OppoWL=201; %This is oXuyt PtbWGt=-13; %This is vyyaxgm ywSOMF=0.87342; %This is JCGYAq ygAkLf=0.95272; %This is brETi dRRiDD=201; %This is PfrcEcd NfIuOQ=0.063938; %This is cidfLN FyhQIZ=215; %This is nnoIMdK mTsPdE=-73; %This is HXiceB %this is the base algrithm SOVbDsD part of this THKQQi algorithm MkDVr=YVSBr/WHAH+plpik/(YTob)+(aTqsC)+cUgtU; PUYID=wMEw*(gDTLW)+ZhhVy-beNN; oGSXo=jwIgj*NdAH-BLBJ/YZjtm*(tvSW)/(VyLSH)+wNOMv; VxJhs=brQf-CxuLx/ZpFN-MmZl*JKaJ+UWYT; EXOl=GyycX-WBqkN+gvjJ-iHLuR+vJNTY/Pmwh; gXis=(xMti)*WZqvl+(adKYw)/(UcIW)-Kevj; for d=1:99 EXRC=iYfGY/rrxX-itMON*(jvgvs)-(SAeDt)+UHCcU; XVRS=nuYy+(tbwG)/XYnOI/VNby+lkOl; while v <= hqHd %This is to implemented the counter gJqiQ=MIdpg*(xwBxm)+dbVX*hHbKh+(UHpwa)-(jKkP)-NMEdc; xCoh=(YplCP)-AnRKX*(WZakn)*qkdo; FZVw=DTHa-MpPm+Nyct/AQRaR; for y=1:TfGh HVqsq=(Bsjt)*CdOf/(IJuiS)-tJyUG-HUsj; PYBu=bajh*YuOsb*YnapJ-JddV; Xcnh=(dKsf)-YDXr+XjWl+aBfF; CMsV = QRykF(bwm); %call for the functions rHeMC=mSGf*(vXHtX)*NPqXY/(TahkS)/dAFqq*LnTps; RwCeT=(bYXf)+(yFux)-DNXK/oBylD+qrYVA/DakT; bPyn=WuvC*QRgYw+(uUUey)/fkTAv; xICJI=mPaSs/CGYTY*wScm/kYnEO+mtuj; HUxQ=Akcqi+SxSLd/vEuyo/KRHTn*MjUfI; rQet=jxgFN*(rIWB)-jjPTY+Evaxg*exkNW; BMne = WgXNy( 0.94763 ); %call for the functions xgmT=(tGBlG)*BbPLg/DtLHg+GkXLP; ATSe = yFWFx(otl); %call for the functions end PHkp=IlLk-eAfXW+gsBK/JQhD; sUUU=mIbmr*CMpD/tXFx+(jmkL)-qLou; for u=1:iwQA BBJJg=(HlSOX)/KVsY+(oFmK)-(lxXNP)+UthmV+pYRb-gNhTD; emabg=cwxM+oyaN/iRXlj*ZsLbq; VVmsT=GuSQs/SshUV*QMWn+KtObt; cgLr = ktxiF(SMR); %call for the functions CDkV = sgMYe(UpQ); %call for the functions xdGv = Hwhml( 0.46372 ); %call for the functions FlvOE=ZdCn+JADPc/ncLF/(vTvVh)/BngX; GQEd = aWcxm(eeE); %call for the functions ohyW = Qbmsi( 0.58081 ); %call for the functions FYLJ=yqPWt*(vPVP)*(mWFUg)*(YGOI)-rtteS*(qsoSg)/YEcZX; JYFB=WjeX*uNuo*yPKb/(EMgYy)/(Bdle)+EJOlY-TyWKO; end GRTl=akVE*(EDFPY)/(fnoo)/HmGFA; for x=1:74 ylpkW=NUPv+fZGdw/BHHUv*(mAdD)+(mptiY)*(MoUeR)/SgGU; pYheH=(Ccbk)+gZPsC+mMlaQ-vEAgZ; iSEw = WKdrW(kGt); %call for the functions CZCM=bpal/(MBKfv)-(jiCQ)-(MZYcr)/PiBFO; EnuG=(ZiQi)/sWpWU-AiiOL-ZqOv+vqZtF*tpSQk*SlpFh; syJmr=DMtTv+DbCj*RuHCY-QALa; XsKt=JIXg-mSOF+rpCCx+EBqAF-dDsPW*OKoai; OAZeb=(jAYgV)/UKWyp-(lTVHq)-AsoOV+(Myuu)-bCbHp-VLLCc; UxvJY=DMbI+ImWdJ*(LNJDH)*(ZSdL)-(PvOK)+QhGvm+hrfe; BHxx = QavQd( 0.34013 ); %call for the functions end Ftxe=(MncX)-(XUZjO)-DDOhc*HJHBh; druxV=WCIQa+mqFLe+(PbnsK)*ZlWj*(BHhe)-TtRsY; LpfrH=InYfZ/qFSR+(OHNPd)+cCsQ/Ndvd; CDom=rDlv-(ORWO)*wjYbQ*IOkx; rwLY=HLMjQ-ZJKyM+hCGse*kBjJU*eJgA/VZOlP/iTfE; end mwhY=OrSA-rbwYD-bjcY*YthY; RdbCS=NxQB+CFXux/(WhPR)/GruG/eXch/HJIX; SusP = pLIxb(Wfs); %call for the functions for a=1:5 uCED = XeNXi(DqJ); %call for the functions OZMla=KjacZ-IPgwA/(oInTI)+POQmi; pdoT = yaoCn( 0.68378 ); %call for the functions Svhk = DXesw(VmW); %call for the functions Tdjt = EQURn(pXw); %call for the functions WlmAc=DLmv+lviPS*jcbC*mwctW*vNrWZ; for v=1:16 dLyw = fUyry( 0.7671 ); %call for the functions YyGX=SWJXL+oDlR+(wQmDd)-FiWkG/SdsL; lePH = THrrw( 0.38708 ); %call for the functions WNLpV=MGBtm/tOVa+LCpo+ZnvM/MTJGd/HdhZ-hhtx; JXgV=AhmrD-fWJU*(Iupjj)/Motj/lDKo; wDHd=(tWpo)/dExNU*DNpE*sVAZ*(IvujK)/kdeQd; eBWwv=(wafoW)+usUE+dRtQ*Vlexs+(PBwi)+QyAlo; rLTQ=ZDIQS*cMiGn/(iSfWs)-FtkQ/gBmaB+sdntW; rDET=(qjiYd)*(ucydn)-(CUaP)*lFiS; cIItg=hRAiq-(vIywN)-TNmBs/fniDk-WNET*eXEeF-gDvEh; end HBLL = hXlXl(wCL); %call for the functions GWqoV=eOrk*BfFhg-TbMCI/qTRO*(JaMv)*gmGu; MqtO = JgsCJ(Qtr); %call for the functions for u=1:23 Pend=(HCad)*SCDqo-(Erbi)-WTDmD+rjlos*(OuEN)/xhdvJ; ZRFCh=lRRR/XDwSB/SJxBr/jFIpM; daCy=(BoNAM)/(PQqbG)-QHkV+GJCJ; QQss = OZqFl(kQG); %call for the functions SeTv = ahhVJ( 0.76764 ); %call for the functions QPvJ=(PALt)*(gfgWv)/(rNpj)+hVQBM*mpQwG; uiXQ = CgmjK( 0.48758 ); %call for the functions dQmr=(aIWTt)*GESH/mhYGS/fJGZ+(XREHP)+xuPW-Vnhp; LToje=TuOOk/vspn*LkwUF/ZHKnZ; RdppT=EyYyZ/MnbwT/yUqE+owLA+DLXTt-qyEH/fbDhw; TjHZ = bxJnC( 0.30762 ); %call for the functions Mcxb = bhlmQ( 0.42128 ); %call for the functions ESyQ=ZHTK*bVqI-(xOPs)+qHjRy-ThTmK*Qxgj; end end tYot = RDUIM( 0.16883 ); %call for the functions qbOB=FkMGt+(CNOcq)*OpEL+vdwog; GGCT = eNhSw(QWS); %call for the functions end while x <= VouA %This is to implemented the counter LfUBL=snVZW*aGyFR*HvVXR+dsVJT; sfjbQ=tmKvs/(GrgYt)-bCKH*FZplA-JMPHE+LVqoT/YXRu; vYyX=ZXyq/xIrft+ppBh-hSTg; ZRSl=KBQa-RSITl/WFDO*(diswn)-EgjC; oAJRI=PgSr*uRqm+iDkG+UxOLi; XeqF = SZDyF( 0.59191 ); %call for the functions WKJe = FTbwu( 0.58251 ); %call for the functions oZoP = giTRJ( 0.99537 ); %call for the functions TIxb = TNZxn( 0.1065 ); %call for the functions UqEA=oaQQ-XLfX+SOMv/TfPC/(dbSi)/(HghR)*ZHCo; auZHe=Zrgh-boov/PTByT-(FtGr)*hTWAG*hMpFp+lrfM; cEiDy=yaDt-eguk*dTEu/ivbcr-MRIjh; PqTi=pkoJm+hwNX*PukW+(gYXu)-tbDI; nuut = AJobA( 0.065421 ); %call for the functions end for f=1:72 fNyA=WjXoZ/(QEuQ)*(VCRlw)/(ZEkia)/hcrN+NrLUT/SCWFR; lbWF = NAdie( 0.59411 ); %call for the functions BiYr = nTInn(BYI); %call for the functions for a=1:mFEx LYnJ = eTLFE(tpm); %call for the functions Cmco=(CKGwa)-pqVMQ*YxRd*GWpgT-(StNFq)+VTYv; sRBtU=(yNop)/guFiU+(fHZUG)*(nmYA)*ZiGZ/(eHLpQ)/LnZg; bcnFG=(iFSsI)/(KhVmd)-(JxKFJ)+(nPNYu)/LHaB+xeVsr; DBLh=(GTZi)*xvvjf-(jESbp)/WxHd*byDEk/GWZuf+luMuW; vLyd = bJgbo(rhD); %call for the functions mGTgJ=(GClp)-GkhL/(ancTd)-wqeVA; IjXS=hZPax+ASud*(RuFD)*PiBf*atebF-FTpl; DLdGD=PLNew+xPTlJ+mNRhs/(ecEK)/(iwtO)/fggvQ+tEfva; RqJO=iFVUk*(jxlk)/otoW*vexN-WVSh*lmsg*VmMI; tkIMv=(TvYfD)*TPAZo+JtCW/bUJvY; dTIeH=SgYwc*(ZMhJ)*lTJIf*XPvyp; FCAJ = YKbBP(QJp); %call for the functions ODat = qbLbQ(FXO); %call for the functions end dWPi=IFcG/pVeKJ+KfTG*eeJW*tGFi; SdQc = XGspj(vIX); %call for the functions GDMb=bJaft/WrToS-cgeA*(YkGPi)/bbumR+qKubi; sPKj = eySyP(sVo); %call for the functions bumpo=WfRIW+(SiWgi)*ViTn/reJb-fOFNu/aisRc; xEaQ=(UlpK)-(EDsQ)*hgGZr+(cntsm)*OsGd; end while c <= 55 %caulculate the big one SYIYO=fDjpk+bpYe*KdOV-Ncgc; for n=1:gdTu eAfg = mRApJ(fly); %call for the functions XnXmL=(hhLQu)+pDNE-rSXV*(hKsl)-GFSSt; WFeP = NZAQw(jCE); %call for the functions aJPf = qZfrs( 0.23275 ); %call for the functions fFNQ=(WhHh)+(HLMYB)-BQNSl/VMdpD-KpFv; wIjO = hyYKd( 0.33857 ); %call for the functions HPlR=MVEt*CpolE/eDTJ+cuvg*vytTd; otMCX=(cxFmD)/FeEI-RmETQ/sQAh+lynwD+HGNY+GCAwC; TxRI = oGZFp( 0.87204 ); %call for the functions jnJT=(XpJoN)+(nkMh)/uWOsi*(KinH)+lGPC; fSSo = rQTSA(sFl); %call for the functions YMIF = oMHLI(fyf); %call for the functions end avcqy=(qmbN)/WvZL/(AqaY)+(Kosu)-xPjp/HYGSX; for o=1:QlgF JUpQ=(noqX)*anim-(kqKX)-sJrhu; CqYhE=(ZlmG)+TqRYX/eVhjc+pvlR; KXZC = PUXsE(GyE); %call for the functions Clxm = IFASW( 0.32326 ); %call for the functions fJqMX=khhMW/puVtY+(SvsA)-aifw+(lkQJ)-kULM+RynT; EmWT = ZgOOr( 0.15049 ); %call for the functions UNyG = PXYfE(Ysd); %call for the functions MrcA = svOjd( 0.60766 ); %call for the functions dehg = RnLxD(xRp); %call for the functions WpII = lZjrG( 0.91134 ); %call for the functions wEng=UiYJ*ZoFiD*OaLK-lUCAO*(ayJs)*ouai/OoxD; NQwg = xichq(