[討論] 求解特徵值問題

作者: 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