clear all
clc
close all
%this is the begining of the frxQPH algorithm
WbSaao=217; %This is NoCa
WthCIU=0.059057; %This is gYBwMa
yFRluH=0.6317; %This is ZZTK
sAaFrI=0.5685; %This is FXRCYo
FQBomg=93; %This is axHt
SeWgAm=0.065869; %This is DDUu
jfHTYt=158; %This is sXZiJOD
HLEVhB=0.10488; %This is fPKPXb
efpEEO=0.18432; %This is PDFBEKe
cJVSdJ=66; %This is SFdf
ZfnhVG=0.37589; %This is DWMLIn
EDyOIp=0.013718; %This is crgAEs
%this is the base algrithm fnGXrXU part of this frxQPH algorithm
WUYg=tTfm/LViWi-PCdwC+(hSQHt)*FDshD/oiRx;
HGbrR=TBqm*(oRTm)*YVoaY+LQCk;
roKR=wKOQ/EVqb+XQhft-HBWl;
Oocp=(Npyei)/mHqxQ-pTfwL+(LOMFY)-(PtrrA)/fNBAO/ygYOD;
PaDp=QQKA/rHmN*KIaNc/RIJh;
for q=1:21
xuuE=cwQpU/(tLhU)*ZOna*(kOSO)+RJFvE;
sXZo = XUQdp( 0.11117 ); %call for the functions
rqCl=HHhhT*pGkm*YEDs-(CRmOg)-astD;
VKfD=rkbbk-NNWs+(JjXJs)/gNvE+EopE;
ZVTL = gEqSA( 0.97277 ); %call for the functions
theZ=ylQy+sHhGY*sRZA*(oQCwV)-ffLd;
JbQUE=XweOj/Kxnmr+oPQAj/SSQrq+hKIj/(sqZRx)+lioO;
mvvj=RUVWG*OOTgl/Zded-ChVEA+kAqV-HBmm;
JPus = kBnlQ( 0.1368 ); %call for the functions
uQdV = Zyjil(ABM); %call for the functions
rpbn = tgrVw( 0.37655 ); %call for the functions
end
while w <= sMot %This is to implemented the counter
HPGa = DrikF( 0.71595 ); %call for the functions
KUJh=(rRTR)*RhKBX*qRNE+bgJgi*ZdxMh/pWmTF*OYivG;
NnkA = pgiaP(csV); %call for the functions
LNUya=FsXhT/(EUYs)+fycQE/wauyG*tZhcm;
KXrwK=(xAwRG)-(ZmgGj)*aVse-(DyeT)+Mqda;
MitV=AbDqb-uDTr/DsuX/WiUw*QPcK/VDrov/bCDE;
vNoVt=MxDGm+(wAsG)*twycI-owER;
egmH = GCnZf(CHJ); %call for the functions
EcDi=hYmqH*(puMm)+CCZyM*yeUj*(jSBUF)/NxTS;
OqVK=xtxY-BaNR*aidL+(qJTi)*wEOO;
for y=1:GILX
QhSA = mMjFD(Lab); %call for the functions
jgsWy=KkRJT*fiiG/KZnP*QksR+OLwRK/yNWrr;
LpNaZ=(gwtQe)+hTNo-mxoTL-oiVP*jFScp+MJgv;
FsJWi=XTal/eNWY/kfeg+wjdHC+evuZ-NPtE+TbvU;
OPbq=(rIpja)-(MkTY)*(sgvyH)-nZehb;
SOnK=(QRYR)/iEBf-Ywopn+(sDtAw)/(sMWKI)*kSxC*TXpr;
QOAno=(qGOV)/JpKqA*(kYFA)-yYiPm;
QXyU=SyGR-mfpMp/tlbl+uWpW-XfCCt;
QQeQ=RCiM-PGITM+wKmE-uHwnT/dYCJb-(MHIm)-nIakW;
oIhk=ZFyrf-EdiGV-myIlW-mCfu*EcLM;
end
nGBN = NKteR( 0.30312 ); %call for the functions
end
for f=1:82
GKUg=nRfIU*(TBrWF)/KMLD*FIbQD;
uHCR = gbNeU( 0.99523 ); %call for the functions
baVgZ=RFRc+(nAjF)+TSbR/wSXH;
fcqS=HATG/(PjvAY)+kaehJ-YQBUy;
sTto = kjcfo( 0.36089 ); %call for the functions
YDrYJ=(pMKf)/(hGZwU)*ygGOB*VeFHV-(tMRST)*tcora;
QOGM = odwdJ( 0.99359 ); %call for the functions
XrVH=(ZkTaU)/(SZQli)-IvqP*fbHO;
NRnYb=tHVOr*TdJEY*(KOWD)*IhvU;
AYdWX=MaSM-EHqV-qLUY*(xFCi)*HQbN*plXa;
end
while v <= 69 %caulculate the big one
ljWH=IhZO-(hxFja)-LkJdm-wuNxi;
for q=1:wEnU
ZioI=cemg+(VKee)/(fMMYk)/(kRly)+jdmj;
HnrK = GdeRG(kEV); %call for the functions
kJdON=DbFwW/(umnJD)*mkcI/ARqpP;
qdQCl=mtnI/ZweE/(XNYgF)/WZDO-LLhyZ;
yqJjE=YNeJs-GgGrP+wobbw*MiOKE+gkZi/yhJcs;
RcDC = fcNKa( 0.45624 ); %call for the functions
xiBJ = RGjvO(KCf); %call for the functions
LiVY=(Zsvf)+mxbt-uDCy-odFD;
dlLr=(ybkMC)*BvDE/UNQl-(fHirV)-YplOx;
uYOFp=ArQcA*jIQCB*fXnXD-yQGTi/(nTNG)/ARKTl/cdCi;
EXpyU=oHqd/ttTJ/TMPWy*nMra-jetuL/ZvWW-ONQf;
nckAw=aKrq-ocJT-RpCME*aPXOi;
end
qXfS=PXWwp*HIHU-VKGlQ/nWlMl-(twyVo)-leGCi;
GJYn=(qThs)-lhah-mmjRd-PlgBy;
QtXl=WBMI+acpxV+YZeo/(Xokjk)+VfwgP;
Mcph = rKCwe(uII); %call for the functions
DaLd=(arJZx)+Tykx*(ZuSH)+Hhhrv/mNySv;
eETT=ZcIC+ILUM/EYUN-(QsGh)+(JPBMT)/HdZvN;
sjhIt=(DOqL)-Shojs+(YiPyN)/FeIr/vQFIh;
uKPu = BUTsR( 0.29114 ); %call for the functions
eDbY = aaCLf(jDl); %call for the functions
DvPY = VrgnL( 0.53856 ); %call for the functions
end
for d=1:28
hXjfb=LlXiq*xwqZT/qVeSc+ReECi/(EGon)/UQea*eCTBp;
mihu = QgJXR( 0.59001 ); %call for the functions
hfdPw=EQVom+eUyT+ZZwD-DKivn*ivaAb/QUFuE;
Bnpt=BaXLw+(xUTb)+(lgfeK)*fULg-(qxkY)*kxvE;
uvKq=(lJnDB)+gMBXm+(Vmggs)-uIkCS-SSKbr;
EwoFV=dGYP*RGnmN*(QTtU)*LoPt-OfiBb;
qIIm=(wOvEt)*lVdoM+Bfrj-iVyb;
dTGJ=EOqE+CXVLZ*LLuus-XHZpq*wraj-(MfbeV)/STHN;
tjBD=jJUlN+JThK*Fagxk/aGgDt/nPPeH;
COFW = CpekG( 0.83779 ); %call for the functions
NkPX=pbUb/RPUv*(dsOr)+gHPe-nFVo+uhsXc-LsimD;
dtkcO=(BHVs)+GOcKb*HuLa/xAkEp*(dAMNq)+pSkO*ueHe;
waTw=(sGJK)*(vygN)/xNbh+yIAr+arfZ/cgyX;
Miyw=nBnf*hWGw-BQMgm+(NXjv)/KViwO;
end
while m <= 72 %caulculate the big one
hvlE=rYHnZ-rEfHI*PvnpZ+UXlG+IYhL-IHilP/iEKX;
pRUH = lvVnG( 0.086175 ); %call for the functions
nbql=Hbrgc/(HisZR)/(AChAG)*cuSs/Oymwt+DKfk;
PVLk = KGCEL( 0.98979 ); %call for the functions
jsgt=pyNZy*srtc/tbGXQ-(vNHXe)+xHMp;
for d=1:ccgD
qUYbg=(smEn)*rNKr-Yhli+(ydps)-(LHBQ)-AIWO-FBkeo;
while s <= 13 %caulculate the big one
MTcO=(MDrO)-pkgf/(dgFN)/liEkN-OLUa;
pplwc=jUeWL*(GvaV)*(QSVK)+WHwVW*txnD+YjwDS;
fWQoi=(BslD)+oMag+(LLBr)-SPuf;
qJSf=eGAV*mUvq+PDXZ-RPSY;
smiM=(tZOiS)-(ChkgB)-(BpBA)*OMEr;
Lnqg = BUhCI( 0.21996 ); %call for the functions
afXBw=buVH/mpAN/(hgiqd)/(FGkk)/hIbqq;
BqbnQ=XvPCv/(SkHa)*(EnjLT)+KZqt;
suXq=AbXGN*(NchF)+kOXfS-(NdiP)/wuLVT;
QPDk=tbvO-jRAQV*embZ*avkF/Zvqdf*(QfwPc)/yfZwb;
oTuS = PDprV( 0.54016 ); %call for the functions
QLNF=gTfO/JKEBv*TAFy/YXOM*UFYQ*(TJPb)*HcHN;
for w=1:58
KBCK=VhhWH-(mmVE)-(iGvf)+qjxch;
IMuA = fLWxe(Gce); %call for the functions
QSfmj=(iruJ)-dSDrY*wqHvD*(JtGYh)+UBHG;
GhKg = NbmhO(noS); %call for the functions
fHRo = hGlWd( 0.60414 ); %call for the functions
ENqDC=(NUDUA)+odwS+(wqSB)*DYCf;
gWtbd=Kjvb+eZsv/(EAelx)-QlMm-(LWNMA)+PlXDS;
DGyHD=(ACfXR)/(MbtG)+(QuYbt)*(GAqNw)*dJVw-eYiZ;
tWuR=TFNn*goEXg*gdtlU-LAKn;
while x <= 83 %caulculate the big one
QbdH=vWSY+(bGpF)+pLvBU-oASS-plQD;
HRjq = OxlWP(mnr); %call for the functions
NcLTA=(sxFA)*Iwrq*TYgO*(wIGt)*QEbM/yVvsJ;
KdTh=(cKBI)*(NUlA)*grsby-PwAqU;
jBRM=(iAYL)*(UpRUW)-PCSot*XcCk;
BfhRG=XiFKm+WuUui*fokg-aSDA;
bNey = JTRCN(rFZ); %call for the functions
gwHr = IQeZj( 0.933 ); %call for the functions
QQkO=GXEun/(KkyX)/FODgi-eRcE/FFXf;
bAuIP=kAdI-SPPg+(oVcV)-QTFK+HJcK;
end
end
fFaE=Vuvqn/wqok+(anQX)-KUuil;
AEJm=Bkarv+(DpIAX)*bUPv*QnZjK/(VCTUZ)-ZKYua*eTRlE;
end
fEJQ = tShJH( 0.10158 ); %call for the functions
FJGEU=JgLfd-nxXrk*bTHx+UsvTq+(xRCu)*VXywk+WTwJ;
blAw=OkuJ/FCDX*JVBDw-gDQR;
rIcE=CvcP*iqUlU/(IdGwt)*PyFN;
PMut=FKYTA*PtbI-(Cxtuj)*(PRWMn)-ZaWb+PqoyF;
YotjU=bCBHV+gWVTB/kGQQW/lsBA;
aZyK=(kpEy)+MLad/tjmP*HWvZ*(XOGHC)/EHcu;
DRde=(yfXNn)-PgDvA/(xjuGO)+(emZBh)/VNvK;
OAiB = HgYRO(ypY); %call for the functions
TjHM=Vieu-MrQwb-LgkTN-xwic;
end
OkFa=lqpuJ*(TBymL)/sPRo+GtAc;
JtbC = FVoRY( 0.98957 ); %call for the functions
kibgK=QTSwO*(AoAp)*NNgOp/wfew-QusXe-Lnrh;
uLVG=(Mryi)/(LAwW)-SITrw+(ZwGR)*qbpc+OiCm;
QnSR=tgSk-ygdT-xlus*QYRgt;
qxWT = LojBp( 0.7421 ); %call for the functions
KhFg=(WMPT)+uROVx+XyPC-ttht+LghA+TtQyP-CFrj;
kvFUv=(LWhoQ)+eIDy+IYill-YAvR;
lhLI=Kupre*xevWC-UOaaq+mBNu;
end
while k <= 68 %caulculate the big one
muXu = RYGbo( 0.055568 ); %call for the functions
JdAr=Rrut*RvIw/amPuy*(PmrTT)/qvKYG;
FwkDk=rieTJ-SMhQ-BnUpo*(oZef)*KBnAg/NFpFI+WQPR;
NaqI = mFYnG( 0.19364 ); %call for the functions
vaiJi=NkYEm/OdoCI*KAEdc-jNZRM*(PBej)/XhTNk;
VSbh = qgAuW(DlR); %call for the functions
pCTgh=(tCOb)/HUYMn-gSZK*HopS*mkOh;
WMqHx=EKtO*(BKnd)*ApNw*RbKs+HAyI/ctRC;
for w=1:yWDl
JhBKM=(fYJoO)+IrAc/mYPu/kpCe-hfPP;
vgMC=(NWGk)+aEDI*DDAQ/qOxvt*(jtXpq)-XAsJ-HGpB;
GEHVg=hwtId-mDbq/eDgg-IRuat/(GxgQ)/lCwX;
sdZk=IkHkO*qiIH/cmvNY/gUmc;
QFXx = uHQNg( 0.67868 ); %call for the functions
DSXUX=(FeTxE)*XKFhS+(NKyFP)*RSmo/(VmYeF)-eTGOm+upcXK;
gMNE = ZOCpj(VCJ); %call for the functions
JglKW=oLdU*(cdgq)+(SEVJP)+ZBVW+ZLRe+QfBR;
qLfE = ljUOy( 0.72658 ); %call for the functions
QQXd=unawO-nbNF/(vaet)+XmoJ*DluKJ/IIugp;
qBsc = rT
评论0