MATLAB Answers

how do you create a model for a sum of exponentials?

15 views (last 30 days)
carenar on 25 Jun 2015
Answered: Image Analyst on 8 Sep 2020
I am attempting to write a code with my input data and model it as a sum of exponentials. Depending on the curve, the algorithm needs to determine if the fit is robust enough (probably compare an error value) or if it is not to add another exponential term(ie a*exp(b*t)). Thus creating an equation of exponentials to properly fit the data the best it can. I have a formula to calculate the sum of squares error. I was going to use that value to determine if another term should be added.
Anyone know some sources that may be of use to help create the algorithm?

Accepted Answer

the cyclist
the cyclist on 2 Jul 2015
If you have the Statistics and Machine Learning Toolbox, you can fit this type of model straightforwardly using the nlinfit function.

More Answers (2)

Arturo Gonzalez
Arturo Gonzalez on 8 Sep 2020
Per this answer, you can do it with the following matlab code
clear all;
% get data
dx = 0.001;
x = (dx:dx:1.5)';
y = -1 + 5*exp(0.5*x) + 4*exp(-3*x) + 2*exp(-2*x);
% calculate n integrals of y and n-1 powers of x
n = 3;
iy = zeros(length(x), n);
xp = zeros(length(x), n+1);
iy(:,1) = cumtrapz(x, y);
xp(:,1) = x;
for ii=2:1:n
iy(:, ii) = cumtrapz(x, iy(:, ii-1));
xp(:, ii) = xp(:, ii-1) .* x;
xp(:, n+1) = ones(size(x));
% get exponentials lambdas
Y = [iy, xp];
A = pinv(Y)*y;
Ahat = [A(1:n)'; [eye(n-1), zeros(n-1, 1)]];
lambdas = eig(Ahat);
% get exponentials multipliers
X = [ones(size(x)), exp(lambdas'.*x)];
P = pinv(X)*y;
% show estimate
y_est = X*P;
plot(x, y); hold on;
plot(x, y_est, 'r--');

Image Analyst
Image Analyst on 8 Sep 2020
The attached demo fits any number of Gaussians to a signal. The demo uses 6 as an example, but you can change it to however many Gaussians you want.

Community Treasure Hunt

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

Start Hunting!

Translated by