help efun matlab,MATLAB：使用ode45时可能有两个事件值吗？

我不确定你是否可以按照自己的意愿去做，但是可以做很多事情。 首先，这听起来像首次通过时间 (又名第一次)的某种数值计算。 如果这些“粒子”是随机的， ode45停下来，不要使用ode45 ，而要使用ode45方法。

据我所知，你可以有多less个事件函数没有限制 – 或者说事件函数的维数(类似于ODE函数的维数) – 而且它们的数量并不与你拥有多less个ODE方程有关。 事件函数接收当前时间和当前状态向量。 您可以创build每个事件的任何或所有元素。 你是正确的，更多的事件function和更复杂的事件将减缓整合。 性能还取决于检测事件的频率。 如果你的每一个粒子到达“屋顶”，就像你所说的那样，只触发一个事件，那就不会太糟糕。

在实现方面，这里只是一个简单的例子，基于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。但是，它很简单，输出数组是矩形的而不是不齐的。

我并不清楚你的意思是什么“联结在一起”，需要粒子“停止”。 如果您不仅需要logging事件数据，例如更改系统参数或以其他方式更改微分方程，则需要在每个事件之后终止并重新开始集成。 看看ballode例子(在Matlab命令窗口中键入edit ballode )来看看一些build议，以提高效率。