← Back to Math 010

5.2 Diagonalization: Class Notes

Definition

A square matrix A is said to be diagonalizable if it is similar to some diagonal matrix; that is, if there exists an invertible matrix P such that P−1AP is diagonal. In this case, the matrix P is said to diagonalize A.

Theorem 5.2.1

If A is an n×n matrix, then A is diagonalizable if and only if A has n linearly independent eigenvectors.

#5

Find a matrix P that diagonalizes A, and check your work by computing P−1AP.

A=[106−1]

Solution

det(λI−A)=0 ⇒ |λ−10−6λ+1|=0 ⇒ λ2−1=0 ⇒ λ=±1

Find an eigenvector for each eigenvalue.

λ1=−1: [−20−60] → [1000] ⇒ x1=0, so p→1=[01]t

λ2=1: [00−62] ⇒ 3x1=x2, so p→2=[13]t

Form the matrix whose columns are the eigenvectors, P=[p→1|p→2]:

P=[0113]

Find P−1 (here det(P)=−1):

P−1=1det(P)[3−1−10]=[−3110]

Find P−1AP:

P−1AP=[−3110][106−1][0113]=[3−110][0113]=[−1001]

This is a diagonal matrix.

Answer: P=[0113],P−1AP=[−1001]

Notes

Note

Listing the eigenvectors in the other order, P=[1031] would also diagonalize A, with P−1AP=[100−1]. So P is not unique.

What makes this work?

P−1AP=D ⇒ AP=PD

Here,

[106−1][0113]=[0113][−1001]

Column by column:

Ap→1=[106−1][01]=[0−1]=−1[01]

Ap→2=[106−1][13]=[13]=1[13]

So [Ap→1|Ap→2]=[λ1p→1|λ2p→2]: corresponding columns satisfy Ax→=λx→.

Proof of Theorem 5.2.1

(⇒) If A is diagonalizable, then there is an invertible matrix P and a diagonal matrix D such that P−1AP=D (by the definition). Then AP=PD. Let

P=[p→1p→2⋯p→n]andD=[λ10⋱0λn].

Then

AP=A[p→1p→2⋯p→n]=[Ap→1Ap→2⋯Ap→n]

and

PD=[p→1p→2⋯p→n][λ10⋱0λn]=[λ1p→1λ2p→2⋯λnp→n]

Matching columns, AP=PD gives Ap→i=λip→i, 1≤i≤n.

Since P is invertible, {p→i} is a linearly independent set (Equivalent Statements), so each p→i≠0→. So the n column vectors p→i are linearly independent eigenvectors of A.

(⇐) Suppose A has linearly independent eigenvectors p→1,p→2,…,p→n with corresponding eigenvalues λ1,λ2,…,λn. Let P=[p→1p→2⋯p→n] and D=[λ10⋱0λn]. Then

AP=[Ap→1Ap→2⋯Ap→n]=[λ1p→1λ2p→2⋯λnp→n]=PD

{p→i} is an independent set, so P is invertible. Then AP=PD ⇒ P−1AP=D, so A is diagonalizable.

Tip: The statement is about the total number of eigenvectors, not eigenvalues.

#17

Compute the matrix A10 by first diagonalizing A.

A=[032−1]

Solution

The process: (1) find the eigenvalues λi, (2) find the eigenvectors p→i, (3) form P and P−1.

|λ−3−2λ+1|=0 ⇒ λ2+λ−6=0 ⇒ λ=−3, 2

λ1=−3: [−3−3−2−2] → [1100] ⇒ p→1=[1−1]t

λ2=2: [2−3−23] → [1−3200] ⇒ p→2=[32]t

P=[13−12],det(P)=5 ⇒ P−1=[25−351515]

P−1AP=D=[λ100λ2]=[−3002]

Note that if P−1AP=D, then A=PDP−1. Then the inner factors cancel in pairs:

Ak=(PDP−1)(PDP−1)⋯(PDP−1)=PDkP−1

Recall

D=[d10⋱0dn] ⇒ Dk=[d1k0⋱0dnk]

Here D10=[(−3)1000210]=[31000210], so

A10=[13−12][31000210][25−351515]

Answer: A10=[24,234−34,815−23,21035,839]

Theorem 5.2.3

If k is a positive integer, λ is an eigenvalue of a matrix A, and 𝐱 is a corresponding eigenvector, then λk is an eigenvalue of Ak and 𝐱 is a corresponding eigenvector.

(Not a proof.) Consider Ax→=λx→. Then

