[讨论] 求解特征值问题

楼主: poringi (Leoric/小奈奈)   2014-06-10 16:43:07
小弟想用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
谢谢:)
作者: ejialan (eji)   2014-06-11 10:45:00
令x=w^2 det(c)会是x的n次多项式 4X4是4次多项式有公式解5次以上多项式没有公式 要数值解可以用roots

Links booklink

Contact Us: admin [ a t ] ucptt.com