Eigenvalues and Eigenvectors of Symbolic Matrix
6 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
I have a symbolic matrix of which I want to get Eigenvalues and Eigenvectors. I want Eigenvalues and Eigenvectors in symbolic form.
syms E t
H = [E -t -t -t -t 0 0 0 0;-t E 0 0 0 -t -t 0 0;-t 0 E 0 0 0 0 -t -t;-t 0 0 E 0 -t 0 -t 0; -t 0 0 0 E 0 -t 0 -t; 0 -t 0 -t 0 E 0 0 0; 0 -t 0 0 -t 0 E 0 0; 0 0 -t -t 0 0 0 E 0; 0 0 -t 0 -t 0 0 0 E];
eig(H);
0 commentaires
Réponse acceptée
Stephan
le 4 Déc 2019
syms E t
H = [E -t -t -t -t 0 0 0 0;-t E 0 0 0 -t -t 0 0;-t 0 E 0 0 0 0 -t -t;...
-t 0 0 E 0 -t 0 -t 0; -t 0 0 0 E 0 -t 0 -t; 0 -t 0 -t 0 E 0 0 0;...
0 -t 0 0 -t 0 E 0 0; 0 0 -t -t 0 0 0 E 0; 0 0 -t 0 -t 0 0 0 E];
[V,D] = eig(H)
gives:
V =
[ 0, -1, -1, 0, 0, 0, 0, 2, 2]
[ -1, 0, 0, (E + 2^(1/2)*t)/(2*t) - E/(2*t), (E + 2^(1/2)*t)/(2*t) - E/(2*t), (E - 2^(1/2)*t)/(2*t) - E/(2*t), (E - 2^(1/2)*t)/(2*t) - E/(2*t), E/(2*t) - (E - 2*2^(1/2)*t)/(2*t), E/(2*t) - (E + 2*2^(1/2)*t)/(2*t)]
[ -1, 0, 0, E/(2*t) - (E + 2^(1/2)*t)/(2*t), E/(2*t) - (E + 2^(1/2)*t)/(2*t), E/(2*t) - (E - 2^(1/2)*t)/(2*t), E/(2*t) - (E - 2^(1/2)*t)/(2*t), E/(2*t) - (E - 2*2^(1/2)*t)/(2*t), E/(2*t) - (E + 2*2^(1/2)*t)/(2*t)]
[ 1, 0, 0, E/(2*t) - (E + 2^(1/2)*t)/(2*t), (E + 2^(1/2)*t)/(2*t) - E/(2*t), E/(2*t) - (E - 2^(1/2)*t)/(2*t), (E - 2^(1/2)*t)/(2*t) - E/(2*t), E/(2*t) - (E - 2*2^(1/2)*t)/(2*t), E/(2*t) - (E + 2*2^(1/2)*t)/(2*t)]
[ 1, 0, 0, (E + 2^(1/2)*t)/(2*t) - E/(2*t), E/(2*t) - (E + 2^(1/2)*t)/(2*t), (E - 2^(1/2)*t)/(2*t) - E/(2*t), E/(2*t) - (E - 2^(1/2)*t)/(2*t), E/(2*t) - (E - 2*2^(1/2)*t)/(2*t), E/(2*t) - (E + 2*2^(1/2)*t)/(2*t)]
[ 0, 0, 1, 0, -1, 0, -1, 1, 1]
[ 0, 1, 0, -1, 0, -1, 0, 1, 1]
[ 0, 1, 0, 1, 0, 1, 0, 1, 1]
[ 0, 0, 1, 0, 1, 0, 1, 1, 1]
D =
[ E, 0, 0, 0, 0, 0, 0, 0, 0]
[ 0, E, 0, 0, 0, 0, 0, 0, 0]
[ 0, 0, E, 0, 0, 0, 0, 0, 0]
[ 0, 0, 0, E + 2^(1/2)*t, 0, 0, 0, 0, 0]
[ 0, 0, 0, 0, E + 2^(1/2)*t, 0, 0, 0, 0]
[ 0, 0, 0, 0, 0, E - 2^(1/2)*t, 0, 0, 0]
[ 0, 0, 0, 0, 0, 0, E - 2^(1/2)*t, 0, 0]
[ 0, 0, 0, 0, 0, 0, 0, E - 2*2^(1/2)*t, 0]
[ 0, 0, 0, 0, 0, 0, 0, 0, E + 2*2^(1/2)*t]
2 commentaires
Tutku Ilgin Özcan
le 20 Oct 2022
When I try this, it does not satisfy H*V = V*D. Do you have any idea why it is so?
Torsten
le 20 Oct 2022
syms E t
H = [E -t -t -t -t 0 0 0 0;-t E 0 0 0 -t -t 0 0;-t 0 E 0 0 0 0 -t -t;...
-t 0 0 E 0 -t 0 -t 0; -t 0 0 0 E 0 -t 0 -t; 0 -t 0 -t 0 E 0 0 0;...
0 -t 0 0 -t 0 E 0 0; 0 0 -t -t 0 0 0 E 0; 0 0 -t 0 -t 0 0 0 E];
[V,D] = eig(H)
simplify(H*V-V*D)
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Linear Algebra dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!