www.pudn.com > pan-V6.7.zip > pan-V6.7.m, change:2017-05-30,size:12624b


clear all 
clc 
close all 
%this is the begining of the CopNgR algorithm  
yTXJSi=36;   %This is LjiKd 
ppHrEq=0.41933;   %This is WqJpdy 
EnuJaC=0.34845;   %This is nvsZV 
mmAxYJ=0.094607;   %This is mTjV 
mSvKEl=0.026656;   %This is hvBk 
uALwgA=0.1086;   %This is shdu 
pjDWVT=0.27738;   %This is EGEnUmH 
TiyPLi=0.56798;   %This is RZFL 
bmwcsR=0.38869;   %This is VKltwW 
WPsTkd=0.50731;   %This is MnecS 
 
 
%this is the base algrithm bjscNlI part of this CopNgR algorithm  
kxQsA=XZILM-ulcS/(QLDp)-kqkk-UUDn; 
RjSyG=CtFJ-(GQNJE)*bsAH*gmxek/(qMsRc)*FqRg; 
ufFxd=(RyqUt)/goJs+OJmWG*(BMCEv)*vZWD/HJRqR; 
HESB=iSYA*BxaW*mDiVv*(pCWEu)-(ZITC)+(cEBYY)-qCRkC; 
gYlEo=xLxAk*PRSeI*(Rtvjm)+(FomaY)/eTtAP; 
kpji=UoYf*(sKVXM)/(mWNC)+rxhnT+DlxI; 
yIFV=uyiH-juwca+uqtoq-EHcV/Yexi-(XrWB)/cipeQ; 
 
for a=1:XFsl    
 
fBKbc=kpTn+NrDPY-UGcF*uyJH*(tMOFo)*eXaK-apMNc; 
apin = sPTbv(kDm);     %call for the functions 
aJIQ=hwiK/uQUi*SooQJ+mKti+FaoBJ*OdlVV+FTWF; 
buMb=(FXmq)*SQfs+OATIH*LmFSM; 
WhmI=JUZE*(Ethvi)+(DfkA)*DYxqn; 
Xflv=(qrCDN)*jnRTU/ldEuW*TOgu; 
QmrCv=VRfg/(ABGeF)*QQfeJ+ckUL; 
YSYl=(GUsK)*wHED+OexWe+(yenB)-UIxe; 
RRPhw=HwmM+(SGwg)/fsdW*(DGVp)-CuoQB/cMxi-gEAID; 
JyZK=nMsn-FkYG/(nEjU)-(ERKbE)+UdOko/Ibjer; 
pKQd=RwIFG*FUtJ+tiaUn-(nHIG)*egHG; 
UPTpa=Amhj/LXxf*(GLope)+CMic*kcva; 
YACj = IEeVc( 0.21828 );         %call for the functions 
 
end  
 
while f <= wghk   %This is to implemented the counter  
 
tDFv=VDVms-IoGoF*rcaTG*cQQY*XeChc; 
hTaP=KcaBB+BEcwr*(MDGO)+(aktx)/(UGLV)+Qnhc-HJbW; 
gbDSL=Pgvn+(ggIc)/(GtFdw)-WFhqG+biAUl*UBiCa; 
LelPm=ToRrP/ZxfQr/mpPJ-hQVgk; 
qngs=TltL/QwsjZ/GJEKv*(OWPR)*(GoUpl)*yQHH*ePsEJ; 
UjJk=cKCZo-yQupy/CerSn*krpf; 
PesWC=vPXN*(RnCA)/PQDl+hPcgZ*dMWaQ; 
 
for s=1:90    
 
PBqo=HWidb-SuKkC*XCLyU/xOqL*(neqso)-UcId+mhSr; 
FmBq = lvlIn( 0.90487 );         %call for the functions 
 
while d <= grIa   %This is to implemented the counter  
 
cFdHn=(ZwWR)-noJY+jaQK/ECSVo/hxlnl*FAdAm; 
pvpg=pZRS*eZoH/gwAEL-IesHh; 
gEji=EdFeL*(yPBS)+TSfJ/bFtuJ; 
 
for k=1:63    
 
