Finding the eigenvectors of a 3x3 matrix - help please

  • #1

Homework Statement


Determine the eigenvalues and eigenvectors of the matric, A:

[tex]A=\left[\begin{array}{ccc}1 & 1 & 0\\ 1 & -2 & 0\\ 0 & 0 & 1\end{array}[/tex]

Homework Equations


I think i understand what is going on. I have found the matrix equation to be:

[tex]\left(1-\lambda\right)\left[\left(-2-\lambda)(1-\lambda)-1\right]=0[/tex]

So:

[tex]\lambda_{1}=1[/tex]
[tex]\lambda_{2}=\frac{-1-\sqrt{13}}{2}[/tex]
[tex]\lambda_{3}=\frac{-1-\sqrt{13}}{2}[/tex]

The Attempt at a Solution


I have gotten some solutions (but this being my first attempt at Tex i think would take me an age to write out!) but im confused as to which is the best method.

In lectures, our lecturer seemed to say that when you get to

[tex]V_{1}=\left[\begin{array}{ccc}x_{1}\\x_{2}\\x_{3}\end{array}[/tex]

You just set

[tex]x_{3}=1[/tex]

To make life easier and go from there. But i dont see how or why?

I always tried to use the normalisation condition that the square of the three components of the vectors equals one and find another relation from there (or sub in something else)

Could some one please help?

EDIT: Sorry for my failed attempt at Tex - only took half an hour!
 

Answers and Replies

  • #2
radou
Homework Helper
3,120
7
Hint: plug the eigenvalues into the equation (A - lambda*I)x = 0 and see what happens.
 
  • #3
vela
Staff Emeritus
Science Advisor
Homework Helper
Education Advisor
15,230
1,831

Homework Statement


Determine the eigenvalues and eigenvectors of the matric, A:

[tex]A=\left[\begin{array}{ccc}1 & 1 & 0\\ 1 & -2 & 0\\ 0 & 0 & 1\end{array}[/tex]

Homework Equations


I think i understand what is going on. I have found the matrix equation to be:

[tex]\left(1-\lambda\right)\left[\left(-2-\lambda)(1-\lambda)-1\right]=0[/tex]

So:

[tex]\lambda_{1}=1[/tex]
[tex]\lambda_{2}=\frac{-1-\sqrt{13}}{2}[/tex]
[tex]\lambda_{3}=\frac{-1-\sqrt{13}}{2}[/tex]

The Attempt at a Solution


I have gotten some solutions (but this being my first attempt at Tex i think would take me an age to write out!) but im confused as to which is the best method.

In lectures, our lecturer seemed to say that when you get to

[tex]V_{1}=\left[\begin{array}{ccc}x_{1}\\x_{2}\\x_{3}\end{array}[/tex]

You just set

[tex]x_{3}=1[/tex]

To make life easier and go from there. But i dont see how or why?

I always tried to use the normalisation condition that the square of the three components of the vectors equals one and find another relation from there (or sub in something else)

Could some one please help?

EDIT: Sorry for my failed attempt at Tex - only took half an hour!
If a vector v is an eigenvector of A, then any non-zero multiple of v is also an eigenvector of A for the same eigenvalue. In particular, unless x3=0, you can find some multiple of v where x3=1. So your lecturer is saying just set x3=1, just to keep things simple, and then find the corresponding x1 and x2. You can then choose to normalize the resulting eigenvector if you want.

(If you're having trouble actually solving for an eigenvector, do what radou suggested.)
 
  • #4
If a vector v is an eigenvector of A, then any non-zero multiple of v is also an eigenvector of A for the same eigenvalue. In particular, unless x3=0, you can find some multiple of v where x3=1. So your lecturer is saying just set x3=1, just to keep things simple, and then find the corresponding x1 and x2. You can then choose to normalize the resulting eigenvector if you want.

(If you're having trouble actually solving for an eigenvector, do what radou suggested.)
THANK YOU!

That was exactly what i have been looking for and have spent the better half of the past 4 hours trying to find. In text books i see "and now we set this to one" and no explanation as to why. It all makes sense now!

Im not having any troubles generally with finding eigenvectors - but for future reference then, which "method" is better / reccomended? The (in my eyes) "standard" method that radou posted or by exploiting the fact that you said? Surely setting x3 as one makes the whole process quicker and easier?

Thanks again
 
  • #5
vela
Staff Emeritus
Science Advisor
Homework Helper
Education Advisor
15,230
1,831
It's not an either-or choice. You do what radou said to get the equations the eigenvector must satisfy. Those equations will have more than one solution. As long as you find a non-zero solution to the equations, you have an eigenvector. Sometimes you can just see what a solution is by inspection, and you can just write down what the answer is. More generally, though, you can manipulate the equation so that some of the variables are expressed in term of the rest. For example, x1 and x2 might be written in terms of x3. Then to get one particular solution, you just set x3=1 to make the arithmetic easy and calculate the values for x1 and x2.
 
  • #6
Thank you for your continued help.

I know see that was a somewhat silly question.

I have been going through the original question trying the various methods avaliable.

Using http://www.arndt-bruenner.de/mathe/scripts/engl_eigenwert.htm" [Broken] I know that:

Code:
Characteristic polynomial:
    x^3 - 4x + 3

Real eigenvalues:
    {-2.302775637731995, 1, 1.3027756377319946}

Eigenvector of eigenvalue -2.302775637731995: 
    (-0.28978414868843005, 0.9570920264890529, 0)

Eigenvector of eigenvalue 1: 
    (0, 0, 1)

Eigenvector of eigenvalue 1.3027756377319946: 
    (0.9570920264890529, 0.28978414868843, 0)

First of all, i have found the eigenvalues (in agreement with the above). I know tried to calculate the eigenvectors.

I first used my original method:

[tex]
[\begin{array}{ccc}1 & 1 & 0\\ 1 & -2 & 0\\ 0 & 0 & 1\end{array}]

[\begin{array}{ccc}x_{1}\\x_{2}\\x_{3}\end{array}]

=\lambda_{1}

[\begin{array}{ccc}x_{1}\\x_{2}\\x_{3}\end{array}]
[/tex]

I then got some messy equations involving horrible fractions of [tex]x_{1}[/tex] and [tex]x_{2}[/tex]. I was able to use the normalsiation condition to separate them both and eventually got eignvectors that match the link i posted.

I then tried doing the same as above, but also setting [tex]x_{3}=1[/tex], which gave me the following (normalised) eigenvector:

[tex]
\frac{1}{\sqrt{11}}[\begin{array}{ccc}-3\\-1\\1\end{array}]

[/tex]

Which does not seem to be in agreement at all?

Finally i tried radou's form and once again set [tex]x_{3}=1[/tex]. This worked great getting me straight to [tex]v_{1}[/tex]. However, for [tex]v_{2}[/tex] [tex]x_{3}[/tex] is actually equal to 0 so the whole thing fell apart. Then what?

Where did i go wrong?

Thanks a lot
 
Last edited by a moderator:
  • #7
vela
Staff Emeritus
Science Advisor
Homework Helper
Education Advisor
15,230
1,831
For λ2, you get

[tex](A-\lambda_2 I)\vec{x} = \begin{bmatrix}1-\frac{-1-\sqrt{13}}{2} & 1 & 0 \\ 1 & -2-\frac{-1-\sqrt{13}}{2} & 0 \\ 0 & 0 & 1-\frac{-1-\sqrt{13}}{2}\end{bmatrix}\begin{bmatrix}x_1\\x_2\\x_3\end{bmatrix} = 0[/tex]

which gives you the equations

[tex]\frac{3+\sqrt{13}}{2}x_1 + x_2 = 0[/tex]

[tex]x_1 + \frac{-3+\sqrt{13}}{2}x_2 = 0[/tex]

[tex]{3+\sqrt{13}\over 2}x_3 = 0[/tex]

From the last equation, you get x3=0. The first two equations are multiples of each others, so you only need to use one of them. So take, say, the second one and solve for x1. If you now set x2=1, you find

[tex]x_1 = \frac{3-\sqrt{13}}{2}x_2 = \frac{3-\sqrt{13}}{2}[/tex]

So an eigenvector for λ2 is

[tex]V_2 = \begin{bmatrix}x_1\\x_2\\x_3\end{bmatrix} = \begin{bmatrix}\frac{3-\sqrt{13}}{2} \\ 1 \\ 0\end{bmatrix}[/tex]

Now, if you want to, you can normalize it.
 
  • #8
Thanks again vela for your continued support - i really appreciate it.

Ive had another go this morning on this problem and i see last night what i was doing was inserting x2=1 prematurely (in the original vector before i had actually expanded out the equations).

Everything works well with the method i was originally using:


[tex]

[\begin{bmatrix}1 & 1 & 0\\ 1 & -2 & 0\\ 0 & 0 & 1\end{bmatrix}]

[\begin{bmatrix}x_{1}\\x_{2}\\x_{3}\end{bmatrix}]

=\lambda_{1}

[\begin{bmatrix}x_{1}\\x_{2}\\x_{3}\end{bmatrix}]

[/tex]

I then tried again with your rearraged version:

[tex]
(A-\lambda_2 I)\vec{x} = \begin{bmatrix}1-\frac{-1-\sqrt{13}}{2} & 1 & 0 \\ 1 & -2-\frac{-1-\sqrt{13}}{2} & 0 \\ 0 & 0 & 1-\frac{-1-\sqrt{13}}{2}\end{bmatrix}\begin{bmatrix}x_1\\x_2\\x_3\end{bmatrix} = 0
[/tex]

Which works great for lamba2 and 3 - but for lamba1 (which equals 1) the whole right column of the matrix becomes zero. So i just get 0 = 0 and not the (0 0 1) vector that i got with the other method? Any ideas how to get around this?

Thanks again, i just want to get both methods straight in my head so i fully understand whats going on.

EDIT: Thinking about in, in my other method i just get x3 = x3, and i just said x3 = 1 - but it could be anything - so i guess this goes to what you originally said about any non zero multiple of the vector?)
 
Last edited:
  • #9
vela
Staff Emeritus
Science Advisor
Homework Helper
Education Advisor
15,230
1,831
Which works great for lamba2 and 3 - but for lamba1 (which equals 1) the whole right column of the matrix becomes zero. So i just get 0 = 0 and not the (0 0 1) vector that i got with the other method? Any ideas how to get around this?

Thanks again, i just want to get both methods straight in my head so i fully understand whats going on.

EDIT: Thinking about in, in my other method i just get x3 = x3, and i just said x3 = 1 - but it could be anything - so i guess this goes to what you originally said about any non zero multiple of the vector?)
Right. You're x3=x3 equation tells you nothing about x3, so it can be anything. In the other method, x3 just drops completely out of the equations, so again, there are no constraints on x3, so it can be anything.
 
  • #10
Perfect. So i could just have said x3 = 5, but then when i normalised it id get "one over the square root of 5 squared = 1/5)" and so just back to the (0 0 1) vector. That makes sense!

Thanks alot. This was driving me mad but ive got it all straight in my head now.

Thanks again.
 
  • #11
One last question. The next part of the problem reads:

"Find a matrix R such that RAR[tex]^{T}[/tex] = D, where D is a diagonal matrix (i.e.,
a matrix whose non-zero entries are all on the main diagonal)."

The answers say its just a matrix of the eigenvectors ive just found - why?

Thanks
 
  • #12
lanedance
Homework Helper
3,304
2
so for a handwaving explanation

composing a the matrix R of the eigenvctors of A and taking D = RART in effect represents a change of basis, from whatever you started with to the basis of eignevectors.

When an eignvector is multiplied by the matrix, the result will be the corresponding eigenvalue tiems the eigenvector.

In the new basisi teh first eigenvector is represented by (1,0,..,0)T, hence why the diagonalised matrix must have the eignevalues on the diagonal

otherwise work the math to understand why
 

Related Threads on Finding the eigenvectors of a 3x3 matrix - help please

  • Last Post
Replies
3
Views
1K
  • Last Post
Replies
1
Views
3K
  • Last Post
Replies
4
Views
692
Replies
8
Views
2K
  • Last Post
Replies
3
Views
9K
  • Last Post
Replies
3
Views
2K
  • Last Post
Replies
1
Views
6K
Replies
2
Views
5K
Replies
26
Views
3K
Top