24小时热门版块排行榜    

查看: 1111  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见