Eigenvalues and Eigenvectors of Symbolic Matrix
6 views (last 30 days)
Show older comments
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 Comments
Accepted Answer
Stephan
on 4 Dec 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 Comments
Tutku Ilgin Özcan
on 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
on 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)
More Answers (0)
See Also
Categories
Find more on Linear Algebra in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!