How to fit experimental data on only one ordinary differential equations (out of multiple equations) with multiple unknown parametrs

9 次查看（过去 30 天）

Learner 2014-6-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/132238-how-to-fit-experimental-data-on-only-one-ordinary-differential-equations-out-of-multiple-equations

评论： Johan Sebastian Diaz Tovar 2019-11-21

Hi,

I am trying to estimate 3 unknown parameters that exist in 3 ordinary differential equations (ODEs) by fitting it to experimental data of only one ODE.

My equations are:

dy1/dt = (c(1)*y1.^2) + (c(2)*(y2)) - (c(3)*y3);

dy2/dt = (c(3)*y3) + (c(4)*y1);

dy3/dt = (c(1)*y1.^2) - (c(5)*y2) + (c(3)*y3);

My experimental 'x-axis data (time axis)' and 'y-axis data' for dy3/dt is as follows:

xdata = [0 1 2 3 4 5 6 7 8 9];

ydata = [3 17 9 4 8 2 1 9 2 4].

Starting guess of unknown parameters: c(1), c(2), c(3) i.e.

c0 = [1, 1, 1]; where c(4) and c(5) are known = 2 and 1.5.

Thanks in advance.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Learner 2014-6-4

Please anyone help me.

Is it possible to solve this type of problems?

请先登录，再进行评论。

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

采纳的回答

Star Strider 2014-6-4

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/132238-how-to-fit-experimental-data-on-only-one-ordinary-differential-equations-out-of-multiple-equations#answer_139484

Definitely possible! The easiest way to describe this is to quote Monod kinetics and curve fitting.

The strategy is to put your differential equation and the ODE solver call inside the objective function you want to fit. It integrates the equations with the parameters passed to it, and returns the solved equation as the output of the objective functions. You can specify that c(4) and c(5) are known inside the objective function file, and only fit c(1)...c(3). Specify the initial conditions inside the objective function, or you can also fit them as parameters if you wish. (I’ve done all those.)

With a vector output, you can use nlinfit, lsqcurvefit, or with an extra step, fminsearch to do the fit.

14 个评论
显示 12更早的评论隐藏 12更早的评论

Star Strider 2014-6-5

My pleasure!

Probably the easiest way is to save the estimated parameters (a save command works well here), copy your ODE and solver system and statements to your main script workspace, then with the estimated parameters in the workspace, solve your differential equation at the values of your independent variable, xdata. Your y3 variable should now be in your workspace, so plot it against xdata. (I suggest you plot your ydata as points rather than a line.)
The easiest way is to print them from inside your objective function. So the first line (just above your y0 assignment) might simply to put one line as c, because that will print the values of c to the Command Window. To be a bit more formal, use fprintf(1,'\\tc(1) = %.6f\ c(2) = %.6f\ c(3) = %.6f\',c) I don’t suggest plotting during the fitting process because it slows the process significantly. If you want to write the c values serially to a mat-file, you can load the file and write code to see the fitting process later without slowing the actual fitting process itself. There may also be a display option for lsqcurvefit, so you can use the optimset function and set 'Display', 'iter-detailed' to see everything. This wouldn’t save anything, just prints to the Command Window.

StarSign1997 2019-4-11

I copied and pasted the code here but an error occurs saying function value and YDATA are not equal in size. I am confused.

Johan Sebastian Diaz Tovar 2019-11-21

StarSign1997 when copy the code make sure to take the transpose of Sdata, if you don't do that, you'll get that error. So: Sdata = Sdata'; Doing this you solve that.

regards.

请先登录，再进行评论。

类别

AI and Statistics Curve Fitting Toolbox Interpolation

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

Community Treasure Hunt

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

Start Hunting!

Translated by