Question

我随着时间的推移测量两种不稳定同位素（在同一样本中）的辐射强度。辐射具有两种不同的能量，给我一个衰减的二阶指数公式 - 类似于指数衰减的公式，但有两个项，这些是相同的但是对于不同的起始强度和不同的半衰期不同同位素。

f(t)=(I_0,1)*e^(-lambda_1*t) + (I_0,2)*e^(-lambda_2*t) %(eventually with a constant term as well)

有许多方法可以将数据拟合到此功能，但我如何获得不确定性（例如，以标准偏差的形式）拟合变量（半衰期和初始强度）？

Answer 1

使用NONLINEARMODEL类，可以得到拟合变量的标准误差，如下所示：

fitFunc = @(b,t) b(1)*exp(-b(2)*t) + b(3)*exp(-b(4)*t);

%# leave off the semicolon for nicely formatted output
mdl = fitnlm(tVector, energies, fitFunc, initialGuessForBeta)

%# SE of coefficients
mdl.Coefficients.SE

在MATLAB中拟合二阶指数函数

1 个答案: