Use ode values without storing

Question

olabaz 2017-9-6

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/355469-use-ode-values-without-storing

编辑： Stephen23 2017-9-7

Hi, I have a system of ODEs that I can solve using ode23. To solve I do:

[t, S] = ode23(@odefunc,tSpan,S0)

and then what I am really interested in is a variable:

D(t) = 2*pi*(S(1)(t) + S(2)(t) + S(3)(t) + ...)

and then I plot D(t) vs t.

Doing it this way, however, I am forced to store all S(1)(t), all S(2)(t), etc. when I really only need their sum for each time. This takes a lot of memory since I have 280k*3 equations. How can I go about just calculating D at each step and saving that in an array to then plot D(t) vs t.

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

请先登录，再进行评论。

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

Answer 1

Stephen23 2017-9-7

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/355469-use-ode-values-without-storing#answer_280575

编辑：Stephen23 2017-9-7

Write your own output function, which is called by ode45 after each successful step:

https://www.mathworks.com/help/matlab/ref/odeset.html#input_argument_d0e708834

Inside the function:

Store the sum value in a persistent variable.
Add a special case that returns the persistent variable (e.g. the third input = 'getval').

To use it:

Use the init call to empty the persistent variable.
On each call append the sum value onto the persistent variable.
After ode45 has finished call the output function with your special third input to return the persistent variable.

Easy! You have stored only the values that you want, and get them in a vector at the end.

2 个评论
显示无隐藏无

olabaz 2017-9-7

Thanks. This is what I was looking for. I wasn't able to get the output working so I just made a global variable and appended to it.

Stephen23 2017-9-7

编辑：Stephen23 2017-9-7

@olabaz: I am glad that it works. If you upload your code I will have a look. I would recommend avoiding global if you can.

请先登录，再进行评论。

Answer 2

Jacob Ward 2017-9-6

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/355469-use-ode-values-without-storing#answer_280531

在 MATLAB Online 中打开

You could create a function that calls ode23, finds the sum of S, and then returns that sum. Something like this:

function [t,D] = odeSum()
    ...
    [t, S] = ode23(@odefunc,tSpan,S0)
    D(t) = 2*pi*(S(1)(t) + S(2)(t) + S(3)(t) + ...)
end

Since the function throws away any variables that aren't output variables, you wouldn't be saving all the useless S data and would be left with only D when you come out of the function.

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

Jacob Ward 2017-9-6

I see. Well, I could be wrong, but I think in order to do what you are saying you'd have to go in and change the ode23 code itself so that it only outputs what you want. Maybe someone else will have a better suggestion.