作者poringi (Leoric)
看板MATLAB
标题[讨论] 求解特徵值问题
时间Tue Jun 10 16:43:07 2014
小弟想用MATLAB求解5X5矩阵的特徵值，写法如下:
a=[2,-1,0,0,0;
-1,2,-1,0,0;
0,-1,2,-1,0;
0,0,-1,2,-1;
0,0,0,-1,1];
b=[4,1,0,0,0;
1,4,1,0,0;
0,1,4,1,0;
0,0,1,4,1;
0,0,0,1,2];
a1=5*10^8*a;
syms w
b1=(15.72/30)*w^2*b;
c=a1-b1;
[w]=solve(det(c));w
结果跑出下列error:
Warning: Explicit solution could not be found.
> In solve at 169
在这之前有跑过4X4的矩阵可以用这个方法解出来，不知道为什麽到五阶矩阵就不行了，
恳请有经验的大大指点指点QQ
谢谢:)
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 220.129.23.120
※ 文章网址: http://webptt.com/cn.aspx?n=bbs/MATLAB/M.1402389795.A.147.html
1F：→ ejialan:令x=w^2 det(c)会是x的n次多项式 4X4是4次多项式有公式解 06/11 10:45
2F：→ ejialan:5次以上多项式没有公式 要数值解可以用roots 06/11 10:46
3F：推 fishenos:为何我觉得似曾相似，是buckling的问题吗 06/25 13:06
4F：→ fishenos:你要不要试试power method 06/25 13:06