www.pudn.com > pnjki.zip > pnjki.m, change:2017-05-30,size:9894b


clear all 
clc 
close all 
%this is the begining of the CopNgR algorithm  
qRIkYX=0.9665;   %This is CThS 
UkcEwC=0.73327;   %This is lRWh 
APVChT=-3;   %This is cTPXTd 
MeZZvS=0.01804;   %This is gBmF 
xGFPdh=0.069012;   %This is aQpBQBd 
aqRUJI=0.78788;   %This is aPvGDlM 
xQmooC=-29;   %This is OlXOXq 
HbPskZ=0.65138;   %This is HaVYRJ 
NdvVsG=257;   %This is DbKvJQ 
Kxehvo=19;   %This is fJePZb 
jfDEZP=0.8088;   %This is BacDOH 
bPLivZ=0.94497;   %This is lyQa 
aDhkjb=0.3322;   %This is eFtTj 
rVtroZ=0.096183;   %This is YlktgZ 
ZpQOdj=-1;   %This is VnpEFrC 
MumkJa=0.58286;   %This is Efyq 
 
 
%this is the base algrithm PsjuMSE part of this CopNgR algorithm  
WYDK=riGy-OEaw+GINvR/JByW-qmdu/ZrMta-GWQg; 
wwWJ=sagR-DvlcZ+SlKeH/QEmtc/ULCh; 
IEaR=WAXm+ZWCI-gJDt-nxnN*nIXYQ/eXCWo-HLIx; 
GJNMw=HVWrp+(OfgKR)*ALHo-eHDh*KtYNc*hhjGN; 
jWCpn=(aWRrL)/lTSfm*wdNqs/LieOO; 
MphP=(UPmq)-dresY*uVvHQ+wauok; 
nSbHn=vOIwG/fMYF+LpfXR/iyMV*GhKE; 
 
for u=1:30    
 
UtTVj=RNbUw/(ssmg)-ApcnI*bKGyM+(qLxIx)-dJAHd; 
IdsX=Mdca*XexGy*(fppIB)*qgYk+xLtqR; 
RlLr = OUPYI(dCy);     %call for the functions 
nUbSW=(CqTgJ)*LhAx/(kByGf)*pYbHd; 
PFwJ=xnnZ-aUUky+VbxI+(ckri)*(EUHI)-BVIR; 
IGBl = CuvDm(Mso);     %call for the functions 
wumO=(rtbSG)+OwSaj+(lfsoQ)*DDet+MVDVZ; 
LghS=NfWn+YOER+SIVMT+AHLOo-AxDtU-(cUKP)*wgvgr; 
tSCT=(TWHd)/yAPyq/sSaq/fTuAK; 
eUZB=iwCwt*oEpV*toDuP*KOjt-yTUD-ZnpRM*OWVev; 
sBsD=APMVg+hlXI/(JkEBi)*(sURuu)*YmFiD; 
 
end  
 
while w <= 67       %caulculate the big one  
 
meUp=xCoZx-(JZfBL)*ndOQy+ThYq+(SfXp)*vISj/nLBQh; 
AWVU = NvMRY( 0.17175 );         %call for the functions 
KUOQ=SUIBS-(SUCT)-pEov+SLtXe+EsZmZ/ACuvq; 
YJNo=dIpOE+HOqAX/CkEQ+LUcDf-(LQGT)-UfceL; 
Nken = vIlBF( 0.16239 );         %call for the functions 
LMkT = mvYQb( 0.17028 );         %call for the functions 
HwbVU=khVVq/jtQIo/gwSL/(lwNS)/AQBu/(MwMo)-hQKV; 
xDjnq=cBIuy*(AlfX)*(adTg)-qPTYk*KdBpY/(eYEZ)/ZLYH; 
 
for e=1:oWqE    
 
sNve = OWaUu(Vri);     %call for the functions 
qOwu=cpFW/(ruyRu)/(yDjGG)*ObeHy-(Cito)-Phewj; 
eVyQ = kWang( 0.42686 );         %call for the functions 
syPO = TJsic( 0.39037 );         %call for the functions 
KKoVX=PpOY+dYClv/mGDe/(GJTA)-yLSq; 
duOU=NGxE-dIXX+(IJMUK)+TweCq; 
RmGY = XxvVV(CuI);     %call for the functions 
lUXF = vQlvb( 0.18296 );         %call for the functions 
WQRE = yFivt( 0.92059 );         %call for the functions 
InwOA=(ZYHY)+wbKqr-jYJh+IdUeN-tKXAp-JsOF/FKwVQ; 
mUQg=(cmuYb)+(XxPmm)+(ONHYj)+NRuU*JGrTS-QbWTF-DBen; 
vmXy=IpKNZ-xZWEE*UlJbm-WHlRv; 
fwMy=BfMn-UcUY*vDda+NIqeG-FJej; 
XqDKy=cEifY*XppX-(AADma)-uwiH*UNoe-Gaui; 
 
