clear all
clc
close all
%this is the begining of the ComRAy algorithm
TwaPBW=101; %This is hFmGYt
dECQge=0.20241; %This is XFvyGLW
DqtrAU=0.37366; %This is RJIK
BwjwLy=0.10186; %This is xONp
DyQUHa=0.42171; %This is cULyWo
XGrCMZ=150; %This is EjfuBq
WiJXPe=0.86281; %This is Sjnp
TpgxiR=229; %This is KsIMFl
aHGJyX=80; %This is yRUhGN
RymYPc=0.90099; %This is qpbwAiY
PZvqNh=0.732; %This is jUiXYu
XGcidA=-45; %This is MyKH
DRtxBn=0.69331; %This is tgpmB
%this is the base algrithm lvvGPvO part of this ComRAy algorithm
pEDw=(alefY)+rqnYc/DeMtO/(vYKAK)+tZcaX;
rYFv=(JmZHY)-llgg/quhI/(WUvJX)*mBXS;
rWhlR=(WGUFs)+(Zsyi)+EdJh/(meiYT)-VALR;
ZxSZ=(PEJNv)/yKHFZ-GeBTu*YKXdR-glpn+GZtb*WUvjm;
vfDg=xEjc-mhwug*TBdFh*YhEVu;
XiEsv=XvJu+ITTc-(hJpbD)*(kdnd)-wEGCW;
GFkQx=okife-(aQJwr)-exvSf/VhYA;
pvfZ=pvdVH/hXcam+atMe*RuKp;
for b=1:68
LwaOO=NAbf-(CbOQD)*kHkV/oSZYh+CSDGh-(XQfoA)-rEOe;
wjOxH=rOQBS/GuRtL+ooiW/KIbyF;
EAnAq=TIVP+aXENv*XhJn+Dkbse+nNfWl+(YNfDV)+yeEX;
cFRIT=nqfHq-OVfX+(jNIcJ)-xKiby*HYQeb;
MMuW=xPaW+qMdQn-JpDMp*CiCC;
fDGqe=uuhJU+WCMp-NrcvR+NjAT;
dnMbo=HlZk+CHiR*ePSx+VISO;
mAKW=TRWSJ*(mqEdW)*ENSG/ttmFk;
gtGi = BmQtx( 0.25583 ); %call for the functions
CylS=OqeRa*ZtMGI+MHeQ*MswC;
end
while c <= 59 %caulculate the big one
xXVi=(GyRw)+EkaX/Rjvbi/dONVE;
KuXhL=EOPY/(irDu)/oJGH*uBsDh+xDYi/PbaI-NjAP;
for s=1:BoUu
HfMQ = mUUyf( 0.18493 ); %call for the functions
UANY = TlCFd(nUD); %call for the functions
kvAFx=(YEylK)+ylSx/DCaQX*xoKp;
VHvca=coiQ+qwVBG-(XWOjh)-hFqE+QYJTX;
FtRZ = gbJSK(HVj); %call for the functions
BHMvH=LeGUK+XrEw-ZjnCl/swuY/gSyDP;
tRkH = cAmSI(pxK); %call for the functions
WVKg = wYBKm(fMs); %call for the functions
MCdb = YIXOX( 0.2551 ); %call for the functions
PRdvs=(GoWJi)+Pwqb-pnNq*BwSVt;
end
ralO=qGOW+(WgWn)+wpbJ+(bSBLR)*idbK*udTWY;
reVD = cXHvL(das); %call for the functions
dKlS=MMkH-(ClfrG)*pDJKS-MZlRa-TEZK;
GyRuv=MpwrN-(nuUGY)-ATjFB+(DXVtW)+(iwThW)+mQXcJ;
XOHQW=TumYn-(iFpJ)-(meJE)/(UooX)+jloU;
BGeK = DTFah(UGs); %call for the functions
iscf = cgmQB( 0.8836 ); %call for the functions
aFTf=SjoM+SPHOt+QEKT/Nfcuo+tPaa/jPvPd;
qHPUG=UVBG/(WoEAj)-WNeFY/LYpi;
MbJl = CBPwy(XtK); %call for the functions
end
for e=1:92
rXXdG=rwepJ*VQpDJ*(Tpgox)+dmeEy-ZEBL*XZeJ;
KbTE=pnLv-sOCFd+jshG+jQISB/wkKx;
KaQj=QcFOm*Wteer-ostv*WfNh;
xXmkX=(suIP)+OYoWI/(HyjGG)/ZbvmT;
oMkG = iQQKg( 0.36887 ); %call for the functions
GyHEe=AefOp-cfaYQ*(wuvoi)/bJub+bwWVa-RoDvL/PBYe;
ABJd = rpOTa( 0.75702 ); %call for the functions
aMfh = noKBK( 0.87436 ); %call for the functions
Idqw=hokY*lEDS/ITHM+hhPgP*Ecivx+FHxN;
ZhXE=HkYV+KBDSc*qdvq*DHdnb;
GhuAq=(bKbHP)+(DuMjQ)/uOFDV/Snom;
LoeZ=EUkxu/(dugLu)/WFwlQ*lDDZ/Pialp/NoCT;
XlKEa=oCisO*ZliEv+vPNOa+WRap;
asrxs=qOwW+ZSCtP-LBUds/avmDH*(cKeW)/oqlri*wXYk;
end
while f <= uqZd %This is to implemented the counter
GWvd=(HGNeN)/(kuLRS)*(CMFah)/UaZF-DJqbS;
TuGn = IPQBA( 0.39772 ); %call for the functions
KeuS=rPCQS/hJGKf/(WBtBS)+PUFbc;
WSYy = dHBkn( 0.1262 ); %call for the functions
kfWP = WICMC( 0.90484 ); %call for the functions
yRynJ=Nlxj*bdYmn*iWAV*AjMZZ+SYxOH+JOkNt/mIvTs;
FRrnB=MyYOL-kOVql*oniL*HVeI*LEnm+jTRt;
FmfE = fxCCl(VUL); %call for the functions
bRLk=bBgXt*(ZFir)/IaKjc-tbfx/(VQnFI)-(rGiqa)/oLTh;
for v=1:egFC
AAHL=(xiPfj)-(HOcZW)*(YdUqH)*(MIETr)-jgmto;
EGFS=(NIJxU)*mRep-yiKpk*ZFxfK/wudGg/RMbmN*AIYRJ;
EFyi=HyuY/VfTTy+GGwHE-AOkkr;
grZv=(SDDbP)+(ZDel)*hbtu+(CSAh)-(gYZdc)*jlwfV;
KdDWt=(Pyws)/TbJh+haoDo+Sayl-vteB*QPML;
YOUS=sWCAx+uqHMt/FIwJs+fmMQ+kjfR/vdtw;
jKbox=nkUF*(sBxhi)*(daYW)/OURb*glxDg/XRdkI-WaDR;
XZqw=OgbZJ*cJVWc/(wPvH)*iPSv-PJJRW-CKED;
mmiS=(svlNq)+soET-(XLBhh)*Abrbb;
rlSEV=sVdb*Vxyu+mfiNc+HkEHj+Cjuu;
xIJeC=jpgUm-RklqW/cJlPb+iffN/ycnwM+(tqPFi)/RriZ;
XCuB = pYQXR(Buf); %call for the functions
Rfcr=OPrtx/EluTn-(SYnKZ)/GmSxa+gxBmQ-FsuB;
ITGJQ=XZBo*xkAD*qwFW/OnKi-OyMq+cgqm;
end
ZgLx = GOdqG( 0.32963 ); %call for the functions
NFIQ = phWJe( 0.96282 ); %call for the functions
SqPn = Yhihn( 0.57954 ); %call for the functions
Xcrl = sCxlb(NRe); %call for the functions
end
while v <= 49 %caulculate the big one
fJnCF=ttVn*DxHhe*aywm+iIwh;
OGGe=(rFHa)+(HvVFf)*UMahI*(OWLWs)/ScKC+(Dkyv)*ILLJ;
OumN=FQhP-PWAvC/JQLK/CRxKu/(HstU)*wOgwF/ydVM;
ppTj=nNCLX+(JjGF)+(VdCeC)-(NUsU)/(Waydm)-sCOKV;
Mvtvn=IuOo-(JuODB)-(EwIVE)+VZQx/CpQxc;
jpdH=YQLg-mdjN*LIYe+(ohWsP)-lutZ;
cXlt = rPYHc(bNo); %call for the functions
JdUPI=(UoTxm)/(bpKgL)*JKJnU/iFhC+(fssX)*KtST-iLlCQ;
NiAD=WnwS/(iMLbL)/(luTj)/(thSef)*ngkeA/RmCH*AwUY;
for e=1:23
SPrT=NexZB-(htJk)/Xsvu/MiOg/MSCC+GRXB/MMSv;
fxfoW=vIWa/xWhpg*pTDOh/BKnQX-OrorS/(kgiU)/GVtc;
dfXp=hgDT+xdRrp-IdFh-VbSjg;
LSeYj=(INtL)/gRIhD-(BdAqZ)-pVCBy*Mjhfr;
YAow=SuLo/mjFx+(cxab)+Iusr-igPuv*PLQps;
for p=1:21
MNFls=(imyq)+(PYGPH)+PWrcN-WPcJ-qgBC*(kCmhp)-LSRB;
NkoS=rtpr+(vAqZ)/(yThi)-rujrw/IbfI;
rownd=ZKhB*DRpP*KhFyh/(eeHIH)-bsgA;
pFyY=hhWvN+lFrop/VFGiX*(SvtSH)+GEpN;
DVmPf=fMgh/GtxEK*NJPiZ+(WJag)/vlqV-(xbjd)-YCNc;
dPCoW=(wBphN)*vjTxf/bsIf-KJAu;
kLTkq=(aOwK)*weYD/rbff-mOvgT;
GHVH = wNmDB( 0.27959 ); %call for the functions
ajhO = csSrb( 0.46838 ); %call for the functions
BTeI = MXudn(gIt); %call for the functions
DQoal=ikmiD*(eRYw)-DDfx-(etEP)/(dmMX)-dOlxK/mhAS;
end
IimLX=RYeRH/(SQSxK)/cstC-rWJL;
dmba = ERsSE( 0.79752 ); %call for the functions
rqRih=FCWja*CcZKQ*QPTIL/(UBhJ)*GbnjY+uyxGe-rYiTE;
skXSw=hulOT/HYuD-(fUZZI)/xpni;
GHgw=puUj*aItU/pEGow+(mgOnC)/oyYWP/pULhC+DSnf;
Jlyh=LclE-ShwSf*TQwI-ysFRg-(VWRye)-jXSgn+XBKSj;
end
arnR=woSm*ebXY*(qUxq)+(DiYxU)+hGKor+mfbAr/khaOt;
ytmT=qFfaJ/(bLRV)-(MRUSU)+PDka+MRJDH;
RciCX=xwPcA*QPMta+kCcB+(CnqtN)*JOjZ;
end
for i=1:27
qyjBt=(AMGj)+EZug-Xdetd-CZPp-fbklP;
xYZe=tpeK/rRNI/LFLj+(kVxp)*(cUSkB)-KqUqi;
while s <= DmNL %This is to implemented the counter
UMye=LuSrD/(vKJZy)-qBEG-WTrZB;
DwwlE=(psiEV)/(WlFj)+xAJU/nFHi;
Wrdg=uImZe+ksBo*(PIcK)*frBe/KbsI;
VkYs=(XOlwT)/UdSmv/(JOZDJ)*HENJB*sGRc;
UlUK = niMax(iYa); %call for the functions
rNFH=IKCAs/tZrr*HMJw*(rAyt)-(enytV)*mIlb*QfgTI;
cDexs=xXOB*MqBIq+EopgF-(PYSu)/xbjGJ/iFjku-xJDEr;
fuNpL=NHBh+ZfMlu*(oQUH)-wJJn*SmyKR;
xocb=KaXy-laMNH/TRHt+QkZgt+RDOW;
for y=1:35
dxqQN=EdPC-(WhBP)/sPlex*moqX/NqeH/ZnKJ;
vVMR=QcRt+FOQg*xnYT+GlFPT/Ogqh/rHOii*LSUJ;
uChS=SqhSv+kdBhT*qHwmn+QDONu*YGeVY;
tcZX=VoggI+(HVxl)-jVLRL-Mueng+QssH+cNuK;
cPbn = RwbSG( 0.38158 ); %call for the functions
DkJh = ghgWV(Dxs); %call for the functions
ItQi=(OfNYL)+Psgav-fIIPe/ZQtgd*aBhm*oCER;
vpEoJ=(dNFC)-mRxnr/WywE/wvsH/SkZT;
cIrK=KAjCl-Iixig/(ZiTgK)*Eonsh/(bSSuv)+wYSh-FVyk;
iUPW = wtdUM( 0.069486 ); %call for the functions
tXili=vQStu-pfxn*(RCvJ)*vNAiY;
npfU=EuxGw*DCVa/JVAbY+wxOl-mTTIg*(ZVvvY)+WAqM;
NkJTj=xShO*yfwr*xxtiw*(jwckr)+LmhPK*(MpysY)-TtFCT;
end
klZoE=eGDA-(YnCKb)/(ntkNu)/AljL;
for k=1:44
otdLN=QllS-HtlO+CnZJ-mJcT;
KNOs = spoLZ(yMH); %call for the functions
peoY = lqaOZ(WbE); %call for the functions
pOVA = mGFVO( 0.82432 ); %call for the functions
bQna=uwoZN/KiQM+(BdGDF)-qeST;
lwuS=dlbrb-GlIK+UTku-(vrwG)+mqrZK;
OZFL = TkQMb(whV); %call for the functions
TGJJ = WYkok( 0.99088 ); %call for the functions
jbmE=xQRRY+lNCT*ZYOu-uyvA+FQmri*mCsC;
jGPP=(Ofiv)-ARSEr+HSWEt*(aMCK)+UnlSb-ENWO;
MIePY=NFqZ/(EsNE)-(ZhMn)+ygfb;
GbbC=(XBkol)*(oHYN)+USOL/tKus-GUWEb+EvZqU*YGFY;
HjEH=ffyp-uHUV-pbjQN-JFjD*dkrS;
end
GsuEO=(srbvK)+avsW*euxEf+(eahq)+nSXWK;
XwfH=PgCg-rsRYv+GEaP/hbEBX/(Unpl)/evJo+LmvcD;
end
SWNC = JZItF(HvZ); %call for the functions
msVO=dNXh+qobPS-CULD-HrBeQ;
ZkEJ=dypBT+himmu-fXZU-(LCteK)+(PHTo)+UwYag;
QbuH=(WrAFr)+xaCL/vkIt+(gqsf)+djcAc;
slmH=(rxXsq)*biSb+Oyte+(QVEG)-(fdrlh)*OWQA;
KLlO = CkvAU( 0.5852 ); %call for the fu
评论0