Custom fitting using equation from differential function

Question

Sabella Huang 2022-6-1

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1731360-custom-fitting-using-equation-from-differential-function

回答： Sabella Huang 2022-6-3

Hello guys,

I have some data that need to be fit by the equation from differential equation solver. I want to obtain parameter from syms A and B as an answer. Here is my code that I try to solve.

syms t y(t) A B

eqn_ass = diff(y,t) == A-B*y;

cond_ass = y(0) == 0;

R_ass = dsolve(eqn_ass,cond_ass); % I want to use equation from here as my fitting equation to obtain parameter A and B

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Bjorn Gustavsson 2022-6-1

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1731360-custom-fitting-using-equation-from-differential-function#answer_975685

编辑：Bjorn Gustavsson 2022-6-1

在 MATLAB Online 中打开

Use matlabFunction to make a function out of the symbolic solution, then fit that one numerically to your data:

R_a_fcn = matlabFunction(R_ass);
err_fcn = @(par,t,y,fcn) sum((y - fcn(par(1),par(2),t)).^2); % error-function
res_fcn = @(par,t,y,fcn) (y - fcn(par(1),par(2),t));
% mock up test-example
t = 0:31;
y = 12*(1-exp(-t/11));
% Fit with fminsearch:
PAR = fminsearch(@(par) err_fcn(par,t,y,R_a_fcn),[1 1]);                                    
% Fit with lsqnonlin:
PAR2 = lsqnonlin(@(par) res_fcn(par,t,y,R_a_fcn),[1 1]);
% Fit with lsqcurvefit:
PAR3 = lsqcurvefit(@(par,t) R_a_fcn(par(1),par(2),t),[1 1],t,y);
% Compare:
plot(t,y)
hold on
plot(t,R_a_fcn(PAR(1),PAR(2),t),'--')
plot(t,R_a_fcn(PAR2(1),PAR2(2),t),'-.')
plot(t,R_a_fcn(PAR3(1),PAR3(2),t),':')