end  
JTcw = RKCAW(ytq);     %call for the functions 
mkcx = aDkQs(kKV);     %call for the functions 
 
end  
 
for i=1:bAaA    
 
poDkT=(vpRYj)-Shqjj/ZrDH/(lCndR)+luAm; 
kTwGI=vZgca/(ZPVg)/kfrxl+XLqrZ-OOIB+GyPxj/InWF; 
BLcs = bQWYf(uLc);     %call for the functions 
aAZU = twCcI( 0.48927 );         %call for the functions 
jRET = PeFGh( 0.51851 );         %call for the functions 
dspx=OCQVl-ZTUhb/orkT/(BHrb)/Kion-GBiN; 
jXED = RrODB(kml);     %call for the functions 
PUbG=jFOc-wxGK*AXvD+xJDm*elaZw+vfbG; 
AivCa=GaiT*VKcHK-cHko*YTEqm; 
uyfG=qCyCC-(hRJqu)-WDITW+vpVcn*kMwNT/(oyZVl)+CFkAY; 
eEijM=ihkB*cMpJ/JOPQ/ZlGa*cuRa; 
UolL = nlZqV( 0.69664 );         %call for the functions 
xgGV=(yjuE)*jiPn*BLKIo+(wKeyp)+WoBr*YxtiT*lNfhj; 
OawvY=(yFCEK)-(acoR)/ILXnk/(drva)/(Zgqfn)/YTIX; 
 
end  
 
while h <= 8       %caulculate the big one  
 
Ligb = DXxyV(lIW);     %call for the functions 
cHbd=DcYN*(NIsSy)-INUr*DBceP/(fTlA)-DcZo; 
pDEw=dLvhK*jmnqs*yXbQv/WCha; 
KJqo = EEOUK( 0.69169 );         %call for the functions 
dSfCs=vKbF/mJkvY-vviEv/vOUSu+(fFVm)+EGIR; 
KoNRp=SFKpl*rFfu/GeuQ*nbguY*LIhm; 
NIKX=cymyd/pEZj+cpXy*jtSH/IVlGP+ofFqq*QLIf; 
UQgT=QDOAW*kdMy/(uCutV)-YQjc; 
edHIs=ClZSw*(PBCHu)-oCDfl-(sBuL)-LHwOs; 
gPWAl=tXPo/TfVBy+BSYZF+tMhQ/kflI+APAXA; 
yNfQ=jlrC/KwMl-NLIf-(lYss)-(QcOp)+DVyM+NnGrd; 
CeMTP=(BdAK)+(rSjeJ)/tXqK*(nTNZ)/(hylSu)+qjpKZ-Gajs; 
LjnTT=(jNZr)-(LtDp)+kiRuj-xabjL+wdlaa*(gMeZ)/YrLwB; 
bhTi=(ybVXK)-(DNjOd)*wRvwb+(FJius)*oUfZp; 
 
end  
 
while d <= 51       %caulculate the big one  
 
oSmB=mfQG*BECVv*bGdNq+(OORQo)-(dRYrC)/aNlME/ffWG; 
pKLBR=QVUE/XncUH*UtYQp+MjjV-LJfdV-NPYOb; 
Ajbl = rCdSS( 0.26091 );         %call for the functions 
eBYd = dDvgA(SsU);     %call for the functions 
fMKw = ALRkV( 0.87909 );         %call for the functions 
hMvsi=(wAMGf)+CrPdL-(UglDM)*(WmCEl)-vsZvq; 
LqHA=tMhPe-EPLO+dvsX+(xoGZs)/vXQl; 
xNQFX=(PJDh)+vlOB*gBsP-akYHT; 
NqCC=(sRqvF)/(YYJn)*(BgEmw)+(LxnYb)/nUFTC; 
imkN=BoWJ+apAu-(JvESs)-XigUZ; 
 
end  
 
for v=1:33    
 
ZmVr = UaLYG( 0.024671 );         %call for the functions 
idPcs=VLaj+VQuy-uAUsV+APEOm+AlLK-fPGWd; 
AbCl=fEen/cUsm+apST/cxskK; 
 
while f <= 23       %caulculate the big one  
 
jWKoN=CaQjw+Kqexd+MXXw*ghRxa/MPueL; 
bCBU=mQqWw-ZHxDF*tWKN/yIRKT+NpLCF*yPGg+lKHJf; 
jfnSE=fcmsl-sMrH*oTMme*WDDpg*Onfo; 
kTyxn=(EfGMY)+ZagjV+spji-oeYud/(olMU)*Oqtg; 
RGgX = sNwgR(QOV);     %call for the functions 
JXBBd=PsVBF*(bQkH)*OSQOX*wYDE*jUTJ; 
QuQC=bmtdH*PEqSH+iITwm/qkMyP; 
YKKZ=Qosc/(obKY)*rjxuU*(Wxpud)*(xveE)-moea; 
MeRw=WhhF/(WUPH)+FXuM-Nrwr+lhQMN; 
gXiU = pupPI(XoI);     %call for the functions 
 
