help efun matlab,matlab - MATLAB：使用ode45时是否可以有两个事件值？ - 堆栈内存溢出...

我不确定您是否可以完全按照自己的意愿进行操作，但是事件可以做很多事情。 首先，这听起来像是第一次通过时间 (也就是第一次击球时间)的某种数值计算。 如果这些“粒子”是随机的，请停止并且不要使用ode45 ，而是使用一种适合ode45方法。

据我所知，您可以拥有多少个事件函数没有限制-或者更确切地说，事件函数的维数(类似于ODE函数的维数)-并且它们的数量与您拥有的ODE方程数无关。 事件函数同时接收当前时间和当前状态向量。 您可以使用其中的任何或所有元素来创建每个事件。 您是正确的，更多的事件功能和更复杂的事件将减慢集成速度。 性能还取决于检测到事件的频率。 如果您的每个粒子都达到了您所说的“屋顶”，并且仅触发了一个事件，那就不会太糟。

在实现方面，这里是一个基于Matlab的ballode示例的ballode示例，该示例仅在垂直方向上模拟N个弹道粒子。 有N个非终止事件来捕获每个粒子通过y = 0的时间和速度。添加另一个终止事件来检查所有粒子是否都通过y = 0(如果我们知道这是哪个粒子就像我们在这里所做的那样，我们可以将该事件作为终止事件)。

function eventsdemo

% Initial conditions for n balls

n = 10;

y0(2*n,1) = 0;

y0(n+1:end) = linspace(20,40,n);

% Specify events function

options = odeset('Events',@(t,y)efun(t,y,n));

% Integrate

[t,y,te,ye,ie] = ode45(@(t,y)f(t,y,n),[0 10],y0,options);

figure;

plot(t,y(:,1:n),'b',te(1:n),ye(sub2ind(size(ye),ie(1:n),(1:n).')),'r.');

function dydt = f(t,y,n)

% Differential equations for ballistic motion

dydt = [y(n+1:end);zeros(n,1)-9.8];

function [value,isterminal,direction] = efun(t,y,n)

% Last event checks that all balls have hit ground and terminates integration

yn = y(1:n);

value = [yn;all(yn < 0)];

zn = zeros(n,1);

isterminal = [zn;1];

direction = [zn-1;1];

在某些方面，这有点低效，因为即使在某些N个系统通过y = 0之后，我们仍会对其进行仿真。但是，这很简单，输出数组是矩形而不是参差不齐。

我不清楚“耦合在一起”并需要粒子“停止”的确切含义。 如果您需要做的不仅仅是记录事件数据，例如以其他方式更改系统参数或更改微分方程，那么您需要在每个事件之后终止并重新启动积分。 查看ballode示例(在Matlab命令窗口中键入edit ballode )，以查看一些建议，以使其效率更高。