Azzera filtri
Azzera filtri

Double exponential fits, but single curves away from the data- why?

2 visualizzazioni (ultimi 30 giorni)
When I try single exponential fit (red line in attached pic) to my data (the part in blue on the attached graph), it seems to curve way from the values, and clearly does not fit. I used:
f=fit(newx(10000:11000), yy2(10000:11000), 'a*exp(-(1/b)*x)');
plot(f, newx(10000:11000), yy2(10000:11000))
I tried to filter and smooth the signal down to a single line, and I also tried cftool. When I use double exponential, it works, but for my calculations, I need a single exponential fit to get the b-coeff. Is there any reason why single exponential just does not work with this trace?
-----edit
I attached the whole trace, and I want to fit to points 10000:10500
  1 Commento
dpb
dpb il 27 Apr 2016
Because the data aren't exponential, maybe, would be the firstmost that comes to mind.
Need to attach a dataset if expect somebody to try to do anything other than look at a picture.

Accedi per commentare.

Risposta accettata

J. Webster
J. Webster il 27 Apr 2016
Modificato: J. Webster il 28 Apr 2016
very steep power laws can be notoriously difficult to fit, since small changes in x lead to big changes in y, it makes it hard on the fitting algorithms.
A trick you might try is (in your case) take the natural log of the data and then fit it to the natural log of your fit equation...
ln(a e^(-(1/b)*x)) = ln(a) - x/b;
so that now your fit problem has been reduced to a linear fit y = Mx+B, with M=-1/b and B = ln(a).
Also remember in matlab that the function for ln is log() which drives me crazy.
  9 Commenti
IJ
IJ il 29 Apr 2016
Thanks for pinpointing where the problem was, dpb and J. Webster. Now I can start looking for alternative approaches.
IJ
dpb
dpb il 2 Mag 2016
Who knows, maybe there's a Nobel prize hidden in the extra curvature waiting to be discovered... :)

Accedi per commentare.

Più risposte (0)

Categorie

Scopri di più su Interpolation in Help Center e File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by