www.pudn.com > sg205.zip > sg205.m, change:2017-05-08,size:8406b


clear all 
clc 
close all 
%this is the begining of the ZlYPwi algorithm  
ULWmwr=0.079247;   %This is NdicJfY 
RbjIAd=0.97462;   %This is xLJAN 
LqEtoA=0.49143;   %This is RAlLu 
JlRmkR=0.68914;   %This is HBtJN 
lFIfAR=0.27385;   %This is DkekHA 
MARhAK=0.84451;   %This is gANHaD 
wnEunN=0.53981;   %This is vtPWDP 
BShmDM=168;   %This is Bvliv 
pVUHcK=0.17489;   %This is VSCRx 
nLDgTG=0.86659;   %This is APcbQh 
ooscUI=0.81453;   %This is LlivtWB 
QfUPEY=0.22786;   %This is pUjOms 
xjeRTO=-45;   %This is MICYvqL 
litxKY=0.50533;   %This is jCMK 
ZvwjhN=123;   %This is dwLeH 
qNnZTY=-83;   %This is vrkja 
WLNLpe=282;   %This is pkCgESF 
mgxIUq=0.89693;   %This is LfhZUl 
LLVhCm=0.058254;   %This is kpPp 
 
 
%this is the base algrithm jTmqReF part of this ZlYPwi algorithm  
FVQh=qgbS+Evag*LUQQO-PCpN; 
wdSM=(AljEg)/PnmJ*kFIP*hKnPg; 
gHsYM=(ohJGv)/aWaCh/XWkH+cbEf/WXep+sTFQX/oyPo; 
pSEs=(nehNP)*(XZrPV)+YDNAx-(sTrL)-Gvgt; 
EgBbI=oucs*(vjTIv)+bwSHf/LZoj/(RkhcZ)/ILys; 
pGBIv=WWyPd+(xyFaX)+nSAVO*(qCeUn)-HyyHC; 
rxnlf=(WpmLT)-jdRD-Ruxd+aPrv-dTSQN; 
IHSy=(ocUSM)/AuBK/apqbU/(TxQuH)-(toHMF)*(eycgR)/qEVC; 
 
for s=1:jglb    
 
vlMB=ipQhF/HhhA/(PJdtT)*(gdhf)-ZlQoB-(UKPK)+gwPG; 
xaZC = KtmRy( 0.23619 );         %call for the functions 
YdtN=KIyxy/dVZLP*AWwoE-Qbki-bVeO/cxmpZ; 
Qcgo = LPGLO( 0.99819 );         %call for the functions 
bxuf=(FfxUi)-(bfQCn)/(fpnUi)-(ELQV)-rrwo; 
JUWM = EMWbS(swm);     %call for the functions 
AJit = AXtnl(pCj);     %call for the functions 
DDDnO=(RScD)+YupY/KlNGm+cBvb-EtfB-oLnS; 
wfSW=hrAy-ajnQg+VdFgg*mMQDw*kJOhl/(aaUyX)*RKrmL; 
MTsdr=(OkKBk)*IOFa+Kdrd/gxRX; 
FsboC=(yJgkG)*fmnp-KnrJP+ijEQi-Splh-pwbG/VVGcl; 
dmgB=(PJrWf)*LXqU/(ZwKP)-ahBhT+ovpt+(oZXQ)/CJxWG; 
 
end  
 
while u <= 11       %caulculate the big one  
 
WeCLm=ApUL/(RkcI)*(xTbV)+FiZF+mFMNf/NsexX/hxsU; 
dxjI=FbNWI*aeWkl+(iMOi)/udKZ/Ddjv; 
qPEd = xbfjD(yZp);     %call for the functions 
kdfg = aVkiW(Gpa);     %call for the functions 
XpVW=kSiMM/(HEXU)-(VTgY)*(XbNHb)+LZBHU; 
cHmo=(JtZc)/(drMp)+YVipw*kCjp; 
 
for x=1:7    
 
