24小时热门版块排行榜    

查看: 6614  |  回复: 7

yuxi1989

银虫 (小有名气)

[求助] 为什么对于对称方阵,其特征值和奇异值是一样的? 已有2人参与

RT。主要是看一个matlab代码,做PCA。它是把协方差矩阵作SVD分解,然后把U作为特征向量,S作为特征值了。正常的PCA应该是协方差矩阵求特征值和特征向量。网上查了一下,有说法“对于对称方阵,其特征值分解和奇异值分解是一样的?”。想了下,作者可能是根据这个来写的,因为向量的协方差矩阵是对称的,而且是方的。可能matlab的SVD比eig快,作者节约时间!但是我还是不理解“为什么对于对称方阵,其特征值和奇异值是一样的?”特征值的话,总感觉是平方的关系,特征向量组成的矩阵和U的关系?这个我就更不明朗了。所以前来求助。谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

iamsad

至尊木虫 (著名写手)

【答案】应助回帖


感谢参与,应助指数 +1
yuxi1989: 金币+1, ★★★很有帮助, 谢谢!但是怎么判定正定呢?貌似只有根据特征值了。成鸡生蛋、蛋生鸡的问题了。 2014-07-17 19:19:23
对于对称方阵,其特征值和奇异值未必一样! 因为奇异值一定非负,而特征值可以是负的。
只有symmetric positive semi-definite的矩阵,其特征值和奇异值才一样。
2楼2014-07-15 21:36:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

math2000

铁杆木虫 (职业作家)

【答案】应助回帖


感谢参与,应助指数 +1
yuxi1989: 金币+1, ★★★很有帮助, 谢谢!但貌似真的如上楼所说,需要正定的条件。 2014-07-17 19:20:28
设A是对称方阵(考虑实矩阵,复的一样),则A^T(A的转置)=A。
若r是A的特征值,则r的平方是A的平方的特征值,又因为A是对称矩阵,所以r的平方是A乘A的转置的特征值,由A的奇异值定义知 r是A的奇异值。
反之,若r是A的奇异值,即r^2是AA^T的特征值,由上知,r^2是A^2的特征值,即r是A的特征值,所以楼主的结论是对的
3楼2014-07-15 22:11:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hy东舟

银虫 (初入文坛)

对称矩阵特征值与奇异值差正负号,特征向量和奇异向量张成空间一样。从谱分解来看,对称矩阵可以酉对角化,A=Q*B*Q^t,B对角线特征值,只要将负号加到Q的列上,使B对角线大于0,得到的就是奇异值分解
4楼2014-07-16 10:31:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuxi1989

银虫 (小有名气)

引用回帖:
4楼: Originally posted by hy东舟 at 2014-07-16 10:31:24
对称矩阵特征值与奇异值差正负号,特征向量和奇异向量张成空间一样。从谱分解来看,对称矩阵可以酉对角化,A=Q*B*Q^t,B对角线特征值,只要将负号加到Q的列上,使B对角线大于0,得到的就是奇异值分解

非常感谢!(抱歉,不知道怎么你这个帖子不能给金币了)。但是我还是有一个问题,就是如我所说。通过SVD(奇异值分解),对实对称方阵A=U*S*V(S为奇异值)。其U是不是就是特征向量呢?因为PCA关键是要用协方差矩阵的特征向量构建变换矩阵。谢谢!
5楼2014-07-17 19:26:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hy东舟

银虫 (初入文坛)

引用回帖:
5楼: Originally posted by yuxi1989 at 2014-07-17 19:26:14
非常感谢!(抱歉,不知道怎么你这个帖子不能给金币了)。但是我还是有一个问题,就是如我所说。通过SVD(奇异值分解),对实对称方阵A=U*S*V(S为奇异值)。其U是不是就是特征向量呢?因为PCA关键是要用协方差矩阵 ...

就我所知,理论上U不一定是特征向量的,但matlab中svd算法当矩阵对称且中小规模时,本质上是就是利用QR算法算谱分解,所以计算的是特征向量,你可以实验下。另外svd一般不会比eigs快,毕竟svd要计算两组正交基底U,V,eig只需要一组特征向量基底
6楼2014-07-18 10:59:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuxi1989

银虫 (小有名气)

引用回帖:
6楼: Originally posted by hy东舟 at 2014-07-18 10:59:46
就我所知,理论上U不一定是特征向量的,但matlab中svd算法当矩阵对称且中小规模时,本质上是就是利用QR算法算谱分解,所以计算的是特征向量,你可以实验下。另外svd一般不会比eigs快,毕竟svd要计算两组正交基底U, ...

谢谢!你的线性代数基础真好,得多向您学习。我再去看看那个代码,想想那个作者为什么要这么写!
7楼2014-07-19 17:42:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

一壶漂泊2015

新虫 (初入文坛)

楼主,你好,我也遇到了和你相同的疑问,请问你解决了吗?
8楼2016-05-11 19:44:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yuxi1989 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +15 低垂的野花 2026-07-31 20/1000 2026-08-01 14:10 by gltch
[基金申请] 2026年国自然面上资助率 +15 布布和一二 2026-07-30 18/900 2026-08-01 13:56 by gltch
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 8/400 2026-08-01 13:53 by y4XuShGyllvT
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:51 by y4XuShGyllvT
[基金申请] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:51 by y4XuShGyllvT
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 jzg3xk1po8 2026-07-31 7/350 2026-08-01 13:49 by y4XuShGyllvT
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 jzg3xk1po8 2026-07-31 10/500 2026-08-01 13:49 by y4XuShGyllvT
[基金申请] 小木虫看见有人已经知道结果了 +16 1234567wang 2026-07-25 20/1000 2026-08-01 08:57 by type2type
[基金申请] 微信指数没变化,科研之友没阅读 +14 wangze12014 2026-07-28 18/900 2026-07-31 18:35 by lzchaihuo
[教师之家] 基础研究怎么拉横向,学校到款任务越来越多,难以完成 拉横向,都有哪些途径啊 +8 锦衣卫寒战 2026-07-28 8/400 2026-07-31 18:13 by beefly
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
[考博] 申请2027年材料/化学类博士 10+3 考研老狗? 2026-07-25 5/250 2026-07-30 19:33 by 白色木鱼
[基金申请] 今年的WR进展到哪一步了? +5 wsgjhwz 2026-07-25 9/450 2026-07-30 12:59 by 可淡不可忘
[有机交流] 产物和副产物价值比较 30+4 小汤02 2026-07-27 5/250 2026-07-30 12:44 by czyzsu
[基金申请] 时间戳他又来了 +14 晓晓爱翠翠 2026-07-26 17/850 2026-07-29 16:30 by 呵呵二哥
[有机交流] 这个自发加氧反应的机理是什么? 50+4 YUAN2273 2026-07-25 7/350 2026-07-29 12:27 by xiaobenpu
[基金申请] 这种情况还有戏吗 +5 drbart 2026-07-27 11/550 2026-07-29 05:50 by drbart
[基金申请] 同事接到电话了,我却没有 +4 1234567wang 2026-07-27 4/200 2026-07-28 10:11 by GOODLUCKER
[基金申请] 准备明年的基金了 +5 Tide man 2026-07-27 5/250 2026-07-28 09:25 by newfuzzy1
[分析] 水蒸气蒸馏法乳化现象 +3 月光娇娇 2026-07-26 3/150 2026-07-26 23:25 by caiyun
信息提示
请填处理意见