Fortran Coder

查看: 2947|回复: 2
打印 上一主题 下一主题

[通用算法] 矩阵的负二分之一次方计算

[复制链接]

58

帖子

21

主题

0

精华

专家

F 币
283 元
贡献
173 点
跳转到指定楼层
楼主
发表于 2022-4-14 01:21:16 | 只看该作者 回帖奖励 |倒序浏览 |阅读模式
大家好,

我想向大家请教一个问题。假设一个复数矩阵A有逆矩阵。

请问有什么办法可以计算这个矩阵A的负二分之一次方吗?就是A^{-\frac{1}{2}},比如LAPCAK或者BLAS库函数里的哪个子程序可以做这样的计算吗?

我知道ZGETRF和ZGETRI可以联合使用,计算矩阵A的逆矩阵,但要如何计算A^{-\frac{1}{2}}

谢谢。
分享到:  微信微信
收藏收藏 点赞点赞 点踩点踩

798

帖子

2

主题

0

精华

大宗师

F 币
3793 元
贡献
2268 点
沙发
发表于 2022-4-14 09:01:05 | 只看该作者
1、矩阵分解;2、对特征值倒数开方;3、矩阵恢复。
https://blog.csdn.net/DoctorCuiLab/article/details/83048600

58

帖子

21

主题

0

精华

专家

F 币
283 元
贡献
173 点
板凳
 楼主| 发表于 2022-4-15 00:49:41 | 只看该作者
非常感谢你的回复和建议。

我的理解是,先通过你提到的方法计算矩阵A的二分之一次方(A^{\frac{1}{2}}),然后再计算其逆矩阵,也就是A^{-\frac{1}{2}}
如果是这样的话,我可以再请教一个问题吗?

如果A是一个2000\times2000的复数矩阵,那在对它进行对角化时,也许会遇到物理内存溢出的问题。请问这种情况下,是否有其它办法计算A^{-\frac{1}{2}}呢?

谢谢,盼复。
您需要登录后才可以回帖 登录 | 极速注册

本版积分规则

捐赠本站|Archiver|关于我们 About Us|小黑屋|Fcode ( 京ICP备18005632-2号 )

GMT+8, 2024-4-28 01:38

Powered by Tencent X3.4

© 2013-2024 Tencent

快速回复 返回顶部 返回列表