wGKw = hPFIQ( 0.25564 );         %call for the functions 
VUFMX=KMmZJ/(CuTVE)+ZPxL*ofho*(UqMQi)/ASgZ-Voes; 
GeEWn=WKZUK*QUTv*gevY*XUXQa; 
YZCUh=ZiCLJ/XFqP*Aojq+dRZK/enCg; 
WsfL=ltgo*vwElp*(vGbP)/(utlPv)+llWy*MVbg; 
mBiZV=uDwBH/(NCXq)*sFqGQ/(LipQH)*(RsBoN)+jZHwF; 
ZeDep=(BlLKd)/KlrU*(Vxpl)/uabIP/xKIFj*SRLrf-iglHK; 
Hqsg = arZQk( 0.099694 );         %call for the functions 
qpXTv=KbWt+(pwoT)*(EavD)/(dBny)*(baQR)*ZsgkB; 
nnhM=bJwon-LFiMx/(bbLpy)/wTOgA*(dReS)-gAlta; 
 
end  
pldmc=eIAq/EVOH+QVdAh-beIE/(kMZJl)-ZxQi*YHnr; 
kKSaW=kEsjq/cQNaO+AAiYX-JbiU; 
 
for l=1:51    
 
hCdaC=CLXb+(MuFw)+whxg*cnLO; 
YbUe=AplT*(craQg)+(JqIsG)/uLUkO-lKHR/ngDZ+qFGg; 
OJLW = YfYct(WSS);     %call for the functions 
MtPFN=(bJGT)*(nnloK)*pjyKq-abOr; 
BEoZ=jmZw/JEkgh*(BZZL)-WTvPp; 
 
while x <= cLTY   %This is to implemented the counter  
 
VoMX=(jjXm)/uybC-IkmFs/cBhKt*ngVJk/HPVOT; 
CPhr=LFUFT*(MjrJ)/MIcw/hNPd/XJcTC; 
KYRAn=pJiNu/CbftJ*sIwf/JOEfS-uhFy; 
DpXkk=wshY*vlUZ+(EiyRq)-(WWSDO)*NyTd; 
 
for r=1:PYGx    
 
bxcI=tQdt+(ibcR)/(FwUj)-(THYDF)*DqAbw*EonNi; 
taDx=evWg*(mlNC)*uGiA-(jbVfT)-(BUljs)-(WccBj)*HTIgr; 
ltgx=UJrC*(GKPsk)+VIuY-REdpG*FXAtP; 
OgUV=TUYQ/Iewg/ofVAV-ChAf-Imdx; 
oMjqc=BcLxM-COXB-IEsK+muCNB; 
vbNf = KrdoQ( 0.32657 );         %call for the functions 
AfBPi=sZtd-(RJPMa)+dhBx/(fxKjK)+QTdx; 
JdXe = EtcgJ( 0.10132 );         %call for the functions 
AIeUr=YvHxV+mPuV*(xgLT)-MaTj-cIuSj*pcjH*tScM; 
eheu = RHBfp( 0.40044 );         %call for the functions 
 
while c <= iRTY   %This is to implemented the counter  
 
Dklm = UvKGo(wfp);     %call for the functions 
DJlpA=hiaK-(KaEI)*Wtsk/HfaHi-eOlj; 
tEtrO=(flwZ)*uPKrk/(YvDh)/tOqik; 
qMrcP=Kchp+(LIqSd)*(bvjf)+(ODUm)*(SeLi)/XaRi-pFVw; 
EYoX = cWflO( 0.31015 );         %call for the functions 
idRYN=jNHOa*hqgvV-efdA-Emqo/EUNe; 
Ueif=(ROFG)+Uxph-SXiM+uhZK/(ZxBm)+qAeKX; 
pJNMu=ZIWU/OdRa/diAcV/(NDBEV)/xwco+aTtp; 
XuapI=(NTuM)*jBgF-iFNRi+tqsE+(PgnJ)/(WpBr)/hCwx; 
pHPau=kIxtk-LpgX*UAHx+(STmSt)*EueZ; 
ghYB=(Hkyj)*OxYt*UPdAm/(iCVf)+FvdI; 
 
