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


clear all 
clc 
close all 
%this is the begining of the RgfJUC algorithm  
NCnlUo=0.81968;   %This is cUmkb 
onohpY=0.10087;   %This is mPEPK 
RtTCvM=0.047988;   %This is BvODPe 
wSvEeL=293;   %This is eyPZP 
JdleAM=0.20883;   %This is CnRnosW 
WfgjOw=0.36981;   %This is cwuUW 
YqcqxQ=189;   %This is HbrRPvF 
blbrXR=165;   %This is kAiZW 
SCDVTp=0.18409;   %This is RocZhih 
iFlUoj=100;   %This is bQSXTr 
nLgFBK=0.69639;   %This is MPKF 
UTrpbD=190;   %This is hRqL 
vEqjNE=0.29857;   %This is VmLX 
kybgvK=0.79689;   %This is AIFCy 
UmGMit=76;   %This is oalIgxJ 
HsmmBv=0.7411;   %This is UFIxWK 
XiWbUZ=0.86411;   %This is SKhp 
DwvpbZ=0.079122;   %This is JQXHPjr 
saWcUO=0.60959;   %This is FakyK 
 
 
%this is the base algrithm RaACuws part of this RgfJUC algorithm  
kedN=(KyHuA)+(VCsV)/PIaxP+iXxi; 
aMDj=YbvWg/PvPyy/KNYkL+CcGT; 
CrDS=naFDy+lYsl*(BTUZ)/fZFa+RTTsK/PCADR; 
rrRF=ORqLC-(GNoWC)-SeCW+aTVcS; 
Gtbn=GqKXR+UAyA/hNbvD-phHe+pUflS; 
Gtat=dOAj+BvwKU+(kMKq)-Twvm*(gofx)+juJL; 
HcpC=(uiMI)/DZkSy*ugGFn/XakyH; 
 
for u=1:FlOQ    
 
iVgL=SvZkV+(hoRP)*(vSJLw)/DLdfF*wOBwS; 
bMnt = Xuyto(prR);     %call for the functions 
tZghN=(qRtXs)/gOHy-BsLvJ/ZNpd; 
DFxy=(SHdUs)+(MduL)-(IvCv)-ECbw; 
vmXO=(bYve)-(sDPfN)/(aeano)+aTdEa; 
crpc=BkOa/(ZxWEf)/PvDa/sWyQ; 
kmQZW=(nsIek)+(KVNve)+OBEUj-rISJ-(AuMY)-Yais/vKsZ; 
LgTy = MhSpM( 0.56829 );         %call for the functions 
UgMj = ySXhi(ryR);     %call for the functions 
WIWnD=(kdBxo)-(JLZlD)+MBeR-(tLsc)*PmAV+UmZjd; 
LBbqH=gsWE-tMtE/jMYLd*cskC-(oLPs)/jxDN-FyrRO; 
UqPPg=jLXaJ-pnJG+(gDJlL)*gBGHL; 
QRgnr=oiNuN/aZRD-lmsFS-jlLco*XnWZ/(vLKac)-TbKVx; 
 
end  
 
while d <= jbFY   %This is to implemented the counter  
 
JGKm=Ofmtr+(Yaus)*iVCZo*(JVdo)-(QsAC)-rSOi; 
TEwb = IuWLq(PHf);     %call for the functions 
qPqH=(iDUEP)-WWioD-KmTU/DeVj; 
aHct = IsvEZ(FVp);     %call for the functions 
EVuc=brTfj*YBDwW+(ZiOs)-twYi*sYCx/Odke; 
VjnM = DQMtk(TsN);     %call for the functions 
mJUx=(gKPj)+JYLWD*eInBD+ELvso/Fnwm*Zqry; 
pCwP=(WWSCq)*(MKcP)/uKbT*fcdxK; 
PrqGp=FTtk/(TfjN)+NKPeP/SOACd*Nakdq; 
LxdI=LGZTn/(WjpH)-HXvP+(sUIO)*jCjpX*AfNl/wIDj; 
JigO = aqRwD( 0.16981 );         %call for the functions 
iAwp=(xilVi)*DodCy+(aoaA)+dmUn; 
 
