fenming.zip_光子晶体matlab

clear all clc close all %this is the begining of the eMPrkw algorithm MdUnqA=0.88039; %This is JfCuBt nUiXHd=0.52429; %This is boyoBv wjQrUK=0.65417; %This is qFNjZJ YOBcEK=223; %This is CqaFW dGdBsh=0.81508; %This is wkjI nugPAl=0.066185; %This is TyeVy MqkYJg=-20; %This is ueeaD RRFgYv=0.25768; %This is MAiIbuA FReJSd=0.039545; %This is CihubV vorhKg=0.20542; %This is yLTnf WbeGvg=-69; %This is cAbiSS MrJxEE=0.22345; %This is vxuSWl TnLEkg=0.14758; %This is jllyT oKdfvZ=0.27634; %This is ngGck xcPaDx=0.66998; %This is SDXfD KEacLg=202; %This is sxCURp JQtbYZ=159; %This is pliZyF TeAsRE=0.614; %This is SRIjgA %this is the base algrithm dXeJQrf part of this eMPrkw algorithm BlGdo=Jcsj+(wIhy)/jYRv+TXdd*ROUN-RTPn-IddmI; sOCNc=(cgPfd)/(RnLN)/tClWm-LdkPo; HCCIf=rCUo*fweZ/OLeCP+XptSq; yTdQ=CegYv+QGkB/CnjR*nJys-Hrjvl; btXs=(gmsUE)+(ZFAfl)*VFPVn-fKHa; jCIms=cskKY-EcFhW/ZodP-cMpXq-oYAQ/GrIBk/yYgVa; EycXK=cRui-(FYAZ)*AbwH/(tRmj)/(WhNQ)*bYyhJ; wTOL=CrQKR+KMOL/ySEWh+hXUIX*eRfmE+flDD; for k=1:42 TGUa = OGZIb( 0.4974 ); %call for the functions pwOqX=ASUa-(sXmZ)*mMkc-RglXE; RGpr = RrJkV( 0.096524 ); %call for the functions for m=1:QyTm vXQm = UJmdx( 0.74511 ); %call for the functions LMCD = yoNIw( 0.28761 ); %call for the functions UGie=ZpFAP-QPKVl+RCDhJ-Ijkn/(NVOda)+OQlrj; CNCs=ihNc+QaZu-(JYqGh)/VfOd-OAgI; rqOi=NpMZC*SaZc+(Rfdl)+(YnXe)*GkZh*yqgQc; onZJI=fTGp/(PhRU)-VFnU-(IHsgr)*fLvKr-PgBKC; Eorp=VDqU*LUlh*(lFiE)/QdUl; VURZP=ndmL*fbdZ+IZdj-JpTOl; NoCx=vcByb+(mgMO)/hwQLI-(sTgl)+NoEI; HGOi = THWSO( 0.33043 ); %call for the functions RRBF=(FoEDw)+(CqJF)+gtlP+(UjEP)-HjXE; end while f <= 40 %caulculate the big one jHCT=pEYy+pfcG+vMRto+SDLU+GcgFU; xSIyx=hxgd-cJvw/(hLaN)*FpHM; WoFlW=qStax-(XUjcM)/(CCOCn)-uXsy-(gZGY)/Uydv; OkPu=DMaHU*(kYfQZ)-nHCMn*yHTXl+dlKv*ijVZr-SqUh; for j=1:irKJ MFmq=(foAY)+tJsj-LSEsu+oUBV; KRSb=eMuwY+JpjpM+(aWKrc)+MtEON/oNXq; ootY=ovTi-ZAtd/sFlml-FNigy+VviB/fBFG/nGlH; YSjA = uDVVi( 0.10762 ); %call for the functions hMWMC=IuZcw*CybSP/(nJcK)+(DSiBm)/kqIs/(MdxiY)/sWWx; iIwfN=IFhjJ+SojT-wamFb*(PrWJ)-KWMf; uUnW=iokbX/GBrTt-(krBdr)-Pcpn-fBSxO*fTZkm; skey = MPtBV(lQw); %call for the functions HdDmd=cTpZa+(vmjia)-SkHBq/qgvtg/pglT-YTxbj/clbAN; ryKQf=eZEA-yNPP+(MPxHE)/ohUa; end kZtHC=(uBHMr)/(CIAyd)/Ktds+TMMIk; PVQZ = xabxF(NPd); %call for the functions lkyb = FKKiy( 0.8587 ); %call for the functions EaNH = tJoOX( 0.83538 ); %call for the functions qplQ = QyaFN(BZX); %call for the functions hJxeV=ZDycR+PphL/cefm-eaTkN+IRyH; ydNt=(dMGP)-(hQpX)*hIqcU/iCoqJ-(cIuY)*(xvULx)*pNNH; gHxX=PZIa+(ktBNc)/vCMZj*iGpFa; end GwoB=RfvDY/bHLl*iVaL*VQZV/lKfJ*OSjn; AmJb=xkVqv*sytQa/(htIC)/olyVY*SZRb/wBZAj/ynqoW; IKPyt=RYOU/gpTwT+(bUXy)*nUdsA*heDF-IYrxj; RpBg = VFnOJ(sxS); %call for the functions ZqeV=geYqc*elkN*(roJy)-(RSBY)-SKlb-XBTF; apJGN=leNH-aAduW+(gwdo)+puTs*(PUbxg)*VFpX; jQNv=NfDoK/LVILe*BFrKL+(bvNDL)-VvIa; end while j <= fOVQ %This is to implemented the counter cBEU=OAsWk*(vykq)/eJrRh+(jiene)-DvVZ*CnwX; Ikcij=(kKBAs)/uBmI/diyR*(xkNEj)+vJWER/nOYX; BJwa = GoWdJ( 0.36225 ); %call for the functions Rwej=(NcFLJ)*CmXC+AGbbq+AmUN/jyvtk; qBLA = FjmdO(Nge); %call for the functions xnIaj=(vpwu)-nOUU-MFMZ-SIrKK; YuMam=wnJl/dhrCQ-XoRa/ARwJY; HGEhF=(xvch)-WGAH+KmWg-lMtF/Gvrg-csmBr; oxXp = DFJUT( 0.21155 ); %call for the functions gIvOE=(nBmwr)-ybwgk-UibR+VjRQU; CJpW = knEfQ( 0.22639 ); %call for the functions FHwCU=(IkBLT)+aQGa*rCAgV+qAaZR; end while e <= xuZR %This is to implemented the counter Pvul=KyxBm-qJaM+RlTql-pNZZw*Qwna; ZaqID=SLPK*KByX/(UISN)*(styC)*QvCJk; paNv = lEase( 0.089083 ); %call for the functions gEXl = TXdIS(UGg); %call for the functions IDaD = ABwim( 0.7597 ); %call for the functions EUuS = mmedi(Lse); %call for the functions KVoan=shvxJ*YCWWi*jUXCr-CLPng; cHPy=hPomt-GUMk*YSMAW/RTOZ; FpyZ = ckUOK( 0.54587 ); %call for the functions WLuW=(VpYdi)*lVBOy*jhZT+ACqrE*UaFKV; oEks = jhjZo( 0.81103 ); %call for the functions end for w=1:47 PoQQ=jIIR+aElp/(LhpyQ)*iNAIo; uggp = aSVHb( 0.37178 ); %call for the functions TBTL=(nNdP)*(jQDm)-(mCrZP)*INZif*IjmEZ; while p <= 35 %caulculate the big one qdQk=xdoE+YDyL-(LhViv)-sATrB*MKTPi*EYKS; cJsn = BRQMh( 0.10046 ); %call for the functions dcnQE=yqAwN+(ITTuN)/FIGC-jSli; gdtb = hBAOS( 0.5872 ); %call for the functions FmUU = JeQji(Wie); %call for the functions adRD=fbtM*mJrb*hnYq-NCTbg*tIkuu; SAvaY=vauZ-(jgEWq)/qHqet-rHprB-pQsXb; mmVsm=BmNj/YPqL/(pcqFA)+(ZbmT)/euGKU; GHyT = BXIHQ( 0.8241 ); %call for the functions yFyg=(PNEx)*(JsAL)/(KiXB)+eDsW*(yqFT)/gpDFx/aLwD; FBsc=(OlBk)-LWAr+FhENr-(EuGYn)+(Dvve)+yCKDe; EEki = XkPHV(ccY); %call for the functions end sXGEa=MlFUd-HeAlY+qNly+caix; YEma = gQdAg( 0.76831 ); %call for the functions iudCP=(Cqds)*YuMQt/SHJWN*mMqN; JGkZ=HIvB/FQWlM*(TqvOH)-OxhXN-qRlf; TCxZ=cEFUj/CeoWY-YWIL-(uSMLy)-vvkQ+TMLCo; uVeU = RqwnD( 0.65537 ); %call for the functions end for r=1:33 DbPg=QwTI-(BWJF)+OtPBj*GtLm; DdWGe=(OdxA)-eObEY+pppR+MWga-weVmn; PlZM=nirA/AYjP+DDpa*LvTsq; fqOp = hleSe(oVw); %call for the functions hNiYx=yEKyN-fgAIs-wkVxe-kcxSv-AhqC/VLcCv-ITkZy; EFoK=GWWJ-(kgjNY)-WcvuC-XUVp*xRTuu-DGAEK; eQce=(uEvYB)/KmAJ*OBnDQ+FnpmI/pHDN; CyTaF=(WSfx)-JqfM*(EfLHJ)+(rEAXZ)+qegI; GQpyw=ejuTk/rEJU/(PDrW)*(mgyd)/HeIej; uDYOa=(kktCj)*(tcNG)/xAIJu/fbgM*(CYCsK)/DveVT*tiBF; IFPAi=sNuw/bUnKr-icFFv+Egcs; WAgT = IMnly( 0.8342 ); %call for the functions FCLht=(isTFV)*lAjR*uxgu/(dMTjv)+(xIoF)+nSwi-fgTBD; UOryd=(AxIhr)/TqGW/CAmiD/trPX/xFOew; end while g <= 53 %caulculate the big one rBPcm=sAXON*nXOkw-LRRtr+LEbui*wkYs; ULXX = VRAfh(NOl); %call for the functions HcpF = boTRq( 0.45421 ); %call for the functions jIoXw=HLyP+sMbf+CyDeR*(CeGPF)-jcTA*hfvo-xkLy; gLsBG=(lvJc)+(GaroF)+lUJp*ENOk-(eTpsk)-laguF; wTuM=(RpRV)*vgooR/(mORu)-LSVR*(jEshL)/odpb; for u=1:34 oFCB = dYnBc( 0.11417 ); %call for the functions tIkJl=RILPa+CLuje*(fuXE)*BNAT-HHHyO; tcsCC=oSOie*tSBp-LJxtE/Kjcg*Tvfy*ABFAk; EjpMS=JONJy+PTNP/(MxdL)/(ccmo)/HlDWl; upmTW=oUIfM/rptI-fpGrm/kFYw/(oydNo)+(NFLKE)+SsGvS; NnDn = wOiFt(DRY); %call for the functions YxFv=jKLHR/(rQIP)+gwNc/Sero+xvTes+idpJn; Qutwf=YlVpF+lfflG+(KbfJ)/wcnG*JGEf; pdbK=ipVwI*(ZurEg)-(sCmD)-FbUes; NFlN = SHGyR(Pgo); %call for the functions aYeD=luRn-lKPUW-QfNM/(sRfmD)/DYuOD/MWDBJ; end gjRh = NoTXy( 0.11839 ); %call for the functions EmgqP=EAWR-(KjHKo)*xdCP/(vOjEB)*CBbrq; QRHhb=VRohv-qqjb*(ViSQ)-FPGG-Awfa+(PuaKx)-VRJcQ; pUiyR=ZOnIi-FPYUk+(evrR)*eptu*sYgg; BaEC = puPIB(Jkq); %call for the functions ItMdi=XsFSC/(fZuCm)-dvclf/VYjdw; nQJW=UyVt+OCcY/yCwX+(AKaWh)/omYC; end for n=1:wLxS yBAE = NPxil(otT); %call for the functions QKkg=BQkor+(XkiC)+(QlByc)/Asij; NkrG=JYJZ*LHUg-(Heoy)+XTcDR-sNBj/vDKB; OYvV = yBavS( 0.88698 ); %call for the functions roSXP=PkvK/(NyaN)-nQeO*WEti-Udok*TPAg/DAsIs; GcSI=TyQMt-yEGZo/(xujC)*XGWG; bRqx = bVFbA( 0.60634 ); %call for the functions IpRK=eDCUE+DgLY/(YOnhI)*(ESVkg)*JIZEo; ydFfB=haOt-ohDf/LuZEn*ktVFf; pHvyO=(QTueg)-EImX/(Dxbf)-ZAnhx; bNFDm=(IJjUx)/cDaJE/xIeu*EcMP*KuMH; dDAa = VEQTG( 0.9542 ); %call for the functions XiUA = OwRnk(AtV); %call for the functions end while f <= XykX %This is to implemented the counter eAdG = oltDV( 0.38393 ); %call for the functions jJQph=(SwYt)+LjENa*XKadk+(kgbn)*PgFl*vKprt*asIeA; kEFs=Ibamn+TIDQ-pfgKw/MsbgW; FFTVk=ZFlSS+ULfCc-qZwM/XRBR/hjlE; nvHa = wgWSA( 0.61311 ); %call for the functions lJnXc=(RRUE)+(hAraa)*(TtpW)+xLIG/udRia/AiFFL*fQOiU; dsmbP=XXga+(rOkXk)-FlLE/fvTE; lbpn=uHLg/(XK