Matrix-Eigen Notes
Matrix-Eigen Notes
Yufei Tao
Department of Computer Science and Engineering
Chinese University of Hong Kong
[email protected]
1 Definitions
Let A be an n × n matrix. If there exist a real value λ and a non-zero n × 1 vector x satisfying
Ax = λx (1)
Example 1. Consider
1 −1
A =
2 4
(A − λI)x = 0
where I is the n × n identity matrix. Introducing B = A − λI, we can re-write the above as
Bx = 0 (2)
Let us consider the above as an equation about x. We know that if det(B) 6= 0, the above equation
has a unique solution x = 0. However, this is not what we want. Remember that our goal is to
find an eigenvector x of A, which needs to be a non-zero vector. Therefore, we must choose λ
appropriately to make det(B) = 0. This provides us a way to find the eigenvalues of A.
B = A − λI
1−λ −1
=
2 4−λ
1
Hence:
det(B) = (1 − λ)(4 − λ) + 2
= λ2 − 5λ + 6.
To make det(B) = 0, we can set λ to λ1 = 3 and λ2 = 2. These are all the eigenvalues of A.
Proof. The characteristic polynomial of A is a polynomial of degree n. Hence, Equation (3) can
have at most n distinct roots of λ.
Example 3. Consider again the matrix A in Example 1. We know from Example 2 that it has
two eigenvalues: λ1 = 3 and λ2 = 2.
x1
Let us first look for the eigenvectors of A for λ1 = 3. Namely, we want to find x = to
x2
satisfy:
(A − λ1 I)x = 0 ⇒
1−3 −1 x1 0
= ⇒
2 4−3 x2 0
−2 −1 x1 0
=
0 0 x2 0
x1
Hence, any satisfying −2x1 − x2 = 0 is a solution to the above system. The set of such
x2
vectors can be represented in a parametric form: x1 = t and x2 = −2t for any t ∈ R. This set—
which we denote as EigenSpace(λ1 )—has dimension 1. Every non-zero vector in EigenSpace(λ1 ) is
an eigenvector corresponding to λ1 .
2
x1
Similarly, to obtain the eigenvectors of A for λ2 = 2, we want x = to satisfy:
x2
(A − λ2 I)x = 0 ⇒
1−2 −1 x1 0
= ⇒
2 4−2 x2 0
−1 −1 x1 0
=
0 0 x2 0
x1
Hence, any satisfying −x1 − x2 = 0 is a solution to the above system. The set of such
x2
vectors can be represented in a parametric form: x1 = t and x2 = −t for any t ∈ R. This set,
denoted as EigenSpace(λ2 ), also has dimension 1. Every non-zero vector in EigenSpace(λ2 ) is an
eigenvector corresponding to λ2 .
Example 4. Consider
4 6 0
A = −3 −5 0
−3 −6 1
Its characteristic equation is:
det(A − λI) = 0 ⇒
4−λ 6 0
−3 −5 − λ 0 = 0 ⇒
−3 −6 1−λ
4−λ 6
(1 − λ) = 0 ⇒
−3 −5 − λ
(1 − λ)((4 − λ)(−5 − λ) + 18) = 0 ⇒
(λ − 1)2 (λ + 2) = 0.
Hence, A has two eigenvalues: λ1 = 1 and λ2 = −2.
x1
To look for the eigenvectors of A for λ1 = 1, we seek x = x2 satisfying
x3
(A − λ1 I)x = 0 ⇒
3 6 0 x1
−3 −6 0 x2 = 0 ⇒
−3 −6 0 x3
3 6 0 x1
0 0 0 x2 = 0 ⇒
0 0 0 x3
x1
Hence, any x2 satisfying 3x1 + 6x2 = 0 is a solution to the above system. The set of
x3
such vectors can be represented in a parametric form: x1 = 2u, x2 = −u, and x3 = v for any
3
(u, v) ∈ R2 . The set is denoted as EigenSpace(λ1 ), and has dimension 2. Every non-zero vector in
EigenSpace(λ1 ) is an eigenvector corresponding to λ1 .
x1
To look for the eigenvectors of A for λ2 = −2, we seek x = x2 satisfying
x3
(A − λ2 I)x = 0 ⇒
6 6 0 x1
−3 −3 0 x2 = 0 ⇒
−3 −6 3 x3
−1 −1 0 x1
0 1 −1 x2 = 0 ⇒
0 0 0 x3
x1
Hence, any x2 satisfying
x3
−x1 − x2 = 0
x2 − x3 = 0
is a solution to the above system. The set of such vectors can be represented in a parametric
form: x1 = −t, x2 = t, and x3 = t for any t ∈ R. This set is denoted as EigenSpace(λ2 ), and has
dimension 1. Every non-zero vector in EigenSpace(λ2 ) is an eigenvector corresponding to λ2 .
The vector space EigenSpace(λ) is referred to as the eigenspace of the eigenvalue λ. The
dimension of EigenSpace(λ) is referred to as the geometric multiplicity of λ.
Example 5.
4
We state the next result without proof:
In all the examples we have seen, the geometric multiplicity has always been the same as its
algebraic multiplicity. However, this is not always true, as you can see from the example below.
Example 6. Consider
−1 1 0
A = −4 3 0
1 0 2
−1 − λ 1 0
−4 3−λ 0 = 0 ⇒
1 0 2−λ
(λ − 1)2 (λ − 2) = 0.
Hence, A has only 2 eigenvalues: λ1 = 1 and λ2 = 2. They have algebraic multiplicities 2 and 1,
respectively.
Let us now focus on the eigenvalue λ1 = 1.We will see that its geometric multiplicity is 1. To
x1
find its corresponding eigenvectors x = x2 , we solve the equation:
x3
(A − λ1 I)x = 0 ⇒
−2 1 0 x1
−4 2 0 x2 = 0 ⇒
1 0 1 x3
−2 1 0 x1
0 1/2 1 x2 = 0 ⇒
0 0 0 x3
−2x1 + x2 = 0
(1/2)x2 + x3 = 0
−t
Therefore, the set of solutions to the above problem is the set EigenSpace(λ1 ) of vectors −2t
t
for all t ∈ R. This set has dimension 1.