Clear Filters
Clear Filters

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

5 views (last 30 days)
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?
I attached the whole trace, and I want to fit to points 10000:10500
  1 Comment
dpb on 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.

Sign in to comment.

Accepted Answer

J. Webster
J. Webster on 27 Apr 2016
Edited: J. Webster on 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.
IJ on 29 Apr 2016
Thanks for pinpointing where the problem was, dpb and J. Webster. Now I can start looking for alternative approaches.
dpb on 2 May 2016
Who knows, maybe there's a Nobel prize hidden in the extra curvature waiting to be discovered... :)

Sign in to comment.

More Answers (0)

Community Treasure Hunt

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

Start Hunting!