首先需要说明的是,数据拟合并不是Mathematica软件的长项,代码及结果如下:
data={{1,12.08},{2,18.6},{4,16.84},{6,15.24},{15,7.6},{24,4.08}};
nlm=NonlinearModelFit[data,a*(Exp[-b*x]-Exp[-c*x]),{a,b,c},x];
Normal[nlm] (*拟合方程*)
Show[ListPlot[data],Plot[nlm[x],{x,0,25}],Frame->True] (*数据拟合图*)
拟合方程为:24.618 (-E^(-0.926056 x) + E^(-0.0777131 x))
RSquared = 0.996983
AdjustedRSquared = 0.993967
AIC = 21.3752
BIC = 20.5422
另外,建议使用方程y=a*[e^(-b*x)-e^(-c*x^2)]进行拟合
拟合方程为:22.3838 (E^(-0.0701378 x) - E^(-0.926438 x^2))
RSquared = 0.999582
AdjustedRSquared = 0.999163
AIC = 9.5221
BIC = 8.68913