stop Ode45 when y is less than a value

22 次查看（过去 30 天）
Lazaros Christoforidis
评论： Lazaros Christoforidis 2020-3-28
%So thats my ode call and I want it to stop (and save datas % obviously before t=13s) when the value u(1)<=minu
%Where minu is calcuated before.
u0=[vf1m(jj) 0]';
tspan=[0 13]';
Opt = odeset('Events', @myEvent);
[t,u]= ode45(@(t,u) ok(t,u,p2,m,maxu,maxT), tspan, u0,Opt);
%I tried something like that, nothing worked.
function [value, isterminal, direction] = myEvent(t, u, minu)
value = u(1)<=minu
isterminal = 1; % Stop the integration
direction = 0;
%Eroor Message
Not enough input arguments.
Error in myEvent (line 2)
value = u(1)<=minu
Error in odeevents (line 28)
eventValue = feval(eventFcn,t0,y0,eventArgs{:});
Error in ode45 (line 148)
odeevents(FcnHandlesUsed,odeFcn,t0,y0,options,varargin);
Error in skr (line 217)
[t,u]= ode45(@(t,u) ok(t,u,p2,m,maxu,maxT), tspan, u0,Opt);

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

采纳的回答

Torsten
Torsten 2020-3-28
As defined, ode45 expects myEvent to have 2 input parameters, but you use 3.
Use
Opt = odeset('Events',@(t,u)myEvent(t,u,minu));
And I suggest you set
value = u(1) - minu;
  5 个评论
显示 3更早的评论
Steven Lord
Steven Lord 2020-3-28
You want your event function to be continuous (as the one Torsten suggested is) rather than discontinuous (like yours that uses <= is.)
The ODE solver will detect when the event function passes through zero.
Lazaros Christoforidis
Lazaros Christoforidis 2020-3-28
Ohh okay, noww it makes sense
Ευχαριστω, υγεια

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by