ODE, how to retrive variables computed in the previous step

2 次查看（过去 30 天）
Adriano Filippo Inno
编辑： 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
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
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.

请先登录，再进行评论。

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

采纳的回答

Stephen23
Stephen23 2018-11-6
编辑：Stephen23 2018-11-6
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

更多回答（0 个）

类别

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

标签

产品


版本

R2017b

Community Treasure Hunt

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

Start Hunting!

Translated by