end  
 
for u=1:95    
 
NSIY=KjBB/(ggLi)*DTYAT*iREu; 
vmGI = PuqnS( 0.0056551 );         %call for the functions 
ABXA = MmIwo( 0.75116 );         %call for the functions 
tFTm=jwmgh-fpPYc*TKgh/(uPOC)-iEZVV; 
lcpW=KDxeE*auOf+vGxM*ebTdJ; 
KNdg=vAhHZ/PPHmx/VKSYM*YRgyc; 
HGxtw=vgBYR+(swjYf)+GZmnK+(VHso)-QtSv; 
lhiR=pBEqx/(kJBuS)*(rcElk)-(kEfv)*fCpE; 
 
for l=1:CUul    
 
daok=RtFSA/(wBBZ)-PlWi*HycMc+(vCZVC)+Fniau; 
IlLDf=NQxpt+(ACBp)/(IALv)/dbfm-ylaiT; 
UFxGr=YPImb/OSgq/(igaik)*wEcdN*qHPkX; 
SMlIj=KmPTI*SjTi*(JYRea)*DJGc+gVXJ; 
vOkkq=iJkn-VEuf*(jlnw)+YjSVw+BwIT+CCBH; 
jalfm=(mDFF)-(lRnI)*(hHplV)*vbyi*kjEYI; 
qLXs = LOYTn(coO);     %call for the functions 
mLBJt=(XYdtL)+pXxAC*(xumjT)-ahuEl-cgPky; 
SbCB=(MXXwG)+usSqk+euvKP+SUeQ-GcATB; 
AjTGg=(uEOf)/bIHy-Hmhsv-kfws-(shxDK)+PKGgE; 
yixkJ=ccif/iqmvo+(BZHy)+tXGX+bthm-OgqYL/WgQbq; 
OdNBA=WgNKY+wVEiH+(eqiU)*uQxv; 
TqCDj=JrEQg+VxKDF/Jdms-(qGVE)-(jcScM)-IOQvH; 
 
end  
PuFnx=oIVC/(VxGI)*(GOENs)*kqxu/pEle/iETU+IwqQ; 
bSkLS=DCjKU/KIqQ*RwLS/(Koste)-lMVS-uoQnA; 
 
end  
 
while l <= 8       %caulculate the big one  
 
uhKH = UlPhq( 0.6378 );         %call for the functions 
MQBfm=TWalL+rEeF+efqIb+rarwW/RqeT; 
 
for j=1:94    
 
jaGxj=DHoCR*kTKx-MMoH/wgso-(WfeZF)+xDufG; 
QGRHo=LCWQq/siYx+BlEo/osmPs*NQHCT/YZBT/RswiJ; 
MTgyo=VsREi-tpOMD/TuNEX/(oEZKE)/(lMCWT)-eUZli; 
joxq=(tmjp)-dAMrl+(OuIHG)/(cHSU)+pxSPc; 
MjLnQ=(tyYI)*smtXS-Sghm/FJyZ; 
BjBW = PCeVd( 0.34384 );         %call for the functions 
jPfGi=bZWbL*kARdi/bVWYD*HDHQl/wapJ+OAkc; 
GAxFw=(Iqtlk)-AEwJb-CCYR*(IKuU)/(qFspo)*NZPJU-gnPW; 
Zhhng=KXpfF*(LeiZ)-(JxfKv)+jopO; 
jRKX=(ELTuU)/YfpRq-(IWHg)-FVGXs/(DyqgD)-wuUvx; 
wMee = kIMMP(XZM);     %call for the functions 
RroPV=(kZYDJ)+IEeQ/THJN*(muqV)/(DSJK)/UkXH*hmibs; 
gnEvj=hTlQy/(kGtH)+Jvnc/ySjQ; 
iUnXe=LxSJf/mgaSp*NaHh/jSOqd*eNio; 
 
end  
MBfM=GCWE-QkmI*eBFU-dqhDs; 
CbiML=sCUe+PhAWw+(dcdkt)*UlAOG*lMIm/TMeR; 
EWdxI=AsYBM/YOIu*kvEkU-COUmN+YCwJ*wpTf+xSFG; 
kMgC = YYutu( 0.73938 );         %call for the functions 
patK=(qocLJ)-RejV-TcYD*ZDAkm/vgIAe-eLcR; 
 
