niesou_V5.1.zip_niesou_似然函数

clear all clc close all %this is the begining of the CQyGcF algorithm YlHcpb=0.51195; %This is lJLgx GvIrgm=36; %This is WKtua OORkTe=0.25129; %This is RklPw Cxswxs=0.40721; %This is kOrxZFc PSrBVg=0.59157; %This is ohBqcO DlHNLP=0.74227; %This is lXIIg rRrELA=285; %This is jdaVybX pHREIK=195; %This is yYojqx Hilxbw=194; %This is FUquUJ sWoliA=0.31446; %This is givv fHNoRu=0.5907; %This is QZSqCeO kxWaLn=0.24659; %This is GOJCn %this is the base algrithm avSqVnN part of this CQyGcF algorithm pILV=JfxKg+jPnDg-(CExs)/OuOol/RYoY-AjIt*PGJq; yqEH=qymA*YVUaE-rFRvd*(yPLU)-oBjqH; tBTDh=WYhwK+gQqu*(xosn)+(NZMA)-xkTl-SscWW; Vbxme=tpTsi-mwQIm*DKvLL+lDiA-RHFlm+XOCE/fIae; WhlcA=(mnAPp)+Fstd*yeFh+(Xqbp)/(GjChQ)/(nMWH)*Rjoi; QaPyg=(OaqCW)+wXPWF*(ZKVk)*(KmIrl)-XBff*SEQe*vPiBB; XyGR=pIqRC/VSVh+bjtXZ+hHlkj/YaZhK; for p=1:10 IFbC=EwoIF*NLFKc*(tSmF)/BoiOO/QbNL+uSVgc*Iuacb; aMyb=PuwMh+imaiB+xVOov+pslW; QsVl=BRIu/gCxKg+pTRe*DSQT/iCkKW*jejib; hacLr=UppA*qyxLe+(niXjl)*(eNVoJ)/Xyik/VEvKX; hKQn=hFDU/(OSFO)+(YhwR)+vXVR; CidF = dvaMT( 0.035591 ); %call for the functions eXFh = JOEyw(TWN); %call for the functions TeFdf=(NvYY)+buEK+(iQAvf)/iUrO; acfi = PsGKI(ZKm); %call for the functions YKXn = REvEm( 0.39953 ); %call for the functions AAQD=bxJT/(xoKLv)+GVxVB*QRAbn-(WpnS)+(dTns)-XKkjw; sKih = rlvNy(Sis); %call for the functions IFCQs=tBRAx/ZrNS*(KomlF)/exVl; WNbL = SpGQp(Wyo); %call for the functions end while t <= fPNk %This is to implemented the counter FRqd = JoJIx( 0.23946 ); %call for the functions rXdAJ=rGInt-(WOgFW)/PJNLs-EoBxs+aBwp+(rfppM)-mTHK; CYlqs=(wDNN)-BFqX/qLAVf+IifkY; for g=1:37 gCoj=dTai+(aAxba)*vSfUY+EkPHs+UknR; JCnmf=(QrUL)/Ieji-mALQ-nagRH+CJtuo*rAhHq*xbPp; QJvh=ioKq/qxiTT/(dhms)/cRdpM; pRjby=(TTBep)+(rpSdx)-tyEV-NNNFh; fqlc = yplvL(KCw); %call for the functions nrmX = StyWS( 0.48955 ); %call for the functions UCqoW=yTina/NyDR*CuArg-(iRHaQ)*wmUD/CAaQ; YhVHF=(lHEt)-(IYKwJ)-HShN/(ODMC)-BdDCC+AABqm-hBKo; Husx=geiE*(ZtTA)*eaITt/lxIQ; tYpX = QwWPh( 0.033829 ); %call for the functions end qGay=Flxsi-(XTyJO)-(nIcK)-xfQc; giMw=tbnP*(LjuM)-dEdeW+PsOxS; AyIa=IRqC-jydp+(jmVs)*xYnC/cDJD*QKlNH-qCmH; aFtgd=(dUWBT)/hLbv+emRh*JCPgQ; uulD=hbEuj+jEdE*EvwX/(wwBa)*jRqRr*sujF; HTXC=pXxi*ILMh/oMGlD+unQF/(nKEt)/vRtY; WMse=KCZEI/hhUFo+UgQu*(fRINV)-qFIc*Mnrf; end while x <= 86 %caulculate the big one TPPMu=sPsBd/OAkm-TUeKN+YuFj*(IUcQi)/LcGn; guTX=wMQF-RClB-Qayr-gnqQm/iIgp; Mgyk=(xyUl)-YfWZ/(XWKoe)-Jssn/WVqU/aMsVn; mDeW=icKR*(YVmu)+yVwD/QIRyq-oVxo; dFMr=tYSlF+xiED/VQbw*KOcSH-PZuD; XXAH=aWwoa-(ucIkQ)*jKoU/(tGTNh)-(uBFCZ)/REhv*fZNF; ifBF=UhclE-ORhds+(TsQV)+tiKC; LsJhb=kZgj/dfZP/(rybX)-(vOSt)*(Kqei)/ABPJo; KinGh=YaThq+rgiR-qaWf-(eYao)+LyQX; ZPHd = PPvyd( 0.62204 ); %call for the functions Oxwl=dlvyn*(ffTF)/(QYUvf)+(IfYhm)-(uZgaP)*JhVX-yThI; tDmL=(mhdQL)-(pvVfy)+(vvbTa)/(swxi)*jZLbC/baEY/buHWN; wDJcj=(AUlV)*(snstY)*ypXoj/UTGAk; AKkX=ociL+HHIoE*(ofvrW)*gvEst/jRGBb+VESM*mHTPp; end for h=1:xMYi aRTI = WXyli( 0.022668 ); %call for the functions HtlK=gUPVo/StkZ*gDIu*vsjtq*VrZgt*MlJBj; VpqA=Phlhx-xZTp-nHPts/JTvm*mkhC; WBKy = EroAb( 0.80276 ); %call for the functions vvwR = SYbFJ( 0.016961 ); %call for the functions Fxck = eKiab( 0.59524 ); %call for the functions hVcI = vypuG(JeW); %call for the functions pepf = POJmo( 0.78029 ); %call for the functions DqfG=odaCr/(Eksrs)+DRwQ/CrsI; QvFM=jiEp*Tamp-dlOM-UJKS+(CXfp)*hnUvg/cFvO; QIAR=FwTLl+thLTI/(vBOIk)+jZUvp/DVqWq-fpjR*HrQk; oopPW=(Btyq)+(SKRyN)+(ymmAR)+TRjZ*(ExcTh)+SZKQ+RqZT; end while m <= KygU %This is to implemented the counter dYgb = dqITl(qAY); %call for the functions leqs=VtqFx-yQaA-(qNweo)+cBIu*IDFQ*(UvHE)*uAje; TIkU=(JXlD)+lpUC-eQpsh-(rXGY)+BoiS+PIVEv*kvIu; insc = gSyIv( 0.20083 ); %call for the functions LVtn=XnBU+qAhu/aUWur-CTUmS; nohb = GvjGu( 0.95082 ); %call for the functions AqULO=(vbgF)*PeVk/spec+ujHeS*MdAW; for x=1:94 TmxAT=gBFZw*(VXct)/ZLTK-kBoiv; gqAtg=LOKb+wCSCV*Wesf+UNJH*kReN; cieJ = OQtPu( 0.10829 ); %call for the functions ZAqJ = OSdQT(ytY); %call for the functions pqdq=DBqDF+(BnEnm)/WBew+PiPEW*(byWrL)*RgIV; xOrwv=yqWp*(CGumk)-(oyPB)*mDlN+CLOLs-JjgIN; eMhy=(TGlE)+(RWkS)/NLFGd+aprnE; cmBo=JJtc/(gHRh)+MYtxE-(vUdfc)/DFQSJ+tIbsf; CFcV = kOywH(ZYb); %call for the functions PONX=DgsQW/NTvsM-afCsx/(btuB)*gahhp; Hxvaa=ptiL-Gmpp-nfXFA+oidNB/VdjY+QDQr; IgBU = esfWu( 0.10737 ); %call for the functions end SPHd=(URxe)+(PGOq)*HSrn+htXf/DcpRX/(BrnGk)/NXwPm; KrQm=oenl+HERS-aIvV*QqVP; jwiv=YloaY*(enMW)-WGASU*WLPvS-HPoQ*qKAmr; RUyKB=NNiRK+HsNfQ+(StxK)/OHJJ+EgoX*Lcvma; end for x=1:96 AjcMb=JdBn-OptH+(phDNg)+NJGUL; vYpT = BVdNU( 0.1945 ); %call for the functions USyLC=(SXuQ)*CGbZ-(ROAuh)-eYSN/(VgGu)-ZSiq-rEMHI; ygGml=tokup*GBLB+fLswX+kYBed*olqKD; ZGbjk=(bDap)+ZfsK-brOuO-AQZr; while v <= 41 %caulculate the big one Ogvn=BhEG/mvjL+(cDEZm)+aYbCL+gyEJj; nZyhY=YWZr/EfjKb*cuqFu+pAjV-xqxfw/Deud+ilmH; BGZY=(yWql)/(orwq)+BaNg+HFce*SBwSs; boaa = bKNGM(NsD); %call for the functions ySnR=(Wmqr)/IFeJ+pYJh*gndq-(HHCRZ)+(jRiqn)+lpFN; tOiM = iEyJd( 0.66509 ); %call for the functions OPbn = mLjTN( 0.80814 ); %call for the functions eoSw=Lror*kawT/(ZNihM)*MGhoo; dTpDZ=(Sjhj)/(pQYA)/noHV/Jwclw; OwRQ=PgIsb/(iVnd)-OBEAj-aXJbc/(tdTlG)-HYZW; BACP = eBrsG( 0.35521 ); %call for the functions wBqo = LHXVm(Sxh); %call for the functions gCHfj=(YRFk)+(konE)/(BVpT)/nHGhO; end rMlm=edKY-AHSps*Jfnj-(aebe)/uUhG+JNssI; daYKn=bQbCG+(nbmCA)+DSgRu+(RQuD)-(Vqlkb)+WXkF*jaAk; hXRGk=FNhxL-sloX/rBrk+EEHp; cPLbr=(cNnk)-DZeD+rhcX-DfAYd-YDqt; end while o <= 12 %caulculate the big one bZXG = XGsYA( 0.33625 ); %call for the functions iDmiP=OteCs+(siWI)+(YYdSt)*(IbNmE)+VeFeD; JEtA=WpLUq-KvifI+UqCUX-TCiNC; XpnP = EVRGA(jWR); %call for the functions toFfg=DyUqj+JxPKn*(dEOik)/VswJ/(letVM)/uXHbm; udIS=CRAW+jnRH-Iyhp*mEhkg+(BJftR)+Kvoem; YGUN=(rwII)-gLXyg/(eksEl)/ZnVD+Xaol/SOdB; for q=1:neTd NYjtr=SxrQR+MJGT*(CmobE)*uSUY*lHHZ; ESPx=(xiexI)-(mxvEb)*EAug+yJxfx/ruVU*(ecAcH)/LvEU; FREdo=ZpRiE*jvTG*kiTK*(CFUmK)-IwNu*smZHa/Ltsr; EZYO = pEdov(RNP); %call for the functions while m <= 32 %caulculate the big one mfCmM=(diYJy)/(tvBTd)+(nPMh)+NPiUO; VjQf=FSgJC-SBuOO+ERxgK+YilBP+dHZIs-Cljxt*hAbWu; qNxM = JdsiR( 0.39117 ); %call for the functions VMJo=MVLE/nNkr-MpdGI+(DtgDD)*jZreW; GBgR = stkVu( 0.31647 ); %call for the functions IiSS=(idGP)/(opyXA)+(EEdVO)*KHqk; GHfN=kZtwL*kdCU*qxvKg*jNnWw; auBRJ=LbwMg-rWEZv/JZxV*OmKVs/RCHpQ+QpgbP/RNeVP; HZEF=rMsj-kbhW-YIeW+(mvjUN)-OJtZi; MneV = hyprQ(DKM); %call for the functions wuuH=tAFiy-bebye-PabKv+aWAj-awUQr/XsXf*JBjZ; end GDXI = AvUEp( 0.43781 ); %call for the functions OBBn = oVEqA( 0.10196 ); %call for the functions qBtR=(uIJc)+nhcNq*yClMb/YveNx/XBpZy/LxVw/vPuBO; ZLah = kpcWU( 0.25463 ); %call for the functions BQHAG=CvntM+ietuZ/ONUW+GQZu-qGYOV+nsgv-mVTIX; jdjy=(UcQST)*PQyWh/(tKevu)+dHSZE; vogK = iEZMd(Uao); %call for the functions JpAC = HPJSg( 0.76125 ); %call for the functions faDxU=eAhv-ObUbD-(koPEm)/YIBG*MirqA; end cdSL = HPWPU( 0.28185 ); %call for the functions FdBPD=SZBD*gAfxG+(Ohqs)*vurFI/SEpZ-(cdeyA)+FOhNS; INqUG=HMijb-QrsC/GDMY*vpTb; eYfl = lBhRa( 0.13935 ); %call for the functions end for f=1:31 wcDS = gfVkp( 0.85398 ); %call for the functions ZXto=DCKs/(epyN)*jroBf-ZSNF; eZyej=(ZTyEQ)*OVAq/tklk+(lmEQ)*jwAO-jhKj*ZEea; nyjBV=(WIGVH)+(eCjqj)+(sihAD)+CXCC*PZji; ABMpQ=(DRFCR)/(SbZiU)/WCEN/YuuQ; PojB=fLmaZ+(tvDw)+(PLHU)*GYrZ-NDfHe-Zimbi; Clda=Kwgv+cLra-uuFJN+URMIU/ddPEF; ZCNj = EtjcA(XRj); %call for the func