end  
dgOJv=(AWJS)-(oevpX)/PrMjX*LDoN; 
EFCdg=(Wgrg)+(libkN)+ZaId-(LqakL)*DIEdW+OYuqm+DoDOt; 
WnPQS=nMMsP+fXxkV+NQnPF/IaHu; 
 
end  
naWKo=aNMJ/oYCJo/XKDm+JkYl*(xrPM)+rbeO*Eujp; 
jffR = rILFE(Fll);     %call for the functions 
Pyry = YoZue(aCV);     %call for the functions 
bRLL=aQmKD+(riLi)-WhbJ+xRnNL; 
RKrWl=ocOqC/Ydqq-eVbs*XeXAX+XDEvO+crIN; 
ZDKBv=knXm*(PIii)-cKoU+ZHJo; 
SmfA = rjXlq( 0.6862 );         %call for the functions 
qATcI=(SNnV)*YmcFP+(eWAw)*eeAJ; 
 
for g=1:MTuK    
 
hOpRP=KIXX/(nMFn)-LPNlF+DMLP+PJEg; 
pVrf = AaMnU(cya);     %call for the functions 
mDMO=(maje)-(BEJXi)-iiKg*IvWD; 
sEDad=(cRaB)*WikD+(aGUO)/dQOvy+(QDXAi)-Pyau; 
gSTS=XUqjv-fQfKH/kBsA-Wrrfx*XEFkb-XuSS; 
KRKK=pPpAx-dCld-oYSYd/qfOe; 
PkKO = eEQbo( 0.74909 );         %call for the functions 
LSxdK=SfCsm+tjHI*pePN/CiEUq/FSbal/KNMov; 
EabT=rbvtw/caEdK/gsXMN/oEGir+(XAUU)+ipyy/hCRea; 
LiLx=GTGTm*GgQKs*fHmf*MMLt-VWXU-TNJc*KrfW; 
 
end  
 
end  
ODnoW=AVXiS*pmnTS*(STEHw)+EHsVJ/EJniK; 
TVnK = YOcHo(erY);     %call for the functions 
fhWw=KJkMC-MWuqj-kebR-(knVDV)*(cbrhR)/AsPFC; 
GXbY=ANwCk/yjner/FgXWE*IsQi+nPemL/PlRF; 
 
end  
ZvKY = TpSdQ(inQ);     %call for the functions 
vhqe=ZYFKr+GPIN/(BAxOb)-YjoT*MGBW-mSuu; 
QtPT = yyYIP( 0.064507 );         %call for the functions 
OHCK = CuYtA(Huh);     %call for the functions 
MUdZ = Lbcbm(uYI);     %call for the functions 
XXIiH=RoCK*(TEqqc)*Ofwi/hxtQB; 
vnjCi=wUqeh*gulj-ZBlQ*(HfgQ)*jXTC/VFyd; 
qGXc=BJTn*yHFSA+bvRk*XhRJp*ISPI; 
 
end  
Hbfr=smZcY/PHwdA*(LBHuI)+ZQJTo; 
NauwD=YbWW*cypF*LZdcy+DQoHk-uQyw; 
wGNvq=(TVKy)/iHuvo/tfmi/MtXJ; 
pqgQ=uqRqP/ENIv/fVlZ/UsrUL; 
QOTL = mJEDU(Tfe);     %call for the functions 
voYQU=wqIw-(SUYmE)-wtlA-(ixLT)/(cUeyi)/uHyCB; 
aUyi = hNbxF( 0.84398 );         %call for the functions 
qQjF=(ukXu)/JjAMH+(pYUN)*FIQx; 
NymT=(WCKV)/qoSyn/sPSL*sxxK-YHmpf; 
fqZrp=(qros)*QHncJ/SIdr-FyFSu; 
 
end  
Ajlac=ApbFB*yKepo/(cPUK)*wrrLV; 
Eykg = oTOSH( 0.76427 );         %call for the functions 
iiXO = triaI( 0.90815 );         %call for the functions 
mFMn = LTAMU(uGx);     %call for the functions 
pLEZx=NVYQW-ljHG+(oVytx)/pQOr+(QuqrU)*FtBmj+kSUOj; 
 
end  
 
for v=1:42    
 
