Least Square Minimization (Levenberg-Marquant method) of damped oscillation curves

12 次查看（过去 30 天）
Louise
Louise 2019-12-2
编辑： Matt J 2019-12-2
Hi,
My goal is to fit my experimental data (attached) with the following equation with Levenberg Marquant method :
A*exp(-c*t)*sin(2*pi*f*t+phi), where A is the amplitude, t is time, c is the damping coefficient and f the frequency and phi the phase coefficient.
As my skills are weak in least square minimzation in Matlab, thanks in advance for your help,
Louise.

请先登录，再回答此问题。

采纳的回答

Matt J
Matt J 2019-12-2
You can use lsqcurvefit with the 'levenberg-marquardt' Algorithm setting
https://www.mathworks.com/help/optim/ug/lsqcurvefit.html#buuhcjo-options
  4 个评论
显示 2更早的评论
Matt J
Matt J 2019-12-2
编辑：Matt J 2019-12-2
This works a bit better. I don't know if I trust the model enough to expect a better fit.
fun = @(x,t) x(1)*exp(-x(2)*t).*sin(2*pi*x(3)*t+x(4))+x(5);
x0 = [max(S),1,0.1,1,mean(S)];
options = optimoptions('lsqcurvefit','Algorithm','levenberg-marquardt');
lb = [0,0,0,-1,-inf];
ub = [inf,inf,inf,1,+inf];
x = lsqcurvefit(fun,x0,t,S,lb,ub,options)
plot(t,S,'ko',t,fun(x,t),'b-')
legend('Data','Fitted exponential')
title('Data and Fitted Curve')

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Signal Processing Toolbox 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by