Finding the eigenvectors of a 3x3 matrix - help please

In summary: But i get an error saying "invalid column: 2". Do i need to normalize the vector before doing this?In summary, setting x3=1 in the equation A-lambda*Ix=0 yields a non-zero eigenvector of the eigenvalue 1.3027756377319946.
  • #1
Major_Disaster
14
0

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 I am 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 don't 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!
 
Physics news on Phys.org
  • #2
Hint: plug the eigenvalues into the equation (A - lambda*I)x = 0 and see what happens.
 
  • #3
Major_Disaster said:

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 I am 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 don't 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
vela said:
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 textbooks 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
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" 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
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 what's 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
Major_Disaster said:
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 what's 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 I've 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 I've just found - why?

Thanks
 
  • #12
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
 

FAQ: Finding the eigenvectors of a 3x3 matrix - help please

1. What is an eigenvector?

An eigenvector is a vector that remains in the same direction after being transformed by a matrix. In other words, it is a special vector that only gets scaled (stretched or compressed), but does not change direction when multiplied by a matrix.

2. Why is it important to find the eigenvectors of a 3x3 matrix?

Finding the eigenvectors of a 3x3 matrix allows us to understand the transformations that the matrix performs on a vector. This is useful in various fields such as physics, engineering, and computer science.

3. How do you find the eigenvectors of a 3x3 matrix?

To find the eigenvectors of a 3x3 matrix, we first need to find the eigenvalues of the matrix. We can do this by solving the characteristic equation. Then, we substitute each eigenvalue into the matrix and solve for the corresponding eigenvector by setting up a system of equations.

4. What is the relationship between eigenvectors and eigenvalues?

Eigenvectors and eigenvalues are closely related. Eigenvalues are the scalar values that get multiplied with the eigenvectors when a matrix is applied to them. Eigenvectors, on the other hand, are the vectors that remain in the same direction after being transformed by the matrix.

5. Can a 3x3 matrix have more than 3 eigenvectors?

No, a 3x3 matrix can have at most 3 eigenvectors. This is because the number of eigenvectors is equal to the number of eigenvalues, and a 3x3 matrix can have at most 3 distinct eigenvalues.

Similar threads

Replies
12
Views
2K
Replies
9
Views
446
Replies
7
Views
2K
Replies
2
Views
621
Replies
6
Views
697
Replies
11
Views
2K
Replies
19
Views
3K
Back
Top