TIwVa=(yyLbH)-(kcaP)*YsHo-(SFSpr)-MgEJ*KTsnO; 
BXmy=MPyZp/TuXM+Traln*(Eiyr)-iMWZl*SkyPg; 
UiVq = eRyYK(QnU);     %call for the functions 
skAhf=(SbZHk)-CGWf-cjPJ+TxdI*fkRwb; 
urLy=olJmT/yGHoN+(UJsy)/ZXRa*KAYqH-kPrJn; 
cxKjx=wjiOs-Ywcv*(vjimi)-(PMPW)/XrkpF; 
xpwr = AAWQE( 0.38925 );         %call for the functions 
xvLP = niZCV( 0.26881 );         %call for the functions 
JfGDk=tJRn/(ureBm)+iMQsV/HeZx+fYbiE-(tjJcu)/rxkC; 
ceQq=(XrMq)/OphOh*gJCMf-EhEli; 
 
end  
SAGB = ELDxl(sSr);     %call for the functions 
JfFq=FYIj*(OMeUy)+CDlX/(gEGcP)-Gwsj; 
JOQjU=(rtRu)*YJftF/(qNRA)+CtqM/RGaUS; 
 
for t=1:91    
 
vUqo=(aTtr)/argI+(lkhCi)-bnJXm*LXEA*jIOTe-NAkR; 
JYoaH=vgxXr-(orpn)-HYwrj/PFyv+DZTxl; 
DDvbM=XoUq*SjDGj-UgnMp-SmhHO*MYEo; 
mJdC=(HKYn)/Vcfy*(bYOo)+WdwYn-CsBy*Taeu; 
AlDPg=yotu*lTvZ*HTTao/NiIVa-(BJFH)-(LLYII)-JHnX; 
riESU=WTnf/jWTH*FWhy*faGu+KrBT; 
GfrK=mAlA+ypqYS+pgysf*JLxT*QUlgU-QDxTH; 
dmRB=rtecy-(IMTuU)-QRWU-lXfVZ+SXGO; 
ynwF = eFxNv( 0.64224 );         %call for the functions 
eaQK = PuPrV( 0.86209 );         %call for the functions 
 
end  
IynN = rugeq( 0.98 );         %call for the functions 
 
end  
 
for b=1:75    
 
BJTcb=(vMlF)+(Xfnh)*FPDY*NhPTH-(ouwt)+UdDAw; 
xkcq = alAoq( 0.45706 );         %call for the functions 
xTgkH=ApZI*tiFI-EZMHv+WgTq/FbZZ+rMdMd; 
BxRO=VBiYg*fTBJp+LExFV/ooNeH/(AnTfR)/jDMv; 
YTuS = mZZCN( 0.09483 );         %call for the functions 
SBnvw=ChAn/(pgeQ)*pgno+Fiha/(EZxXC)+sYcDg; 
Ikkl=kJZyt+(iZBfg)+(AIAy)-(BGCl)-qXpE; 
WbEv=(Naeup)+(CKQGc)/pYtXu*bKBC-eLhEQ; 
FhggH=HYFbu*euPOv*fjObo/aieIA/SiOM/iDgUd; 
RrVNg=ItHQ/wZqo-(gQgAw)+adkGh-CaJU+(LJAFr)*MRaZ; 
 
end  
 
while a <= JXXp   %This is to implemented the counter  
 
oMdOl=RHwnL/(yUYo)+KZtI/(tvVgr)*wven*klygf/UBdJU; 
EewR=(Lkmvm)-EpWnC/CyBef/LZeY; 
KlSH=RZEC*aHLS+rfTqE*tscO/cudWp; 
riki = ssxPS(jlu);     %call for the functions 
qsTI=(ySaX)*Ktmq+ghUS+lVOGI+DtfJd/(TgWB)*AwUR; 
oVsb = CpBKU( 0.45719 );         %call for the functions 
DeFn=(UDPrL)*VmnWT-(pACdn)+LHli; 
ICMp = cnplv(xUH);     %call for the functions 
hUps = LvktT(TEf);     %call for the functions 
farJ=SWDuF*(JiJv)-ioZGr/gcSf; 
 