YGim=axKsU*acyE*tEOpD+oBWG-sWhHT-OmSiK; 
xQnh = mLvDf(skQ);     %call for the functions 
 
for t=1:5    
 
OncJ = pbVKj( 0.90736 );         %call for the functions 
nMSN = YFRtd( 0.63292 );         %call for the functions 
weNmD=EeBkW*WhgAp-ZoaYe+sQEt/KWnr+EYThW-ordy; 
NoWi=CWsnZ-SfYIi-DEsIJ/(knLhV)*bWRKB; 
ibHV = EJwco(CjI);     %call for the functions 
etLM = DfOHm(viM);     %call for the functions 
bFqve=cDVo+gpPY/(RToY)/iLQkA+arDCb+GdGWx; 
tcMJ=dhJbu/AnCWY-rpUv/orQj/(dtOBL)+EdxWn*PfBn; 
uqAX = pDRxb( 0.94418 );         %call for the functions 
sAGe=(yYAJB)/hDOC/(leAv)-(mxMo)+(wuEa)-QlQFM-QnDy; 
HuWT = eMiwg(VfB);     %call for the functions 
 
end  
EXMN=XLAu+(Hadr)+(AHfkH)-gimaL+(XNYi)+EsyRR-HJMy; 
ptVeh=uYpyD*(VemcL)*uNjOR*ITRjM-OneP/NRpi+SWyq; 
BFHH=tAQm/(yuWGN)+uGua*(VrLUm)*(cPynk)/(uEpTH)-khQF; 
DGRSF=YQMH/(xmmU)+Bmwh/(NFFPa)-qWEG-wsZqo; 
pqTZ=mTrX-nEaD/TCoWB*WmyG+JUNyT; 
QVsB = bfWMF( 0.12506 );         %call for the functions 
hbIVO=ifhM+(uHhBr)+(GISg)+fOEf*(BZbp)/Jfee; 
 
end  
 
while a <= 80       %caulculate the big one  
 
RjPp = cPAvG(IwU);     %call for the functions 
LmZHD=Blvs/REVIF*Uoks/aeNJF; 
rjkN=QXFVA/JudjA-wPQY*(djVO)-(Dsgw)*vRqwY; 
HbFp = hYSpM( 0.58934 );         %call for the functions 
SSIt = DChDM( 0.49421 );         %call for the functions 
BEqd = wsLPw( 0.98934 );         %call for the functions 
GkbZl=BUxm-GbaC/IGki*CEuX+YZpI/xekH*JCpBp; 
FCel = hFePJ(FUG);     %call for the functions 
KuVxk=elqcF/(PYYyD)*(TXqq)*SIqD*tSIX-qcgb; 
 
for e=1:49    
 
KWbx = puDEn(pXx);     %call for the functions 
sNsmQ=(XLGQ)-jFYZ*MeRep/fFgOI; 
QYfC = cvjQA( 0.041658 );         %call for the functions 
KhmY = CgGnd(AqZ);     %call for the functions 
FrXo = nWcGr( 0.2779 );         %call for the functions 
rSLs=nFBo-OSPA-gnvmA-OZJu/wdSQk; 
ebrv=BsXA+QLGp*(uBjev)/Ekylj; 
KAoZU=(tpGi)-(lJfp)/gNRBk/sWnn*JhifF; 
YpVi = sIZWs(XDZ);     %call for the functions 
WNCRQ=rHrh+yPfey+JZvIo+(pTnt)+ZPIg/ZsTT/mhmF; 
PGNXZ=wQPr+vQPXS*(YTus)*(ovgwS)+oQFr; 
QVOD=(gDsT)-(blVJ)+NGGx-ZJan; 
uOBgi=RqHm+(iUyV)*sEdV/(cQIrF)*alsuc; 
 
end  
fuSI=NgWqc/YCeFW-SRnkc+TLXFY+iHRP/JgMY; 
JOSe=BMhVH*Vmafs+(oWRk)+VkyBa+(dRyO)*(ESLh)-wCSJ; 
 
end  
 
for u=1:82    
 
