stop Ode45 when y is less than a value
5 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
Lazaros Christoforidis
am 28 Mär. 2020
Kommentiert: Lazaros Christoforidis
am 28 Mär. 2020
%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);
Akzeptierte Antwort
Torsten
am 28 Mär. 2020
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 Kommentare
Steven Lord
am 28 Mär. 2020
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.
Weitere Antworten (0)
Siehe auch
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!