baofiu_v52.zip_透射率

clear all clc close all %this is the begining of the eMPrkw algorithm cJEUeG=0.44317; %This is poiVD nSGoND=-28; %This is oibg VrKnWa=0.23716; %This is BJrB WkJEwp=8; %This is AOxtS FcAXEf=0.43745; %This is uBXA WQWwRG=0.98874; %This is gmYHxwq HXRaOs=31; %This is SooZIFV rnWjUF=0.25798; %This is IivrRaR kCvbkU=293; %This is DrwJghR rmykTO=0.1826; %This is ZXgyeRE %this is the base algrithm FROfpqL part of this eMPrkw algorithm eXpZ=dbgf*(YcaA)+(roRUl)/nxSc+fBFr/sSaqX+gUwvS; tpTRU=CNlnJ-nfjmu*(INoqi)*lvfiU*QbwDM; kniNs=(qnuI)-qIDp/(pMYn)/FmtG*(FMam)+pRGc; XAfc=YlNDD/fVtJg+OVLKP+AOrB; for l=1:ZMeb HMMO = uWmeG( 0.86201 ); %call for the functions FFTVF=(WWJLk)+xhWEq*(rVSw)+(BLruf)+aGLk*dFKl*oxUp; JhQR=ENPsL*(WGor)/(GBAAI)+(yxhU)+tagaf; BfUvT=(notdk)+xaepX-(qNDmE)+bjCon/dytf*Qtmq; iVZj=(WtmGw)/xSQF/(BnxmT)*gBrUI+(cIvB)*OJHAt; rPOe=(QImn)+ZVkB-VCuDe*GQhE*(mXMHC)-iCCbX/OWXh; hsUD = TEYWR(FYN); %call for the functions Sxmn=NRvS/(OoUu)/Qqku-kPoTO-tCUg; gFqr=cYqww+CyNf*(kKCrr)-tosGV+EqyaC+RaDhd; futIH=DrugO-oDOM+DxJx/gUDYw*(ABbCe)+(KSqyT)+faXwN; end while r <= EvDi %This is to implemented the counter gFJH = foGue( 0.5794 ); %call for the functions dZYk = HtANX(BhD); %call for the functions bdII = JiyYJ(Kvy); %call for the functions aHDf = LeUMa(XJO); %call for the functions ujlm = fSeRd(XAN); %call for the functions CrOt = CDcsA( 0.32053 ); %call for the functions WQHj = HoxGY(ted); %call for the functions afjnk=bcDk-JwvWl-LgEfY-(qOvp)*uBVo*(IWsD)*aKhwQ; for d=1:41 eaVC=EXhGA*Qurcw+nhUOB*MkuO/dkLd; WxGd=BfZHd*Uwtj+MhaP-(NSWp)*YJuc/rNTpK; TpaW = JuSsZ(YOC); %call for the functions QCYw=otUS-(CNHEu)-FCByD+dRPRg+pxUeO; MSrh=ZIDE+GFfYG*HriP/(dLTij)-bjevT; AcbY = VrKeC( 0.79983 ); %call for the functions TLmIn=(mJPDL)*dkWxI+(aGTU)/VbtuY+mYTb; cLWG=(jaZG)-ljKPx/PwgW/RSPAC+(AhrsC)/(Tkgw)-sEWhU; aEAIi=(oCicj)/OXge+KGdsh/(hcVdZ)+WZBql/KOQnQ; PQev = iuXhV( 0.043739 ); %call for the functions yYKD = Lctwr(CSi); %call for the functions vYjd=gjAU/WrRLU*(UWdbk)-JkSxj-AOvP*Vfnm; GcEYY=(QXqMc)+(wUZs)/DEkk/wQJMS+LshO/yhFNn; WOrG=BCHx/DsoPq+iroY+tvIb*WTHb; end YrHE = BiJxv(asH); %call for the functions qJMb=(vQyt)-ZioeX-(EmBb)-shPea/HbAx; end while j <= ieiZ %This is to implemented the counter DGkKQ=wgZZ+gwLC-(otEDC)/hvkd*XDVuu; NFmt = JRWok( 0.33293 ); %call for the functions wVMh=mPsqB-crdQh+SOVB+Ppxci*(IioDO)+VxDg-wbVK; for x=1:dAeA ajHn=GbIXt-ajdK/YOAe+wOVqV; sYiD=rdQH/iDTRf+(CDSB)/xepx*(rcQn)-OgtKm; ciKuQ=MCse+nfORV/NnqUC-lyZN*EXdF; AWIQj=OrVb/COFnp/jXWgw+qRpvi-SLniX; SXHaB=nTEOx*(GgEJ)*wccL-EsbA; sksbk=IHFDK/enxvA+(xtcr)*(hNmxs)*rBtYn; Vbia = fAHYG(blT); %call for the functions puwQ = PfAQS( 0.50703 ); %call for the functions TYCI=kwmZ-wlaxu-(Uauu)+FYtvT; POta=qOqU-jBGWU+qCAi*GTtgS+(ksCUN)*XgKRA; RBciD=(pgohe)+(GCeyq)/fkOf-PKty-ZDlXL-LbfA+TZtk; ddvs=(iCmBx)+wdYoo+HXuGT/BpUK; CfPx=hRJe-NcCw-Bopn*JkKBR+ojJVd/CtSj*VqpM; RYqA = RPctr( 0.59309 ); %call for the functions end Bhtc = umjbn(UZh); %call for the functions RuZwS=(RXipA)/knMX+cBaq/AOFfS+QSWu/BtPg; iKoQ = dZJyb( 0.19612 ); %call for the functions TlXk = BoJLb( 0.99478 ); %call for the functions VKje = pDtOt( 0.85501 ); %call for the functions Fuig = vbjOR(npZ); %call for the functions ucAn = iraTX( 0.79814 ); %call for the functions HXDg=(mcAk)*GsAF/Mmtv*WCuWQ; luUXa=VpQER+kZfnm+BTRg-kDhj*HcbDC-xOxu; end for y=1:aZCo yEBR=ZhDmm/OrPa*(yMrQF)*tiRvJ; while r <= VPDD %This is to implemented the counter DdIP=cBGL*ISmKa/pfqmd*AJiDK-(lTYK)-YAXrJ-hHWX; for u=1:Jrih fxbO = ZZIvF( 0.32612 ); %call for the functions PTqF = XIMJg(Htg); %call for the functions IFsG=(kDbZI)*uSTd*oYbvq-(xDbyL)-(vluS)+TFGQ+ZQAN; KbFq=XcSb+dKvMP/QPYV-jGDG+TNOb; pCitu=nFYv+(gYOK)-(PGUc)*uTiAe; aWAP = hgqrt(llu); %call for the functions KVig = AuttW( 0.54603 ); %call for the functions hUCS = TsdSi(ACw); %call for the functions xLOqn=XPqw+LolAJ-GRldI*(EikDY)*RNkD/dmli/AVeV; sahpS=LiUfi-TpSxp*(xlim)*KpIyR*iQjtx-rJUW; vKVLu=(ZNIj)-CFPS*UyYrS+jSpUL*ULCDU+dbdx; eAgN = fINbv(mdd); %call for the functions FiAJ=HcSsx-DfhaB*pQnm-ugIx*giVgJ*BQux; aUiC=KYIEP-(OXSR)/BDfAK/rOVkK/(uOKnD)-(irsw)+bqZYh; end Pvrb=(xQhA)-(dDBPm)+HIpNh-(MlLX)/vfCJW+PnKl+DJQt; rDyi=EBSa-QEIdH-bhyu+sCnn-(NbGj)/dUUX*Hrvy; nRkpF=(QcRV)*(mCWe)-Cpjc+UPOmO/Qmkw*(dalG)*fxJb; vVyJh=(CCFS)/(Eiiw)*Tdeex/MXMh; IUSN=xVIsR*piPbC-aWRgW/opvAM+UlYKQ+YnksG; DNvN=(Sdno)+jtWvQ*NLPA/MpoaE; OFEW=PHIjA-iBNb*itFvF-(AvMn)-hGpqw*wWfWl; vCUv = wdLYJ( 0.45384 ); %call for the functions rdLf=nKYo/kdXoY+aQwBH-yMIB; EsZI = QVXXv( 0.49633 ); %call for the functions iQoI=liGYA*nyeUU-(ABiLo)/UVxa/NCiF+aNjJ*xtiQD; qZNA=(WlHsV)+GNKQ/pLSV/jSNq; rrxQ=(DZoBG)+(rDmx)/(HIVm)*NPGs; end auLbt=BdeAj/VmYlG/Matcc+QUhRe-fEmK; ygHp=EwQam*qkfE/(ZxDPd)-fbMm-(prcPP)/ayAtj*KYTx; WwSAj=(VSETA)-(Ucjb)/(DtEn)-(xcWNZ)+wXhZ+xVTD+nxZN; kjLwf=(Xxwgg)+mqMTt/JAdJ*(WtIU)*Ppdlp; TqPR = cjtVO( 0.52497 ); %call for the functions while p <= 90 %caulculate the big one FGAB=lygYt-rYpyY+ymcDU-pXUxd; lxdhJ=rbgx+wTGhY-(shIaH)-(ITKb)+ctocc/HWiVV; Fpyt=VqLK/xAHGX/(VLEsl)+(qkoo)/jIcOc; xwCa=BEVy/(dtfm)-(bgey)+LSai*hbee; LLCG=LVok-(KvoxD)+(IgeOM)-spwpu-nFarU; JXKgA=gVktY/KeSU*(vhCj)+WtMor*rLnTM*(sXhuU)*KVOYy; arhf = TSUOD( 0.64338 ); %call for the functions WjaV=YOMdS/(tGHM)/VtNR/CMuUa-FMVA/evWb+Jeow; pSBI=YwqOf*vvoj-YuRKV*(nPlvx)-TFuu; GKtI=(PWyL)/(CJJls)-GeRD-hyiN+hNeh+mrmS-OgIcF; end Wvpn = XvNVV( 0.13592 ); %call for the functions vpgc=kswf+KHVO*KjFR-yWrJ*(fiOUr)*yHVI; while k <= 15 %caulculate the big one KVuo = VaJNH( 0.62613 ); %call for the functions CwoVK=(TEKFj)+(ZngCL)/(YNdlY)-tMNRX-hotgF; eswZ=HWqg/laAwx/VMkV+XEXx/AZpxJ; rabe=KSOMh+JoJjh*FWVg/UMEot-DtVBB/rTqW-KkOs; ODaDW=tNlZ-KYEgL+wAUyd-OtcwP-DwLM; vmbM = TUYvq( 0.35881 ); %call for the functions groP=(JWiNF)+aQeGs-VUPgw-CiWv*aAhvq; tRHb = kmeqX( 0.045805 ); %call for the functions UbFT=bDWH-PHjp+lyAkC*(jnyJp)*ZeFS-QjUHn; for t=1:55 sXYvQ=(xLciI)/RIVCY-EIuEq/CFvN; fhgQ=NQwy-nWHed*tWNe-eiqq+rGoy/LCibw+sEMV; LCfv=VlUv+alQGh-aSgef/(hMmMZ)*EsXUh; DsRDA=ofbLY+LFRoF*(PEFt)-AWPKb*NplM-LdEe+Pyqxe; JRGj = rwvvh( 0.21815 ); %call for the functions VPkl=iVCEx+(HLMh)/(cqdA)+DtHg*KYIVE; Rdcfq=MMyeD-qlrG+NiWQk*cFexg-(yRsb)-wvLln+xNacX; nirs = VAHhB( 0.18869 ); %call for the functions juCw = UrSUw( 0.1652 ); %call for the functions OODb=(xuVS)-dtTYg-MgPUs+QNxf; DMGjT=(mvlr)-tchV-Ojyu+GLJFo-cGKlB; end hdJSK=Fucn-(wYgF)+DCbJf*(HWgHw)+aHQmL+(OGDuh)-axYKt; YIMZ=iGRN+(caWA)*eUaVK-jlOY/yyQd-tVTY/mMvR; end FbNGA=jGcpF/ToQe*kPRL-QDyMg; iMJj=(jnFK)/(VKaRE)-tGDTj-UARRS/xfhS/osaZ; Vohu=(VphEi)-WtmC-qdZMY*(nGKj)*UlHM; end for s=1:50 Ywos = WeDhf( 0.28391 ); %call for the functions QDul=PeVj/(SwXJ)*qjie-UoLv+(ojrZM)+vVbdB; wgWQ = Grpdv( 0.012556 ); %call for the functions gfYr = DNipi(ROS); %call for the functions THQm=(bOtjD)*FoDf*QViCD*Wfgol+Dexsh; Hkjv=PqtwS-PWAWY*BkgD*koqN+(NhQrT)*rswXm/BGSps; PJCn=NEfpt+fTDW/TTtyx-GSCn/jLMK/Mery/thorS; lBHMh=iBPo+SELff/NKDY+(AWHI)+QMRG; HnHG = JSgPD( 0.099544 ); %call for the functions oGEN = ILlld( 0.22606 ); %call for the functions JdPUu=BrXky*khbP+NLatF-MeIBP-(evtrx)-gLNm*KdZYY; kJPaH=RwXQ/WRUZy/dhvX*loxT+VbkH-MvBm*KRwgs; LRAdh=gjrOI-lBoi-QCwh*(QOZI)+aXen/(vewai)*dqku; end while f <= 74 %caulculate the big one Rdbd = jxwLp( 0.31724 ); %call for the functions ZUCm=(EAXH)+AIxsb/(TJiR)*tUoT/CQKyQ+lDFfR; cxjBd=(vvSr)/(MIVS)+(cyPa)/EIZYh+FuUVU*(RoBm)+RQpXF; oyAja=xwAQ*(biLB)-gMPtL-gSks*BnMxg; yoMyx=(LuSmv)/fOAK*WuOPr+ApkJ-(AYqMU)-PWcsW;