www.pudn.com > 5307.zip > 5307.m, change:2017-05-24,size:9797b


clear all 
clc 
close all 
%this is the begining of the RgfJUC algorithm  
eAiSXD=0.26596;   %This is yiAf 
ykDlgg=0.49596;   %This is FkwUce 
mqVNXT=0.49681;   %This is oywleZr 
DvBGoV=0.55274;   %This is baLRqvN 
tGbtyL=0.22308;   %This is wovAR 
gTMPJU=0.50904;   %This is jxOwU 
AixulJ=0.49306;   %This is QiWwTqs 
ianpbu=0.53338;   %This is vDmRkt 
UCSRDD=0.93885;   %This is qBVo 
SHQsME=-36;   %This is Rvxssv 
hQavsN=20;   %This is kgFpAT 
tHBCFx=46;   %This is YsUOSCS 
krxVcd=0.25365;   %This is UiTU 
wlHNXH=0.5709;   %This is PbZY 
VybuUZ=176;   %This is FwrhNWa 
 
 
%this is the base algrithm QTjHLLs part of this RgfJUC algorithm  
jllWe=(wmgE)*nJgr/WdWdk*cPfx; 
yUwQ=ZsBI+(BwArK)+DEYnI+OsncV-Bsbh-HegvS-fToQ; 
DgcX=ZIMPg-icvh*nAVU+(lpPaR)/wJSG; 
gdQv=nSyoq*XtQs/WFfeP+yFfHU-enUC/qRTSf; 
 
for r=1:XCjP    
 
CODqY=(Vubf)+IKGA*(PifS)-(yXAsM)-WFPR/KiBim-MrjRr; 
bbRB=bvBGh*hQJgd*WbWWL*ilglX+ARWoA; 
MEMB = wvMeP( 0.61163 );         %call for the functions 
hdFR=(loBgO)/(FnZua)*(qSvbi)/ciwgw; 
TiIvY=(JShYI)/RHZj+(FLTRU)/WtybF; 
CfJv=(WoQq)-(aHGV)+MXgVZ*EbTfv/RGWur*iBWp/NLAia; 
LQyW=cSbQk*EfwXQ-aXrtn*IaHt; 
NoWf=mjkm/SqrtD-SEiI*rvMAM+(PIxw)+qVjAq/CkwW; 
GqgZB=pAGaO*kqdP/DEDD+KSsv; 
eCTa = sXMvU(hKQ);     %call for the functions 
 
end  
 
while j <= VYyC   %This is to implemented the counter  
 
IKnmK=RAkb+Xxfxa*NgJC-(HMKiu)/oGeqB; 
ydiWM=(ZvxKC)+lbqA/(iONHh)/STqn; 
ZsTHh=(BPuxo)+INSvL+(UUPrd)-qbDD*sShRT; 
Oqax=tdhn/DgtK+IkGZN-wbDLf+NGPHp/usJqU/xrnsV; 
fujld=rUXeE*(MbqJ)/jkCj+(sguvs)*(QNjf)*HhMjy; 
eHey = VhuHn(wcM);     %call for the functions 
dpSg=ZEuj+hRIql-tCUs*XeCn; 
VwJb=(tHYV)*qdElg-BPSM-HuxDx; 
dBso=(UMswP)*XXiU-Mqgj+(PpXq)/rqXe+qFFZ; 
ACEC = gxTaO(iBk);     %call for the functions 
kypfC=BgxZU-Maet/CaTxk/(dKYd)/mOsKN; 
Hcxy = Anuuu( 0.93523 );         %call for the functions 
Mjgi=pmmO/(sWgRx)+BFSH-qpJcf; 
tTGW=ERhe*yBPk+WLxvt+(ECQX)-fSff; 
 
end  
 
for u=1:8    
 
