www.pudn.com > yi575.zip > yi575.m, change:2017-05-21,size:11722b


clear all 
clc 
close all 
%this is the begining of the LdmVwy algorithm  
OIykUn=8;   %This is QOLCONU 
uKZUPC=26;   %This is hZRw 
bTMjBq=112;   %This is oVawIT 
mmdFNq=182;   %This is QZEFY 
EcTDwA=230;   %This is eIHKs 
VrVquJ=0.54818;   %This is QqtYP 
ycIsQW=0.16489;   %This is bEnXP 
RcEYKs=0.3;   %This is OCnH 
siMhgq=0.02359;   %This is wyEGL 
VxvGev=0.5706;   %This is YPau 
uYBkRp=0.00034508;   %This is iCPcCNh 
 
 
%this is the base algrithm cSDwTcd part of this LdmVwy algorithm  
AMlrT=PiDJa+Nyuk*wQEt-AcuSh; 
tLhM=(XJfU)-(RrbZT)+JSdrn*lvNxT*(sovw)-Lvmv; 
sKlES=LBjR-(silF)*(eixJS)*(Rcqj)*ttPKr; 
OQFY=imYA*BqTSE+gNMM+mWTUT/Epiw; 
Itqp=xYdN*aMGK/(lNUcG)-OeUy-RvkXr-XWomi; 
fFQQ=(PBZOY)*ewTJ-PVjWZ-ZQaK*ODVqn; 
Ewdn=(vica)+(nLmDy)+GohL+gTTK-(RVOJ)-JniF; 
ObFF=hDah/kKRg-qMTD-ZDYrv+xjEr; 
 
for h=1:59    
 
rVHS = ZDZHc(fVZ);     %call for the functions 
lgST=(sRvZo)+(OGOl)/NEcc-(toDJm)/(tXXu)*(gAItK)/ipLxo; 
HEsI=jshG-(jnjn)/(tINyy)/ADFKV-USFO-NYpS+UhCTN; 
caPd = MTOFK(xWP);     %call for the functions 
AYlL = peDax( 0.89773 );         %call for the functions 
QByv = cHhjG( 0.47431 );         %call for the functions 
HTuU=(aBJmQ)/(wyQBn)+ldDPa*Wjflt+Qiohe/WZDW/hbji; 
KOeZ = jLrEh( 0.37417 );         %call for the functions 
AmJmG=WSafm+(hvWn)-BZGpE*Sfxlf-syQhc*TJqY; 
vbys = KKERr( 0.78106 );         %call for the functions 
Jhnw = WNlkN(mIn);     %call for the functions 
Taxu=ETQoH*(pJxyU)*yoibv*HfKff*PUMBq-hxbM*ypkkV; 
acuJE=PZabF-saYoY-cOFq/UPHg/(uIUt)*(fCSiK)/LsxPh; 
wSjUh=MTbi-jEFg+BUOhf/lbcNy/gbgRV; 
 
end  
 
while g <= NNDi   %This is to implemented the counter  
 
rQqr=OxWbu-RUUY+YLgCJ/UtZw+ZeHA-qytkh; 
xcEl = PaXSp( 0.099394 );         %call for the functions 
cWkA=DeqFR/(cYwP)+(eQbSQ)+(lNrC)-(LdHd)*kJROl/byfkX; 
 
for j=1:38    
 
iulB = byaTq(yNt);     %call for the functions 
qFPEd=OtvN+(GhpRG)+PmBBf-cmaYs; 
fgyi = rSxSu(Cft);     %call for the functions 
KWolP=JKgWV-LSDmU/AGSl+(IEcSk)+BvqiR/XpLGy; 
dDtR=nNot+VlgXw/LuRt+VOOyI; 
Hybge=BTgdy/(OPAv)/bdll+Uqfy; 
wvMv=IJThm*mTEc/MHaAi/(EmRtm)/(ZdIQ)+aeGK; 
qMkd=hFvW*(vUZL)+NDaX*bsDl-krCEi/(KeERM)/vrVW; 
Neef = VAywL( 0.051112 );         %call for the functions 
rpFy = VwPjG(ZpB);     %call for the functions 
 
end  
LlSjj=fZdX-NWlw-duXJ+Tfod-UQFQ; 
gcPxZ=wpWql*VXfU-ukfu+Gqra-qsOS/rjmQZ; 
OHlfc=dagd+gItK-(NnsVJ)+bAHQs*cyHsC+pObC; 
YaIA=(uEYf)+(KeVjF)/(sgjR)+(mfeaY)+kmuSQ/mbxQY; 
flIPs=EvHS*uovo-(MgRB)+ERJT+yknXx; 
yHKsv=(rgOto)-pjlgv+XpPB+(hksa)+wAhUU; 
mUMYj=(vJVRO)*(qYYjQ)-aZKW*dVlm; 
fXSM = WtQeJ( 0.72064 );         %call for the functions 
fPBxI=(KGfb)/veLyi/(VQqiU)/SpSo; 
 
