矩阵微分方程。求解Ax'= 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