24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2807  |  回复: 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的回帖

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的回帖

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的回帖

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的回帖
相关版块跳转 我要订阅楼主 fo3mt2 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上本子正文33页,违规吗?会被低分嘛? +5 1234567wang 2026-05-17 6/300 2026-05-17 21:18 by maxjing
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 xx7gd5zq4e 2026-05-15 6/300 2026-05-17 19:36 by Equinoxhua
[基金申请] 这年头没有找到涵评专家,还有中面上的可能吗 +13 dd921ww 2026-05-12 16/800 2026-05-17 12:38 by 于轩
[考博] 26/27博士推荐 +3 1木头人13949 2026-05-13 3/150 2026-05-17 09:41 by YuY66
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 v9tggjlwd0 2026-05-15 4/200 2026-05-17 08:11 by 11n4dfd8yn
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 xx7gd5zq4e 2026-05-15 7/350 2026-05-17 07:58 by 11n4dfd8yn
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 l7k6xnh0yc 2026-05-14 6/300 2026-05-17 07:31 by 11n4dfd8yn
[硕博家园] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 cjf4bx70cj 2026-05-14 7/350 2026-05-17 06:55 by 11n4dfd8yn
[找工作] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:37 by ue3ir18jc3
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 k37jurhrau 2026-05-16 3/150 2026-05-17 01:25 by ue3ir18jc3
[基金申请] 精华III评审感受-评审感受-评审感受 +16 ferrarichen 2026-05-11 20/1000 2026-05-17 01:10 by 南开小綦
[考博] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 ky2p12rrjj 2026-05-15 4/200 2026-05-17 00:57 by ue3ir18jc3
[公派出国] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +3 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:35 by ue3ir18jc3
[论文投稿] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +5 x0mp7owy2b 2026-05-15 5/250 2026-05-17 00:29 by xiangfeng
[考研] 售SCI一区T0P文章,我:8.O.5.5.1.O.5.4,科目齐全,可+急 +4 x0mp7owy2b 2026-05-15 4/200 2026-05-17 00:25 by ue3ir18jc3
[有机交流] 求助2,4-二氯-5-嘧啶甲醛的合成方法 20+3 光吃不拉 2026-05-14 6/300 2026-05-16 19:46 by Equinoxhua
[有机交流] 求有机合成大神指点三硫酸乙烯酯(CAS:2793408-99-6)的合成路线 30+3 Leekmid 2026-05-13 10/500 2026-05-16 16:37 by czyzsu
[硕博家园] 申请博士 +3 呃?呃 2026-05-15 3/150 2026-05-16 11:01 by a4742549
[基金申请] 青C资助名额大幅增加! +11 西葫芦炒鸡蛋 2026-05-13 15/750 2026-05-15 14:36 by ambravo
[论文投稿] 求助大佬sci投稿哪个好中 +3 江沅188 2026-05-12 4/200 2026-05-13 14:35 by 江沅188
信息提示
请填处理意见