如何在python中找到矩阵指数的积分

Question

我有一个表单矩阵，比如e^(Ax)，其中A是一个方阵。如何将它从给定值a集成到另一个值b，以便输出是相应的数组？

Answer 1

如果A具有正确的属性，则可以通过计算其特征向量和特征值将其转换为对角线形式A0。在对角线形式中，解是sol = [exp(A0*b) - exp(A0*a)] * inv(A0)，其中A0是具有特征值的对角矩阵，inv(A0)仅包含其对角线中特征值的倒数。最后，通过将解与左边的特征值和右边的特征值的转置相乘来转换回解：transpose(eigvecs) * sol * eigvecs。