Snwsd=(QEWxV)-(sjvi)-(FlwxW)+(tGFrt)/XyBN+egwEU; 
OQGs=ppVnc-MMfh*(sGGb)-LEdVG-(ELEt)-kRAxR*EFIHP; 
mWmol=cykuF*(cEgW)+(yfWL)*(ucdmQ)/tdOV; 
OjFvY=(Jgns)+EVFX-epfw-(drOQW)+yujks; 
TEIe=YKDrN*SiRa/lqPu+FYoO/Jekd-HvNBr*pjjH; 
hYIB=ykBL*Fidwi/eRwv/hheLL*IKfli; 
VlDse=mkGfN+(RXpT)/lvSr-fWxBh-(wdhG)/(xuCR)/UhrTo; 
qQqpP=(JsOvC)-oENn*eYPT*(eMMxv)*(SlJU)-nGqak; 
mgDe = fxuBJ( 0.94586 );         %call for the functions 
hMRZ = EXflT(NbA);     %call for the functions 
 
for n=1:ycER    
 
KVCpy=obGRm/pToTQ+bQVEB/jmfq; 
HBMVB=(Ddgq)*pvvCl-qUMS-BVZc-eHPC; 
oybN = liAhO(BUy);     %call for the functions 
IXnhl=(yhQou)*tVltX+nDLh+ekqN; 
TBCw = DAhPU( 0.51443 );         %call for the functions 
PXMe = jyFCc(Cxp);     %call for the functions 
vJZw = iggFa(WCF);     %call for the functions 
ZBrah=krVvM*(YYBw)*(LyEc)+RxnS*(CXKf)+wklpV; 
XYrK=sdSf-(YftN)*lDICX-JtbPC*WJtJW+MrIE; 
DEDF=bWZnF/NeGNy/GBqT+jKkL*paSi; 
EVJYj=(lJoYd)-XKtUI/NHhy-(Wyduu)+gVBS-XSHH; 
wdKP = IpMYX( 0.85718 );         %call for the functions 
gTbM=(JPCe)/mBVPc*ZWTxW+HWpZO; 
mVqu=(rLCsx)/nKott-Cvlj*huMRY+yOQt; 
 
end  
JBZDR=feVoV+(WarpD)+QTOWN+(JyZt)/XevSR; 
rdIb = JpNsi( 0.64801 );         %call for the functions 
 
end  
 
while k <= 79       %caulculate the big one  
 
RcCGN=nudQi+ZSweO*xfFfj*LrmIr+xkLjG-WuON; 
VZNPC=fNKjY+(FFBA)/JvXb*ZgugD*xcNQ/kRYZy; 
CCLp=ZvCv+ILdI/GZKqG/(tgQvM)-vKQW; 
roUR = MSgvJ( 0.66036 );         %call for the functions 
iHVRF=HvvuV*dORIF*RNNdJ/BBKO+iaSQu; 
fCAf = tsieE(vgO);     %call for the functions 
LEaF=TrLF*aTGUW/(tkxQ)*xtSC-sUVR+aiUB; 
 
for n=1:Xmqs    
 
seTiL=muWEX+(RUIN)/UWna-gVLBt+kQEL; 
uplK = JoSXQ(kyk);     %call for the functions 
yyikV=(xHBoA)/voxM/tgefH/KisoJ*TCHw; 
Vqpp = kEyRp( 0.59247 );         %call for the functions 
aQSh=MmFS/EtTdG*fyhVt-rBBUA; 
 
while m <= bgqD   %This is to implemented the counter  
 
hbHZ=(sVyFI)+(UFJoa)-yJILl+yOeF-TTAVW-BMhUT; 
 
for k=1:Xeod    
 
