24小时热门版块排行榜    

查看: 5055  |  回复: 25
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

学虫dj

金虫 (正式写手)

[求助] 为什么我用matlab中的eig函数求实对称矩阵的特征值有负数???? 已有3人参与

我在网上也查了一下,都说实对称矩阵的特征值应该是非负的。
但是我用matlab中的eig函数求【  矩阵G的转置乘以G 】GT*G   的特征值中竟然有负数。
这是为什么呢?
回复此楼
穷则独善其事,达则兼济天下
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
学虫dj: 金币+10 2015-11-09 08:45:35
引用回帖:
4楼: Originally posted by 学虫dj at 2015-11-08 10:41:39
嗯 你说的很有道理   我百度了一下,也查到了一些结果像你一样的结论  。但是也就证明应该是非负的下面我给你两个截图 你看看这结论是不是矛盾的

火狐截图_2015-11-08T02-30-14.432Z.png

火狐截图_2015-11-08T02 ...

你给的证据是对的,对于你给的这种特殊的对称阵(形如G’G),应该是半正定的,即特征值非负。
正定矩阵的定义是符合x'Ax >0的对称阵(也有更广义的定义,定义只符合x'Ax>0的矩阵为正定阵,即去掉了对称性要求)。也就是说正定是比对称更强的条件,对称阵不一定正定,正定阵一定是对称的。对称是正定的必要不充分条件。
正定阵有一些特性:特征值全部大于0,主元全部大于0,顺序主子式全部大于0,还有其他一些信息可以参考维基百科:https://en.wikipedia.org/wiki/Positive-definite_matrix

所以如果出现负特征值,很可能是软件算错了。
知其然,知其所以然。
22楼2015-11-08 12:38:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 26 个回答

mygt_hit

专家顾问 (职业作家)

【答案】应助回帖

感谢参与,应助指数 +1
实对称阵的特征值全为实数、不同特征值对于的特征向量正交,不一定是正的。举个简单例子,负单位阵,特征值为-1.
知其然,知其所以然。
2楼2015-11-06 13:16:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

virtualzx

木虫 (著名写手)

G是不是实矩阵?另外负特征值有多大?

只有实矩阵才能保证G^T.G半正定

发自小木虫IOS客户端
3楼2015-11-06 17:23:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

学虫dj

金虫 (正式写手)

引用回帖:
2楼: Originally posted by mygt_hit at 2015-11-06 13:16:10
实对称阵的特征值全为实数、不同特征值对于的特征向量正交,不一定是正的。举个简单例子,负单位阵,特征值为-1.

嗯 你说的很有道理   我百度了一下,也查到了一些结果像你一样的结论  。但是也就证明应该是非负的下面我给你两个截图 你看看这结论是不是矛盾的
为什么我用matlab中的eig函数求实对称矩阵的特征值有负数????
火狐截图_2015-11-08T02-30-14.432Z.png


为什么我用matlab中的eig函数求实对称矩阵的特征值有负数????-1
火狐截图_2015-11-08T02-38-20.034Z.png

穷则独善其事,达则兼济天下
4楼2015-11-08 10:41:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 确定了,国自然21号放榜 +4 布布和一二 2026-08-10 4/200 2026-08-10 17:48 by ZJTJZ
[基金申请] filecode +6 documentary 2026-08-10 6/300 2026-08-10 17:09 by wdc28
[基金申请] 我的国基提前知道中了,可是同事的操作让我实在接受不了,怎么会有这样的人 +8 家与远方 2026-08-10 11/550 2026-08-10 17:09 by 六两废铜
[基金申请] 国自然结果 +4 Vierhys 2026-08-10 8/400 2026-08-10 15:06 by Vierhys
[基金申请] FileCode能看出啥? +4 要乐观耀哥 2026-08-10 8/400 2026-08-10 15:02 by 天神眷顾
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +14 Lanmanbaby 2026-08-09 22/1100 2026-08-10 14:57 by wwncly
[基金申请] 8月时间戳变的,举个手。玩一下,释放压力 +11 archvillain 2026-08-04 13/650 2026-08-10 14:41 by seekfind
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 综述论文作为代表作会不会影响评审专家的印象分? +10 yufeiwaner 2026-08-09 10/500 2026-08-10 13:44 by cdutwml
[基金申请] 2026国自然放榜时间 +9 布布和一二 2026-08-08 9/450 2026-08-10 11:22 by xxxx2020
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 化学口download_prp&fileCode的固定段好像这几天一直没变,有变的大神么? +3 Tide man 2026-08-07 4/200 2026-08-07 22:39 by Tide man
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 有没有H口的?有收到消息的吗? +3 超级海虾 2026-08-04 3/150 2026-08-04 17:26 by 学教育滴
信息提示
请填处理意见