温度与发射率分离（NEM_TES）

function emi=NEM_TES(L,Lsky,threshold) k=1.3806504*10^(-23); c=2.997924584*10^(8); h=6.62606896*10^(-34); x=[8035.8,8178.3,8301.8,8391.3,8502.8]'/1000;% 修改为对应传感器响应波长 emi_max=0.97; emi=ones(5,1)*emi_max; R=L-(1-emi).*Lsky; B=R./emi; T=1000000*c*h./(k*x.*log((2000000000000000000000000*c^2*h)./(B.*x.^5) + 1)); Tmax=max(T); B=2*c^2*h./(x.^5*10^(-24))./(exp(c*h/k/Tmax./(x/10^6))-1);% 注意单位的换算 emi=R./B; temp=R; while 1 R=L-(1-emi).*Lsky; if abs(R-temp)<threshold emi Tmax break; else temp=R; end B=R./emi; T=1000000*c*h./(k*x.*log((2000000000000000000000000*c^2*h)./(B.*x.^5) + 1)); Tmax=max(T); B=2*c^2*h./(x.^5*10^(-24))./(exp(c*h/k/Tmax./(x/10^6))-1); emi=R./B; end