2006年2月22日星期三

用gnuplot拟合Gaussian曲线

在gnuplot里定义Gaussian函数
f(x)=a*exp(-(x-b)*(x-b)/(c*c))
这个Gaussian函数的宽度 FWHM = 2*sqrt(ln2)*c = 1.665c 建立一个参数列表文件gaussian.par:
#three parameters
a = a0
b = b0
c = c0
其中a0, b0, c0是人为设置的初始值。在gnuplot里运行
fit [xrange] f(x) 'datafile' via 'gaussian.par'
上面[xrange]替换成相应的自变量区间,datafile是数据文件名。

没有评论: