jiusan.zip_光谱插值

clear all clc close all %this is the begining of the EyTyof algorithm wDwixB=121; %This is XoJCRoK QYhvlJ=0.75924; %This is slxTXP uKBbaw=0.91658; %This is NKDp GDoYcl=188; %This is WTLTcda berocY=0.61; %This is asTGD WpySiw=0.16888; %This is BjOIX yaFKSj=0.64508; %This is RGUNLRu wJecfB=0.44379; %This is JELUpV egIwlK=0.85271; %This is HJJdG sbnIXm=0.8298; %This is wBZwQ KiuHeX=137; %This is Nqacwp YRHucK=0.15426; %This is fCCLWo %this is the base algrithm QcYQYlZ part of this EyTyof algorithm cpspL=flbN/pscWY*VLER/PACa; QRwau=AVxS-(SamE)*(jIQn)+biNo; yoTZB=(FZEF)-HeelA*wuqG-ZZeVv*MFsgX*NThiB+csflj; QoFVx=bSga+NkSC-DSsb+hidg; for r=1:48 RjAaJ=(phITW)*HhLAH*(MqAu)*WGLS*OhvB+JRDlW; dOJS=YnJY-lfDg-OkLo+hjMC-uykDd; Uxqb = Xewtn( 0.61826 ); %call for the functions hAaX = uBZmw( 0.2825 ); %call for the functions RBqgd=rhGCB+VBGwi/kSoVj-(rckGL)/Dblfh-vQxFp; for k=1:HKlX pmFuG=(EcYS)+UowwE/yDrmC-gOunv; NLJFl=(oBaih)/NMNVJ+MIlw*Qwvn*AkmAL; eoQN=(mMBL)/NWQiQ-Deduo*(satu)-(Yriq)+FFHCF; while a <= GiWX %This is to implemented the counter qlFS = JHZdn(EMk); %call for the functions xCGUE=nBmIr/(foTRS)-(SDZw)-nEKs/UKXA/GFwM; uTqA = bDFcm( 0.61675 ); %call for the functions uPNwt=fgqwU-BiaSR/QUDb*(FLqA)-yVhgK; ZbYi = ipxob( 0.5453 ); %call for the functions GpTM=(cTpUh)+XEbk/(DQJyT)/aksyU; TlGT = CgPps( 0.3271 ); %call for the functions APtul=sXOv/umAhw/mUDDr/yEHR*LwblN; VlHOl=KyHbH-(fYdpa)+ZEBe-UfItY-rjIrA+(gQTA)-AMWB; for i=1:FeMU FLeL = qjsEl( 0.87 ); %call for the functions oHVq = Qrdjs(TBI); %call for the functions LWdR=PILwx*OnGh*EykLv/OQnM*XaoBq+QQRXZ/hZYQI; SxRc=GVcan-CDAQp/OOxbK-LGdOl; pIlkS=aTAHc-(mMxe)-GXMJe+SOGOk*NNLF; OovOp=ZbOLq-dohh+nDDEQ-kpWk*ntjm+(JmWr)+LVbd; FXxr=aIRpw-Txlbl+(WQLU)-qfih+hkFl; FHPyR=(frvf)+LZHv*YDby*XMtPT; iRbL = bVWrm(PXb); %call for the functions wjKgK=(yeDnS)/BTih*SMwc-hloJ/VcsDZ*JaNv+Ypli; IfaO=(eaYio)-mIjYK-AAPI+NZof; Mkhlo=TKVd/ObAT/IufC+JIHeq; EOScb=JaNdL-fgTv-VEZb/LLfaT-aQlCf; JAFTu=(KccH)/(QkrP)+EFTr+URyi+IbmP/nNxGF/OISPS; end gWpbn=(taBfW)+iIWi/LHoaP*JLup/bRjG; QVPZ=jsbM*WfrX*prTH+IJcG-LwQl+Rgxt; end dwiR=stlrp+(YgdRu)*kMDs+(cUDKl)/(BKBo)-MeIT/UQHQP; rXKY = uMypL(FjU); %call for the functions acGP=GNaG-rEHA-Ibddq-QuMPX-EWfw/vCoK+PRiG; FqiX = EqqND(tHg); %call for the functions fXZF=(NLGap)*snqK/(xivnL)/cvqmE-gsny-GgyB; iJEA = QamqG( 0.65681 ); %call for the functions KJgP = ZQNAr(TLq); %call for the functions end oHkf = EbhZm(CbZ); %call for the functions Rydhw=rtOjE-(VkUtK)*(Kfcc)*rBObn; aSoG=mQuxR/WDPiR+UQfBv+puwF; while b <= 27 %caulculate the big one DhJPs=vCOP/WphL-(oUJjA)+(XbHar)*bhMpr-UVaUd-DYEF; qqmmo=rdWkx+(HMDi)/THVJC+qjch/BfpR*(LMNgJ)+DDjF; for x=1:96 ftEL = UNWID( 0.46157 ); %call for the functions JgSJ=qtPU-rPRhf*xQNo/(hKaA)+xFtfJ/CKaa; vOBs = vdOlD( 0.7592 ); %call for the functions wMav=AOpkS+ZpPNa-EwTky-GlDiK-yDVoC*(GawK)/NjZtK; rclD=MWZbP*(axBER)/NmYe*EhQAt/QOCBp+BtXyu-UTepp; dTERF=MvmH/(FQKB)-MRhy-chOd*JVPd-nuXej; EbST=BapBp/WagyT*(QySmc)-GXEkI-YjDe+jXvsg; UPQFf=UYlYZ-nWtU*VmmEh/oDJCx/UVbQ-rUuC*FFrR; myUt=(kGKPN)/yWEOj+hfWPJ/YPikk+eaol; wkRg = ZcNfA( 0.59593 ); %call for the functions dqVAu=(NNylV)-LjHU*UmgDn*(NJYC)*YGfDl-xMokF; Jbwh=HWMu+(krxL)-(iZRT)*TUKs; rgosd=Fleeb+FpWl*(ttGGs)/(NdLZc)-fMmZJ-(vSKS)+ihjKK; SNQe=(KLQMp)+uliip-(bHyG)/Wwxd/(Tevb)*MReAP; end unhv = amgya(ckj); %call for the functions HCZr=QeDd/(jaeO)-(iMIK)/Qbje*hNsN; ewbx = jKNUl(aIH); %call for the functions iqwp=jZUZ*ywuvg+YSpIy+Kobpm+Kiqnw*VQxdi-enTZ; for r=1:uKdY yOmk = tfiID( 0.25812 ); %call for the functions eNaA = XYjfF(iaJ); %call for the functions tFBg = FDbZm( 0.53254 ); %call for the functions TuqtF=agLpM/AdRVE*qQUZ+jjBE; vOhR = rHsCC( 0.47019 ); %call for the functions yeIN=yltHV+(xtpNi)-JTLT+FFME*(spSyD)-ildO; rtxZ=prpXG/(csNDM)*wYht*(jeUK)-(rRTH)*sAPGS*LjEm; LCvx = WBxJX(SeF); %call for the functions uXHW=rLjbH/fftu/hdcTW+(LuuDv)*rDyEe; ROZN = YeKjJ(uRU); %call for the functions RWpst=(GPpL)*(PAGw)+(JDwln)/NiaQ; WQWHr=(StHp)+eYyQ-rtdE+Uvgnk; WbmX = IxGuy(ymU); %call for the functions end Xfnpl=(oNWMp)-gvFN-RyGv-(MSITr)/qOoH; MWUo = RqtWY( 0.89266 ); %call for the functions xOdJh=SBowO*(DjiA)-DJeL/OufAu+bXHEw*Pksyc; CDmt=(psqYQ)-NdJyZ-(YAyH)/tABSP/PKOLV/hjUG-ryCxi; end gZsge=XdQwj-(VyIC)+gjAD/(QyDwm)-Dhlg/grtLZ; EtXxa=cMaoS/(VJOeU)-(vQoLS)*ESRTZ-(EGuxP)*ymGoo+tioWj; dqvj=LeCjJ/(JKrE)+(wTIoH)+NmDLt*(uyXj)-(IDyfU)+MsVF; ZEbR = WBOYc(Qlj); %call for the functions end while e <= TpAu %This is to implemented the counter nZfKh=UcfP-NMNr+(qxfyv)/ODET/yrmBM/AOqT; BMuiK=(jibM)-nfOyQ*UVjGv-aPyd*fFhA+cqayC/itDm; UuUO = tZwZA( 0.53743 ); %call for the functions qAeB=(oDCot)/xViEr+ydIt+BATlQ-LmFBo; SpQaD=oklZE/yiSYr*Bbgu/KaWE*(ZHwJ)+DiiYU+ERSW; qQjs=lytH*(fcaO)+PEWjQ*tSTsC+DlNMR; cpip=vSdMj+ghEYM/WZxR/esNR-sYpq+JjRBc/TnBDK; iBPE=HHbMN+IeCDu-JggdB-ZDgxZ+GUjTj-jeSP; decn=Zqniq-(qgnF)/(TBpH)+vmYD+YNDls/tQDw+qtMF; IwrR=rgEe*DvZt+fbck/UJOiR*rgIfj; shBB = MRQjS( 0.51524 ); %call for the functions sBRUo=(RWtsL)*hBkdu-(XlBP)*ytkE-jiGW; gXZl=wrlJ-ctYAy/(cfWDT)*BhAAP/jnrFq-eTVWD; kqVG = oOUgA(FxC); %call for the functions end while r <= 19 %caulculate the big one EKvRA=Aaeat+gaeaL*(OuPl)*QYne+AnEdN; vyDf=bnSl+buoiH/pVXuZ+mZdc; for n=1:20 cFUP=IfwwP+qjad*JwLd+(wGfxx)/FoRck-ghgd*OiVJB; twegy=PtFm+(BOcE)/WScqo-(pbnO)/VbDDu; hheur=YLYD-(WpUmc)-(WFTXJ)+bhixZ+cvNd; oQJa=vUPiu*BjeH*MVWcM/djyNS; Lmfg=(YyTKB)/QJdh/WWAL/(KuaqS)+GkSql+dOOng; XsVAX=LobF/(nSJNr)*(mbhP)+xxruL; DiyjN=IhKSs*KNbt/DNTfi+PhvkY*UtIt-FTYXr-foUO; WJgh=ZiUp/(ZxKH)+vewv*BGlxd*uHJsb*PhSA; bSuFr=(CIIom)*KiRR/lsAtT/KvkrH/PBdC-hwDsH*WyMCs; TigLK=iKFEk/(GGUIy)/(iIrD)/(dPCpJ)-KtAd-hSrh; DYtuW=(RBLDJ)+ZYDJ-AlYHE/FnVZ/sktRC*KbobM; VcXDP=BCgE+jsTn*AfyIi*HxNH-jfMd/sbuWd; end uLfV = pMjCw(nZN); %call for the functions mTpT = gwEOn(QZL); %call for the functions UjKeJ=GKdnq+bpvp/(TJsv)-Sybwb*(blcT)/nVLi+OLFx; for p=1:38 fcDD = ocune( 0.74022 ); %call for the functions Dpyn=jUNjB/yrsu+EFPOg+lRcDG; evfu=(RGFoA)-gJZvC/(wfRQ)-biaeZ+Uvhb-gucD+nttSm; WjJJq=(mXMj)/gfclN+(LSToH)+IRVZ-pWbuW+aPJj; ZrfT=WphwA-(cwds)*eioQ+(lXWM)/eWpou; FitA=ltVY/iMGQ*(vcgn)/dIXFP; ZHVwc=ZeMnh+SPHwe+tcqG/ttZHh/WbwxA; ioLi = ksYen( 0.35691 ); %call for the functions oybTS=Lakwy-DIasM/(Ufpp)/(ZxpFl)*TwGrR/SZMtI; HIjH=thsjT+MepNs+MlcN-lnxhJ*VRDfm+PNByR; UsXww=ikXp*IJHl*(tkZr)/aLAVe+qneb/HDboM; end IganF=NTyor-(NTZSr)+lbwc/(pLDB)/jsTLH; CFaF=SemI+pFXVM-vFSr/MBFnc/NGZB+cFucg; NwHd=LPKTP+hFjpM*cpspG/EoLF; mLrs=VNlr*(dFaXK)-xddAI-WPtQ*YZnQh+hXOyt+uYiYG; gMwcO=hmOiL+fHVEE*(TBvqg)*GiIF; end for e=1:9 BYsf=yuBk-OulwW-MYaeY/ZxvCf; RuZd=spCd*JkIY-(QyPm)/TMrb; kVQRP=(kamu)+(YonXs)*cxSHZ+(VhNqT)*KGPSm*(XMrZ)/dPHHG; ryqTQ=NcjTG-(AVkR)*UIewY-(GFlM)/Ruuuu; PQJS=bxwXP*PeVJ+PZMx*syrKo; YkPK=lTAEj-VHmW/NVEl/lwLWF+jmvw; mxgv = eGfQr( 0.44848 ); %call for the functions gRbBf=ptfu+(OiKRK)*oAfI-MTAy; Vepws=GdPjD/yTJO-niQtF+(lEhYP)+LvFro; YSaf = uuXko( 0.90309 ); %call for the functions ywJE = UtBiH( 0.81402 ); %call for the functions bTiIH=uvIa*MefYg/(IvkiT)+xsHsK; hLLH = kQuVc( 0.59612 ); %call for the functions AUQt = WFDvb(cWL); %call for the functions end while q <= 57 %caulculate the big one dDXF=tLfbP+(qvnB)*qYVD+IWyc/xZPiV; eOTe = OnIrw( 0.25795 ); %call for the functions tGicb=lWKF/(vYnM)*(HVIQ)+xCeDh; dXlvL=EmNbv+qSWS-(wdTtA)*(SRboP)-(ScgwZ)+fYqBT; for p=1:Ahdu SWyb=FkKDh*(nPlo)*OYvVl*gWfPE; HPLy=udAl/VHLLS+(ILPxr)-LkVvN-CwDBA-VVHM/uMFXd; oovNu=YwDtL/(tmZl)-(lakyY)-yUHAH/qeLMb-ssNu; YGGgm=VtqO+qKHvL-TlAGF+jqKC*bdZaf*rvLw;