Processing math: 100%
Skip to main content

Section 29.5 Examples

Example 29.5.1. Scalar-triangularization procedure: a small matrix.

Is matrix

A=[6145121251]

similar to some upper triangular matrix?

We begin by investigating the eigenvalues of A. Compute the characteristic polynomial:

λIA=[λ+61451λ2125λ1],
cA(λ)=det(λIA)=λ33λ23λ1=(λ+1)3.

So A has eigenvalue λ=1, with algebraic multiplicity 3. From our analysis in Subsection 29.4.3, since A has a single eigenvalue, we are confident it can be made similar to a scalar triangular matrix.

Next, compute a basis for the eigenspace E1(A).

(1)IA=[5145131252]rowreduce[101010000]E1(A)=Span{[101]}.

As you can see, the geometric multiplicity of the eigenvalue λ=1 is not equal to the algebraic multiplicity, so A is not diagonalizable. So we continue with the scalar-triangularization procedure (Procedure 29.4.1).

Now compute a basis for the generalized eigensubspace of degree 2:

(IA)2=[131000131]rowreduce[131000000]E21(A)=Span{[101],[310]}.

Note that we have included our basis vector for E1(A) as the first vector in our basis for E21(A).

We are only up to two linearly independent vectors, so we need to keep going by computing a basis for E31(A). But if you compute (IA)3 you will find that it is the zero matrix. What does this mean? Every vector in R3 is in the null space of the zero matrix, so we have E31(A)=R3.

However, we don't want to choose just any basis for R3, as we need to maintain independence from the basis vectors for E21(A) that we already have. That is, we need to choose a vector in R3 that is not already in the span of our existing two generalized eigenvectors (Proposition 18.5.6). But since dimE21(A)=2, we know that it is not possible for all three of the standard basis vectors for R3 to be in E21(A), and so we can just choose any of them that is not. We will choose the third standard basis vector.

We now have our column vectors for a suitable transition matrix P:

p1=[101],p2=[310],p3=[001],P=[130010101].

Finally, let's determine the form matrix P1AP. We already know that Ap1=p1, since p1 is an eigenvector of A corresponding to the eigenvalue λ=1. Next, compute (A(1)I)p2 and you will find that it is equal to p1. Thus, the entry above the diagonal in the second column must be 1. Lastly, compute

[p1p2(A(1)I)p3]=[135011102]rowreduce[102011000].

Thus, the entries above the diagonal in the third column are 2, 1. Putting all this together yields

P1AP=[112011001].
Example 29.5.2. Scalar-triangularization procedure: a larger matrix.

Is matrix

A=[31883250252732431115728209]

similar to some upper triangular matrix?

Again, we start by investigating eigenvalues. If you compute the characteristic polynomial for A, you will find that it factors as cA(λ)=(λ3)5, so that A has a single eigenvalue 3 of multiplicity 5. So from our analysis in Subsection 29.4.3, we are confident that A can be made similar to a scalar triangular matrix.

Next, compute a basis for the eigenspace E3(A).

3IA=[01883280252702431112728206]rowreduce[100132010010011120000000000]E3(A)=Span{[10110],[32102]}

We need five vectors for the columns of the transition matrix P, so this is not enough.

Next, compute a basis for E23(A).

(3IA)2=[0022100442006630088400442]rowreduce[00111200000000000000000000]E23(A)=Span{[10000],[01000],[00110],[00102]}

The above spanning vectors are the ones arrived at from row reducing (3IA)2. But remember that we want to use our already obtained vectors from E3(A), so we just need to choose two of these four that are independent from our eigenvectors. Those first two standard basis vectors will work, so take

E23(A)=Span{[10110],[32102],[10000],[01000]}.

Still not enough vectors, but (3IA)3=0, so we can choose any vector in R5 that is not in the span of the four we already have. The third standard basis vector for R5 will work, and we now have

E33(A)=Span{[10110],[32102],[10000],[01000],[00100]}.

Now place these matrices in the transition matrix

P=[1310002010110011000002000].

To determine the form matrix, we can reduce [PAP] as

[131003931802010062501100133273100003031110200006282]rowreduce[1000030311101000031410010000301000010000320000100003],

obtaining

P1AP=[303111031410030100003200003].