What do the components of eigenvector represent?

11 次查看（过去 30 天）
Jay
Jay 2014-5-6
评论： Jay 2014-5-7
When I use [V,D] = eig(A), V has 2 X 2 matrix. What do the components of V represent? V(1,1) represents the angle between x-axis and principal direction, and V(2,1) represents the angle between y-axis and principal direction?? I would appreciate any help. Thanks,
Jay

请先登录，再回答此问题。

回答（1 个）

Youssef Khmou
Youssef Khmou 2014-5-6
编辑：Youssef Khmou 2014-5-6
N dimensional matrix is associated with N dimensional canonical base, in this case N=2, you have a plane (x,y), after eigendecomposition you have the diagonal matrix D which contains the spectra of the matrix A and the columns of V are the associated eigenvectors V(:,1)= V1 ex +V2 ey such as V(1,1) and V(2,1) are the x and y components of the first eigenvector .
V(1,1)= ||V1|| cos(theta)
V(2,1)= ||V1|| sin(theta)
  8 个评论
显示 6更早的评论
Jay
Jay 2014-5-7
B =
0.2696 0.4800
0.4800 -1.0203
>> [V,D]=eig(B)
V =
-0.3145 -0.9493
0.9493 -0.3145
D =
-1.1793 0
0 0.4287
Jay
Jay 2014-5-7
I am trying to draw principal directions on a plot.
% these codes are in for loop
[V,D] = eig(e2_j_k);
if (D(1,1) > D(2,2))
eigen_j_k = D(1,1);
p= [V(1,1);V(2,1)];
if (V(2,1) < 0)
theta_j_k = -acosd(dot([1 0],p));
else
theta_j_k = acosd(dot([1 0],p));
end
else
eigen_j_k = D(2,2);
p= [V(1,2);V(2,2)];
if (V(2,2) < 0)
theta_j_k = -acosd(dot([1 0],p));
else
theta_j_k = acosd(dot([1 0],p));
end
eval(sprintf('eigen_%d_%d = eigen_j_k',j,k));
eval(sprintf('theta_%d_%d = theta_j_k',j,k));
% I already know xc and zc
xc2(j,k) = xc(j,k) + 0.3*cosd(theta(j,k));
zc2(j,k) = zc(j,k) + 0.3*sind(theta(j,k));
figure(6);
plot([xc(j,k),xc2(j,k)],[zc(j,k),zc2(j,k)]);
hold on;
The code above is that I am trying to draw the directions of principal strain. I don't get any reasonable results. Do you find anything wrong or would you have other ways to draw? Any ideas? I would appreciate any help.
Jay

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Loops and Conditional Statements 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by