ODE, how to retrive variables computed in the previous step

Question

Adriano Filippo Inno 2018-11-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/428282-ode-how-to-retrive-variables-computed-in-the-previous-step

编辑： Torsten 2018-11-6

I need to solve a nonlinear system inside the ODE at each time step in order to construct the current state derivatives. Since the system is quite involved I wish I can use, as initial condition to solve it, the solution of the system computed at the previous step. How can I do that? I also want to know if there is a faster procedure to get all the solutions of the system after all, instead of re-evaluating the ode function (not with global or persistent variables because I want just the accepted steps).

2 个评论
显示无隐藏无

Jan 2018-11-6

What exactly is "the previous step"? Which ODE integrator are you using? Even the one-step solvers evaluate the function to be integrated several times for each step. Afterwards the step size controller can reject a step. In consequence "previous step" is not clear.

"faster procedure to get all the solutions" - faster than what and what do you mean by "all the solutions"? The integrator replies one solution only.

Torsten 2018-11-6

编辑：Torsten 2018-11-6

ode15i is a solver for systems that are implicit in the state derivatives. So no nonlinear systems have to be solved by the user during the solution process.

Maybe it's worth changing the solver.

Best wishes

Torsten.

请先登录，再进行评论。

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

Answer 1

Stephen23 2018-11-6

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/428282-ode-how-to-retrive-variables-computed-in-the-previous-step#answer_345456

编辑：Stephen23 2018-11-6

在 MATLAB Online 中打开

The simplest solution is: use nested functions.

Write the input odefun (usually the first input) as a nested function. In its parent workspace you can initialize whatever variables you want, e.g.:

function dydt = main()
prev = [];
dydt = ode45(@myode,...)
%
    function dydt = myode(t,y)
    ... = prev % whatever using previous value
    ...
    prev = dydt;
    end
end