我正在编写一个考虑速度的数值算法。我在 scipy/numpy (scipy.linalg.expm2, scipy.linalg.expm) 中遇到了两个矩阵指数函数。但是我有一个矩阵,我事先知道它是对角线的。这些 scipy 函数在运行之前是否检查矩阵是否为对角线?显然,对角矩阵的求幂算法可以快得多,我只是想确保它们正在做一些聪明的事情 - 如果不是,是否有一种简单的方法可以做到这一点?
最佳答案
如果一个矩阵是对角矩阵,那么它的指数可以通过对主对角线上的每个元素取幂来获得,所以你可以通过以下方式计算它:
np.diag(np.exp(np.diag(a)))
关于python - Python 中的对角矩阵指数,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/16081797/