ODE45 change y at event

Question

bstechel 2019-10-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/483590-ode45-change-y-at-event

编辑： Corentin Dewez-Raszka 2020-9-24

Can i influence the 'y' during an ode45 at a certain event?

I want to let the ODE45 stop at a certain event, manipulate the result at that timepoint and let it run further using the manipulated result as initial condition.A basic example is this:

spawnday=250;
function [value,isterminal,direction] = half_y_at_event(t, y, spawnday)
  value = [t-spawnday];  
  isterminal = [1]; 
  direction  = [];   
end
options=odeset('Events', @half_y_at_event, 'AbsTol', 1e-8, 'RelTol', 1e-8);
[t1, y1, te, ye, ie] = ode45(@dget_y, [0 500], [y_0], options, [pars]);
ye=ye/2;
[t2, y2] = ode45(@dget_y, [0 500], [ye], options, [pars]);

Is there a possibility to combine both ODE45 into one?

Mergin the matrices later is not an option.

2 个评论
显示无隐藏无

darova 2019-10-4

在 MATLAB Online 中打开

Using ode45 two times?

[t1, y1, te, ye, ie] = ode45(@dget_y, [0 500], [y_0], options, [pars]);
ye = ye/2;
[t2, y2] = ode45(@dget_y, [0 500]+t1(end), ye(:,end));

bstechel 2019-10-14

Using 2 ODE45 Is not convenient here since the event will occur multiple times over the total timespan of the ODE45 i'm interested in.

请先登录，再进行评论。

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

Answer 1

Steven Lord 2019-10-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/483590-ode45-change-y-at-event#answer_396291

在 MATLAB Online 中打开

You can use the ballode example as a model. That has an events function that stops the solver each time the bouncing ball touches the ground. It then changes the initial condition (based on the solution the solver gave when the event occurred) and starts the solver again with the new initial condition.

edit ballode