blackbody.m

来自「来自澳大利亚Qeensland大学的计算机视觉Matlab工具箱。 This 」· M 代码 · 共 25 行

%BLACKBODY	Compute blackbody emission spectrum
%
% 	E = BLACKBODY(lambda, T)
%
%	 Return blackbody radiation in (W/m^3) given lambda in (m) and 
%	temperature in (K).
%
%  	If lambda is a column vector, then E is a column vector whose 
%	elements correspond to to those in lambda.
%
%  	e.g.	l = [380:10:700]'*1e-9;	% visible spectrum
%	 	e = blackbody(l, 6500);	% solar spectrum
%	 	plot(l, e)
%
%
%	Copyright (c) Peter Corke, 1999  Machine Vision Toolbox for Matlab

%			  Peter Corke, 3/95

function e = blackbody(lam, T)
	C1 = 5.9634e-17;
	C2 = 1.4387e-2;
	lam = lam(:);
	e = 2 * C1 ./ (lam.^5 .* (exp(C2/T./lam) - 1));