我一直在努力寻找一种方法来访问在每个求解器成功时间步长返回的求解矢量,以进行运行中的其他计算。将我的问题简化为:对于时间跨度[0,.001]
我求解一个微分方程Xdot =A*X+B*U
(状态空间方程)使用
sol=solve_ivp(conv,tspan,X0)
在每个sol.t
点,计算X
的解。我想在运行期间访问该解决方案。sol.y
在模拟结束时返回整个数组,但是我希望在每个时间点都可以使用X
的解决方案,以便在运行期间在函数conv
中使用它。
如何包装函数conv
?
def fun(t,y):
print(t)
return conv(t,y)
每次调用该函数时,都会打印t
的值。