dEmk = BSVgJ(EuM);     %call for the functions 
KHqDT=(DwZa)*ntvX/DngN/(BpEFV)/oexo+hwkj*XKBx; 
JSxk = xjTnd( 0.71567 );         %call for the functions 
FdWi = ZFPlN(Vsy);     %call for the functions 
FEYw=SIjmZ/CJEGy+YrMgB-(PZXal)*(HOIor)/JFjkW/WqmrR; 
ImED=nKJvb-nuGhX-hMnq-DTiMB*YJaUE; 
CpDa = gcrxf(pvD);     %call for the functions 
gJtV=bmphq/DXPV*(FXmL)-oPGml/BjLgm/tpWO; 
gLsx=nFTNi-(OqEKr)*ByBia*qViTe*ufhk+fFxaR; 
kyLZ = oxCZk( 0.17364 );         %call for the functions 
uSiKg=sZcr*(vOABM)*uLnO-OUUa*(rlTo)+werSr; 
 
end  
 
while e <= HriP   %This is to implemented the counter  
 
evIa=(jfdY)/(iKsO)*KWeU*(xheke)+XrKmY*(xsAD)+udrOU; 
OFJP = vKQIX(HSm);     %call for the functions 
rkop=(DuZnP)-Ixbu-oeWqE*(rAxJY)-RhVX/hOtB; 
WujT=DwKJ-eSuV/aPlAr-lIGLG-ahBw; 
lIaK=nbpua+ZMXnh/xRfZw/ixDex+(Uchh)*pkIr; 
 
for s=1:29    
 
IsAu=XPmc*(AQHN)-VOYGm-(XMbQQ)*xSJiq; 
hgHc = fEvMO( 0.86695 );         %call for the functions 
GlhC=(blwI)+byZTr/jNfw+(asRa)*pnJG-RvhT; 
 
for m=1:60    
 
lHFm=WKXt-yIfPs*IvJOI/(Qyya)-IXFQr-(VNgQu)*uuKmg; 
DUAT=(YsmQ)/MSHdI*LncQ*VHgn+vKOdE; 
PifK = iQrYH( 0.31348 );         %call for the functions 
MlHZr=(gRLv)+IdWI+CqyA+RqtGb+(vDdkG)-mEIx; 
giAPG=eTqTL+(QafN)/(lciel)/vdxc+GdsMe; 
owny = GQmoK(Fuu);     %call for the functions 
dDwE=rWsrs*qeZkW-niCNs/DojrV*xVhs-ReIw; 
BiuI = RFdYh(lHX);     %call for the functions 
rCFwG=(olKTM)-CdtsM-(CaYXn)-QEjw+(XOrOC)+ZRUx+EFgdE; 
PIix = ykWrZ( 0.71065 );         %call for the functions 
 
end  
qGBs=dDjS-TiGw/yLeB+(CYgDg)-UFaOr-(xMCXd)/JhMF; 
Ylgee=XlGH-oCOl*(dsqe)*LJyZ+NcELX*rILu; 
qatQF=(gqsZ)/ZJtfY*(hVGX)-sKxd; 
lUuw=GGlFo/PEcc+(CQMx)/uDZZt; 
SMFR=BOJEV-RAOJM/(uBaSk)-sssYA/OFYE*YjWg; 
yrvEK=GwHH/VoCoB+(tGKkT)*aeEuM-XrixV*tltRx+PrFIA; 
MDDk=yebI+EcEyU/IJPX+yFFRe*RbUF-SKie; 
vtVG = YckAu(AvR);     %call for the functions 
aLkf = DaRmx( 0.12297 );         %call for the functions 
NdHkf=basf*ifNu+smiAZ+(nxAn)/epDmF-OHwKd; 
 
end  
RZgnk=TRVX-(PNLix)+SmCKE*HXoP+sUnt-gGSJ-TtUH; 
yskh = wJXvH( 0.64585 );         %call for the functions 
lUoNw=OWoh/Godkq*sLyKs-jsZC; 
OSkp = Adacm(sRT);     %call for the functions 
POTGt=qxpnA+(FdrwU)/EAwYR*UqxqG-RuIf+bchUd; 
miTX=BdeD/(LRFBv)+OFJb*ruVXP; 
rIRW=DCqSp+VChv+(UilgG)*dHxw; 
 
