How to vectorize or speed up ODE solver with vector ?

4 次查看（过去 30 天）

Ole 2016-1-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/265890-how-to-vectorize-or-speed-up-ode-solver-with-vector

评论： Daniel Armyr 2016-2-3

I would like to vectorize or speed up ODE solver with vector E? This is extremely slow following this routine.

nTime     = length(timespan);
y_list    = zeros(nTime, length(y0));
y_list(1, :) = y0;
for iTime = 2:nTime
  [t_, y_] = ode45(@yourFcn, timespan(iTime-1:iTime), y0, options, E(iTime));
  y0       = y_(end, :);
  y_list(iTime, :) = y0;
end

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

Daniel Armyr 2016-2-3

If I understand it correctly, you are calling ode45 in a loop and letting ode45 just compute a single step for each loop iteration.

As far as I can tell, your code does nothing that it could not do by a single call to ode45 and then extracting the data you want from the return vectors of ode45.

请先登录，再进行评论。

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