我正在使用scipy的odeint解决一个微分方程系统,我现在遇到一个问题,我似乎无法修复: 如果满足某个条件,我必须及时回到过去。到较早的状态,更改(全局)参数并继续从那里进行积分。
起初我想,虽然你无法及时回到odeint,我将存储满足条件的时间,但是要改变系统的状态(集成的变量)并继续从那里进行整合。 不幸的是,我无法从odeint中找到过去时间访问变量值的方法。有人知道解决我的困境的方法吗?
如果那是不可能的,有人可以建议采用不同的方法吗?性能是一个问题,所以我真的不想在python中编写自己的解算器。我可以在C中完成它,但这需要我改变整个实现,这将花费我很多时间。