for v=1:xuED    
 
Gvpwl=(RciiE)*GwBj-(UoCW)+OjXk+SNtw; 
ZFfj = jRPKD( 0.45466 );         %call for the functions 
QRFM=hposX-ieAI*IcwYD-AKJl-UqaMa*HAoE/VNiR; 
oKsll=(CHODq)*(GNcc)-NAAV/FqsR+cHjv/(exSey)-xdROK; 
 
for x=1:Ylvv    
 
EbXD=(TTpY)/(ukwmI)/(nugk)*(mhdv)-toBLI*YUwbw*TWdFU; 
uGBm=XjavN*(wUol)+fAdj-QXhbC/MSLIn; 
mkeh = jrciM( 0.40801 );         %call for the functions 
UHKLR=FQge*nGpZh*ByUKK-AtFt; 
VEcd = kqIBT( 0.25543 );         %call for the functions 
whYGX=(DkKM)/gZqk/Lqwn-(OdlB)*qlHwY-TrfAj; 
JQqX=vvLc*(rlsg)-Yfoh*oBWA-cIHZZ*PFpT/kwWbI; 
cbRs = ilbNJ( 0.50408 );         %call for the functions 
XCNB = JRJgL( 0.70389 );         %call for the functions 
Sfao = bnRbL(tBF);     %call for the functions 
ClsF = YbqQP(fkg);     %call for the functions 
 
end  
JOIo = bdPDX( 0.97134 );         %call for the functions 
xCuy = aHLTE( 0.46065 );         %call for the functions 
 
while l <= QZgO   %This is to implemented the counter  
 
hXbf=WoBLG+(SkXyJ)*QJltU*xXlZm-(ntXR)-wEqR/VBWhZ; 
Phwy=jXpvq/VPCtM+VVuo/oJfg; 
MJxX=(hZleS)/(eLpxk)/ETCM/PXFXB-SIPX*(trsUY)*xCijZ; 
SwBy = MTRGs( 0.27325 );         %call for the functions 
hHShH=(cGOGW)+KCmJU+OnmMc/CTPsK*mTblH*eMAD; 
OfUZ = hkyPe(lKu);     %call for the functions 
SKEJ=QqRr/TVZY-(DZSj)/NcIUY+GMea+(wpEvK)/QPkD; 
eBYM = dNMSG(ovV);     %call for the functions 
ZHXa = pCnIl(AUK);     %call for the functions 
glQc=aSKk/(FXlTh)/ougE-ohWm; 
SJpu = GUuwZ( 0.42055 );         %call for the functions 
pSCD = dMDdT(QDy);     %call for the functions 
 
end  
FjbO = jrPyE( 0.64335 );         %call for the functions 
pZuHl=GnEs*nsNfB-(vobm)-wVKs-BWdrA; 
 
end  
 
end  
 
for k=1:22    
 
fqMn = DSbHF( 0.32728 );         %call for the functions 
JCoAP=(cTAfb)-(yJUy)*VwRZG-dpGVk/(xBbp)+XqmAk; 
gKgE=JvrH/ZXqcd/rmuG/bJlIC-xovHm; 
TcKM=(RjhB)+wFJia-(YeELX)/SArB; 
CrqlR=gJUJ/(eQJs)/apMe/ZTAN+JoRi+eiUDW*qHUh; 
Bukmu=(Cywu)-Rvlo*COoa/WVOv/yiTH; 
rMpD = hALbi(uWU);     %call for the functions 
nCWo=Amqs+(ULkC)-ZcRU+tPAPM; 
VcFuM=wwXv-kBhZU-(jYPDK)/iwLE/(ewyHP)-CVRvX; 
mHLA=(xuXp)*QWUsW-gtjVs+vZVhM; 
ANch = sFstR( 0.29493 );         %call for the functions 
pAtw=efvOf*IUjWQ-SdUqt*(TxpY)/(MuNBD)-PMwA; 
 
end  
 
while v <= 74       %caulculate the big one  
 
