对分子云中团块进行高斯拟合的时候有若干参数，但是不少对最终的拟合结果影响不是很大。最近根据计算结果，我忽然意识到高斯拟合中的权重函数对拟合结果有比较大的影响。

想当然地，一个显然的参数是THRESH，也就是多强的信号被认为是团块，而小于此值的团块认为是虚假的而被忽略。但是实验发现对结果影响更大的其实是另外一个参数，WWIDTH，高斯拟合时的权重函数的宽度。我一直以为这个宽度表示的是权重函数的实际宽度，但实际情况是，这表示的是权重函数的实际宽度和模型函数的实际宽度的比值。

关于这个参数，cupidgcsetinit.cgen里有相关的这样一些语句
static double wwidth; /* Ratio of weight function FWHM to observed FWHM */

/* Get the ratio of the weighting function FWHM to the observed FWHM. */ wwidth = cupidConfigD( config, "WWIDTH", 2.0 );

beta = 0.5*wwidth*sqrt( -log( wmin )/ log( 2.0 ) );

d = ( xx[ iax ] - cupidGC.x_max[ iax ] )/( fwhm_obs[ cupidGC.dax[ iax ] ]*wwidth );


我一度将这个值取得很大，这样的结果是计算十分缓慢，并且很多团块被漏掉。至于这个值应该是多少，我现在觉得，按照程序的默认值取2.0应该就可以了。


转载本文请联系原作者获取授权，同时请注明本文来自钱磊科学网博客。
链接地址：https://blog.sciencenet.cn/blog-117333-471099.html

上一篇：C/C++语言学习笔记（四）gsl的随机数产生
下一篇：IDL的颜色表