Eigen Problem
Eigen Problem
Michael Haag
Justin Romberg
This work is produced by The Connexions Project and licensed under the
Creative Commons Attribution License ∗
Abstract
This module denes eigenvalues and eigenvectors and explains a method of nding
them given a matrix. These ideas are presented, along with many examples, in hopes
of leading up to an understanding of the Fourier Series.
http://cnx.rice.edu/content/m10736/latest/
Connexions module: m10736 2
(a)
(b)
http://cnx.rice.edu/content/m10736/latest/
Connexions module: m10736 3
1.2 Examples
Exercise 1:
From inspection and understanding of eigenvectors, nd the two eigenvectors, v1
and v2 , of
3 0
A=
0 −1
Also, what are the corresponding eigenvalues, λ1 and λ2 ? Do not worry if you are
having problems seeing these values from the information given so far, we will look
at more rigorous ways to nd these values soon.
Solution:
The eigenvectors you found should be:
1
v1 =
0
0
v2 =
1
And the corresponding eigenvalues are
λ1 = 3
λ2 = −1
Exercise 2:
Show that these two vectors,
1
v1 =
1
1
v2 =
−1
3 −1
are eigenvectors of A, where A = . Also, nd the corresponding
−1 3
eigenvalues.
Solution:
In order to prove that these two vectors are eigenvectors, we will show that these
statements meet the requirements stated in the denition (Denition 1).
3 −1 1 2
Av1 = =
−1 3 1 2
3 −1 1 4
Av2 = =
−1 3 −1 −4
These results show us that A only scales the two vectors ( i.e. changes their length)
and thus it proves that Equation 2 holds true for the following two eigenvalues that
you were asked to nd:
λ1 = 2
λ2 = 4
If you need more convincing, then one could also easily graph the vectors and their
corresponding product with A to see that the results are merely scaled versions of
our original vectors, v1 and v2 .
http://cnx.rice.edu/content/m10736/latest/
Connexions module: m10736 4
0 ... ... 1
So, A − λI is just a new matrix.
Example 1:
Given the following matrix, A, then we can nd our new matrix, A − λI .
a11 a12
A=
a21 a22
a11 − λ a12
A − λI =
a21 a22 − λ
If (A − λI) v = 0 for some v 6= 0, then A − λI is not invertible. This means:
det (A − λI) = 0
This determinant (shown directly above) turns out to be a polynomial expression (of order
n). Look at the examples below to see what this means.
Example 2:
Starting with matrix A (shown below), we will nd the polynomial expression,
where our eigenvalues will be the dependent variable.
3 −1
A=
−1 3
3−λ −1
A − λI =
−1 3−λ
2 2
det (A − λI) = (3 − λ) − (−1) = λ2 − 6λ + 8
λ = {2, 4}
http://cnx.rice.edu/content/m10736/latest/
Connexions module: m10736 5
Example 3:
Starting with matrix A (shown below), we will nd the polynomial expression,
where our eigenvalues will be the dependent variable.
a11 a12
A=
a21 a22
a11 − λ a12
A − λI =
a21 a22 − λ
det (A − λI) = λ2 − (a11 + a22 ) λ − a21 a12 + a11 a22
If you have not already noticed it, calculating the eigenvalues is equivalent to calculating
the roots of
det (A − λI) = cn λn + cn−1 λn−1 + · · · + c1 λ + c0 = 0
conclusion: Therefore, by simply using calculus to solve for the roots of our
polynomial we can easily nd the eigenvalues of our matrix.
vn λn vn
set of n equations with n unknowns. Simply solve the n equations to nd the eigenvectors.
3 Main Point
Say the eigenvectors of A, {v1 , v2 , . . . , vn }, span C n , meaning {v1 , v2 , . . . , vn } are linearly
2
x = α1 v1 + α2 v2 + · · · + αn vn (3)
where {α1 , α2 , . . . , αn } ∈ C . All that we are doing is rewriting x in terms of eigenvectors of
A. Then,
Ax = A (α1 v1 + α2 v2 + · · · + αn vn )
Ax = α1 Av1 + α2 Av2 + · · · + αn Avn
Ax = α1 λ1 v1 + α2 λ2 v2 + · · · + αn λn vn = b
Therefore we can write, X
x= (αi vi )
i
and this leads us to the following depicted system:
where in Figure 4 we have, X
b= (αi λi vi )
i
Main Point: By breaking up a vector, x, into a combination of eigenvectors,
the calculation of Ax is broken into "easy to swallow" pieces.
2 http://cnx.rice.edu/content/m10734/latest/#span_sec
3 http://cnx.rice.edu/content/m10734/latest/#lin_ind
http://cnx.rice.edu/content/m10736/latest/
Connexions module: m10736 6
Figure 4: Depiction of system where we break our vector, x, into a sum of its eigen-
vectors.
4 Practice Problem
Exercise 3:
For the following matrix, A and vector, x, solve for their product. Try solving it
using two dierent methods: directly and using eigenvectors.
3 −1
A=
−1 3
5
x=
3
Solution:
Direct Method (use basic matrix multiplication)
3 −1 5 12
Ax = =
−1 3 3 4
Eigenvectors (use the eigenvectors and eigenvalues we found earlier for this same
matrix)
1
v1 =
1
1
v2 =
−1
λ1 = 2
λ2 = 4
As shown in Equation 3, we want to represent x as a sum of its scaled eigenvectors.
For this case, we have:
x = 4v1 + v2
5 1 1
x= =4 +
3 1 −1
Ax = A (4v1 + v2 ) = λi (4v1 + v2 )
Therefore, we have
1 1 12
Ax = 4 × 2 +4 =
1 −1 4
Notice that this method using eigenvectors required no matrix multiplication. This
may have seemed more complicated here, but just imagine A being really big, or
even just a few dimensions larger!
http://cnx.rice.edu/content/m10736/latest/