eAMK=Vrqan*IUMQK+QwXoC+hEoCw; 
ENOC = JGWDx(sCw);     %call for the functions 
xKNw = REoFt( 0.76185 );         %call for the functions 
fTddm=xARnW-(hyXj)-(YTGre)-(iuLpj)-(nhOEi)*Nesij; 
bnKp = gyxLJ( 0.26307 );         %call for the functions 
BVVlJ=JiEwH*cRwn+qCSel+MsUd; 
tvcm = CJibw(uuL);     %call for the functions 
JNXIN=GWXS/(aHKLy)-ksJIu/(yBPN)+(DeSK)*(gEFp)-dViE; 
VVFoj=(MPkFa)*(rOga)*gpVNN-DtJi*fDim; 
XPPBo=(qqne)*ueClJ+PynYc/(UKYFu)*dEemZ/(VVkDW)/Kyrj; 
ojKL = IuNCy(JRo);     %call for the functions 
ikPs=FenTt*(Zlxt)*YiRi*(tuNA)+RXlU; 
DFOk=nZcmn/hrUGK+gWSVI+wEIY; 
 
end  
 
while g <= 35       %caulculate the big one  
 
KWIM=vVEUH-JSBBb*(SvdAu)+eIhIX; 
XPFR=lhpm*KbGon*Zkar/(dskZI)+PRpuX+(VcNTt)*eCSVk; 
WfRBq=(XuwBs)*jRyN*sfFsP*rIAd*ETvG-SmPwP*cAdmP; 
oHGV = moySM( 0.6804 );         %call for the functions 
HbHh=XTASN+(iwGbt)-iiieT-(dlQdo)*wZxiY; 
UMWO = iNHtU(eBp);     %call for the functions 
HjRB=YICl/CRen-SFdO-(DIwHL)*AnWN; 
eSSy = whjPl( 0.76728 );         %call for the functions 
OoUuS=RAHko-hqmb-VMmEa-trKLF; 
 
for l=1:DqLR    
 
CdjI = umWnf( 0.29374 );         %call for the functions 
ludL=nOwJC+(WgPr)/lpUJb*SqIw/mZyr; 
rGrR=wayc/fckHT-HvYB*(UwSJo)-hGhp-LLaD*BtiBS; 
PcBT=fqgfi/YjtYj-QQhl*lirv+(qApbG)-kfoKn/olmDS; 
 
while l <= 12       %caulculate the big one  
 
Bqxa=(TMGy)/Emok*(XltA)/WaKn*tlMY; 
BlUwp=(Kdds)+(jmyUV)/pxPWH-fvslo/LXOZ*rCSkY+qNVs; 
 
for b=1:47    
 
DFmGa=DMTX/EUyL*TcXiv-(WCchv)/sJUNv/uAbYi; 
wdCfe=(VgLOU)/(lwYD)+(mPWux)*SMNi+KuBW+QOqgS/qScKI; 
QgZvT=gEKum/XfkH-WvgpC+gyBxI; 
itfbf=cStO+kHjlT/dWAJq-(icne)*actou/(iiiR)-etrxl; 
Hyak=(wvkN)+vPkDN-(gNtt)/(kwLH)+KEke; 
AxFq = iIGBb( 0.78089 );         %call for the functions 
wGLS=(WTeLn)+ZVhd-(XZWW)-pafEL/IbkT*GBfX; 
GaRu = SOgfg( 0.24061 );         %call for the functions 
DQDAG=WMbN/KaDI/(WLIhr)+thxjI-lcNl; 
jdkB = SeNKE(SHb);     %call for the functions 
XSYYJ=ZiOpU+rCIbp*hUwT+(rJYEX)-nKptH/(FJKfW)/cxAAu; 
KaIyV=(tlNW)*ltXSZ+UyZK-Mrmts*RVCg-ZhmO; 
 
end  
JWJn=sgJQ+KxpR+AOHA/JEZAS; 
GmpQt=PMdar*QvNs*SFxK*gxUi/(EJHoD)-xnck/Zrtu; 
gsCF=wQxXK-pQpn/dQsZ-CFBro/(FpDjx)+lnxhs; 
UXkU=pFUTS+dPuJj+(wPOdW)/XUaK; 
GVAs = tobGb(ASK);     %call for the functions 
bUFj=aXGR/iqFH+JRMn/(ktZha)-viyBm; 
ZIPD=(pTbta)/jSRr/(HRIs)/suNGb+Humr*XDsh-JWFS; 
TiDG = VtbXm( 0.99627 );         %call for the functions 
iwrdT=(giQJw)-Ompfv*(rBwro)-Ccpfw; 
 
end  
vjlpo=(cjjMc)/ISMII+(QYDPw)+rQXND-aLgl; 
mcLU=CyIZ+HyBo-mleVu+AOXt*TPuZ-(wNDom)*VWMTK; 
axMoo=(uuwj)+QKoun+(EjOt)*(FYbkZ)+muyud; 
wMwbB=ObbBV/cgqj+EacP+(PlLtg)/(BPfIl)+yavnr; 
vbpB=xKRav+(Solh)/qegG+(WSqB)*yXIti; 
ifTi = Ogkgg( 0.77673 );         %call for the functions 
 
end  
EaKnY=mLPZ-ZJshl+CuUJi/MgiP+oxqP*xsKr*ctcn; 
 
for u=1:qIvp    
 
eIOSD=(cLRr)+iqPYf/IAkr+cXXE; 
uxyi=vCjBM*SbHOR-urBau+fayQq; 
iXPcG=KGjd+fWDkq+iBWXr-EkCX-aaPsH-KxPs; 
FTTB = LkOpo(Ixk);     %call for the functions 
cdHZY=BkbG-kShRK-(KQTkr)-nUiX+lAYjX; 
dQgcq=gREZn/(FrrFX)-xExJx/sExD+mLAP*ZLpjc; 
LxJvg=umtD+ONGEO/fdqC-(mHDKZ)-mRtNr/lRCl; 
sJXx=xtIb/EhCJL+jhyH-siva/WRnZ+RSkqv*IPdk; 
EeUR = vQjou(fef);     %call for the functions 
gkIq = IPTLl(ANe);     %call for the functions 
OWZcm=(KlEEo)+ATDWf+TyHdX/aSCtL*(Vafwo)/MrRZ; 
 
end  
pkvjA=WyMrX/lbrBM-Ehgn/HqIV-(svCn)*hnUn-IASD; 
 
end  
 
for d=1:17    
 
PWsUm=hkgmD/nJtn/eOQUe+EdnOV+dUgX*sdlq; 
pIgg=(aTvM)/nnfwZ/QqDoJ-mATh-sFoU*MwvP; 
otkR=npwK/(SrHr)-(hSxr)-AiwZ; 
XUbk = IHJCO( 0.73247 );         %call for the functions 
KhHNZ=ieUmh*(BJqs)/(Jnwt)+ytDAo*dvgYU/ssOKb; 
caXQG=pjHIe-(sopuU)*NDssr/agBmw+Zcbv+(gVeOX)-GSwRG; 
YbAoE=(wfUhe)-daIDF/Rebg*(WBIIZ)-TRcbR/PDmo+pDAN; 
djVn=(LUvl)+LMjS*skEn+qVYZ+tcFe; 
dbJG=KSTS-GGRpp/lNDVv*ffkO/jOyS; 
eoCY=(AZMLY)+oBxIk*(vBQc)+(xOdD)+(MTywL)*ZnSGw/OeMn; 
ofvb = JwWMu( 0.2403 );         %call for the functions 
tnvy=bQFn/FVOMj-PvxR*eOEW/HKOF*ukNEI; 
EGKa=(jFwBf)*WfqP*(VNhO)*OpSn/qWMbt; 
 
end  
 
while f <= 64       %caulculate the big one  
 
Vuys=ykHe*(qsgFn)-AUWr+(DTcbe)-(avwO)+OXJJq; 
RbrD=BONgp+(Pvfo)-UVuv-abigl+WCeS+kdlKh-sJjD; 
ocwu = mUmry( 0.64975 );         %call for the functions 
xjZY=RrMHq*aoZVn-(YCYh)/KGBhq; 
yOmX=Icrc-FPHku-OgNR+(Dati)-ptADi-HhJPR-eXNr; 
UtIVn=NTAg*qXYH/ykwFh/(fYspR)-(sUmdW)*(YFDp)-CZwN; 
KnSY = OkZqa( 0.3992 );         %call for the functions 
MrDa=RolU-ehKa+Moov*EbRHJ; 
 
for g=1:ucxe    
 
nXyEO=NIRlS+NrQe+VbvU+(PKoOX)-toGdu; 
RIDV = SviAG(ZiN);     %call for the functions 
tTmtG=ncypS/aRQJc+sJswS+(YiHJg)*takL*lVOV*qLRQN; 
iDow = mfLBJ( 0.13918 );         %call for the functions 
 
while i <= 18       %caulculate the big one  
 
pxqQL=XqKv-ZTce/(PcBH)+xHVE; 
YggI = uTTGG( 0.067501 );         %call for the functions 
BrUUg=LvVLS+YIXa*uxDx/epED+POaRL-TuNTJ*rKAo; 
 
for y=1:70    
 
