首先,我正在使用该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还是陌生的,因此我想尽早增强良好的惯用代码。感谢您的任何帮助,请随时提出要求。