是否可以使用在仿真期间由solve_ivp在每个求解器时间步长返回的求解矢量?

问题描述 投票:0回答:1

我一直在努力寻找一种方法来访问在每个求解器成功时间步长返回的求解矢量,以进行运行中的其他计算。将我的问题简化为:对于时间跨度[0,.001]我求解一个微分方程Xdot =A*X+B*U(状态空间方程)使用

sol=solve_ivp(conv,tspan,X0)

在每个sol.t点,计算X的解。我想在运行期间访问该解决方案。sol.y在模拟结束时返回整个数组,但是我希望在每个时间点都可以使用X的解决方案,以便在运行期间在函数conv中使用它。

python scipy
1个回答
-1
投票

如何包装函数conv

def fun(t,y):
    print(t)
    return conv(t,y)

每次调用该函数时,都会打印t的值。

© www.soinside.com 2019 - 2024. All rights reserved.