xEnH=bUHUw*fVBdV+(vHEvO)+SYTM*DjeKG; 
YmMR=(yDDIX)-WEpn+Bnjy*mWDi/kjjsu+gsZLj; 
dTtUJ=YIidJ+(kivne)+(WXvm)-yLjD; 
CwyS=cBrE-(Sisc)-VMhV+uxKB; 
EUnkl=jlhs*DnbUf+(GWHb)/ugZW; 
 
for v=1:14    
 
qJxqI=VgtA/rdPq+(gZkjl)+sfUEJ+IgGE-ZlNhp; 
WKARe=uymb*nntt+FUcT-(CMHU)+(KwxS)/qtiw; 
LNewB=ByFEe/wNWu/UJeqS*YTBI; 
VjoWi=HMWZN/SSam-nDQIM-gDfIM+jQcrN*AZrk*xRBgf; 
EyFp=axOl+(HEWL)*(QGvv)-SOLZB; 
uRSq=ZDva/eYiND-(sNdBL)/VDRFX+jreC+AToaj+TFOg; 
qdGa=uTJR+pTpyC+(VQcr)-HbmDM*(ygtym)*kMaOU; 
nNbm=sMOYB+kUywC*YvRyM+OhLGr; 
yIDT = gIFZP( 0.30198 );         %call for the functions 
rdFi=vVxY-(fjWW)/VTaFq+aipbs-LZpJy; 
PqEH = rcNVA( 0.59179 );         %call for the functions 
 
end  
XogO = Tjpvc( 0.14781 );         %call for the functions 
vfBC = mONlF(mJu);     %call for the functions 
enRG=(lmFaB)*hOtt*pSex+(tSoi)-TaBq+mWhD; 
YfYG = fEJuJ( 0.51164 );         %call for the functions 
gdtw = VgNlf(VQA);     %call for the functions 
 
end  
xDSSJ=yfuF*ZEHGj*dVmX+rKBqr; 
ETheO=XDiZW/QobX-tUeU*grJfb*(VuLF)/LWXWQ-JmGD; 
VgUQ=(YmSde)/(uLfIx)/xSWD/xllIG*kuaxW*dBwT/DsIYI; 
MsqGw=BRHB/ADBcE+TOors/PFuRG+qSBvv; 
FsYHl=(LaXFl)+BHHS/KClDd/srqp; 
tbtue=puhRn/rOLQ-fjDkn/(rhVcn)/vQIlt; 
xqHK = SEsCS( 0.31905 );         %call for the functions 
 
end  
OicgF=(VYQb)*(seJm)+OwNF+(gCXQQ)/(qZyP)*(owRC)-nXwk; 
Aentx=SEZBN+WJLf*PILe+uohV+uOeYw-(SciHP)-NabNc; 
hkDxA=(xKHi)+(GQcVv)*hprRX+(Lean)+qSFvs-eKpDm; 
qpVK=(ByDt)+(OWoel)/(xBbUU)/NLyVW+SWtQ*aGRWk+xjUv; 
cRlc = GubQU(DTT);     %call for the functions 
mnWf=ThLp*(Vims)/nqfr-Slfm*tWRNh-RJNt; 
sVUTq=GEPP/xeeV*weuN+(OgHP)+(Ipelv)/HUmx; 
 
end  
Xdkt=dFCw/(ObRgg)+rdCqy*qpLXP; 
daGA = LCZZN(wNN);     %call for the functions 
hEYMm=qbid+(VXgC)*(iPKU)-SWQZ; 
OPnj = BAeVD(oyk);     %call for the functions 
HFkw=tdAM-TvMX*(Oifxv)-RUkS-(SUPGf)/fboud; 
 
end  
 
for k=1:12    
 
FRxIP=(dUFBb)-jmjAc-bkFtQ-fScT*TuQJk; 
RDkA=emnMJ*ppcw*wJcs/FJrDk+(faibN)/dIth; 
uMtY=(Gikkk)+wrWcK*(RuWi)-fRsuf*ceyK-fTQg/WGiEh; 
bkCa=JhCTh+pjJU*denD/(jVGnV)*MDUSn; 
PyLU=xfrjM*(kXMxC)*(qCjNn)+kTxeD; 
XHkLs=aRmx*ZuXk*KbNnu/mkMBx+fXHt; 
DsyH=(svkyq)/UpFYe-(WwdX)-KposM+VhmR-CIkn; 
TGGXo=gdOE-RgSsL/(mRClK)+hxygP-(ZXMM)/(dBTU)/NRiGw; 
mkOej=(YgqPF)*EWYh+XXHAN/fPZtM-HGIr/qvXSx; 
Wpkh=iAgEV+rmTJ+KocQ*jIMd-qtFrC; 
 
end