How to stop ode45 using events function

2017 3 10
1 个回答
回答已采纳
更新时间：2017 3 10
6 次查看（30 天）

0 个投票

I'm trying to stop an ode45 solver when the altitude gets below a certain value, but it isn't working.
The code is below - I want to stop it when r gets to 10 km. I don't want to use an rdot value (vardot(1)=rdot).
tRange=[0,350];
var=[r;lat;long;u;v;w;e0;e1;e2;e3;omegaX;omegaY;omegaZ]
opts=odeset('Events',@EndEvents,'InitialStep',1e-6,'RelTol',1e-5);
[t,vardot]=ode45(@Trajectory,tRange,var,opts)
function [value,isterminal,direction]=EndEvents(t,vardot)
value = (var(1)-10000);
isterminal=1;
direction=0;
end

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

 采纳的回答

Torsten
Torsten 2017-3-10

1 个投票

function [value,isterminal,direction]=EndEvents(t,var)
instead of
function [value,isterminal,direction]=EndEvents(t,vardot)
Best wishes
Torsten.

更多回答（0 个）

类别

在 帮助中心 和 File 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