lsqnonlin and true parameters value

Question

Malgorzata Wieteska 2021-3-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/762586-lsqnonlin-and-true-parameters-value

评论： Malgorzata Wieteska 2021-3-4

Hi,

I'm trying to calculate true values of parameters and get an information about how good the found parameters are. For this I want to use lsqnonlin, while finding the parameters in my ODE. I wonder how to get:

a) least-squares optimal estimate of the parameter values;

b) Chi squared error;

c) asymptotic standard error of the parameters;

d) asymptotic standard error of the result of the simulation;

e) correlation matrix of the parameters;

f) convergence history;

g) R-squared coeeficient of multiple determination;

while using lsqnonlin. Any suggestions will be appreciated. Thanks in advance, Malgosia

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

请先登录，再进行评论。

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

Answer 1

Alan Weiss 2021-3-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/762586-lsqnonlin-and-true-parameters-value#answer_638926

在 MATLAB Online 中打开

I can help you with a) and f), but for the rest you will have to look elsewhere.

The examples Fit ODE, Problem-Based and Fit an Ordinary Differential Equation (ODE) show two ways of finding a fit to an ODE" problem-based or solver-based. To look at the convergence history, set the Display option to 'iter' before you call solve in the first example or lsqcurvefit in the second:

opts = optimoptions('lsqnonlin','Display','iter');
[rsol,sumsq] = solve(prob,r0,'Options',opts) % first case
% or
opts = optimoptions('lsqcurvefit','Display','iter');
[xbest,resnorm,residual] = lsqcurvefit(@fitlorenzfn,x0,tspan,a,lb,ub,opts);