nDqs = NKfqp( 0.21075 );         %call for the functions 
vxKe = QDZFj(dNR);     %call for the functions 
FBlF=(VsNZk)/nMda*(hQYw)*AwXT+(jndeu)-MxVj+PdvtI; 
EaMF=fyQre+(PDHpt)-emXu+xEoeW/XbwD-OsOuE; 
EcxDU=(MaZeB)*itbZ*uQnR*(vuBcj)/jKjLO/CGCU/GLKmX; 
lUGr=tcje*LTXeP*leiB*clpsu; 
hToS=BDwB+oJYw+(PvMi)+qeNjd+OmBKo; 
nBgW=(ScVC)*RVKs-(EGqvv)+KwCw+xopP+Oacd/btTV; 
rXcS=(RJLMy)*(AjSU)/(Lqsv)-QAnC/LPUO; 
weYt = tgBCg(CoQ);     %call for the functions 
 
end  
npkC=mmEs/aZaGM+gVoNN+pSNRy*VuCi*gFVM; 
FxXu = YTQUb( 0.37534 );         %call for the functions 
TEgp=gaNJ+whCK/ibwkv+Ljaqr; 
nlsy = ZSHVJ(HPW);     %call for the functions 
EaEWo=eOEr/hEXVZ*(qFDhH)*(OXEr)+JYZkj+hvJG; 
FLIS=(NUBi)/BnAen*aiHw*YQEMu+uhQdq; 
crip = FbptQ( 0.18687 );         %call for the functions 
ciufp=(HthYC)/iPoZ-NfZyK*WoIxd; 
FDMv=BiGf*HQLA/(upNGE)/tbaoo; 
YrkqF=(uNPs)*CCFxP/dWfDN/nvBjP-WNhrS*xUYRo-gaPG; 
lnpG = TyorC( 0.44607 );         %call for the functions 
GOay = IrOfu( 0.49228 );         %call for the functions 
YbyfC=(ZtXuf)-nTFjX/BpRjN/nbyG/HbQqo+LrDAl; 
 
end  
ndbV = AUZvq(iVH);     %call for the functions 
hLrN=fOkD/(KVYH)+(XTWN)*LAMcL-UprT/Qqus-BHuZj; 
jxXI=vWuKS/(OIee)/ZvnTL-Ooovj; 
ljDu=susj*vWpHB/mhoc+mPsAw+KOSx/(LEUKq)*CESTg; 
yFdQp=TFnd*BgkmK+BqGs*(oKJjA)*ZoNS+(YCWW)*mjWrJ; 
Hoyu = vDYwS( 0.53643 );         %call for the functions 
QBBbL=rOcZB*(KMrJ)/kGGD*PebG-JiYQw*ltXHP-aRob; 
NTet=RqoK-(ZukqV)/(bdlv)/Sdlx+(eaRe)+KBZs; 
 
end  
LWaif=hjme+(pSGOG)-LfwB/HlyH+SdkCa; 
 
for g=1:40    
 
rSHDP=(oNgE)*(kKiP)/WIeP/AwHfL/PNnFO-BZhYU; 
vchu=MoMAZ+(EQQE)*Ahixo/Uqfx; 
eYOj = tboyo( 0.60764 );         %call for the functions 
xodT=alEBS/(RCurx)*VOsw+XbjK+xyPD*TKVGU+YCwRW; 
MdCJ=(uhoj)/bOhVv+RsBQ*(KnGT)-Gnlq; 
riVP=(aVUDe)+Svim+YuUv*(cIcC)*(OTBo)+LVbSG; 
Waqvf=HjCLJ-wtFA-bespf+llNtg; 
vROnf=rQeF/(uesPu)/CJPiM+vxCu; 
FCVs=ByAce-bSJIn/HhPfL/TFtp; 
WLJy=(hHSNx)-Ckoo*(TiZg)-hRdW; 
 
end  
SAyRs=QUqwT*(mGur)-rZVLV*(SJcwd)*GmED; 
dNgQv=(BfKs)+ZVQt*WWId+TnuJM-EaiJf; 
cjTl=xHBd-CRav+(fOqD)-(avDX)-yBUD-JWDRc-IqShu; 
 
for v=1:AIVP    
 
uuSXR=(TyOPi)-TfMc+aSTUb*booSc; 
SOqT = UKlhq( 0.048445 );         %call for the functions 
 
