数学建模社区-数学中国

标题: 求助高手:高次矩阵求解特征值和特征向量 [打印本页]

作者: oybs98    时间: 2007-5-22 19:58
标题: 求助高手:高次矩阵求解特征值和特征向量
<p>简单的、低次的教科书上有求解方法,也比较容易,但如果要计算高次,如10次或以上的就不知道怎么求解了。如果按教科书的方法求解起来非常困难,有谁知道怎么解吗?或有这个编程软件。如有请发到我邮箱:<a href="mailtoybs98@sina.com">oybs98@sina.com</a> 非常感谢!!!急!!!</p>
作者: hzhangjun1103    时间: 2007-5-23 19:37
如果你是在VC环境上实现,可以去下载我上传的《VC++和BC++数值分析类库》随书头文件和连接库试试,使用方法当然就是这本书了。如果需要电子版,请站内联系!
作者: oybs98    时间: 2007-5-24 17:26
<p>谢谢!我想在matlab里来计算。</p><p>本这个在matlab里利用矩阵符号来计算比较容易,但还是没弄出来。郁闷啊!!!</p>
作者: madio    时间: 2007-5-25 00:17
要是在matlab中直接使用eig函数就可以!你在matlab中查找Eigenvalues就可以找到相关的帮助!
作者: oybs98    时间: 2007-5-25 17:03
<p>谢谢管理员,我用的就是eigensys(A),但运行结果不对啊,成了复数。比如三阶的矩阵:</p><p>&gt;&gt;A=sym('[-1,1,1;1,1,2;1,1/2,1]');[V,E]=eigensys(A)</p><p>九阶的矩阵:&gt;&gt;A=sym('[1,2,1,1/2,1/3,3,1/3,2,2;1/2,1,2,3,3,5,3,3,3;1,1/2,1,3,1,2,1/2,2,2;2,1/3,1/3,1,1/2,3, 1/3,1,2;3,1/3,1,2,1,4,1/2,3,2;1/3,1/5,1/2,1/3,1/4,1,1/3,2,1;3,1/3,2,3,2,3,1,3,1;1/2,1/3,1/2, 1,1/3,1/2,1/3,1,2;1/2,1/3,1/2,1/2,1/2,1,1,1/2,1]')<br/>[V,E]=eigensys(A)</p><p>运行结果吓死人,太长了!!!!这个精度设置又不知道在哪?结果也对不上。郁闷</p>
作者: 水木年华zzu    时间: 2008-1-16 14:56
<p>a =</p><p>   -1.0000    1.0000    1.0000<br/>    1.0000    1.0000    2.0000<br/>    1.0000    0.5000    1.0000</p><p>&gt;&gt; eig(a)</p><p>ans =</p><p>    2.6075<br/>   -1.4778<br/>   -0.1298</p><p>&gt;&gt; </p>




欢迎光临 数学建模社区-数学中国 (http://www.madio.net/) Powered by Discuz! X2.5