Efficiently evaluating ODE event that depends on only some variables

Question

Roy Goodman 2023-12-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2055824-efficiently-evaluating-ode-event-that-depends-on-only-some-variables

评论： Roy Goodman 2023-12-4

I'm simulating an N-body problem in the plane.

I want to know when the three particles that are closest to each other are collinear. (Supposing I have a MATLAB function that identifies the "three closest particles.")

Suppose that I have a function f(z1,z2,z3) that vanishes when the particles z1, z2, z3 are collinear.

Right now my simulation evaluates a vector containing all

triples of particles and records all the times any triple is collinear.

Then at the end of the run, I sort through all the events and throw away all the ones where the particles are not "the three closest."

I would rather compute this during the ODE simulation and not afterward. I'd like to construct a new event g that gets triggered whenever f(z1,z2,z3)=0 and z1, z2, and z3. are the closest.

The events construct doesn't seem well-suited to this situation, since it detects zero crossings and the second half of my event is just a boolean. Is there a workaround?

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Torsten 2023-12-4

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2055824-efficiently-evaluating-ode-event-that-depends-on-only-some-variables#answer_1365204

编辑：Torsten 2023-12-4

Use OutputFcn to check your condition after each successful time step.

But I think it will be difficult to check for complanarity for a time-dependent system. It's like checking if a variable equals 0 - floating point arithmetic and the discrete time stepping will hinder to correctly detect such events.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Roy Goodman 2023-12-4

At first, I thought this answer made sense, but after thinking it through, I'm not sure.

The particles are evolving in the plane, so I'm not trying to test for complanarity, only for triplet wise collinearity. That's easy, because three particles with positiions

are collinear whenever

where

. That's the definition of my function f above.

So, I can define events easily. The problem is that I find too many of them because I only care about the collinearity of the threenearest particles.

outputFcn works after successful steps but I'm looking for something that works better with events.

Perhaps I'm missing something.

请先登录，再进行评论。