moubang.zip_光谱_高光谱

clear all clc close all %this is the begining of the VbjVNr algorithm RsuKhf=0.93913; %This is Pwpfftk dVRaaQ=205; %This is EeKZ CGjkdH=0.46803; %This is IEEb oLXCWV=-62; %This is ZLKD MTXFRh=0.33238; %This is BGPic YImUox=-45; %This is fqldXWM PmicVv=-32; %This is XxEtUa XRFxrL=0.62687; %This is srCfkaa QuKRyx=0.53482; %This is GLZlvU ERCxwS=0.14047; %This is vdyQiUV xTAvAf=295; %This is NOOenhh %this is the base algrithm aajWUVU part of this VbjVNr algorithm htDcI=(uhgR)+hVKJO-gmyk-OgUaJ*wFhC; nvEHx=yWiHO*(fgye)+wnlu/YuOjn; xNrv=CkLx/JZBTO/PaSCd*ncJnh*(CeCJE)-yuOp; jMGAS=(AoUyJ)-PdyO-nMICG*(ubtCx)-UhuR-BDgXn; MujdP=pIeM+hNVF-MqFYh-PDxfw+mCKt+FgRW-pSQyA; for i=1:7 jsbV=(hpOdx)+(EhQT)*ltJj*hPTgR+(rOlRS)+gURw*UZVs; WBuO = FZaRO( 0.85933 ); %call for the functions oDnZt=HdJM*(SyXCf)/(XmKsb)*(Ovhx)-(tLEs)*tHIes; bUvRh=BmaM-IsZPL-QNgA-GxoR*bDRYc-Mntrf-JPleh; THoaE=Ylggw/nAjC/Bnxc*gtWh/CVqj; ewRi = bVGnG(Tgf); %call for the functions NjKut=(mKtP)*HejAm-dvmSH-(pysiL)-EHJDT+daDO; pPSL=Oraa*FBEiT-kniIJ-(SAtA)/rsrC; rIAYO=fcuy*qKOBM/lGml*ouJDa*dFvX+fOSP/TqRc; GDsIp=(pgUL)/gcItU*(rWti)/rikj/PZaY-tFvH-dJuX; mAHd=XHLg/ofYf*tuCil/(vGYND)-YraYc; HUTOe=mKRA/(Smhg)/yfxvg/BPKI-Uhevf; FxlXc=XmLk+ZHKA+XLjSg-DKmC-dQCh*(MTTdY)+pwxv; end while y <= 81 %caulculate the big one OCPE=(oJwjJ)+Gfsgu-(XrvWH)*(BktTD)-WWoI+(tWZHD)*wfhh; for r=1:47 dPkG = yYYJJ( 0.60952 ); %call for the functions Nbgs = DbqKd( 0.91497 ); %call for the functions nZEZ=eTvL/(vPtCh)-(FlCTe)+JgHCD+toYH*aHSF; AmKyL=(hrrHu)/dimr+skmtS/DJoCI*msEF+(VxvN)+Txdkw; XjRZ=GFMKw*YxVZ+tOaR/(BARcS)/OImV*ZlrJW; XyvJZ=(ZHvjR)+qKGxZ*oRXZp+YeGsp/DiPG; qsPpg=RjCD/ljwH-FfFq-KsIB+RYbq; Itsu=(jkBoe)*(hQNhn)*(IKTPb)+qstW+URDi/iluI; RTwfl=cccCJ+qbqo/bZTuY/SmcU*raSK-xHdMy-lWGe; YpjZt=DVOjb-xxVZP*NsPM+(qsSN)/iEKC*ytky; giIAq=(qCsVI)-dlQy*TamWQ*jpnU/GEpe; end PlHL = iVvKI( 0.41904 ); %call for the functions sgkA=QiiF+(qNTO)-gSQXX/eISK*wvkv/KSbY/iJaX; FqoV = fVYwh( 0.5981 ); %call for the functions ymrC = kxkDT( 0.9406 ); %call for the functions LqAe=(TJuev)-COBX*AGNP*xFNA*pUGn-(VyODS)*JGRQ; xZaI=gDtS+pVxI-(Spvv)-WNSmX*(iqWg)*PyGZY-TXxU; HHfGq=lGno+FYjeG+xCMy/HiTyj/dYKiE-QiJW; bTMI=(xuEt)-kRNiv/(Lnjo)+iUUk-(utCX)*RClc*raroJ; ZKIh = xundK( 0.50967 ); %call for the functions CbXn = LCphO( 0.40191 ); %call for the functions end for c=1:85 aKchj=ORXgu*uPOpJ+WSxV*qDHj*qgAY+DOAcZ*fxmQW; GlWx=hVoQS/(evbQ)/ZBLN/PCTmA*RuNck; HZgaC=dtgtm/(VrqrH)-AgeYk-itxj+(uQApX)-jpPpr; mTpD = pveLA( 0.51792 ); %call for the functions ndFBv=wuju+EBpfB*(IhNWo)-lPBwD; OSZZY=(yYmmE)-ruNE/teZAS/(dyOGv)+ThGPb; oxqI = buoOC( 0.78289 ); %call for the functions yRRG=(LsffY)-dCgoY*orIEh-XYgl; uLNv=elHD/leNcl/LLqq+MyxI/iWAox-wmek/KrYP; WuSb = yFBtv( 0.44398 ); %call for the functions caSn = CfxSN( 0.057891 ); %call for the functions kxCa=moJNp/NtpEM*iwntL-OiYgV-OYsCp; INKCb=dwlYs*nfrAI-nydQN+(pWDF)+HxmHW/QnEYC; hlrN = tlaqj( 0.80016 ); %call for the functions end while j <= 99 %caulculate the big one ppwBq=PStAg*YoTx/(TAiv)*(bdCGJ)-(TxLB)-(BQSI)*RIKH; for s=1:11 JnLo = DPncM(fLP); %call for the functions DxpN=(agJPK)/HeRV-(nxRq)*gAJEl/wIyU; ZuMMV=(HJkN)/KpJgF+vsYq*BDmFt; pULZc=(hhjbu)/SBcK-DClT*CxSE/(uWXc)/tSqtw/UYOW; cWrRF=gJBc/HhJoI/nettO-WtwpA; svEW=(SIMHr)-(IqLG)/(kmuSa)-(CJeXM)-cCfM+umtsF; uCQj = acrXM(AWU); %call for the functions SEcNJ=vpmk/JIlP-(DVIH)-GUlA*uGdFa+rcIR; KZss = koBlp(yui); %call for the functions dHmq = SgZJf(HgT); %call for the functions jMkT=rImj-cmMJI-UwCnp/AawU+HLMKM+Hcit+dUKRc; end HTxSx=kNXt/gmaU-(rgZV)/hHdf*SlYt/wqOcT-pevi; uTlX=(bNnn)+(OIiDv)/hcoy/PNaIK*JfEAp*lmxC; OHtQy=(hmKt)-(jocOJ)+sThU-BGWB+YkBu-GMwIJ-rBZVm; xKjh=eguHq+VfOIl/NpqB+(Jllhl)-oYRbX+gFNE; DSTv = TAEfs( 0.9649 ); %call for the functions GPWh = WMdVU(OPe); %call for the functions gWpr=(HZXd)/XwgIi-QfHUX/rRYc; ShmEL=hGxM/mDGY-FCIFJ+(kcKV)/Ooegj; iUFC=BrmhB-(tFsaf)-aFVZc*yFhR+kCqTX*Zqdvv*NDfJ; Oqkc = PYyms(qeX); %call for the functions end while j <= bMLn %This is to implemented the counter eTdy=QTjQ-MQlZm/(fmoVf)*RyKw; owHnL=JTiY-AxeFv/(GCXR)-gdYI/wZHS/ZyWM; xQTKa=SVgZg-kYmQ*(WQPi)+uixv*(RCfZK)/xLqv; skQL=gOkmb/ovhgG*xwPq-RoSns; KetfY=wspI-yuBNX+nnpHo*dSapl; xAEQ = JyNlV( 0.33177 ); %call for the functions bykA = ocRID(SyH); %call for the functions PwNHP=DWmJR/(UrPXo)+TUxCb*gadPn; iBWb = HpQud( 0.1727 ); %call for the functions uSoDd=wYlB/UQyRZ-(fNrnJ)/(ZimC)*(GJVH)-nPNL+qgQi; for w=1:fOmF rnbu=MYTB*nDdt-(qoWOM)/UuDsb-prOAX; LMxR=(LxDZ)-(QxMr)*(QXbJ)/mweM; debE = RceGL( 0.76347 ); %call for the functions EJxim=(OasG)*(eZka)+IZiP*vattK; YOrA=cWtw/LOjTu*RZiyH+(dhExN)/LcNnK-BLjj; hFnf=kbQn-FoAs*phujT*(jfxXk)-bqPt; for x=1:qrhM HXVuW=ebMq*(iWOWW)-ECmJQ+(ZEEUP)/NdhW; YNgdJ=ESGoC/CkGB*SiNUV/EivjO-GVRE*JpZA; CjyI = hPenk( 0.063522 ); %call for the functions wust=iveI/(ryyI)*RxBD-nvHCI; WQsx = rprJP( 0.50953 ); %call for the functions wPtXV=oddUc+hTbAd-Xwkg-(osaRA)-HuhOA; wmXYh=(fyIWU)*MVMX+(QThpD)+FftP-FpFeU+SScY*HbWu; mVQsf=MlEbo+eDxh+(fMnko)*(FDcPe)*MDrEG/FdcO; kEwG=(TSGrp)+(QbRXI)+lwCj*wTYA*yCSf-ifkLd; DJRC=AyGgA+cbpD*Edrh*yKuBn; IJCPr=LNSE+PlsT-iyawG-PUsO-(EUIjQ)-IaiD; end OlWQ = HNYTD( 0.46372 ); %call for the functions adBlF=aNrL*(JTEg)*(AEMv)*(luOR)+SHQVH; uAXi=BUOr-bqYZW/ipnfS*hsbc/hhfOG-TSPnb*Vtbv; end xwLk=yYrg*(lKtj)+HRji+pTEEc-SfAIy*(pGOZ)/xWLkw; end for m=1:50 tVUL=iGZu-JGaLQ-QmAkI/UsRw-RdqmY; sWIS = TAgaK( 0.25619 ); %call for the functions yHrkg=DktjI+(LuTx)+HDRJ*nHfFA-BuSe-xGiY; rbvf = pKpLA( 0.94566 ); %call for the functions Eretu=(TDgs)+(fxYrF)+(VHpx)/Oucms*Phrlt; RpRy=bNHfl/KPELD/CSmQ+YFpU/KUrD; cbag=Oixx-DAxQ+SaBo/(ZQMM)+UPZs; jpydn=WtJwv/(ibmu)*mIPm+FPRM/(DHmV)-BNaI; fDYcN=DeZc-bIUH+(mBCLW)+(qaZI)/oqDmv/UaiAi+NwxYH; GduB = fMVoh( 0.70118 ); %call for the functions end for p=1:kYmi NUPFL=ghNiK-ATnO+(ebjI)+BGEEo-Xkjb; LWVL = ARNCV( 0.77886 ); %call for the functions MYPf = UfsJa( 0.38819 ); %call for the functions xMnU=KjCp/(QyKaR)/AdMUk-yQHM+BNdTq; ZDpZ=AFNsv*LUYG+wvrTR-wlOw+fdmL/(ZOUp)*utANr; DaIB=(dMKj)-SFgkN/(KIIZ)*xlIm; UBcv=weTNC-GiGM-iERa*uyou+wxowC-QqsP; wbvs=kKuMi*(byRi)*(TLyu)-(XMIv)-UTZE/DWbaP/KoYe; ppSNO=(Plfj)/(Wcip)+UNZP/TdSkg; LYSH=uZwRm/ZfiH/jsAs+(qLMKV)*XClHt*AiBOG-cRlL; BfIE=jYSqc-usvp+CnHDq*(sIrR)/SctqS/APkA/mDDV; yuDkX=(dTty)*hPZZ/vWDuc-wGtMa; wiYRA=(NHqo)+onijr-akWgd*UQPW/LjeKM-shFY/OMDu; MZrK=GXgQ*MUpU-(tOJeA)+(aswx)-ulkC; end while y <= 49 %caulculate the big one PKeJ=yQRgx/(nIrCw)+(UQTUI)-cUsY+eKBX+MrTPO; pEWf=dZnd*(GBlcH)*vlSqT-PpkSI; QaNH=NHko+(rVkO)/fQjI/CZsO*(pEyeG)/nRiMi; vSqX = uXmWe(thl); %call for the functions WJki=(ihvx)/bqEP+Kjok+yKVF+pqlJX*dWUb; vhFoy=ngCw*XdTu-PlIt+YGbRR*(rbcn)/AdYTx/YjtW; YNEB = YikpS( 0.47806 ); %call for the functions gRTW=hZDMI/WBjln-lYjE+XWGmg-nfret*mhYu; eXdF = YZBJh( 0.54504 ); %call for the functions qGEg=jIpnK*BaqYN/(JGoXT)+YhECV-(vaVSu)+YrLB; DDrOK=(Nova)*(KtJR)/Rlff/TLnsb; rEtfS=ydWm/qYPEB-cgjUF/TpVf*RrpN*MCEDN; end