→ ejialan:令x=w^2 det(c)會是x的n次多項式 4X4是4次多項式有公式解 06/11 10:45
→ ejialan:5次以上多項式沒有公式 要數值解可以用roots 06/11 10:46
推 fishenos:為何我覺得似曾相似，是buckling的問題嗎 06/25 13:06
→ fishenos:你要不要試試power method 06/25 13:06
小弟想用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://www.ptt.cc/bbs/MATLAB/M.1402389795.A.147.html