24小时热门版块排行榜    

查看: 3220  |  回复: 4
【悬赏金币】回答本帖问题,作者fo3mt2将赠送您 10 个金币
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

fo3mt2

金虫 (小有名气)

[求助] 复数矩阵的-1/2次方计算 已有1人参与

大家好,

我想向大家请教一个问题。

A是一个n乘n的复数矩阵,如何计算A的负二分之一次方(-1/2)呢?

我看到网上有这样的计算流程。
1.对角化复数矩阵A,得到矩阵A的特征值
2.计算每个特征值的-1/2次方
3.用特征值的-1/2次方组成对角矩阵(非对角线上的元素全是零,对角线上的元素是每个特征值的-1/2次方),这个对角矩阵就是复数矩阵A的负二分之一次方(-1/2)

请问这种计算方法是否正确呢?

如果不对,麻烦大家给些建议,要如何正确计算一个复数矩阵A的-1/2次方吧。

谢谢啦。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fo3mt2

金虫 (小有名气)

引用回帖:
2楼: Originally posted by lq8989 at 2022-08-15 09:29:40
这个方法思路是对的,结果不完全对。我猜测楼主是想计算一个矩阵B满足 (B B) = A^-1。如果A可逆且可对角化,也就是 A = P^{-1} D P(特征分解),其中D为一对角矩阵,对角元为A的特征值且非零,则可以构造B = P^{-1 ...

谢谢你的回复。

我想计算的是一个矩阵B,满足B=A^(-1/2)。若A可逆且可对角化,并且若d的每个对角元也只有一个值的话,你在回复里提到的B=P^{-1}dP,是否就是我想要计算的这个矩阵B呢?

谢谢。
3楼2022-08-15 17:07:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

lq8989

至尊木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
这个方法思路是对的,结果不完全对。我猜测楼主是想计算一个矩阵B满足 (B B) = A^-1。如果A可逆且可对角化,也就是 A = P^{-1} D P(特征分解),其中D为一对角矩阵,对角元为A的特征值且非零,则可以构造B = P^{-1} d P,其中d的对角元是D的对角元的平方根的倒数。由于楼主的矩阵A是在复数域定义的,d的每个对角元都有一到两个可取值,因此可以构造出最多2^n个B。
2楼2022-08-15 09:29:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lq8989

至尊木虫 (著名写手)

【答案】应助回帖

具体要看你如何定义B=A^(-1/2)。就像正实数有正负两个实数平方根一样,满足B^2 = A^(-1)的B矩阵可以多达2^n个。如果不在意特征值的取值,那用这个方法构造出的任意一个都满足B^2 = A^(-1)。如果对特征值有要求,比如假设A是半正定的(也就是特征值全部非负),你也可以要求B半正定,也就是d的对角元全部取D的算术平方根(非负平方根);如果A是非正定的(也就是特征值有负数或复数),你也可以选择d的对角元全部取D的复平方根主值。
4楼2022-08-16 15:01:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

flankway

新虫 (初入文坛)

只能是半正定矩阵才可以开平方。你这个问题,我们分解一下,先开平方,然后求逆。如果能开平方,还能求逆,那么这个矩阵一定是正定矩阵。正定矩阵可以对角化,即可以存在酉矩阵,在有等价意义下变成一个对角矩阵,其中特征值都大于零,可以开平方。然后你把那个酉矩阵带上,就是原矩阵的平方根。平方根求完了,你再求它的逆矩阵。

总的来说,正定矩阵就像是正数,半正定矩阵就是大于等于零的数,他们和数的性质有点像。负数是不能开平方的,同样的,特征值有负的,也不能开平法。
加油!
5楼2022-09-08 23:42:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[职场人生] 40+中年人失业后,还能找到好工作吗? +3 xqy2009 2026-10-13 3/150 2026-10-14 20:17 by 43358052
[有机交流] OpenAI又攻克数学难题了,离AI做有机合成还有多远? +8 wangyl_123 2026-10-10 13/650 2026-10-14 15:12 by zuesli
[基金申请] 申请科研项目纯靠实力能成功吗? 200+4 David1982 2026-10-13 7/350 2026-10-14 14:24 by 13108017953
[硕博家园] 极限拉扯,如何抉择 +10 Isak0320 2026-10-12 14/700 2026-10-14 11:39 by thegarfield
[教师之家] 课题组会帮你通过考核么? +4 好运来~9527 2026-10-10 7/350 2026-10-13 08:41 by xli1984
[基金申请] 入职后第一次管课题组的专利缴费,对账才发现官费比别人多掏了一倍 +6 13108017953 2026-10-08 6/300 2026-10-12 20:07 by w3w3y3
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:34 by XQNyDCa1Tk5n
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:26 by XQNyDCa1Tk5n
[找工作] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 22:26 by XQNyDCa1Tk5n
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 9/450 2026-10-11 22:15 by XQNyDCa1Tk5n
[教师之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 22:06 by XQNyDCa1Tk5n
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 (EPI+-1)(金币-50) +5 qR2AdCWTVtkw 2026-10-09 10/500 2026-10-11 14:25 by VOMexwnDzOFZ
[教师之家] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 (EPI+-1)(金币-50) +4 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 14:05 by VOMexwnDzOFZ
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:51 by VOMexwnDzOFZ
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 qR2AdCWTVtkw 2026-10-09 8/400 2026-10-11 13:45 by VOMexwnDzOFZ
[找工作] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 qR2AdCWTVtkw 2026-10-09 7/350 2026-10-11 10:54 by VOMexwnDzOFZ
[考博] 招收博士生一名 +3 mandolin 2026-10-10 3/150 2026-10-11 10:21 by 不是山谷!
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 qR2AdCWTVtkw 2026-10-09 4/200 2026-10-11 03:14 by M0zrYVSAnxpE
[有机交流] 各位大佬,路线设计求助 +6 zapen 2026-10-08 8/400 2026-10-10 15:59 by zapen
[有机交流] 我应该选哪条合成路线让实验成功率高一些? +5 jing816y 2026-10-08 9/450 2026-10-10 09:05 by zapen
信息提示
请填处理意见