具有外部输入的动力系统的Python数值模拟

Python的Dopri5集成商需要dy / dt = f(t,y)。我想用零阶保持输入(逐段常数)模拟连续系统。实际上,根据某些外部输入,我们可以更改f,因为f是时间的函数。但是,在这种情况下出现的问题是,积分器使用了一些自适应步长,当您要求它查找输出直到某个时间t时,它也可能会使用f(t,y)的某些将来值,这可能是错误的,因为在那之后我们的输入以及f(t,y)都会改变。

有什么建议可以克服这个问题或其他选择吗?

谢谢

评论