2014-01-23 1154 views
2

我坐在兩個不穩定同位素(同一樣本)上的隨時間變化的輻射強度測量值。輻射具有兩種不同的能量,給我帶來了衰減的二階指數公式 - 類似於指數衰減的公式,但是具有兩個項,這些是相同的,但是對於不同的起始強度和不同的半衰期同位素。在MATLAB中擬合二階指數函數

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

有許多方法,以適應數據到這個功能,但我如何才能不確定性
(例如在標準偏差的形式)進行擬合變量(半衰期和初始強度)?

+0

你曾經使用的gnuplot(http://www.gnuplot.info/)?它返回非常好的不確定性;在http://people.duke.edu/~hpgavin/gnuplot.html中介紹了擬合程序 –

回答

0

使用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