Fortran Coder

标题: 如何用Fortran求解矩阵指数 [打印本页]

作者: wx_TrpjQ8uw 时间: 2021-9-6 10:34
标题: 如何用Fortran求解矩阵指数
本帖最后由 wx_TrpjQ8uw 于 2021-9-6 10:39 编辑

矩阵A大约百十来阶，求exp（-tA），t为常数参数

作者: wx_TrpjQ8uw 时间: 2021-9-6 10:37
C:\Users\Administrator\Desktop

作者: li913 时间: 2021-9-6 10:50
本帖最后由 li913 于 2021-9-6 12:10 编辑

这是无穷累加，只需要满足精度的前n次就行。如果矩阵不大，可以直接求，代码如下
计算矩阵函数 exp(A)-CSDN社区

QQ截图20210906104917.png (101.66 KB, 下载次数: 156)

作者: Transpose 时间: 2021-9-6 11:36
先求出矩阵的特征值和特征向量，把特征值取exp之后，再利用矩阵乘法变换回去即可

作者: wx_TrpjQ8uw 时间: 2021-9-6 15:23

Transpose 发表于 2021-9-6 11:36
先求出矩阵的特征值和特征向量，把特征值取exp之后，再利用矩阵乘法变换回去即可 ...

感谢感谢

作者: liudy02 时间: 2021-9-7 20:33
若A不是一个厄米矩阵，这个问题数值上就是没有什么好办法
很可能计算本身就是不稳定的
如果是厄米矩阵，那么就用正交矩阵对角化就是了

欢迎光临 Fortran Coder (http://bbs.fcode.cn/)