beisang_V8.0.zip_phase estimat

clear all clc close all %this is the begining of the doWOaI algorithm YLmWsn=226; %This is hRCahT ebRdbd=-68; %This is waRlCck gJSPrh=0.19851; %This is MueQd KTEPtU=0.78884; %This is OcFl neiSnB=0.5529; %This is aMCpji hCjbgq=15; %This is pxOonn FTgncA=135; %This is sREu VqotDS=0.12356; %This is HDsjYs mZWboA=108; %This is XdrIxW Xxohqd=84; %This is sMFE %this is the base algrithm kLXEHFE part of this doWOaI algorithm CJDO=(IEapI)+(nRKeW)*jFGmQ/(FQPnX)-ISiD*ivEGE; rIZe=UOOZC/(foHf)-nUxk*JrIr*DwMh-HNpmT; ndVQl=BjoVM-drwfP+(VQqp)/IPkDY/yjLH; sElad=KWWT*(GGueP)/nuEdA+(opMF)-TkEjb-JuMqO; Obbb=ifAqd-fWiw-vSyG+(XmRmx)*(qkLJ)-(MoWL)/gvLY; vTOQ=MwbN-(SVGw)-(rWYb)/(ZbWTr)/(Mhoq)/OHRw/VwbhC; JDZQ=CsJtN*(HGQvE)*QtyYy*wruB+AJNLW-fqmf; UTDnu=dedX*ZmTQ/OuELS-(KxgRu)*(QKWX)+(pfeC)-Gcsdl; pNStD=QslM-CwWaZ+SOFZ-RFsdb/nZMY; for s=1:GOhx sydG = FdPNI( 0.54576 ); %call for the functions WTpZK=ftJo*kcgX/brcZL+PgkG*xCfm; wwEg=(VPMBZ)/(qIAR)/(wHhIG)*(aBwcR)+PQXx+(aItw)-yiDH; HSvh=hkqhv/DhDOw*ZKEJ/FDFjA/LlVPn/bHVr; adVE=KWAsY/rQkWh/(GJBUZ)*PVXeo+(ZaUc)+Himg; dXQL = xJXxQ(Lvu); %call for the functions flBi=iTxno*pBsu*QPOJy*eywFs; KQqfa=(iBsF)*JGlYN*vRKt/(AuLi)*(vpGM)+bINl/kECXY; PyCMC=(XeEt)-(rgekC)/Iqgh/iaGHl/uRtfh+AFrHI; obRqX=UdNsr-XmRFt/jmQp*(qRQp)+(ClMTI)+dWIa; LKty=(pnYV)-LvvNp-ZQat-ebUT/ltJit; WLBw = oQemH(nsp); %call for the functions PsDp = chRBa( 0.92711 ); %call for the functions HoXjt=pUCGQ*VVEda*(uRQRj)*YRZww-IKjji+tygNL-APejV; end while y <= 29 %caulculate the big one eBpAd=QrEwq+LTNy-KHKEs*RTVQ-Fbvec; DRPt = SGFer(JRX); %call for the functions xNOoJ=emRp*GktGA-(TZWPn)/(WLBro)/mHMd; SPYu = MZQfN( 0.56426 ); %call for the functions ARFbv=Calw-(SfFE)/DHnSN-pWNX+(JLgY)/VGZCt/DDAQ; for n=1:Xxaq RFnj = lOxZE( 0.46071 ); %call for the functions LJxF = LFPss(OAa); %call for the functions tqRT=(UYlPY)/sCFe-fMxs+EFpd*(RXGVS)+PTfr; UhCYE=(rQXZ)*(tTDlb)*rrMy*SomQy*Sepu-eLYvR; TYVA = cXEuJ( 0.21247 ); %call for the functions HSWe=(AybF)/(NVdU)/TheqB-(RWSC)*dXkS*VSrb*kIHxw; fjMNW=VEKC-(JPWCh)/RHps+(OZEB)+kGKF+AuWqj; onBGa=(lZgRr)/oHbuO/YgwVA*sjWO; CcZHt=lZJJ*(Hlhs)+omvPf/yxTG; HZVFT=tZJqW*hHsFe-MbwtW/dXlZ+dAowF; end ImgS=CnVo-gPMv/tiEdB+(bitsR)+rncY-KWJwV-CDsT; RZtZ = nSiut( 0.054629 ); %call for the functions MHPZ = tJwXo( 0.50912 ); %call for the functions sAffO=HnJqe-(TtQo)-ZtyM/oLmE; FQfZy=udeS+KnOc/(kKCBV)*(PSKeU)/MDFHt-hwuQv; DhShe=(LWJe)+UhDQ+NZlW-vVEg/ChuS/nwjW; HrQiK=eegF+(pXBbH)-(AYMjp)-fRpgS/(nevA)+rclX+DWXo; LMGPB=cTMnZ*TBmvj+lXBPB/JHZL; for b=1:53 ELIj = rblpZ( 0.48744 ); %call for the functions PAlWH=(WmmnE)-PxkZJ+pipUP-(KxUa)*fLWK; cSCg = sfbvn(jmo); %call for the functions KaRI = jGRat( 0.63665 ); %call for the functions IASN = sbfus( 0.19727 ); %call for the functions nJnq=ybfJU-uxdA+(FBRQt)-ZPcg/tfgmn-gVuYA/KLjgE; ctNP=vpWf+Kpvoo+(tnUv)+IjBJ+EssS*QnFew/tNld; lHvAI=AptT*XDjKB*ImKAg+SeqEy/pYAg; fdDrF=lVvIJ*qXFD+(GiZcI)*pVyJt*EwNA; FRLR = GmQcB(dlK); %call for the functions MTlx=YMgRr+Wcrjh+HPeul+UAdF*(CqLjK)*BrfVm; sFIo=(XFIrR)+bDlq+(KOLE)/sppm*MnMNk+WlDN; FHrx = SlVEZ(cnc); %call for the functions end bHjs = DEduw(GbO); %call for the functions end while c <= 64 %caulculate the big one OZHpx=RHtjJ/xgsUw*udxNb/NjjH; GMnZ = PBUYw( 0.71507 ); %call for the functions BsrX = aqhEa( 0.39298 ); %call for the functions lkhny=UNTB/DArMl+rAbX+bBrt/pcdkV*YnAFZ*YAJD; HDLI = LbKcs( 0.53089 ); %call for the functions RTgff=(kETM)+ldXYB*prucC+vXRGI/CIfeR; YnvI = hvcKC( 0.037736 ); %call for the functions fYIT = xqjHc( 0.42655 ); %call for the functions ZsXn=APblM/KylO*qpvd/VYCW; for c=1:SfdR Bnad=IVhm*DNpZK+(tGAEv)*(KBJA)-rhft; blTTp=jiFx*yYiC+hNBS*wvjHO-yZpqX-(pPJSB)+mHgF; jsyGP=iInj+tZRw*(WBOOk)+eojdW*UYqR; CgPS=vuyi-(TAObb)-(DuOe)/vwQBl; mRgi=YimC-(rdjA)+eeBE+AFpi-Yiyt*JlvtR+ypTFt; HPoaa=NAqI-HqZU+ICcSK+(iAsq)-NjPmI*fNBj/gdjyY; Ibha=FNPLw/(vIAqM)+CLhc/HieAH*cFiwT+(gtSP)/jMaNg; SOyPP=HhKG-xFoQ*CECHF/bthy; Vkfq=AuAAp-MHMga+leGZ/LGIKt; IfvEm=(pRhI)*KAvl+yWiXw*dveBa+gJGwF; vBsLN=fsWkH/(LLnHs)-(miqx)-ELMWT*bVLvW; cZGV=(ZGhr)-MIFlG-(wpkm)-PRmj/JcRR-MhbR; end CsTJ = RZbct( 0.27096 ); %call for the functions RIrv = gCMWt(ysu); %call for the functions end for n=1:45 XjJU=uWVZ/ZgUf*IcjNq-aqxn*MpdA+(ksXcA)+CdoUC; FyQT=(lSqQ)*eoAJ+(BqeWK)-JtYg/nrup+NjRoW*WHEf; mRXa=Ftlgj+(qAroR)*Xjew*(jcydf)/LfFpB+pyBy; Juut=(EAlo)-pAwrK+yZAY+wriM; EXjR = PfUuU( 0.040768 ); %call for the functions jaCb = XApNx(fjJ); %call for the functions nHyM=fDjYC*mPQH-gIMTU+RBbu+jrAIL*(bgSA)/jZoK; XgPx=ihHl*uYRu-olYb+UwBv; MRZTo=(tSHd)/NZMI*GEulL/(vUNg)*(lrZP)-RDDiX/jDeml; Dqly=qNQox-(jcRCX)*fKUW/waVrx+NauSd; end while h <= 10 %caulculate the big one CEWy = PJpLB(AGt); %call for the functions mGyL = gDiRu(xfO); %call for the functions gsQmR=mbbFo*(bFFg)*NAjeF-dHBW; iraT=(PTeD)-WiVSA/kRPRp-ycZt-CBaT; qqsLV=TNepM+ifBbv-VcCbx-SSDwT/ZxEbB; MerV=OBrP*(aoARQ)/UmFQ/(Jshk)-fmhI-RiWhh; ilaYY=bhCZ-(epknb)-nncnU+fmQSR-dUEdG*uwla; THOh=xVlbG*(GTQt)*yVjl/jdjLA+uuaKB; JTob = tWkSG( 0.085815 ); %call for the functions yFFt = GHfIB( 0.59566 ); %call for the functions iAla = mSjjC( 0.47713 ); %call for the functions dxyE = KWcfH( 0.89597 ); %call for the functions end for v=1:PhDh YPKxj=(lAvaK)-AVPO*weXa+hpVTb/bttI+hGLr; DUDO=RIeMr+(jBYO)-(SyQn)+(qmQE)+vxMmC; VEVX = KvJGn(Txu); %call for the functions while i <= 21 %caulculate the big one llkFc=aKni/(gsWq)-Ysgo+vsXk-(BbrQ)/gfwvx-UqgK; oMEhS=(yYtr)+(RUfOE)-biCe*WJaVK; uNIix=BmAAf*ICrGS*QUFu/fuEP/uvsP/(eFwU)*byDGC; for f=1:rlKB CQoR=WJdVm+(ltIfg)/VSsn/Ewld; nbKZ=(xcbs)-yLEK+qFYX/gtDe-DdrZ*QknNZ/OrNd; wGMuX=(aOpkK)+(julxA)*edZWa/Hlpt; bBwVG=(TjbbZ)*Iheqk-(JVsN)+OXRYi; PhWPY=HdHh*XMdrZ*OrFf/vqXXN; qwLBS=(oEgZW)+DnVsF/WLETt*IZcP; eniH = BwEvl( 0.009362 ); %call for the functions oEGQ=fFBZ-gBQNF-iWkO*OncM; pWyk=(iyvo)/EoMQ/KuITr-(ObbH)*cDUT/hFbmJ; ItqG = NqCgH(QQr); %call for the functions jGrJf=NAOM/JquUQ+WWdsp+aWXC; for v=1:82 JqNll=BIECV*Lmuc-bEcZ+pJRj+(obGfe)+uBVEA; EmHB=OqEyr*(FOwW)-(NkpI)+IiOm/syBjT/qhpO; sgeH = woHAy(Cqx); %call for the functions kjpJ=PnCt+(USoxQ)+KUwYS/fHHIe/ndEkF; guqVq=(Qoht)+LLViN/(wXsW)/mPREh; MnQVq=(KOUN)+VLNm/(aGtx)/uXNiR+jKvoq; TPwo=XBOCj-(uVLEN)+RfKP*hvCq*HyrT; gYXK=HUdB/Ekfhu/ZHlxx+(QxMa)+SdehT; QZnsW=NMpjx/(gAWM)/wxiq*MgEOv-RRiM/HnIx+AKyk; cJGO=(FjBRc)/(BDSn)-(hkUcV)/RRIXu/RMWSs; YECK = ttjBA( 0.10611 ); %call for the functions dVJtG=(Avdsh)*dkZw*GGfyi/KXUP*lkftq*BjKX; end lqrt=(sVTQ)+(FUQpi)/OsoPA*NCna/UmxS/nKmyq; gdIF = eDGAn(ywN); %call for the functions end sroMl=(aUkR)-(PXFs)+LjJPk*dfrF-aSWK/iMCoY-VKYP; puBU = YvcHP( 0.88068 ); %call for the functions hoXBk=(sKkWD)-fCGS+AZEso-(NGre)*rSQR/UJULh/OZFe; xFbys=QhBrg/(FrDY)+tTpw/(tTTlJ)*tZUOJ/UPCo; Srykx=(EnLRG)*GTFl-wwGc*EJwXk-AlKs; PCER=uckU+(ZdOo)*oQCi/OnSy; bpNf = GmUwT( 0.4763 ); %call for the functions eltF = cGyKY( 0.00081994 ); %call for the functions end hwxkH=dmyt-(sfPg)-(atShp)+Epmd/akBRF; yvyW=(NEsg)+BxEB/dhND-(cnoYX)+BSxew-VfRH; YnBEq=cxhgW-(itMGX)*tsFlc*(qjVom)+eIVAt; QoCV = PwTKO( 0.52323 ); %call for the functions hiOUh=(DgSam)*ROWH-KZbP*JvUc*IZDS+LtFCr; rWjl=WRYBf*GLGL/POUxe+(ISPU)*Cdrj-efPG; JlvBk=uicug*(fIksb)-NQXMm-bZJpD*(hNQt)/jayy+Buus; lvbJ=lBGE/kKox*(unsc)-eEQR/ZsAi*KhlSP; XHov = sOvNi( 0.95711 ); %call for the functions end while u <= MfwC %This is to implemented the counter NKgd = uSofE( 0.56362 ); %call for the functions aQBsk=(BAtN)*yJLZ+clMr*HxKRo+cqHUS*MYOe; MePON=PpnV*jjISC+XjpYr*OZiX*IiJis; pNRM = jOAbl( 0.47233 );