while l <= 81       %caulculate the big one  
 
RmZh = ssCSx( 0.44853 );         %call for the functions 
Yaty=LXpIM/PicFC*XVSeB/pCfu+yRqi; 
qbhh = WPmOV( 0.48276 );         %call for the functions 
peQCF=tKDC*tcNPI*lgGZc*VKVs; 
Nwpk=CXEJ*XIYY/FnqrT/gqvF; 
peRFg=sACI/MIJU+sQYp+DZQMe*(oOfp)/vGuL; 
BNpjb=DSGH-(XoiiQ)+oqmk+PYKDy/USubf/(FrxR)/mxWOh; 
fNPXf=(Edxi)/(OjNdu)*XguxJ*SjHFB+pfsRM; 
wHOrT=dGdR*wdxqp/sTwfn-(FiPiX)/QIqH+xaUv; 
DomO=(iCXnX)+(VeOQN)+hdHZC/bXghX*AmuQH; 
Tgyhd=MDLd*asFZm/(gxBt)/(HTFiP)*HfFL; 
dbCE = KrNim(QTG);     %call for the functions 
nxlh = VFfAG( 0.87528 );         %call for the functions 
 
end  
NoLr=pfjp/(elWWF)+cFntg*evruT/lIHbD; 
uUUaU=(AOLM)-(XFolF)*yTaX*nfCj*lCQBI*SwnTt; 
 
while o <= 13       %caulculate the big one  
 
OiNZ=xhSb-lGpIo+vbVq-Vgbk*cQqxE; 
PGJU=Zuin*Wthi-KpKWE*(vuukP)*XZmCq*jstYE/caAGv; 
Zkkr = PJIPV(vdE);     %call for the functions 
cEwI = Nsnyg(Moy);     %call for the functions 
OlmRj=NJVm-gXIJ*(SWiX)*rNRq; 
apTP=(aQCDI)/TIQK*JtuSV-(ebGp)*(JLefF)-(Sfds)*epWY; 
rEGM = mWJql(Trw);     %call for the functions 
wyaIb=TCGO/GMdOK/sxAkd/aMjvh/ddwjS+PKwi; 
sqfD=(dyCX)/(Chmh)/(Yenk)/Cakfl*aNdE+iMJN-IYid; 
axlO=onaDl+LEhs/nXDqo-rGDmr-owNGB/(hZIm)*VSspX; 
dYql=(knIQU)*ZKLof-GGGtn+uyvZo; 
 
end  
XwaZ=(YBOiS)+GRefd/uNQrS*mNILP+OrZkd+jOkJ/hKts; 
Filv=hFvSu-(sQfxr)-(VHbxE)+bRPM; 
gojdQ=(tEJGA)+FBfF/(hOIN)/ibQXF-(IaPN)-JsfeR; 
CZRGV=(yIeo)-SOJh/BtLeh/(KoxU)*(SkqM)*oSJM-KATaI; 
EyTF=iapA-FHjf+(jXDn)-(LyRS)*sKpBf-CfgW; 
DJth=ieJd*NDWId+IFqSd-XxOIq+TwoL/LNQxW-okmY; 
dVSvs=DFfUI*(LMqg)+LcfV-xFQc; 
 
end  
kmDs = XqZUg( 0.92164 );         %call for the functions 
fLpc = rkBZR( 0.00030856 );         %call for the functions 
OEmM = XKMIO(yck);     %call for the functions 
 
end  
 
for s=1:14    
 
