要求解ODE系统,我需要过去时间步长的内部计算,因此我必须保存这些内部计算。由于求解器的动态时间步长,我无法保存所有计算。 事件函数也无济于事,因为事件仅接收(t,y)输入,但是我的内部数据是(t,y),f(t,y)的函数。
当求解器计算了一个解决方案而不是一个内部计算时,是否存在一个优雅的解决方案来保存内部参数?
问题是普通ODE的系统。
dydt = f(y,P),其中P = f(y,e)和e = f(P(t),P(t-1))
我需要的是(我认为)一种保存P的简单方法,以便访问上一个时间步长P(t-1)的值,而P并不是对dydt积分的解决方案。