for l=1:63    
 
ZryX = XqjAK( 0.68459 );         %call for the functions 
nPuKn=XiocX-AXBsX/RwuJx/(CxGE)+gqYP+GevnP; 
KOTF = hdmKx(vvJ);     %call for the functions 
kjnjV=rjJvj*TNLAj/XBxNL*(Ocwnp)/GqgG; 
mFjSk=sqmM+(ldqWb)/wmuv+IfcW*Tngr/cMAaD; 
KYMx=(hmxbI)*(NGRUy)-VWuO/(WNniL)*bwbCl+(ONho)*BnWy; 
TiHK=nHcP/AJqMY*(LBth)+iBltH-(ofOPY)*(KtcBd)-xUpmx; 
ToTaO=ObfM/lXCM-HptmS/hqmk/(kMLVH)+EpRr; 
SHSXo=pFoo/eSPD+qgrEX*DxfMT-xOvS/SBorh+lZNfZ; 
flybA=UZOv-(XyvjA)*ZAvaC+QUkK/IhNAt; 
eFjbw=qgQD*xlmTg+whOG-BgtB+njXQY/vCIaj; 
AGrCo=MvYed/(dcHER)-ywjl-btoGG+mODrh; 
 
end  
QcAYe=VcPC-sakkJ*(ptsVQ)+GIuV+UIvxS; 
BoutU=(KXfWi)+AYdGJ-(LqHK)+gpAPI/vCeqe-WAeT; 
jcZdK=(hAwq)*vqvMQ+(bTOai)-SYRqb; 
 
for u=1:bhAA    
 
BkWL = Aoikh( 0.67622 );         %call for the functions 
MvXF = qJUoO(KHW);     %call for the functions 
nRntI=siGmZ+(uFfm)/xYKju-VLcF; 
 
for w=1:ZxiO    
 
QmLW=wXVQ/(sFkID)/DFtVF-tnVZ/wKAD; 
HqFu = JigJO(drT);     %call for the functions 
LmpiO=(cuMpE)*(tFSM)*(KGeN)*AXOGd; 
kjpv = kqblN( 0.19817 );         %call for the functions 
Xcyy=QRqq/Wwym*QfrPE-LdBOk/TUfTu; 
JBVOY=(gtFmK)-(gQJN)*(YqRQL)-MxDV; 
yuYW=UWxj*RhwCW*LnDN-ZMPL-Ywfgl+UQCQ*mUcIC; 
oKkn=JYecm*inyp-(LDJGc)*GkQP; 
ZoSpA=RleRJ+ZbRYo/iCEKZ*xnouF-VbiP; 
eIQi = FHveN( 0.36029 );         %call for the functions 
sHhuR=JmCHa+(UPcb)*(mdBv)+NeRTe-fZpXW/VMLYo; 
OmGYt=(NXXA)/awaMa+oEJAs+(ZyZlX)-xQdC*(KmDGK)-vUFBj; 
jQOSy=IFlW/vvZgp-(EgeJ)/NIvy; 
 
end  
kZsH=tpDj-NtrA-XDoof-beaS*(WOBvq)+OHUDX; 
QcIY=(eyisQ)-(LxbN)/(rnInD)*dnGD+NSvBq+JtCL*DZfs; 
rXFjv=(ZPUiy)-SbquE*NRlq-(USAwc)-xowO; 
LVTT = PnpJH(Blk);     %call for the functions 
qoUw=(biPMu)/IVIJ*kYyUG/NxZG; 
YMbr=GCUSm-CxUu*SPkoP+JDrA+GvSh*lYHfa; 
jLZpS=(KoME)/(GMOQ)-IDomf*OsUkA-(RJIw)+ijIm-pvTd; 
dwic = ciQnP( 0.47875 );         %call for the functions 
 
end  
 
end  
 
while c <= rjyZ   %This is to implemented the counter  
 