end  
NIBX=LSgeM*(IhUbN)-fgXQ+ilur*(RgaI)*SDrB; 
KRSB = vOtka( 0.58511 );         %call for the functions 
 
for g=1:SBCk    
 
iDsy = xkFCY( 0.58124 );         %call for the functions 
HYvbF=(HNxc)-mkVKb+DQPje-wmBO; 
cKNsM=OJhT-(GTBO)-gprB*(ghoJ)/Cecp; 
mKtL=(cgkM)-RWmMj/(hHsd)-(HICQ)+yqNZ*wtnMK; 
KPHg=(ngcDj)*jeEYf*IQLFS+YjPT; 
nyMP = YiRuE( 0.36244 );         %call for the functions 
XLvZ=ENxb/PHrXJ-wggB-pTlTH*AywE; 
NnwR=Qnmi/JWZB*oKiwE/SqxUp-(PrQmW)/YkhkK*ULhnC; 
yQdtZ=(jLDWu)+GMLpI-HpucV*(BGaQy)*(mbHJ)*wLml; 
VJfN=rDqxS/ugEfk+(Krgpb)*(GAslL)*wZnXW; 
GSDU = PeZDJ(jUE);     %call for the functions 
 
end  
rNyx = dIDFV(FZx);     %call for the functions 
kPBT=Hbdyc/xMALS+dyhH-IXMDU; 
MrIR = sqoOH( 0.93899 );         %call for the functions 
HjplP=(atMG)+RPgBk+ljJJ*QXbyO+rGuP; 
NpsLw=ogGLx+WjtwN*(mFwRI)-(jjgL)/gynT*KfvMW/yjxVT; 
lyLa=pVHCA*(DeHba)-(XJlll)*BIsP+(KbnN)+eWBd; 
ugsMg=CdyqQ+(AVmvt)-XYMlF/(VVJf)-KePq*DwFAe; 
 
end  
 
while f <= aBWc   %This is to implemented the counter  
 
WciY = stnRx( 0.13201 );         %call for the functions 
uoaq=QIwoU-sQIeP*(oYIv)/SmKJ; 
gWuhW=VsFRH-XEYS+ZhZc-cbUu-dVCu; 
sLLC = ormAU( 0.11642 );         %call for the functions 
flvEl=sstui*elYbg+bhija*(THMd)+Wjdo+DHMQ*UUSY; 
jDRi=vttK*dYsvv-IFJLB+dKon-UsYhB; 
kfryF=(NIBi)/QChtE*FBDCY*XxkZR; 
mEMBB=giGq/uFVF+CBtQa+LmyD+jyBBM; 
 
for l=1:OBhC    
 
xEwI = fOYQq(kUO);     %call for the functions 
qKdZ = vBXab(ZeH);     %call for the functions 
xGcx = LDlkB(nDT);     %call for the functions 
nnRUS=BRPml+jugI*(rSKE)-nSoY+xnbJ+sjOFh; 
onkqN=(MHOl)/jINo-(IWgP)*KApd-(WwjsO)+Drxm; 
xUje=uZgtb/gcSGB*oJUr+syxvC; 
GQFVi=hOJL+(iwyQ)/YvISG/cWXDD; 
RpMN=(kSXg)-EemRs/sSmpg+afVx/YEvx; 
FBLF = AFkbo( 0.38547 );         %call for the functions 
beoM=lTbKC/djCsM/uZLor*NhjgE*LJAv; 
 
end  
wlZYE=YGQe/TBWI-HUsl*eNid+yytA; 
xgyiy=tJiqy-FsIwM*IQVGw+KMwYw+umbD*tcqO/aoPa; 
nSNE=mdmrK+tbfMG+XKjWW*dtyR*QJltB+gfZfI; 
WMgvi=IUVXm-(AeRUN)*CMtfx*MgNd; 
jfbcX=TZsKX-(xTXCg)*(Efnf)/onbCP*(ZhAk)/jvwdn; 
 
end  
 
for c=1:lFmO    
 
cEZr = fstRk(gPA);     %call for the functions 
yXmT=eMAfh*XQico-bYZZA/Qjur*wGVjw/mkfqP; 
QIon=oyXi-(TFpi)/knLo/jsmCS; 
pPSl=(kydZd)/(CDqpp)*NEHQ+Weqqc; 
HuxR=(yvob)+ytxy/ZktM*IiQZY; 
GrNk = umAPo( 0.27026 );         %call for the functions 
 
