Need help for cross checking of eigen Vector

조회 수: 7 (최근 30일)
moonman
moonman 2011년 11월 11일
I have a matrix
A=[2 1 1;2 3 2;1 1 2]
I have calculated its eigen values and associated eigen vectors by hand. I have cross checked the eigen values by using command eig(A). My eigen values are correct.
But how can i cross check that my manual calculation of Eigen Vectors is correct.
I have got eigen vector [1;2;1] for eigen value=5 I have got eigen vector [-2;1;1] for eigen value=1
Plz tell me the method to cross check eigen vectors

채택된 답변

Jonathan
Jonathan 2011년 11월 11일
This gives a possible set of unit eigenvectors.
[V,D] = eig(A)
If there is an eigenspace of more than one dimension, the vectors in V are not unique. Note also, that since the vectors in V are unit vectors, you need to normalize your vector in order to compare.
From the command above, we can see that the first and third columns of V correspond with the eigenvalue 1, and the second column of V corresponds with the eigenvalue 5. Normalizing [1; 2; 1] and taking the difference with V(:,2) yields the following.
V(:,2) - [1; 2; 1] ./ norm([1; 2; 1])
ans =
1.66533453693773e-016
-2.22044604925031e-016
0
Since the difference is on the order of machine epsilon (this can be checked with eps(V(1,2)) and eps(V(2,2))), we can say that to within machine precision Matlab's answer agrees with your answer for this eigenvalue.
It is possible that your eigenvector for the eigenvalue 1 is in the same eigenspace as given by columns one and three of V. To test this do the following.
coeffs = V(:,[1 3]) \ [-2; 1; 1];
coeffs(1)*V(:,1) + coeffs(2)*V(:,3)
Since the result is the vector [-2; 1; 1], we know that your eigenvector is in the eigenspace determined by the span of V(:,1) and V(:,3). Thus, your eigenvector for eigenvalue 1 is correct.
All of this has been fun (for me at least), but there is a much more direct method of testing an eigenpair for correctness.
>> A * [-2; 1; 1]
ans =
-2
1
1
>> A * [1; 2; 1]
ans =
5
10
5
This shows, directly, that your answers are correct.
  댓글 수: 1
moonman
moonman 2011년 11월 11일
Is the answer for eigen vector unique??
are my ansers correct

댓글을 달려면 로그인하십시오.

추가 답변 (1개)

Honglei Chen
Honglei Chen 2011년 11월 11일
The eigenvectors are not unique, they can differ by a scale. Based on MATLAB's answer, your eigenvector for the eigenvalue of 5 is correct, but your eigenvector for the eigenvalue of 1 is not. Note that you have two eigenvalues that are 1, so you need to perform extra steps when calculating the eigenvector.

카테고리

Help CenterFile Exchange에서 Data Type Identification에 대해 자세히 알아보기

태그

Community Treasure Hunt

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

Start Hunting!

Translated by