setting ode45 to solve at even timesteps

Question

Windell 2013-4-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/73932-setting-ode45-to-solve-at-even-timesteps

Hi I'm attempting to force ode45 to solve a set of odes at at least some constant timestep. I understand that ode45 is a variable step solver so it is okay if it solves at intermediate timesteps. However, I would like to have it solve at at least some regular timestep also. For example, of the timestep is 1 second then I would like the solution to at least contain solutions at seconds 0,1,2,3,4,5 ...

Do any of you know how to make this happen?

thanks, windell

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

Jan 2013-4-29

This is really surprising for tspan=0:10000 (btw, neither :1: not the addition square brackets are required): Actually the output times are expecte to be [0, 1, 2, ...]. So where does 0.9960 come from?

Please do not describe what you have intented by the code, but the code itself. Otherwise we cannot find bugs.

Walter Roberson 2013-4-29

Wouldn't it make sense to interp1() to exact timepoints if your data is sufficiently close ?

请先登录，再进行评论。

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

Answer 1

Jan 2013-4-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/73932-setting-ode45-to-solve-at-even-timesteps#answer_83808

It is explained in the documentation how to force the intergrator to output the trajectory at certain time points: Use a vector with more than 2 elements for tspan.

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

请先登录，再进行评论。

Answer 2

Bjorn Gustavsson 2013-4-29

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/73932-setting-ode45-to-solve-at-even-timesteps#answer_83807

If you give it an array with the times for which you want the solution instead of just the start and stop time in the TSPAN input variable. Read first paragraph of ode45 help.

HTH,