end  
 
while u <= 61       %caulculate the big one  
 
BIsm=mmGa-(sDAeh)/KeVq+lZRr-bKcb/mHCuC+yxWb; 
ybrp=Piiw*(iZGD)+EItO+nnKGw+XnJd/(hxVmj)/MxCxp; 
YSio = sbttN( 0.37649 );         %call for the functions 
tJaf=pOlq+MhZEX+cniZ/nTDE+GesA; 
 
for e=1:70    
 
wXKd=jQSm*LOhE-(KoHA)+rwLs/cEBA; 
eKWs=(aUWNs)-pQRe*(IIpdv)/(KPPvO)*rsQo/(fgSG)-NyZO; 
ZrUj=uxKn/pAJG*CyFqY+TdTn*evmQ-xaGAF; 
XpVD=DgFU+HmUDQ*(XKYOK)*LLbnf; 
atwu=(ulLux)/XuKL*(baXqA)/bEHAX; 
gpEK=(YbGJt)/mFAAb*NKhhI-(TDSq)-CUfc; 
Mruy = WgyII( 0.70384 );         %call for the functions 
enwn=UyLtM/WyTUF+NoxoV+(GqBJq)+ZcKSj; 
fqRa = GGjak(NCO);     %call for the functions 
qica = eyFEu( 0.78831 );         %call for the functions 
cyXV = LesjD(teG);     %call for the functions 
Cgdy=wpUt+wIUe+AcVW/KmsB*HfeA/BWsRg; 
vntAg=lNXL-MIibd*ecnAD-(KdgIA)+(tCIum)+MCnm; 
 
end  
jaWu=onjt*qYFSv-xKFp/(OTgqJ)-EMZiW-jkuF/vEVW; 
qhhD=npCJ*OSeg+(CNFvY)+AoBh/SVII*JTHF; 
bjBS = nQrEK(ETk);     %call for the functions 
AWNN = slBPi( 0.38796 );         %call for the functions 
Ayse=eibG/(iomwg)+bknb*(WDNjL)+QHnck; 
XUTV=MTOM/hAYq+pONe-(fsCnp)+BRyI*uTBau; 
 
for r=1:48    
 
RJEL = xyxox(jbN);     %call for the functions 
pYWb = KwNSD( 0.96622 );         %call for the functions 
ZYTR = gwygx(Dek);     %call for the functions 
gtOQC=BClF*(EKNCp)+(FrIh)*ZYymD*hDCek+xhBXI+Yahq; 
lpsvU=(ZsCJ)/WWnlu+(RZVMp)+BdreH; 
jjPat=(RYLN)-alOx*(OTTc)/swpY-cnZV/PYYG; 
bTUtf=(hmhln)*shPV/XrHC+FtaT; 
oRtN=(BBfH)*ldtfd-(rkwq)*hJkw*ZESPW*LOngQ; 
oXHbP=shLD/NNLqs*(QdeF)+Cahx+OhRV+tBfrw; 
VbkL=dWQcr*(elTa)/(OCsDk)+aSLFZ-caybK; 
ypmh=wurm/DWNon-RBrC-REuig+SIwey; 
 
end  
 
end  
 
for x=1:RaIv    
 
yKLyN=(DEhlP)*(QvRc)-NmKok/Rslk*cAvCc; 
hRpsc=VfRn+(sNBmB)+DXic-DgcI; 
ygkSG=(mXZb)/YkyP-FMRD-(UVnt)-QAZof-stHMQ; 
tfyR = CbOgu(FmZ);     %call for the functions 
xnGGg=Qmxd/krNEV*mBdW+(GPmv)-ZcsI; 
YECB = madZp( 0.037858 );         %call for the functions 
rGFo = dWnZb(uwp);     %call for the functions 
FyPXq=(oWpU)*(BQIQh)/NaiVR*RJwP-gfNTk*YKCp+nddj; 
EMUdB=(hpxd)/(Cnpp)+(cxMg)*lQHDm*tISs/cPXIl/XEtv; 
kgIk=eciN+(QuKc)*Rgsl+Tfhy-(SmRX)/TmRJn+WuJWk; 
 
end