Mat121 Module 5
Mat121 Module 5
Example 1. If A is the identity matrix 𝐼𝑛 , then the only eigenvalue is λ = 1. Every nonzero
0 1
Example 2. Let A = [ 1 2 ]. Then
2 0
1
1 0 1 1 1 1 1
A [ ] = [ 1 2 ] [ ] = [21] = 2 [ ] so that 𝑋1 = [ ] is an
1 2 0 1 1 1
2
1
eigenvector of A associated with eigenvalue λ1 = . Also,
2
−1
1 0 1 1 1 1 1
A [ ] = [ 1 2 ] [ ] = [ 21 ] = − 2 [ ] so that 𝑋2 = [ ] is an eigenvector of A
−1 2 0 −1 −1 −1
2
1
associated with the eigenvalue λ2 = − .
2
0 0 1 0 0 1 0 1 1
Example 3. Let A = [ ]. Then A [ ] = [ ] [ ] = [ ] = 0 [ ] so that 𝑋1 = [ ] is an
0 1 0 0 1 0 0 0 0
0 0 0 0 0 0 0
A[ ] = [ ] [ ] = [ ] = 1 [ ] so that 𝑋2 = [ ] is an eigenvector of A associated with the
1 0 1 1 1 1 1
1 1
Example 4. Let A = [ ]. Find the eigenvalues of A and their associated eigenvectors.
−2 4
Solution:
𝑥1 1 1 𝑥1 𝑥1
Let X = [𝑥 ]. Then [ ] [𝑥 ] = λ [𝑥 ] . Thus
2 −2 4 2 2
𝑥1 + 𝑥2 = λ𝑥1 or (λ − 1) 𝑥1 − 𝑥2 = 0
−2 𝑥1 + 4𝑥2 = λ𝑥2 2 𝑥1 + (λ − 4)𝑥2 = 0
Note that if A is an nxn matrix, then the homogeneous system AX = 0 has a nontrivial
solution if and only if |𝐴| = 0. We have
λ−1 −1
| | = 0. This means that (λ − 1)( λ − 4) + 2 = 0 or
2 λ−4
1 1 𝑥1 𝑥1
If λ = 2 . Then AX = 2X or [ ] [𝑥 ] = 2 [𝑥 ]. Thus
−2 4 2 2
𝑥1 + 𝑥2 = 2𝑥1 or 𝑥1 − 𝑥2 = 0
−2𝑥1 + 4𝑥2 = 2𝑥2 2𝑥1 − 2𝑥2 = 0.
𝑟
associated with the eigenvalue λ = 2 are [ ], where r is any nonzero real number.
𝑟
1 1 𝑥1 𝑥1
If λ = 3 . Then AX = 3X or [ ] [𝑥 ] = 3 [𝑥 ]. Thus
−2 4 2 2
𝑥1 + 𝑥2 = 3𝑥1 or 2𝑥1 − 𝑥2 = 0
−2𝑥1 + 4𝑥2 = 3𝑥2 2𝑥1 − 𝑥2 = 0.
1
Hence, 𝑥1 = 2 𝑥2 where 𝑥2 = any real number r. Therefore, all eigenvectors
𝑟
associated with the eigenvalue λ = 3 are [ 2 ], where r is any nonzero real number.
𝑟
𝑓(λ) = |λ𝐼𝑛 − 𝐴| = 0
Theorem. The eigenvalues of A are the real roots of the characteristic polynomial of A.
1 2 −1
Example. Let A = [1 0 1], find
4 −4 5
a. characteristic polynomial of A
b. eigenvalues of A
c. eigenvectors of A
Solution:
1 0 0 λ 0 0
a. Note that λ𝐼𝑛 = λ | 0 1 0 | = | 0 λ 0 |
0 0 1 0 0 λ
λ−1 −2 1
𝑓(λ) = |λ𝐼𝑛 − 𝐴| = | −1 λ −1 |
−4 4 λ−5
= λ3 − 6λ2 + 11λ − 6.
b. The possible integer roots of 𝑓(λ) = λ3 − 6λ2 + 11λ − 6 are ±1, ±2, ±3, and ±6.
1−1 −2 1 𝑥1 0
[ −1 1 − 1 ] [𝑥2 ] = [0]
−4 4 1 − 5 𝑥3 0
or
0 −2 1 𝑥1 0
[−1 1 ]
−1 2[ 𝑥 ] = [ 0].
−4 4 −4 3 𝑥 0
We have
−2𝑥2 + 𝑥3 = 0
−𝑥1 + 𝑥2 −𝑥3 = 0
−4𝑥1 + 4𝑥2 − 4𝑥3 = 0
1 −𝑟
0 −2 1 ⋮ 0 1 0 2 ⋮ 0 2
[−1 1 −1⋮ 0] → [0 1 − 1⋮ 0] → [ 𝑟 ]
−4 4 −4⋮ 0 0 0 2⋮ 0 2
0 𝑟
−𝑟
2
Hence, all eigenvectors associated with the eigenvalue λ1 = 1 are [ 𝑟 ], where r is any
2
𝑟
nonzero real number.
−1
If r = 2 then 𝑋1 = [ 1 ] is an eigenvector of A associated with λ1 = 1.
2
2−1 −2 1 𝑥1 0
[ −1 2 − 1 ] [𝑥2 ] = [0]
−4 4 2 − 5 𝑥3 0
or
1 −2 1 𝑥1 0
[ −1 2 − 1 ] [𝑥2 ] = [0].
−4 4 − 3 𝑥3 0
We have
𝑥1 − 2𝑥2 + 𝑥3 = 0
−𝑥1 + 2𝑥2 −𝑥3 = 0
−4𝑥1 + 4𝑥2 − 3𝑥3 = 0.
1 𝑟
0 −2 1 ⋮ 0 1 0 2 ⋮ 0 2
[−1 2 −1⋮ 0] → [0 1 − 1⋮ 0] → [ 𝑟 ]
−4 4 −3⋮ 0 0 0 4⋮ 0 4
0 𝑟
−𝑟
2
Hence, all eigenvectors associated with the eigenvalue λ2 = 2 are [ 𝑟 ], where r is any
4
𝑟
nonzero real number.
−2
If r = 4 then 𝑋2 = [ 1 ] is an eigenvector of A associated with λ2 = 2 .
4
3−1 −2 1 𝑥1 0
[ −1 3 𝑥
− 1 ] [ 2 ] = [0]
−4 4 3 − 5 𝑥3 0
or
2 −2 1 𝑥1 0
[ −1 3 𝑥
− 1 ] [ 2 ] = [0].
−4 4 − 2 𝑥3 0
We have
2𝑥1 − 2𝑥2 + 𝑥3 = 0
−𝑥1 + 3𝑥2 −𝑥3 = 0
−4𝑥1 + 4𝑥2 − 2𝑥3 = 0.
1 −𝑟
0 −2 1 ⋮ 0 1 0 2 ⋮ 0 2
[−1 1 −1⋮ 0] → [0 1 − 1⋮ 0] → [ 𝑟 ]
−4 4 −4⋮ 0 0 0 2⋮ 0 2
0 𝑟
−𝑟
4
Hence, all eigenvectors associated with the eigenvalue λ3 = 3 are [ 𝑟 ], where r is any
4
𝑟
nonzero real number.
−1
If r = 4 then 𝑋3 = [ 1 ] is an eigenvector of A associated with λ3 = 3 .
4
Note: The characteristic polynomial of a given matrix may have imaginary roots and it may
even have no real roots.
0 1
Example 2. Let A = [ ]. Then the characteristic polynomial of A is
−1 0
Exercises
A. Find the characteristic polynomial of the following
2 1
1. [ ]
−1 3
1 2 1
2. [ 0 1 2]
−1 3 2
4 −1 3
3. [0 2 1]
0 0 3
2 2 3
1 −1
1. [ ] 4. [1 2 1]
2 4
2 −2 1
1 0 0 2 0 0
2. [−1 3 0] 5. [3 −1 0]
3 2 −2 0 4 3
1 2 3 4
2 −2 3
0 −1 3 2
3. [0 3 −2] 6. [ ]
0 0 3 3
0 −1 2
0 0 0 2
1 1 1 1 2 −1
Example. Let A = [ ] and P = [ ]. Then 𝑃−1 = [ ] and
−2 4 1 2 −1 1
2 −1 1 1 1 1 2 0
B = 𝑃 −1 𝐴𝑃 = [ ][ ] [ ]= [ ]. Thus B is similar to A.
−1 1 −2 4 1 2 0 3
Properties
1. A is similar to A.
2. If B is similar to A, then A is similar to B.
3. If A is similar to B and B is similar to C, then A is similar to C.
1 1
Example 1. Let A = [ ] . From page 2 example 4, the eigenvalues are λ = 2, 3 and
𝑟
−2 4
𝑟
eigenvectors: [ ] , [ 2 ]. Consider λ1 = 2 λ2 = 3. If we let r = 1 and r = 2, then we have the
𝑟 𝑟
following eigenvectors
1 1
𝑋1 = [ ] and 𝑋2 = [ ]
1 2
1 1 2 −1
P=[ ] and 𝑃−1 = [ ]. Thus
1 2 −1 1
2 −1 1 1 1 1 2 0
𝑃−1 𝐴𝑃 = [ ][ ][ ]=[ ] = 𝐷 . On the other hand, if we let λ1 = 3
−1 1 −2 4 1 2 0 3
1 1
λ2 = 2 then 𝑋1 = [ ] and 𝑋2 = [ ] . Thus
2 1
1 1 −1 1
P=[ ] and 𝑃−1 = [ ] . Hence
2 1 2 −1
−1 1 1 1 1 1 3 0
𝑃−1 𝐴𝑃 = [ ][ ][ ]=[ ] = 𝐷.
2 −1 −2 4 2 1 0 2
1 1
Example 2. Let A=[ ] . Then the eigenvalues are λ1 = 1 λ2 = 1. Eigenvectors
0 1
𝑟
associated with λ1 and λ2 are vectors of the form [ ], where r is any nonzero real
0
number. Since A does not have two linearly independent eigenvectors, thus A is not
diagonalizable.
Theorem A matrix is diagonalizable if all the roots of its characteristic polynomial are real
and distinct.
0 0 1
Example 3. Let A = [0 1 2]. Is A diagonalizable?
0 0 1
Solution:
λ 0 − 1
𝑓(λ) = |λ𝐼𝑛 − 𝐴| = | 0 λ − 1 − 2 |
0 0 λ−1
= λ ( λ − 1 ) (λ − 1)
= λ( λ − 1 ) 2
1 0 − 1 𝑥1 0
𝑥
[ 0 0 − 2 ] [ 2 ] = [0].
0 0 0 𝑥3 0
1 0 −1⋮ 0 1 0 0⋮ 0 1 0 0⋮ 0 0
[0 0 −2⋮ 0] → [0 0 1⋮ 0] → [0 0 0⋮ 0] → [𝑟 ]
0 0 0⋮ 0 0 0 0⋮ 0 0 0 1⋮ 0 0
0
Thus all eigenvectors associated with the eigenvalue λ2 = 1 are [𝑟 ], where r is any
0
nonzero real number. The dimension of the solution space of (1𝐼3 − 𝐴)𝑋 = 0 is 1.
There do not exist two linearly independent eigenvectors associated with λ2 = 1 . Hence,
A cannot be diagonalized.
0 0 0
Example 4. Let A = [0 1 0]. Is A diagonalizable?
1 0 1
Solution:
0 0 0 𝑥1 0
[ 0 − 1 0 ] [𝑥2 ] = [0].
−1 0 − 1 𝑥3 0
0 0 0⋮ 0 0 0 0⋮ 0 𝑟
[ 0 1 0 ⋮ 0] → [0 1 0⋮ 0] → [0]
−1 0 −1⋮ 0 1 0 1⋮ 0 −𝑟
𝑟
Thus all eigenvectors associated with the eigenvalue λ1 = 0 are [ 0 ], where r is any
−𝑟
1
nonzero real number. Thus 𝑋1 = [ 0 ] is an eigenvector associated with λ1 = 0.
−1
1 0 0 𝑥1 0
[ 0 0 0 ] [𝑥2 ] = [0].
−1 0 0 𝑥3 0
1 0 0⋮ 0 1 0 0⋮ 0 0
[ 0 0 0⋮ 0] → [0 0 0⋮ 0] → [𝑟 ]
−1 0 0⋮ 0 0 0 0⋮ 0 𝑠
0
Thus all eigenvectors associated with the eigenvalue λ2 = 1 are [𝑟 ], where r and s are
𝑠
0 0
any nonzero real number. Hence we take the eigenvectors 𝑋2 = [1] and 𝑋3 = [0].
0 1
1 0 0
𝑋1 = [ 0 ] , 𝑋2 = [1] and 𝑋3 = [0] are linearly independent. Hence A is diagonalizable.
−1 0 1
2 3 0
Example 5. Let A = [0 1 0].
0 0 2
Solution:
−1 − 3 0 𝑥1 0
[ 0 0 0 ] [𝑥2 ] = [0].
0 0 − 1 𝑥3 0
−1 −3 0 ⋮ 0 1 3 0⋮ 0 −3𝑠
[0 0 0 ⋮ 0] → [0 0 0⋮ 0] → [ 𝑠 ]
0 0 −1⋮ 0 0 0 1⋮ 0 0
−3𝑠
Thus the eigenspace associated with the eigenvalue λ1 = 1 is [ 𝑠 ], where s is any
0
−3
nonzero real number. Tthus 𝑋1 = [ 1 ] is the basis for the eigenspace associated with
0
λ1 = 1.
0 − 3 0 𝑥1 0
[ 0 1 0 ] [𝑥2 ] = [0].
0 0 0 𝑥3 0
0 −3 0⋮ 0 0 0 0⋮ 0 𝑟
[0 1 0⋮ 0] → [0 1 0⋮ 0] → [0]
0 0 0⋮ 0 0 0 0⋮ 0 𝑠
𝑟
Thus the eigenspace associated with the eigenvalue λ2 = 2 is [0], where r and s are any
𝑠
1 0
nonzero real number. Thus 𝑋2 = [0] and 𝑋2 = [0] are the basis for the eigenspace
0 1
associated with λ2 = 2.
𝑋1 𝑋2 𝑋3
−3 1 0⋮ 0 1 0 0⋮ 0
c. Since [ 1 0 0⋮ 0] → [0 1 0⋮ 0] which means 𝑐1 = 𝑐2 = 𝑐3 = 0, then
0 0 1⋮ 0 0 0 1⋮ 0
−3 1 0
𝑋1 = [ 1 ] , 𝑋2 = [0] and 𝑋3 = [0] are linearly independent. Hence A is
0 0 1
𝑋1 𝑋2 𝑋3
−3 1 0 0 1 0
diagonalizable. Then P=[ 1 0 0] and 𝑃−1 = [1 3 0]
0 0 1 0 0 1
λ1 0 0 1 0 0
D = [ 0 λ2 0 ] = [0 2 0] . Hence
0 0 λ3 0 0 2
0 1 0 2 3 0 −3 1 0 1 0 0
𝑃−1 𝐴𝑃 = [1 3 0] [0 1 0] [ 1 0 0 ] = [0 2 0] = 𝐷
0 0 1 0 0 2 0 0 1 0 0 2
Exercises
A. Tell whether the following matrices are diagonalizable or not.
1 2 3
1 4
1. [ ] 4. [0 −1 2]
1 −2
0 0 2
3 1 0
1 0
2. [ ] 5. [0 3 1]
−2 1
0 0 3
1 1 −2
3. [4 0 4]
1 −1 4
3 −2 1
0 −1
1. [ ] 4. [0 2 0]
2 3
0 0 0
4 2 3 1 1 2
2. [ 2 1 2] 5. [0 1 0]
−1 −2 0 0 1 3
1 2 3
3. [0 1 0]
2 1 2
2 1 1 1 2 2 3 4
1 2 1 1 0 2 3 2
1. [ ] 2. [ ]
1 1 2 1 0 0 1 1
0 0 0 1 0 0 0 1
Remarks
1. All roots of the characteristic polynomial of symmetric matrix are real numbers.
2. If A is a symmetric matrix all of whose eigenvalues are distinct, then A is
diagonalizable.
3. If A is a symmetric matrix, the eigenvectors that belong to distinct eigenvalues of A
are orthogonal.
2 −2 1
3 3 3
2 1 −2
Example. Show that A =
3 3 3
is an orthogonal matrix.
1 2 2
[3 3 3]
Solution:
2 2 1
3 3 3
𝑇 −2 1 2
𝐴 =
3 3 3
1 −2 2
[ 3 3 3]
2 2 1 2 −2 1
3 3 3 3 3 3 1 0 0
−2 1 2 2 1 −2
𝐴𝑇 𝐴 = = [0 1 0] = 𝐼3 .
3 3 3 3 3 3
1 −2 2 1 2 2 0 0 1
[3 3 3] [3 3 3]
Theorem. The nxn matrix A is orthogonal if and only if the columns (and rows) of A form
an orthogonal set of vectors in 𝑅 𝑛 .
Theorem. If A is a symmetric nxn matrix., then there exists an orthogonal matrix P such
that
𝑃−1 𝐴𝑃 = D, a diagonal matrix. The eigenvalues of A lie on the main diagonal of D.
0 0 −2
Example 1. Let A = [ 0 −2 0 ] . Then the characteristic polynomial of A is
−2 0 3
𝑓(λ) = (λ + 2)(λ − 4)(λ + 1). Thus the eigenvalues are λ1 = −2, λ2 = 4, λ3 = −1. Hence
0 −1 2
the eigenvectors are 𝑋1 = [1], 𝑋2 = [ 0 ], and 𝑋3 = [0]. Since
0 2 1
𝑋1 ∙ 𝑋2 ∙ 𝑋3 = 0(−1)(2) + 1(0)(0) + 0(2)(1) = 0, { 𝑋1 , 𝑋2 , 𝑋3 } is orthogonal. If we normalize
1 2
0 − √5 √5
{ 𝑋1 , 𝑋2 , 𝑋3 }, we find that T = {[1] , [ 0 ] , [ 0 ]} is an orthonormal set of vectors. The
0 2 1
√5 √5
matrix P such that 𝑃−1 𝐴𝑃 is diagonal is the matrix whose columns are the vectors in T.
1 2 0 0
0 − √5 √5 1 1 2
Thus P = [1 0 0 ]. Since P is orthogonal, 𝑃 −1 = 𝑃𝑇 = [−√5 0 √5]. Hence,
0 2 1 2
0 1
√5 √5 √5 √5
0 0 1 2
1 1 2 0 0 −2 0 − √5 √5 −2 0 0
𝑃−1 𝐴𝑃 = 𝑃𝑇 𝐴𝑃 = [−√5 0 √5] [ 0 −2 0 ] [1 0 0 ] = [ 0 4 0 ] = 𝐷 .
2 1
0 −2 0 3 0 2 1 0 0 −1
√5 √5 √5 √5
0 2 2
.Example 2. Diagonalize matrix A = [2 0 2] by an orthogonal matrix.
2 2 0
Solution:
To find the eigenvectors associated with λ1 and λ2 , we solve the homogeneous linear
system (−2𝐼3 − 𝐴)𝑋 = 0:
−2 − 2 − 2 𝑥1 0
solution space: [−2 − 2 − 2] [𝑥2 ] = [0]
−2 − 2 − 2 𝑥3 0
−1 −1
𝑋1 = [ 1 ] and 𝑋2 = [ 0 ]
0 1
Since 𝑋1 ∙ 𝑋2 = (−1)(−1) + 1(0) + (0)1 ≠ 0, 𝑋1 𝑎𝑛𝑑 𝑋2 are not orthogonal.
We can use the Gram-Schmidt process to obtain an orthonormal basis for the solution
space above .
−1
Let 𝑌1 = 𝑋1 = [ 1 ]
0
and
1
−1 −1 −2
𝑋2 ∙𝑌1 −1(−1)+1(0)+0(1)
𝑌2 = 𝑋2 − ( 𝑌 ∙𝑌 ) 𝑌1 = [ 0 ] − ( (−1)2+12 +02 ) [ 1 ] = [− 1].
1 1 2
1 0
1
−1
Let 𝑌2 = 2𝑌2 = [−1] . The set {𝑌1 , 𝑌2 ∗ } is an orthogonal set of vectors.
∗
2
Normalizing these eigenvectors, we have
1
−
1 −
−1 −1 √6
√2
𝑌 1 𝑌∗ 1 1
𝑍1 = |𝑌1 | = [ 1 ]=[ 1 ] and 𝑍2 = |𝑌2 ∗ | = [−1] = − 6.
1 √2 2 √6 √
0 √2 2 2
0 [ √6 ]
The set {𝑍1 , 𝑍2 } is an orthonormal basis of eigenvectors of A for the solution space above.
4 − 2 − 2 𝑥1 0
[ −2 4 − 2 ] [𝑥2 ] = [0],
−2 − 2 4 𝑥3 0
1
1 1 √3
1 1
is consist of 𝑋3 = [1] and normalizing this vector, we have 𝑍3 = [1] = as an
√3 √3
1 1 1
[√3]
orthonormal basis for (4𝐼3 − 𝐴)𝑋 = 0.
−1 1 −1 1
√2 √2
0 −1
√6 √3
−1 −1
2 0 2 2 √2 −1 1
−2 0 0
Therefore, 𝑃−1 𝐴𝑃 = √6 [2 0 2] 1 = [ 0 −2 0] = 𝐷.
√6 √6 1 √6 √3
1 1 2 2 0 √2 2 1 0 0 4
[√3 √3
√3] [0 √6 √3]
1 2 0 0
2 1 0 0
Example 3. Diagonalize A=[ ] and find an orthogonal matrix P such that
0 0 1 2
0 0 2 1
𝑃−1 𝐴𝑃 is diagonal.
Solution:
The basis for the solution space of (−1𝐼3 − 𝐴)𝑋 = 0 consist of the eigenvectors
1 0
−1 0
𝑋1 = [ ] and 𝑋2 = [ ] which are orthogonal.
0 1
0 −1
Normalizing 𝑋1 and 𝑋2 , we have
1 0
√2 0
1 1
𝑍1 = − and 𝑍2 = as an orthonormal basis of eigenvectors for the solution
√2 √2
0 1
[ 0 ] [− √2]
The basis for the solution space of (3𝐼3 − 𝐴)𝑋 = 0 consist of the eigenvectors
1 0
1 0
𝑋3 = [ ] and 𝑋4 = [ ] which are orthogonal.
0 1
0 1
Since eigenvectors associated with distinct eigenvalues are orthogonal, thus {𝑍1 , 𝑍2 , 𝑍3 , 𝑍4 }
is an orthonormal basis for 𝑅 4 consisting of eigenvectors of A. Hence, matrix
Remark. Suppose that A is an nxn matrix for which we can find an orthogonal matrix P
such that 𝑃−1 𝐴𝑃 is a diagonal matrix D. Thus 𝑃 −1 𝐴𝑃 = 𝐷 or A = PD𝑃−1. Since
𝑃−1 = 𝑃𝑇 , we can write A = PD𝑃𝑇 . Then A𝑇 = (PD𝑃𝑇 )𝑇 = (𝑃𝑇 )𝑇 𝐷 𝑇 𝑃𝑇 = 𝑃𝐷𝑃𝑇 = 𝐴
(𝐷 = 𝐷 𝑇 , since D is a diagonal matrix). Thus A is symmetric.
Exercises
Diagonalize the following matrix and find an orthogonal matrix P such that 𝑃−1 AP is
diagonal.
0 0 0 0
−1 2 2
0 0 0 0
1. [ 2 −1 2 ] 7. [ ]
0 0 0 0
2 2 −1
0 0 1 0
2 2 0 0
0 −1 −1
2 2 0 0
2. [−1 0 −1] 8. [ ]
0 0 2 2
−1 −1 0
0 0 2 2
0 0 0 1
1 1 0
0 0 0 0
3. [1 1 0] 9. [ ]
0 0 0 0
0 0 1
1 0 0 0
1 0 0
4. [0 1 1]
0 1 1
2 1 1
5. [1 2 1]
1 1 2
−3 0 −1
6. [ 0 −2 0 ]
−1 0 −3