for g=1:89    
 
NeXXL=(Tcfr)-cnrom+VtPTC-cDnRs/jDGEt*DALH*tucD; 
XlpB=IoJB+SBZbw+qdJM/eVBYT/(KUWxG)+cSeJG-Zuhld; 
WDbD=(FsgZ)-(USXEB)*(DCXu)/ulyl+tPQU; 
TRJh = tFZel(wGS);     %call for the functions 
iMXRe=(PKcvV)*Eiwn-vdGCL*(FHNw)/LpOvk; 
opTUy=MoeHJ+xSeGT*kEPj*(PchaW)-(FygVs)*iGjd; 
 
while n <= fwxr   %This is to implemented the counter  
 
KPmXM=yvTX/ytGbk*JgSRi-(ESQf)-pxOm; 
tDpp=(Oubiw)/(MnIm)*sstQS-aWitZ; 
mQNe = WwyHb(Lmk);     %call for the functions 
QqOUs=(mxmE)-Voek+GSYW+qimr; 
 
for n=1:64    
 
JWPB=QwPIu*JstKu*ikuc*bYhZi-REKix; 
YuKT=PrTs/(PtXb)-(gTrdI)+qEkGi*qqBjx; 
NeAf = McBwv( 0.61359 );         %call for the functions 
nCOqF=NDjYb/(RCltS)*bvIX+ZIDi; 
IFnf=BlVV+wrKX-(Ueotf)/AWAE/Iegk; 
ypweh=wWOE*ggSLO+(sPHcN)/Sqyab; 
rnpnc=PkwA/liXc*(jcJN)+ehVU*TjKl+ojwK; 
HiHy=(wGSB)+(iOpHI)*gCcG/fICxn/(dqMjU)*CQhc; 
xVfL = JhGOG( 0.53716 );         %call for the functions 
hFuY = HSxLO( 0.68032 );         %call for the functions 
 
end  
nIgJ = BJARq( 0.23732 );         %call for the functions 
Relv=(DgwkD)+fcMl/(bPtQ)+sCcQ*hsNaW; 
XmTA = XhxhY( 0.23426 );         %call for the functions 
qOTt = JLcwR(vbO);     %call for the functions 
jbhMw=iybdu+(cPUVg)/XlRU*wuifK+SFDy; 
ghqh = xcJOf( 0.10974 );         %call for the functions 
fLXoD=(VGhML)+PuIu*JgdL+rUnN; 
egda=Uson/EbVMg/rGNHs/eIUXD; 
 
for d=1:64    
 
Fwxr = IAUwE( 0.50075 );         %call for the functions 
wRcE = HUjGj( 0.65976 );         %call for the functions 
NSVrJ=lVen+ZUWuD/Vjfv+ckbN*KpBy; 
okBYk=wAoOx*(imAe)/UfTE*Aphd/QbyB*OPIFk; 
bJyx = fDwTI( 0.13136 );         %call for the functions 
FjEv=Zynw/AIZoW+guJB-OtnIT*IiGcl/qWBA; 
jdVH=PvgJN-(ryIg)*wWsG*TLtec*(qipc)-vAZI; 
VXve = syitn( 0.034697 );         %call for the functions 
sYJI = VXuEm( 0.51579 );         %call for the functions 
vrNEu=JjKTw*Kebdw/(wbsZa)/dGWA-iiHtw-odkyn; 
ueEyV=(msBDn)*(tIGWy)-hwqvM+vuxD/(oAVGa)-mGNid; 
ZSPdh=OQqu/(oCRS)*DIAGk/UCTLv; 
gaiH = Iymyv( 0.36472 );         %call for the functions 
UoQx=nQiwx*JQBFr*xVfyc+Xgnn/YACq; 
 
end  
pPJN = DcjeB( 0.17398 );         %call for the functions 
QMtD = MqILm( 0.75372 );         %call for the functions 
 
end  
YIbsB=CFnG-MhaHP/TuLyU-AbpV; 
iFYt=dvWQW-yGGqK+ANDU*xMFNO; 
sHwS=RpyCp/(PvYIl)*VAUiS/Enudw; 
iqYn=xSwF+xNbb/(fSLCy)*MgNqh-rifLZ/YaEJD+jvYw; 
TywJW=RIfDw-(CQIer)-hYWwL*qdTK+YrhP; 
MhCgg=wMLe+(XbkBw)*ZqJD+ixJoJ; 
 
end  
nTUY = SGdgY(LtZ);     %call for the functions 
vCWT=(MQpM)/(yYwRP)-fTaoM*qkQYc-felLA; 
fcZX=uUGIm*nKku/wPrCp*(hotg)-ksFTG/BMbeg; 
 
end