在Matlab中用事件函数求解ODEs

Question

我现在正在做一个项目，我已经建立了一个多体机械系统的运动模型。它是一个复杂的非光滑系统，系统的运动并不总是由同一组方程组描述的，这意味着如果我用数值方法求解它，我必须有时停止积分（然后用新的初始条件和方程组重新开始）。

function dy = firstODE(t,y)
dy = zeros(2,1);
dy(1) = y(2);
dy(2) = 1/(A1)*(B1-c*y(1));
end


function dy = secondODE(t,y)
dy = zeros(4,1);
dy(1) = y(2);
dy(2) = -c/A2*y(1);
dy(3) = y(4);
dy(4) = D2/B2 + c*C2/(A2*B2)*y(1);
end

function [value,isterminal,direction] = myEventsFcn(t,y)
value = [y(1)-K1, y(1)+K1, y(1)+K2, y(1)+K1, y(1)-K1];
isterminal = [1, 1, 1, 1, 1];
direction = [-1, -1, -1, 1, 1];
end

我如何告诉ODE-solver在事件发生后要使用什么新的初始条件和方程式？

Answer 1

调用两种不同的解决方案。即。

sol1 = ode45(@(t,x)firstODE(t,x),tspan1,init1);
% set some conditions
sol2 = ode45(@(t,x)secondODE(t,x),tspan2,init2);

如果希望tspan2从前一个中断的位置开始，并到达最初定义的tspan1的末尾，请将tspan2设置为[sol1.xe,tspan1(end)]。