小编Rol*_*wer的帖子

在DifferentialEquations.jl中的n个回调之后,是否有惯用的方式终止集成

首先,我正在使用该DifferentialEquations.jl库,这太棒了!无论如何,我的问题如下:

例如,我有以下微分方程:

function f(du, u, t)
    du[1] = u[3]
    du[2] = u[4]
    du[3] = -u[1] - 2 * u[1] * u[2]
    du[4] = -u[2] - u[1]^2 + u[2]^2
end
Run Code Online (Sandbox Code Playgroud)

我有一个回调,每次轨迹越过y轴都会触发:

function condition(u, t, integrator)
    u[2]
end
Run Code Online (Sandbox Code Playgroud)

但是,我需要积分恰好在3个交叉点后终止。我知道可以通过使用以下效果终止集成:

function affect!(integrator)
    terminate!(integrator)
end
Run Code Online (Sandbox Code Playgroud)

但是在满足终止条件之前,允许对回调数目进行计数的正确方法是什么。此外,是否有一种方法可以将此方法扩展到具有n个不同计数的n个事件?

在研究中,我经常需要查看Poincare映射,然后将第一,第二,第三等返回到该映射,因此我需要一个框架来允许我执行此计数终止。我对Julia还是陌生的,因此我想尽早增强良好的惯用代码。感谢您的任何帮助,请随时提出要求。

differential-equations julia

5
推荐指数
1
解决办法
85
查看次数

标签 统计

differential-equations ×1

julia ×1