如何用Fortran求解矩阵指数

wx_TrpjQ8uw · 发表于 2021-9-6 10:34:11

本帖最后由 wx_TrpjQ8uw 于 2021-9-6 10:39 编辑

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

wx_TrpjQ8uw · 发表于 2021-9-6 10:37:05

C:\Users\Administrator\Desktop

li913 · 发表于 2021-9-6 10:50:37

本帖最后由 li913 于 2021-9-6 12:10 编辑

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

Transpose · 发表于 2021-9-6 11:36:21

先求出矩阵的特征值和特征向量，把特征值取exp之后，再利用矩阵乘法变换回去即可

wx_TrpjQ8uw · 发表于 2021-9-6 15:23:30

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

感谢感谢

liudy02 · 发表于 2021-9-7 20:33:05

若A不是一个厄米矩阵，这个问题数值上就是没有什么好办法
很可能计算本身就是不稳定的
如果是厄米矩阵，那么就用正交矩阵对角化就是了