矩阵微分方程。求解斧头'=Bx+b

Question

我有矩阵微分方程Ax'=Bx+b，其中A和B是N*N的矩阵，b是向量。

我想用python解决它。希望有人能帮助我。

干杯！

Answer 1

如果你的矩阵A是正则的，则传递给odeint的函数是

def odefunc(x,t):
return numpy.linalg.solve(A, B.dot(x)+c)

当然，您也可以计算A的倒数，并用它左乘方程。

B = numpy.linalg.solve(A, B)
c = numpy.linalg.solve(A, c)
odefunc = lambda x,t: B.dot(x)+c