Eigenvalues and Eigenvectors of Symbolic Matrix
조회 수: 6 (최근 30일)
이전 댓글 표시
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
채택된 답변
Stephan
2019년 12월 4일
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
Tutku Ilgin Özcan
2022년 10월 20일
When I try this, it does not satisfy H*V = V*D. Do you have any idea why it is so?
Torsten
2022년 10월 20일
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)
추가 답변 (0개)
참고 항목
카테고리
Help Center 및 File Exchange에서 Linear Algebra에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!