在MATLAB中拟合二阶指数函数

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

使用非线性类，您可以获得拟合变量的标准错误如下：

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