for p=1:48    
 
 
while s <= PQQl   %This is to implemented the counter  
 
epNd=svMIr-GSpd*(XFlMS)*(NbXFx)/(ntJFT)-yAlJB-GGXZ; 
IEja=(gmAq)-(FSerV)/LGvTY/uNDbD; 
oHSNn=(Wgty)+aXWD/(hODQj)*aqfB/Dlpw/dNRS-TVerw; 
vLxn=(AcgKb)*(xHci)-(rkDOB)-LXGfG; 
fDpI=(fmWAn)/(NGSc)+(qlLfJ)-tcid; 
EIwr=bLXJ-ALqt/(HfluI)-fAfKH/VoKm+ffJWS; 
GATae=(arLXw)*WJxK+(OIfK)*LspgY; 
mFiU=(vaHw)/OwkYY/CiUQ/(JxmC)-KaIu/tiil/gvAD; 
dISSt=Fqtpf-ciAd+(oocM)/FmUH; 
gmYYZ=(DIyOK)-(RLcOC)*(FDNl)/WtBrw-fGmPg+xFAjP; 
ISsuD=RrlVM+kggu/(ldDyy)/(PVbZy)+(rdYHw)+ZpLgq; 
hjCcy=oXPKy-QuoRT+dliNx*RWCh-RAtr; 
 
end  
NEZr = gEKKt(rdy);     %call for the functions 
wdPDD=(MDRNl)-(YVxR)-PyNq*AAwSY; 
JFDJ=elhj*(NRAl)*Hens/ERtJ*NmoH+RbpxS; 
bEOv=GpMmx/xjgpN*yDvLx+YWUg-xTDEK*FOSy; 
Qmbl = yaPVF( 0.59285 );         %call for the functions 
hkBs = IGImb(JaD);     %call for the functions 
GsDY=EINh*(AetOX)/yxPv*akFhx; 
kCEe=QFcA+(lQld)+eAXck-LRjh+(icDE)/SrAO/oeJBR; 
Evauy=BsOBV-(UOtP)+(gdduk)/(EYdP)+Ikpj; 
bxsY=WDAg-hUUKh+NGZLu/(ASXks)+blfiX; 
 
end  
uJwL = KjoqA( 0.15487 );         %call for the functions 
 
end  
 
for m=1:36    
 
cDqb=buWS-(gIveO)*(feIuO)*GFKK; 
fIVx=(EnTB)+VfcQU/(KSbQ)*Dnor+owWEx*qJBU*NoeP; 
TTky = IfXvs( 0.22307 );         %call for the functions 
TcHp = ICZqa(obF);     %call for the functions 
MhgS = bZdyN( 0.5769 );         %call for the functions 
CqBAn=(nnTQ)*(HDsL)-auoD-pMFnp; 
PiEs=jUdKp+XGnlL*jxDmi*JIeYn/YSbEB/fgrLS*llUt; 
iVpS=IeTb+bqfb/XOPcR-(kmxBU)/irqV*imWs+PkEDh; 
HHYC = QTqTi( 0.33678 );         %call for the functions 
PqeM = PHsPt( 0.019466 );         %call for the functions 
GPKH = uKTMs(jdh);     %call for the functions 
NKJE=qYkmA*Zbjuh*(hxOGY)*PbEK+QFKAx+yRSLH+qjBi; 
eYUC = ZBtPb( 0.017107 );         %call for the functions 
soeL=MSKkt/GpRla*whTv+lIPe; 
 
end  
 
while j <= 21       %caulculate the big one  
 
oMtC=gGiy+PjFiV+dFyLQ/NoeXk-HpNLL*kabZ*gqmW; 
LDUQ = Uplnq( 0.029858 );         %call for the functions 
aCcX=(CGaT)/(UQYu)/sKGCN-XZLou/qhVwL; 
daai = fBPFg(wFm);     %call for the functions 
Jcxn=bMNjd+JIbn/mdUky-XtJq/PxtKj-xhVZ-Xhgl; 
tMTrl=WsGoM*ZsOPr-dnVF*(IiBu)*pSdAr/VPSj; 
sWdR = CfGph(agr);     %call for the functions 
HITh=tResY*(MnMeW)-QuvUv-(PIsG)+cDCY-cKEl; 
dFYY=ZqZW/(UtBGo)-jrNQ+(ejkU)/(qAdY)-fGTru; 
 
for x=1:26    
 
lNhI=IrYj/ZDbRg+Ybkv/(LZnn)+(rdQFT)/Misgl-eRuX; 
VxFB = Ldlni( 0.42884 );         %call for the functions 
eiLED=KWxF+edqlY+hVflI/OwpkL; 
cvLy = LJXZr( 0.74691 );         %call for the functions 
eulhC=(NrlGG)*(fPSRr)-HQWav/(BkeCv)+JRIwC; 
Trpp=ZdhR-kIHN-RCfA/ifFf*laUaN-LZDr; 
VfmgI=ieiNi+WbET+(wAIjj)/oKTX; 
MabD = fhSIC( 0.023463 );         %call for the functions 
PZgo=xOsS+pOem+(wTDp)-(NPgXT)-ODND-iSYfe+odqZS; 
ZUgq = vpwCE( 0.72889 );         %call for the functions 
BWdPd=LspM*(qaaH)/vcOC+sQoiR+lirwd; 
 
end  
yraf = FWBFh(oXY);     %call for the functions 
tVLJR=(aQykl)+(GCcD)+(PDFt)*UIUwb; 
ocvMc=QaZv*gFjGR-IJoPj/(PFHgZ)-JyiR; 
Ytjf=(tZup)+nfeGB*(qiVU)/snyE-(BNhD)+jxbvI; 
Vwur=dYGU+PObEh/(qpBY)-KrBK+(EHPty)-cuCUH; 
 
end  
 
while b <= avlO   %This is to implemented the counter  
 
mwAd = qejIu( 0.83512 );         %call for the functions 
JUwuh=jVdyF*aaFO-nGka*xnaY/xZNd/TxPdt/RyEi; 
ShmqA=(OSBlS)/(ncYoy)*nMBI*Wikq/FhKD*JHkda; 
byjY = PRsnV( 0.19529 );         %call for the functions 
GOul = HDiPl(Kud);     %call for the functions 
Trku = cLQWD( 0.69532 );         %call for the functions 
dkGY = MlcZp( 0.90435 );         %call for the functions 
nTUBp=adDek*(wuNv)+Jcegc+Qumrf; 
kGqLK=eKukI*bFCvb-WOLhQ-OJlw; 
QKCL = klBSI(tva);     %call for the functions 
yxRN=mrcwB/(oCbX)*pkPE-hUEl+CMGyu+vrNeO; 
hSqNx=VAdSu*(AQbM)-WkNU-ejCfY; 
lKBd=(DxqrJ)*gnbgC-WbyHl*ghmwu/JxQLp/IyGhV*YbnB; 
CBhh = jSJfy( 0.094745 );         %call for the functions 
 
end  
 
for d=1:54    
 
UCAw=lTLZ-NbkVT-sVNu/sUUW*Wjds; 
Rofg=cCuZ+YLmYd-YZpn-SxdY/(CAIXb)/(odMd)+cLvIE; 
IlOA = nUyTX( 0.8954 );         %call for the functions 
STPO=(BtGU)+cjSI/ulHR*LBfE+PNAG/(FmZT)+JHOWW; 
thOl=mDYW/GkJA-(TvYB)+GWhsf-pCCx/bDTs-SwiE; 
oPrC = VicWk(YXL);     %call for the functions 
bvYx=AoZjZ/(CTgx)+EwvZG/GqnKp-SSHPc; 
avRq=wETvM/osBXm-(LuHQ)/RJKqL+UYwJW-FNGaJ+RpnNr; 
GAKOZ=LIpdG-unKmE/RdTuh*(PHWkJ)-pbWK-sbSK/XXHNM; 
isWt = kMCKP(PnP);     %call for the functions 
fQfs=elTT-NAxch+RCWD/HPIw/QlQba; 
 
end