Finding an eigenvector of 3x3 matrix

  • Thread starter enc08
  • Start date
42
0
Hi,

I'm trying to find an eigenvector of a matrix. I know that λ = 1, so my matrix (A - λI) is
[tex][-0.5253, 0.8593, -0.1906; -0.8612, -0.5018, 0.1010; 0.1817, 0.1161, -0.0236][/tex]

And from rows 2 and 3 I get these simultaneous equations

[tex]-0.8612t_{1}-0.5018t_{2}+0.1010t_{3}=0[/tex]
[tex]0.1817t_{1}+0.1161t_{2}+0.0236t_{3}=0[/tex]

I eliminate to find [itex]t_{2} = -4.02t_{3} [/itex] and [itex]t_{1}=-2.23t_{3}[/itex]

Thus the eigenvector is

t=[itex]k [-2.23, -4.02, 1][/itex]

But using an online solver gives the eigenvector as (-0.016, 0.206, 0.978).

Thanks for any pointers.
 
Last edited:

mathman

Science Advisor
7,703
392
Before you do anything else, check to see if it is an eigenvector (At = λt ?).
Also how do know λ = 1?
 
42
0
Hi,

The reason I'm assuming λ = 1 is because A is a rotation matrix (I didn't mention this in the original post), and I'm looking for the equivalent axes eigenvector.
 
42
0
As a sanity check, I have checked the eigenvalues and λ = 1 is one of them.
 
42
0
I realised I typo in my original matrix (matrix term 3,3 in the first post should be negative). I've correct it.

I get a closer answer, but my first term appears to be incorrect:

My answer: (-0.0137, 0.225, 1)
Given answer: (-0.0088, 0.216, 1)
 
32,837
4,554
Since your matrix is a rotation matrix, the entries are sin and cos values. I would try using more precision than the 4 decimal places you show, and see if that makes a difference in your resulting eigenvector.
 
42
0
I've just been given the numbers as is, and told that it's a rotation matrix.

Can you tell if my method for computing the eigenvectors is correct?

My approach is

λ = 1, so (A - λI) is
[tex][-0.5253, 0.8593, -0.1906; -0.8612, -0.5018, 0.1010; 0.1817, 0.1161, -0.0236][/tex]

From rows 2 and 3:

[tex]-0.8612t_{1}-0.5018t_{2}+0.1010t_{3}=0[/tex]

[tex]0.1817t_{1}+0.1161t_{2}-0.0236t_{3}=0[/tex]

Eliminate to find [tex]t_{2}=0.225t_{3} [/tex] and [tex]t_{1}=-0.0137t_{3}[/tex]

Thus the eigenvector is

t=[tex]k (-0.0137, 0.225, 1)[/tex]

But the actual answer is given as (-0.0088, 0.216, 1).
 
32,837
4,554
I've just been given the numbers as is, and told that it's a rotation matrix.

Can you tell if my method for computing the eigenvectors is correct?

My approach is

λ = 1, so (A - λI) is
[-0.5253, 0.8593, -0.1906; -0.8612, -0.5018, 0.1010; 0.1817, 0.1161, -0.0236]
Here's your matrix in a nicer form. To see what I did, right-click the matrix to view the underlying script.
$$A - I =\begin{bmatrix} -0.5253 & 0.8593 & -0.1906 \\ -0.8612 & -0.5018 & 0.1010 \\ 0.1817 & 0.1161 & -0.0236\end{bmatrix}$$

To check your work, I multiplied A - I above times your eigenvector x. That multiplication should produce a zero vector, but what I got was off by a little. My result was approximately < .0099, .0088, .00003>.

When you were row-reducing A - I, if you weren't careful with your arithmetic, you could have introduced a certain amount of imprecision in your results.


From rows 2 and 3:

[tex]-0.8612t_{1}-0.5018t_{2}+0.1010t_{3}=0[/tex]

[tex]0.1817t_{1}+0.1161t_{2}-0.0236t_{3}=0[/tex]

Eliminate to find [tex]t_{2}=0.225t_{3} [/tex] and [tex]t_{1}=-0.0137t_{3}[/tex]

Thus the eigenvector is

t=[tex]k (-0.0137, 0.225, 1)[/tex]

But the actual answer is given as (-0.0088, 0.216, 1).
This is easy to check. Just multiply (A - I) times <-0.0088, 0.216, 1>.
 
42
0
Thanks Mark44. I have used more precise values, and using your technique for testing the eigenvector, I get an answer much closer to a zero vector.

That you've confirmed my method is correct will do, especially since the discrepancies seem only due to rounding errors. I would have much preferred the question to use round numbers, but never mind :)
 

Related Threads for: Finding an eigenvector of 3x3 matrix

  • Last Post
Replies
8
Views
844
  • Last Post
Replies
1
Views
12K
Replies
9
Views
3K
  • Last Post
Replies
1
Views
528
Replies
6
Views
5K
  • Last Post
Replies
7
Views
10K
  • Last Post
Replies
2
Views
3K
  • Last Post
Replies
8
Views
1K

Hot Threads

Top