24小时热门版块排行榜    

查看: 1285  |  回复: 5

ljb1672

铁虫 (初入文坛)

[交流] 计算矩阵特征值和特征向量 已有3人参与

关于计算矩阵特征值的论文较容易查询,也论述得较深入。想求大牛比较下,计算是一起算出特征向量的方法有多少?只知道Jacobi在矩阵迭代时可以几乎同时(在编程代码上最多增加10行代码)求出特征向量,想问下,其他方法是否也可以较方便的方法,或有论文比较过。
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

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

feixiaolin

荣誉版主 (文坛精英)

优秀版主


小木虫: 金币+0.5, 给个红包,谢谢回帖
matlab   eig()一下下

[ 发自小木虫客户端 ]
2楼2015-03-01 22:39:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

鸡蛋_汤圆

铁虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
虫友,你能把语言组织通顺吗?
1.维基百科中关于各种求解矩阵特征值和向量的方法有注解http://en.wikipedia.org/wiki/Eigenvalue_algorithm
2. Matlab中可以指定算法的  http://www.mathworks.com.au/help/matlab/ref/eig.html
(lambda = eig(A,B,algorithm) specifies the generalized eigenvalue algorithm as one of two strings: 'qz', which uses the QZ algorithm, or 'chol', which uses the Cholesky factorization of B)
具体的也可以参见www.mathworks.com/moler里面有一篇关于数值算法的一本书。从数值分析教材上我的理解,eig用的算法跟乘幂法有关,可能再利用了正交变换(householder变换)来处理矩阵。至于其时间复杂度,大概是N的三方。
3楼2015-03-02 11:54:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

鸡蛋_汤圆

铁虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
你也可以上图书馆随便找一本书,无论是中文的“数值分析”还是英文的“matrix analysis”或者“numerical analysis”,相信都会有“矩阵特征值和特征向量”的专门论述和讨论。
在冯果忱、黄明游主编的《数值分析》第四章中讲述矩阵特征值问题有讲到乘幂法、反幂法、对称矩阵的子空间迭代方法、QR方法(QR分解,利器)、对称矩阵的Jacobi方法等等。不同的方法适用不同的目标和处理对象,有不同的处理效率。
至于具体算法我一时记不上来,你可以针对你自己的具体问题参见我提到或者没提到的书籍。
至于matlab中eig函数的具体算法,你可以在matlab命令窗口中help eig 查看,或者edit 函数名 查看函数代码。
4楼2015-03-02 12:04:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ljb1672

铁虫 (初入文坛)

我使用C编写Jacobi方法求解特征值(因为可以并行,在增加很少代码的情况下方便求出特征向量,并且精确性可以相当的高)。现在想问下,其他方法如QR(Matlab使用),Lanczos,Arnoldi,还有其他是否有几乎同时求出特征向量的特点。
5楼2015-03-02 13:27:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dracula0724

铜虫 (初入文坛)


小木虫: 金币+0.5, 给个红包,谢谢回帖
QR与jacobi都不适用于大规模矩阵的特征值计算 lanczos用于对称矩阵 arnoldi用于非对称 两者都是迭代方法 生成一组基向量 将矩阵投影在这组基向量张成的子空间中 然后可以用QR或者Jacobi算法求解投影后的小规模特征值问题。
6楼2015-04-08 15:40:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 ljb1672 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] E0414, 我的本子有没有希望? +7 布布和一二 2026-06-17 9/450 2026-06-21 19:29 by 布布和一二
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 19:18 by q1p9ixelzp
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:58 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 18:00 by q1p9ixelzp
[考研] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 4/200 2026-06-21 17:45 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:40 by q1p9ixelzp
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 9skg9htng1 2026-06-20 3/150 2026-06-21 17:37 by q1p9ixelzp
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 4/200 2026-06-21 17:25 by ky8v3skvft
[考博] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:20 by ky8v3skvft
[博后之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 6/300 2026-06-21 17:05 by ky8v3skvft
[公派出国] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 5/250 2026-06-21 17:00 by ky8v3skvft
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +5 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:57 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +3 vlfdw50o0u 2026-06-20 5/250 2026-06-21 16:45 by ky8v3skvft
[教师之家] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 7/350 2026-06-21 16:40 by ky8v3skvft
[论文投稿] 售T0P一区SCI文章,我:8O5.51.O.54,科目齐全,可+急 +4 vlfdw50o0u 2026-06-20 9/450 2026-06-21 16:37 by ky8v3skvft
[基金申请] 青A35岁以下通知答辩了吗 +3 暨阳一只柴 2026-06-17 3/150 2026-06-21 09:39 by kudofaye
[考博] 【全奖博士/科研助理/博后招生】新加坡南洋理工大学机械与航空航天学院 +3 robohaha 2026-06-15 4/200 2026-06-21 00:36 by lanbing1
[基金申请] 希望面上有个好结果 +7 碧水00 2026-06-16 7/350 2026-06-18 12:18 by wuke100666
[论文投稿] 三区计算机方向期刊推荐 +5 1457340941 2026-06-15 5/250 2026-06-17 13:16 by 会议编辑
[论文投稿] sci论文二审求助 +5 潘倍倍 2026-06-15 5/250 2026-06-16 10:15 by xs74101122
信息提示
请填处理意见