Find Eigenvalues without function "eig"
显示 更早的评论
I performed this by hand and know that the eigenvalues are 10, -1, and 4. But clearly that's not what I got... Any help would be nice.
回答（2 个）
syms x
A = [10 0 0; 1 -3 -7 ; 0 2 6] ;
p = det(A-eye(3)*x) % this is the characteristic equation
val = solve(p) % solve this for roots
The polynomial equation, you get is:
poly = 
Using roots:
p = [-1 13 -26 -40] ;
val = roots(p)
Paul
2020-10-11
roots(poly(A))
Though I’m sure one could come with an A for which this approach will fail due to numerical accuracy issues. In which case converting A to a syms object as in another answer may be preferred.
类别
在 帮助中心 和 File Exchange 中查找有关 Linear Algebra 的更多信息
2020-10-11
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
