如何解决计算上三角矩阵的矩阵指数的最快方法是什么?
我正在从事一个科学的 Python 项目,在该项目中我需要评估许多 exp(M)
形式的矩阵指数函数。目前,我正在使用 scipy.linalg.expm
函数来评估这个表达式,但即使对于我的小测试用例,这一步也需要很长时间。在查看了实现之后,该函数不会寻找或使用输入矩阵中可能存在的大部分结构。我提供给函数的矩阵实际上相当不错,因为矩阵 M 是块对角线,第一个矩阵 A 是上三角矩阵,第二个矩阵 B 是一个对角矩阵。
M = [[A,0],[0,B]]
我可以通过在第一个块之后获取对角线值来加快计算速度,但这仍然给我留下了 A 矩阵,它可以变得相当大。提前知道矩阵A是上三角时,有没有更快的计算exp(A)
的方法?
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。