Generate the spectral power distribution of a blackbody radiator
spec = blackbody( wave, colortemperature, [unitType = 'energy'] )
Purpose:
The formula for blackbody radiators is computed for an array of color
temperatures. The returned values are in W/(m2 nm sr) by default. If
you request photon units (unitType = {'photons','quanta'}], the
returned values will be in quanta/(m2 nm sr).
The color temperatures are in degrees Kelvin.
The typical range of color temperatures for visible light runs from
roughly 2500 to 10,000 K. Perceptual differences between these are
roughly equal in mireds which is inverse to these degree Kelvin
measurements:
Equal steps of temperature difference in Kelvin are not equal steps of
perceptual difference. Rather, they are in equal steps of mireds. A
good formula for approximating equal perceptual difference is
perceptual difference is ~ (1/K1)- (1/K2).
Reference:
http://hyperphysics.phy-astr.gsu.edu/hbase/mod6.html
Examples:
wave = 400:10:700;
blackbody(wave,5000,'energy')
temps = [3000, 5000]; sp = blackbody(wave, temps,'quanta')
Copyright ImagEval Consultants, LLC, 2005.