Question about eigen-vector compute value

Hello all,
I am programing a code for deformation analysis to geodetic applications. I have a question about my results when I use [U, S, V]=svd() because I need to compute the eigevector belonging to the maximum eigen value. My code is:
clear
clc
o=eye(3);
z=zeros(3);
A= [o -o z z;
z -o o z;
z -o z o;
-o z o z;
o z z -o;
z z o -o];
m = ones(size(A,1),1)*1;
Q = diag(m.*m);
Qxx=pinv(A'*inv(Q)*A)';
Qdd=2*Qxx;
[U,S,V]=svd(Qdd);
ev=U1(1:12,1)
My results are:
0.0000
-0.1037
-0.0870
-0.0511
-0.5908
0.4316
0.1379
0.5637
-0.0650
-0.0868
0.1309
-0.2795
the first term is 6.47407614942793e-17 and this is my doubt, maybe it is a math problem with my matrix (Qdd). I don't know,
I would appreciate any help
Thank you very much,

回答 (1 件)

Paul
Paul 2022 年 10 月 7 日

0 投票

"I need to compute the eigevector belonging to the maximum eigen value."
Why is SVD being used? Just use eig?
o=eye(3);
z=zeros(3);
A= [o -o z z;
z -o o z;
z -o z o;
-o z o z;
o z z -o;
z z o -o];
m = ones(size(A,1),1)*1;
Q = diag(m.*m);
Qxx=pinv(A'*inv(Q)*A)';
Qdd=2*Qxx;
[V,D] = eig(Qdd,'vector');
D.'
ans = 1×12
0.5000 -0.0000 0.5000 0.5000 0.0000 0.0000 0.5000 0.5000 0.5000 0.5000 0.5000 0.5000
Lot's of repeated eigenvalues. Which one is considered the maximum eigenvalue?

カテゴリ

ヘルプ センター および File ExchangeLinear Algebra についてさらに検索

製品

リリース

R2022b

質問済み:

2022 年 10 月 7 日

回答済み:

2022 年 10 月 7 日

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by