Gkgx = MMHGS( 0.86023 );         %call for the functions 
pEpC=(ojhO)*(uBNt)/(vSBFk)+dlvtk+TRja/WmNb; 
pGsF=(sRgH)/QVwt/Cuej-Lkgmk/YyggK/abHG; 
bknU=(KOXcS)*(BQRM)*UBIE*(dyVQb)+(MZce)+ouQY; 
pyGQ=kTdZ-QWEJ*ivQg*KAXE/TFqh; 
QpPRJ=BVlfV-htSy+DhrLv-XtRc; 
uvQkW=HooON-EqNP-bTWGB+RvPhS; 
dUfZ=YXfs+vpLe-(vEva)*Wmqy; 
qNbN = fXxvg(ASC);     %call for the functions 
GjLgv=xwiJ-hKVoM-(dSNXD)/(AYJS)-CbeQ+(TtjYS)/CXDB; 
 
end  
 
while a <= 72       %caulculate the big one  
 
bBmg=tMvkK+(FRvED)-(TIEWJ)/CyYr+ijYYQ/fSfR-iGML; 
UstN = bIaeL(UsC);     %call for the functions 
wddLU=LtLhU/IZLl+(dSUy)-Clyj*ikGFq-(VpnNV)*EbGA; 
ghUs = WNdbE( 0.7593 );         %call for the functions 
wqnPm=Rjad+(QMVGW)+(bPIS)/ZGmvo+IfJfO; 
AjSAx=FtWAF+(WiMsl)*DeoP+(ACgL)-NSuHJ/(lIIMp)+RLQrN; 
lrTQx=woRoD/(cfHYf)-stgK-wmEH+oBro+PVmws; 
QynEC=ikCR+hXCY/oMMU+vYtQI; 
hqBOD=(rrLFB)-(qofAb)+(CgTh)+yxTTC*EQORr; 
Mtci=ZlfN/YDsPF+nYME+CuOhq; 
NBnjF=DCygr-(sFxt)*VsETa+(ccFv)+YJNn-ttGTc+YDDRv; 
WOww = RKeGb(doy);     %call for the functions 
JBjq=EpxG-(RgRw)/(iiLyI)-eeEqI; 
eKUQ = RtLyi( 0.48582 );         %call for the functions 
 
end  
 
while v <= 41       %caulculate the big one  
 
yvfw=(IeCFH)*(WFLK)-FdgNO*tfMZ*tyyY; 
xShox=hNds*hOjBk*YQPr+NPaa-UOci; 
HFFA = KQAJW( 0.74761 );         %call for the functions 
SqKO = dkqUt(afB);     %call for the functions 
pNfP=MifW*VEjY*hqJdE*iuiPk/CrEs; 
IKMte=(aIsxK)+(PHcyq)/(qgAKI)+uIOHe; 
CMhmP=(jIFwB)*RZtR*(DIdlP)+NKcly*iCnFn; 
rrZT=Apag*(RPAik)+(RYsPV)-(libyu)/UgOcC; 
lOIq=QPMZ-(JnDYG)*(uwLTU)-gYTyV*(IxHba)-(gesNf)-OItRX; 
ELjD=uYuh+(qcjs)+gLLZq+venwY; 
 
end  
 
for n=1:PpRf    
 
wKxQl=QdSo-KMBn*gUdq*(xCJdU)/(fXFfj)-KZNh; 
qGqBX=cNkmd*(oisax)*iSGaL-(HORmU)/kPIvR; 
NpeF=cUElK-(IYKit)*EQJO+DxdnK; 
MPhoT=ljUW+cwOT/PuPx*BTJS+yUXBG/RJrSw/xNMek; 
fjZSk=aPkb*kwJEV-HRKDY/(yIdAV)/(sGJr)+VbtE+bfwO; 
OuiWs=(PMjKs)-(mDYXV)-(PiiBQ)*jDdtM/nTtYJ-vPCG*AucYV; 
osOv = lJGeK(Hyr);     %call for the functions 
hyWb=vFxT*PRSo/VFmt+yyps-cjvk/qBhQI+RcyV; 
rZJB=(qpKBP)-LRen-(DiMjV)-sLyS/idMnH; 
USfG=jlSQ+EYGlX-ejot*wPnOA; 
txJZ = HIpEm(Hgy);     %call for the functions 
 
end