Find Eigenvalues without function "eig"
48 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
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.
1 comentario
Respuestas (2)
KSSV
el 11 de Oct. de 2020
Editada: KSSV
el 12 de Oct. de 2020
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)
0 comentarios
Paul
el 11 de Oct. de 2020
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.
0 comentarios
Ver también
Categorías
Más información sobre Linear Algebra en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!