how to fit exponantial function matlab ?
3 次查看(过去 30 天)
显示 更早的评论
hi everyone i have this exponantial function and i want fit it so i can get the A,B and C values
can you please help me with this. thanks for advance
2 个评论
Walter Roberson
2017-4-22
Is that "+ C" inside the summation or outside it? Is it part of the exponent of the exp() ?
回答(2 个)
Walter Roberson
2017-4-22
Fitting a sum of exponentials is difficult to do well.
I suggest you start with cftool, request a custom equation, and enter your model. Then click on "Fit Options" and enter whatever constraint information you have available.
The underlying routine for this will be lsqnonlin()
lsqnonlin() does not necessarily find the best fit inside the given bounds: instead it goes for 95% confidence intervals. Once you have an estimate of the parameters, you might be able to use fmincon or fminsearch to refine the estimate.
John D'Errico
2017-4-22
编辑:John D'Errico
2017-4-22
In theory, this is a problem solvable using the curvefitting toolbox, or lsqnonlin (optimization TB), lsqcurvefit(optimization TB), nlinfit(stats TB),etc. So it depends on what toolboxes you have available, although each of those tools has a somewhat similar calling structure. You could also use fminsearch to do the optimization if you have no more than about two terms in the sum, since fminsearch tends to be poor for large problems in general.
In practice, sums of exponentials problems quickly become impossible for more than a couple of exponential terms. You need to have great data, with very low noise. Lots of data is good too, as long as it is low noise. Even then, you will ABSOLUTELY need good starting estimates for the parameters.
The problem is this quickly becomes highly ill-posed. Any noise at all will blow it out of the water, yielding garbage for results. If you think you can solve this with n=5 or so, you are just wasting your time.
22 个评论
Walter Roberson
2017-4-24
My tests for the form a0 + a1 * exp(-x/b1) suggest
a0 = 0.000549262354656014
a1 = -0.000133005477037469
b1 = 22.3758394031684
You can get these values with cftool and a custom equation, but you need to apply some boundary conditions, such as 0..10, -2..10, -10..50
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Interpolation 的更多信息
产品
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!