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!