标签: python matrix differential-equations exponential
我需要求解以下形式的微分方程:dy/dt = H*y,其中H是n x n矩阵,而y是n列向量。 / p>
dy/dt = H*y
H
n x n
y
n
我可以将等式改写为y = y_0 * Exp(H*t)。
y = y_0 * Exp(H*t)
选项是在第一个方程式或第二个形式上运行ODE求解器...天真地我想说矩阵指数在数值上会更密集?
就计算操作而言,这两种方法将如何扩展?