(这里有一些理论公式推导：http://eagle.lamost.org/?p=37798 ，比此文更加详细，可以参看)
前一阵被X.L. Fan师兄问到关于尘埃能谱 (Dust SED) 做图的问题。然后仔细检查了一遍自己本科论文时候编写的数据，发现有错误，于是花了些时间重新编写了一下。虽然拟合尘埃在亚毫米波段的灰体谱是一件很简单的事情，但在我看来，其中涉及到的单位的问题比较典型。我们平常在做SED的时候，常常会搞混 Flux Density，Intensity，Flux 等概念，他们的单位等等。并且平常老师在课程中讲到各种SED的时候，也没有仔细的讲讲观测中这些东西要注意的地方。这点比较遗憾。
我把程序贴出来跟大家共享一下(点击下载压缩包)。如果你也是在做冷尘埃灰体谱(following the gray-body function or the modified Planck distribution function)辐射的，可以试试这个程序。程序是用 Matlab 写的，一共包含一个主程序 dustsed.m 和两个调用的函数: errorbare.m 和 fitblackbody.m, 其中 errorbare 是由Gerry提供的 ( 详见：Gerry-云居 )。拟合的公式是：
(1) 
下面测试一下，输入的波长单位是 micron，我取了60, 100, 122, 158, 170, 450, 850 micron 这几个波段的数据，其流量按照Jy的单位分别是65.52, 114.74, 92.0, 91.5, 75.60, 8.832, 1.393 Jy，然后输入测量的误差，调用函数。拟合出单一成份模型中的谱指数、
和温度T以及
（实际上，我在拟合的过程中用的这个参量包含了尘埃的投影面积和距离的信息，但是把含有β的那个因子独立出去了，单独作为一个参量拟合，所以程序里面的 emissivity 不是严格的）。图如下:
给出的结果是
，=1.160865 ，给出的相对误差为0.208285。可以看出，这个Beta明显偏小了。所以应该尝试用更多的温度来拟合，如果用两个温度成份拟合的话会得到更好的结果
， 
， 相对误差为0.205005，这里基于 Dunne L. et al., 2001 MNRAS 的工作，取 =2(有点跑题了……)。
目前 Matlab程序里面拟合使用的函数是 fminsearch ，因为 
很小的缘故，必须用相对误差才能有效的找到比较好的拟合参数（通常在1E-10以下）。下一步要尝试其他的拟合办法，然后在程序中加入 Chi-square 检验的部分。
什么叫相对误差？我觉得一般表示参数的后验概率不是用的“相对误差”这个概念……
我还没有考虑那么多。这个不是严格的“相对误差”，是norm(x’/x-1)的值。
The extinction efficiency
is called 
in 
, which is not right.
, in local thermodynamic equilibrium, is
The macroscopic emission coefficient or
given by Kirchhoff-Planck relation:
Then what are you trying to say
维基百科里说：“The emissivity of a material (usually written ε or e) is the relative ability of its surface to emit energy by radiation. It is the ratio of energy radiated by a particular material to energy radiated by a black body at the same temperature. A true black body would have an ε = 1 while any real object would have ε < 1. Emissivity is a dimensionless quantity.” 我怀疑严格的物理定义里面，emissivity是没有量纲的量，而天文里面为了吧emission coefficient简称，所以才那么叫。
, in local thermodynamic equilibrium, is
,
given by Kirchhoff-Planck relation:
至少，物理上 emissty 的定义不是Hildebrand1983 所说的。天文上物理量的名称都是自己想咋说咋说。。。
不同意。
An introduction to radiative transfer: methods and applications in astrophysics 第11页。
The physics and chemistry of the interstellar medium
公式5.32
在 Rybicki 和 Lightman 合著的 Radiative Process in Astrophysics 教材的第九页，上面清楚的写了：emission coefficient 是你说的那个j.并且国内的教材也叫J发射系数。
Emissivity then, is a fractional representation of the amount of energy from some material vs. the
energy that would come from a blackbody at the same temperature.
这里说的已经很显然了。
tes.asu.edu/MARS_SURVEYOR/MGSTES/TES_emissivity.html