A2x→=A(λx→)=λAx→=λ2x→

Example

Let A=[324202423].

a. Find the eigenvalues of A.
b. For each eigenvalue λ, find the rank of the matrix λI−A.
c. Is A diagonalizable? Justify your conclusion.

Solution

a. Expand det(λI−A) along the first row.

|λ−3−2−4−2λ−2−4−2λ−3|=0

(λ−3)|λ−2−2λ−3|+2|−2−2−4λ−3|−4|−2λ−4−2|=0(λ−3)(λ2−3λ−4)+2[−2(λ−3)−8]−4(4+4λ)=0(λ−3)(λ−4)(λ+1)−4(λ+1)−16(λ+1)=0(λ+1)[(λ−3)(λ−4)−4−16]=0(λ+1)(λ2−7λ−8)=0(λ+1)2(λ−8)=0

Each term has a factor of (λ+1), which we factor out.

Answer: λ=−1 (algebraic multiplicity 2), λ=8

b. For λ1=8:

[5−2−4−28−2−4−25] → [5−2−402−1000]

Using R2→5R2+2R1 and R3→R3−2R2:

5R2+2R1:(−10,40,−10)+(10,−4,−8)=(0,36,−18)R3−2R2:(−4,−2,5)+(4,−16,4)=(0,−18,9)

Both new rows are multiples of (0,2,−1), which gives the reduced matrix above.

The rank is 2, so the nullity is 1. This will give 1 eigenvector.

For λ2=−1:

[−4−2−4−2−1−2−4−2−4] → [212000000]

The rank is 1, so the nullity is 2. This will give 2 eigenvectors.

Answer: rank(8I−A)=2,rank(−I−A)=1

c. 2 + 1 = 3 total eigenvectors (the total number of basis vectors for the eigenspaces is 3), so

Answer: A is diagonalizable

Find the eigenvectors and P; find P−1AP.

λ1=8: [5−2−4−28−2−4−25] → [5−2−402−1000] → [10−102−1000]

So x1=x3 and 2x2=x3, giving x→1=[212]t and p→1=[212].

λ2=−1: [−4−2−4−2−1−2−4−2−4] → [212000000] ⇒ x1=−12x2−x3

x→2=[−12s−tst]=[−1210]s+[−101]t

Scaling the first vector by −2: p→2=[1−20], p→3=[−101].

P=[21−11−20201]

Answer: P−1AP=[8000−1000−1]

Why this works: PD multiplies each column p→i of P by its eigenvalue, so

PD=[21−11−20201][8000−1000−1]=[8p→1−1p→2−1p→3]

and each column is λp→i=Ap→i, which is Ax→=λx→.

Definitions

If λ0 is an eigenvalue of an n×n matrix A, then the dimension of the eigenspace corresponding to λ0 is called the geometric multiplicity of λ0, and the number of times that λ−λ0 appears as a factor in the characteristic polynomial of A is called the algebraic multiplicity of λ0.

Theorem 5.2.4 Geometric and Algebraic Multiplicity

If A is a square matrix, then:

#14

Find the geometric and algebraic multiplicity of each eigenvalue of the matrix A, and determine whether A is diagonalizable. If A is diagonalizable, then find a matrix P that diagonalizes A, and find P−1AP.

A=[500150015]

Solution

A is lower triangular, so the characteristic equation is (λ−5)3=0. So λ=5 is an eigenvalue with algebraic multiplicity 3. If we can find 3 linearly independent eigenvectors, we can diagonalize A.

5I−A=[000−1000−10]

This has rank 2, so its nullity is 1: it yields only [001] (up to scalar multiples), and the geometric multiplicity is 1.

So λ=5 has algebraic multiplicity 3 and geometric multiplicity 1. There are not enough eigenvectors.

Answer: A is not diagonalizable

Theorem 5.2.2

Similarity

Definition

A transformation of the form A→P−1AP is called a similarity transformation.

Definition

If A and B are square matrices, then we say that 𝐁 is similar to 𝐀 (or that A and B are similar matrices) if there is an invertible matrix P such that B=P−1AP.

Example

Prove that if A and B are similar matrices, then det(A)=det(B).

Proof

Since A and B are similar, there is an invertible P with B=P−1AP. Then

det(B)=det(P−1AP)=det(P−1)det(A)det(P)=1det(P)det(A)det(P)=det(A)

Definition

Any property that is preserved by a similarity transformation is called a similarity invariant and is said to be invariant under similarity.

Similarity invariants

The following are the same for A and P−1AP: