NGC3603中尘埃颗粒.pptVIP

NGC3603中的尘埃颗粒 新模型 内容 拟和: 检验 模型 方法 结果: 温度 尘埃密度 Gas-to-dust ratio TODO 拟和检验 卡方检验: 测量次数50,样本?个体?(BingJiang) Cash检验(cstat): 中心极限定理之一,大数近似,复合泊松-高斯分布. 测量次数不确定(MSX,helpdesk) Ref: Ciao(Yangchen),xspec(Yangchen),cash et al Xspec manual及附录, Xspec作者主页/~kaa/head_stats/title.html Least-square. Goodness. – 没有好的goodness评价方案,拟和时舍弃sigma数据(errorbar). 拟和模型 4PAH+blackbody 考虑到21um附近没有PAH辐射,令c=0,得到3PAH+blackbody 结果都不好,怀疑mathematica的结果,局部最小?且有警告. 换其他软件拟和?编程找到合适的a,b,c值? 4PAH+Blackbody Q~nu 拟和软件 Mathematica: NonlinearFit – ?分段函数 Matlab: cftool, Optimization Toolbox - ? .. GSL:/software/gsl 1stOpt : “世界上最好的拟和软件” 算法最多(UGO)!支持全局优化!不受初值影响!可设定参数范围!支持代码分段! GSL代码 #include t.c int main(void) { size_t np = 3; double par[2] = {1.39, 0.8}; const gsl_multimin_fminimizer_type *T = gsl_multimin_fminimizer_nmsimplex; gsl_multimin_fminimizer *s = NULL; gsl_vector *ss, *x; gsl_multimin_function minex_func; size_t iter = 0, i; int status; double size; /* Initial vertex size vector */ ss = gsl_vector_alloc (np); /* Set all step sizes to 1 */ gsl_vector_set_all (ss, 1.0); /* Starting point */ x = gsl_vector_alloc (np); gsl_vector_set (x, 0, 1.0);/* U */ gsl_vector_set (x, 1, 100);/* v */ gsl_vector_set (x, 2, 20);/* t */ /* Initialize method and iterate */ minex_func.f = my_f; minex_func.n = np; minex_func.params = (void *)par; s = gsl_multimin_fminimizer_alloc (T, np); gsl_multimin_fminimizer_set (s, minex_func, x, ss); do { iter++; status = gsl_multimin_fminimizer_iterate(s); if (status) break; size = gsl_multimin_fminimizer_size (s); status = gsl_multimin_test_size (size, 1e-10); if (status == GSL_SUCCESS) { printf (converged to minimum at\n); }