setting ode45 to solve at even timesteps

3 次查看(过去 30 天)
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 个评论
Jan
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
Walter Roberson 2013-4-29
Wouldn't it make sense to interp1() to exact timepoints if your data is sufficiently close ?

请先登录,再进行评论。

采纳的回答

Jan
Jan 2013-4-29
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.

更多回答(1 个)

Bjorn Gustavsson
Bjorn Gustavsson 2013-4-29
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,

类别

Help CenterFile Exchange 中查找有关 Ordinary Differential Equations 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by