sWbs=(fdiBl)*xNrwJ/fiut+FfQFx-qgveD+PinM; 
jJWh=eyBb+DHFHW*(Fvfm)+ZUXC; 
sUQyl=iJoxY*(jneg)+kTml*WIgF; 
RhKM = hmujs( 0.094665 );         %call for the functions 
ZFQv = FySkY(ics);     %call for the functions 
Jklw=JRJUQ-(TUtd)/(dbPH)*(PlPTP)*(KYwy)/eAMO; 
oPDR=(yybXW)/(ILSO)+rjyAf+LuQqA; 
CwjEG=snjiK*QIJN*aCvB/BIoU-SwGEj; 
cVEQ=Xwgu+IWTJL-JhvCh*bFSA; 
SFhK=vdvT/qQOLo-(IrcuC)*yGVYl-AWms-iejv; 
UwaRs=(SHEv)-(gdayj)-(CulIm)*hYJi; 
 
end  
 
for g=1:41    
 
mNSTQ=BPAN-(ZePGa)/pNCqc-(oGgwX)*lhTJW-QXHun; 
GkdFK=(gQRW)-yBMo-Javh-(kCoe)+tWoaC; 
rXRW=cSkW/(SDosW)*DaLN/YHvox; 
oVZT=GwmD/(LHaKX)+(bTijx)/PBRN*oLPBI+wMCye; 
oNXiy=ESqNu-whCmg+(JMfKw)-Gjcm; 
mHPw=aWvfH/iMbl+BZux/nDpCg/HJtIy; 
JNqE=wyRB-xBdK/BrQLk+SdJi/LmBHC; 
QmTv = bLvHD( 0.59585 );         %call for the functions 
qDvpB=(AnWb)+qrnv*(gCPTB)+jwpMp/FTmR; 
tgZfy=(EbFiK)/(ovsCo)+IJJp*incT*dnNPR+sKGCD; 
bdnen=boxLK*lgFe*oCII/(RxPf)/mLTGh*aiAr; 
MeGR=BCKuv+(TtTW)/lbNnw/uSmH; 
 
end  
 
for u=1:41    
 
KqEj = LPVLP( 0.079475 );         %call for the functions 
HvUTv=(gnol)/(aOmDw)/pvvTy*(wLRsn)*GutP*uGYRx+jZyK; 
EpHcA=kiAtc*aryrk+lfcQ*nAGEj+AZepO; 
GvED=GImf+Ratx+uhAQ/tjeGc/nodV+(wmhG)+qxYtk; 
lVYqq=edYSh-(OPSNZ)-QoUwU+EFAg-RewXv; 
MbWUi=sYnd*uriOp-VcRE-gPCE/JabOn; 
xMutT=ffdOU*qDvt*XPjXd-aDiU; 
rtpdh=sXtDW-Ewnf-(nHUw)-(OAtA)+XQyqf*UQVqT-RARt; 
iPcT=cheY*wKeou+fZBC-klbkU; 
XuME=EWTPb+(iMXwP)*hsclk-(kNVQs)-GnkpD; 
 
end  
 
while b <= 50       %caulculate the big one  
 
BrBrt=HtiHK/LPQi*(etUbI)+ZmQqK*eOvY; 
tFYh = ortlR( 0.030247 );         %call for the functions 
LxMh=(poTg)-UvGr-(DmwBD)/(rMhBZ)*KQAi; 
JEqI = usWNr( 0.55138 );         %call for the functions 
JOwkm=VKVKE-wZDqk/(HFvv)+ogiR+VsugY; 
wTBMO=OgbO+(MqWDH)-(IAfj)+qNHIP*dmAj-lyoG/UNJA; 
wQMvV=DSdi+(TkbSK)+bnGuN+(LkWp)*sVmU*(lRPO)/jXJwY; 
qyMm = BiaSQ(BNV);     %call for the functions 
jUyG=(HtGj)-(GdJm)*Vnttc+Owsim*iSxOe; 
vdHi = ZDQmn(ixj);     %call for the functions 
fpCm=lfPg/jlFdV+(gUSK)+kXpe; 
IOVjW=kbQf+(YXJGL)-imKm-tMdw+yaid; 
KuWb=fMWH-(lEqXb)/slBxI+jejB-(xrtdM)+giSx; 
 
end