Every problem in the chapter beside the book’s answer, a chapter at a time as you scroll.

Chapter 5: Similarity

166 problems · 2 figures

Five.II.1.5 §II.1 source ↗

For

T=(13−2−6)Tˆ=(00−11/2−5)P=(42−32)

check that Tˆ=PTP−1⁠.

Answer

One way to proceed is left to right.

PTP−1=(42−32)(13−2−6)(2/14−2/143/144/14)=(00−7−21)(2/14−2/143/144/14)=(00−11/2−5)

Five.II.1.6 §II.1 source ↗

Example 1.4 shows that the only matrix similar to a zero matrix is itself and that the only matrix similar to the identity is itself.

(a)

Show that the 1×1 matrix whose single entry is 2 is also similar only to itself.

(b)

Is a matrix of the form cI for some scalar c similar only to itself?

(c)

Is a diagonal matrix similar only to itself?

Answer (3 parts)

(a)

Because the matrix (2) is 1×1⁠, the matrices P and P−1 are also 1×1 and so where P=(p) the inverse is P−1=(1/p)⁠. Thus P(2)P−1=(p)(2)(1/p)=(2)⁠.

(b)

Yes: recall that we can bring scalar multiples out of a matrix P(cI)P−1=cPIP−1=cI⁠. By the way, the zero and identity matrices are the special cases c=0 and c=1⁠.

(c)

No, as this example shows.

(1−2−11)(−100−3)(−1−2−1−1)=(−5−421)

Five.II.1.7 §II.1 source ↗

Consider this transformation of ℂ3

t((xyz))=(x−zz2y)

and these bases.

B=⟨(123),(010),(001)⟩D=⟨(100),(110),(101)⟩

We will compute the parts of the arrow diagram to represent the transformation using two similar matrices.

(a)

Draw the arrow diagram, specialized for this case.

(b)

Compute T=RepB,B⁡(t)⁠.

(c)

Compute Tˆ=RepD,D⁡(t)⁠.

(d)

Compute the matrices for the other two sides of the arrow square.

Answer (4 parts)

(a)

ℂwrt B3→Ttℂwrt B3id⁡↓id⁡↓ℂwrt D3→Tˆtℂwrt D3

(b)

For each element of the starting basis B find the effect of the transformation

(123)⟼t(−234)(010)⟼t(002)(001)⟼t(−110)

and represented those outputs with respect to the ending basis B

RepB⁡((−234))=(−2710)RepB⁡((002))=(002)RepB⁡((−110))=(−133)

to get the matrix.

T=RepB,B⁡(t)=(−20−17031023)

(c)

Find the effect of the transformation on the elements of D

(100)⟼t(100)(110)⟼t(102)(101)⟼t(010)

and represented those with respect to the ending basis D

RepD⁡((100))=(100)RepD⁡((102))=(−102)RepD⁡((010))=(−110)

to get the matrix.

Tˆ=RepD,D⁡(t)=(1−1−1001020)

(d)

To go down on the right we need RepB,D⁡(id⁡) so we first compute the effect of the identity map on each element of B⁠, which is no effect, and then represent the results with respect to D⁠.

RepD⁡((123))=(−423)RepD⁡((010))=(−110)RepD⁡((001))=(−101)

So this is P⁠.

P=(−4−1−1210301)

For the other matrix RepD,B⁡(id⁡) we can either find it directly, as we just have with P⁠, or we can do the usual calculation of a matrix inverse.

P−1=(111−2−1−2−3−3−2)

Five.II.1.8 §II.1 source ↗

Consider the transformation t:𝒫2→𝒫2 described by x2↦x+1⁠, x↦x2−1⁠, and 1↦3⁠.

(a)

Find T=RepB,B⁡(t) where B=⟨x2,x,1⟩⁠.

(b)

Find Tˆ=RepD,D⁡(t) where D=⟨1,1+x,1+x+x2⟩⁠.

(c)

Find the matrix P such that Tˆ=PTP−1⁠.

Answer (3 parts)

(a)

Because we describe t with the members of B⁠, finding the matrix representation is easy:

RepB⁡(t(x2))=(011)BRepB⁡(t(x))=(10−1)BRepB⁡(t(1))=(003)B

gives this.

RepB,B⁡(t)(0101001−13)

(b)

We will find t(1)⁠, t(1+x)⁠, and t(1+x+x2⁠, to find how each is represented with respect to D⁠. We are given that t(1)=3⁠, and the other two are easy to see: t(1+x)=x2+2 and t(1+x+x2)=x2+x+3⁠. By eye, we get the representation of each vector

RepD⁡(t(1))=(300)DRepD⁡(t(1+x))=(2−11)DRepD⁡(t(1+x+x2))=(201)D

and thus the representation of the map.

RepD,D⁡(t)=(3220−10011)

(c)

The diagram

Vwrt B→TtVwrt Bid⁡↓Pid⁡↓PVwrt D→TˆtVwrt D

shows that these are P=RepB,D⁡(id⁡) and P−1=RepD,B⁡(id⁡)⁠.

P=(0−11−110100)P−1=(001011111)

Five.II.1.9 §II.1 source ↗

Let T represent t:ℂ2→ℂ2 with respect to B,B⁠.

T=(1−121)B=⟨(10),(11)⟩,D=⟨(20),(0−2)⟩

We will convert to the matrix representing t with respect to D,D⁠.

(a)

Draw the arrow diagram.

(b)

Give the matrix that represents the left and right sides of that diagram, in the direction that we traverse the diagram to make the conversion.

(c)

Find RepD,D⁡(t)⁠.

Answer (3 parts)

(a)

ℂwrt B2→Ttℂwrt B2id⁡↓id⁡↓ℂwrt D2→Tˆtℂwrt D2

(b)

For the right side we find the effect of the identity map, which is no effect,

(10)⟼id⁡(10)(11)⟼id⁡(11)

and represent those with respect to D

RepD⁡((10))=(1/20)RepD⁡((11))=(1/2−1/2)

so we have this.

P=RepB,D⁡(id⁡)=(1/21/20−1/2)

For the matrix on the left we can either compute it directly, as in the prior paragraph, or we can take the inverse.

P−1=RepD,B⁡(id⁡)=1(−1/4)⋅(−1/2−1/201/2)=(220−2)

(c)

As with the prior item we can either compute it directly from the definition or compute it using matrix operations.

PTP−1=(220−2)(1−121)(220−2)=(33−2−1)

Five.II.1.10 §II.1 source ↗

Exhibit a nontrivial similarity relationship by letting t:ℂ2→ℂ2 act in this way,

(12)↦(30)(−11)↦(−12)

picking two bases B,D⁠, and representing t with respect to them, Tˆ=RepB,B⁡(t) and T=RepD,D⁡(t)⁠. Then compute the P and P−1 to change bases from B to D and back again.

Answer

One possible choice of the bases is

B=⟨(12),(−11)⟩D=ℰ2=⟨(10),(01)⟩

(this B comes from the map description). To find the matrix Tˆ=RepB,B⁡(t)⁠, solve the relations

c1(12)+c2(−11)=(30)cˆ1(12)+cˆ2(−11)=(−12)

to get c1=1⁠, c2=−2⁠, cˆ1=1/3 and cˆ2=4/3⁠.

RepB,B⁡(t)=(11/3−24/3)

Finding RepD,D⁡(t) involves a bit more computation. We first find t(e→1)⁠. The relation

c1(12)+c2(−11)=(10)

gives c1=1/3 and c2=−2/3⁠, and so

RepB⁡(e→1)=(1/3−2/3)B

making

RepB⁡(t(e→1))=(11/3−24/3)B,B(1/3−2/3)B=(1/9−14/9)B

and hence t acts on the first basis vector e→1 in this way.

t(e→1)=(1/9)⋅(12)−(14/9)⋅(−11)=(5/3−4/3)

The computation for t(e→2) is similar. The relation

c1(12)+c2(−11)=(01)

gives c1=1/3 and c2=1/3⁠, so

RepB⁡(e→1)=(1/31/3)B

making

RepB⁡(t(e→1))=(11/3−24/3)B,B(1/31/3)B=(4/9−2/9)B

and hence t acts on the second basis vector e→2 in this way.

t(e→2)=(4/9)⋅(12)−(2/9)⋅(−11)=(2/32/3)

Therefore

RepD,D⁡(t)=(5/32/3−4/32/3)

and so this matrix.

P=RepB,D⁡(id⁡)=(1−121)

and this one change the bases.

P−1=(RepB,D⁡(id⁡))−1=(1−121)−1=(1/31/3−2/31/3)

The check of these computations is routine.

(1−121)(11/3−24/3)(1/31/3−2/31/3)=(5/32/3−4/32/3)

Five.II.1.11 §II.1 source ↗

Show that these matrices are not similar.

(104113217)(101011312)

Answer

Gauss’s Method shows that the first matrix represents maps of rank two while the second matrix represents maps of rank three.

Five.II.1.12 §II.1 source ↗

Explain Example 1.4 in terms of maps.

Answer

The only representation of a zero map is a zero matrix, no matter what the pair of bases RepB,D⁡(z)=Z⁠, and so in particular for any single basis B we have RepB,B⁡(z)=Z⁠. The case of the identity is slightly different: the only representation of the identity map, with respect to any B,B⁠, is the identity RepB,B⁡(id⁡)=I⁠. (Remark: of course, we have seen examples where B≠D and RepB,D⁡(id⁡)≠I — in fact, we have seen that any nonsingular matrix is a representation of the identity map with respect to some B,D⁠.)

Five.II.1.13 §II.1 source ↗

[Halmos] Are there two matrices A and B that are similar while A2 and B2 are not similar?

Answer

No. If A=PBP−1 then A2= (PBP−1)(PBP−1)= PB2P−1⁠.

Five.II.1.14 §II.1 source ↗

Prove that if two matrices are similar and one is invertible then so is the other.

Answer

Matrix similarity is a special case of matrix equivalence (if matrices are similar then they are matrix equivalent) and matrix equivalence preserves nonsingularity.

Five.II.1.15 §II.1 source ↗

Show that similarity is an equivalence relation. (The definition given earlier already reflects this, so instead start here with the definition that Tˆ is similar to T if Tˆ=PTP−1⁠.)

Answer

A matrix is similar to itself; take P to be the identity matrix: P=IPI−1=IPI⁠.

If Tˆ is similar to T then Tˆ=PTP−1 and so P−1TˆP=T⁠. Rewrite T=(P−1)Tˆ(P−1)−1 to conclude that T is similar to Tˆ⁠.

For transitivity, if T is similar to S and S is similar to U then T=PSP−1 and S=QUQ−1⁠. Then T=PQUQ−1P−1=(PQ)U(PQ)−1⁠, showing that T is similar to U⁠.

Five.II.1.16 §II.1 source ↗

Consider a matrix representing, with respect to some B,B⁠, reflection across the x⁠-axis in ℝ2⁠. Consider also a matrix representing, with respect to some D,D⁠, reflection across the y⁠-axis. Must they be similar?

Answer

Let fx and fy be the reflection maps (sometimes called ‘flip’s). For any bases B and D⁠, the matrices RepB,B⁡(fx) and RepD,D⁡(fy) are similar. First note that

S=Repℰ2,ℰ2⁡(fx)=(100−1)T=Repℰ2,ℰ2⁡(fy)=(−1001)

are similar because the second matrix is the representation of fx with respect to the basis A=⟨e→2,e→1⟩⁠:

(100−1)=P(−1001)P−1

where P=RepA,ℰ2⁡(id⁡)⁠.

ℂwrt A2→TfxVℂwrt A2id⁡↓Pid⁡↓Pℂwrt ℰ22→Sfxℂwrt ℰ22

Now the conclusion follows from the transitivity part of Exercise 1.15.

We can also finish without relying on that exercise. Write RepB,B⁡(fx)= QTQ−1= QRepℰ2,ℰ2⁡(fx)Q−1 and RepD,D⁡(fy)= RSR−1= RRepℰ2,ℰ2⁡(fy)R−1⁠. By the equation in the first paragraph, the first of these two is RepB,B⁡(fx)= QPRepℰ2,ℰ2⁡(fy)P−1Q−1⁠. Rewriting the second of these two as R−1⋅RepD,D⁡(fy)⋅R= Repℰ2,ℰ2⁡(fy) and substituting gives the desired relationship

RepB,B⁡(fx)=QPRepℰ2,ℰ2⁡(fy)P−1Q−1=QPR−1⋅RepD,D⁡(fy)⋅RP−1Q−1=(QPR−1)⋅RepD,D⁡(fy)⋅(QPR−1)−1

Thus the matrices RepB,B⁡(fx) and RepD,D⁡(fy) are similar.

Five.II.1.17 §II.1 source ↗

Prove that matrix similarity preserves rank and determinant. Does the converse hold?

Answer

We must show that if two matrices are similar then they have the same rank and the same determinant.

First, we’ve shown that the rank of a map equals the rank of any matrix representing that map, so the rank must be the same for all representations. (Said another way, rank is a property of the map since it is the dimension of the range space.)

As for determinants, |PSP−1|= |P|⋅|S|⋅|P−1|= |P|⋅|S|⋅|P|−1= |S|⁠.

The converse of the statement does not hold. For one thing, there are matrices with the same determinant that are not similar. As an example, consider a nonzero matrix with a determinant of zero. It is not similar to the zero matrix because the zero matrix is similar only to itself, but they have they same determinant. The same thinking gives an example for rank.

Five.II.1.18 §II.1 source ↗

Is there a matrix equivalence class with only one matrix similarity class inside? One with infinitely many similarity classes?

Answer

The matrix equivalence class containing all n×n rank zero matrices contains only a single matrix, the zero matrix. Therefore it has as a subset only one similarity class.

In contrast, the matrix equivalence class of 1×1 matrices of rank one consists of those 1×1 matrices (k) where k≠0⁠. For any basis B⁠, the representation of multiplication by the scalar k is RepB,B⁡(tk)=(k)⁠, so each such matrix is alone in its similarity class. So this is a case where a matrix equivalence class splits into infinitely many similarity classes.

Five.II.1.19 §II.1 source ↗

Can two different diagonal matrices be in the same similarity class?

Answer

Yes, these are similar

(1003)(3001)

since, where the first matrix is RepB,B⁡(t) for B=⟨β→1,β→2⟩⁠, the second matrix is RepD,D⁡(t) for D=⟨β→2,β→1⟩⁠.

Five.II.1.20 §II.1 source ↗

Prove that if two matrices are similar then their k⁠-th powers are similar when k>0⁠. What if k≤0⁠?

Answer

The k⁠-th powers are similar because, where each matrix represents the map t⁠, the k⁠-th powers represent tk⁠, the composition of k⁠-many t⁠’s. (For instance, if T=reptB,B then T2=RepB,B⁡(t∘t)⁠.)

Restated more computationally, if T=PSP−1 then T2= (PSP−1)(PSP−1)= PS2P−1⁠. Induction extends that to all powers.

For the k≤0 case, suppose that S is invertible and that T=PSP−1⁠. Note that T is invertible: T−1=(PSP−1)−1=PS−1P−1⁠, and that same equation shows that T−1 is similar to S−1⁠. Other negative powers are now given by the first paragraph.

Five.II.1.21 §II.1 source ↗

Let p(x) be the polynomial cnxn+⋯+c1x+c0⁠. Show that if T is similar to S then p(T)=cnTn+⋯+c1T+c0I is similar to p(S)=cnSn+⋯+c1S+c0I⁠.

Answer

In conceptual terms, both represent p(t) for some transformation t⁠. In computational terms, we have this.

p(T)=cn(PSP−1)n+⋯+c1(PSP−1)+c0I=cnPSnP−1+⋯+c1PSP−1+c0I=PcnSnP−1+⋯+Pc1SP−1+Pc0P−1=P(cnSn+⋯+c1S+c0)P−1

Five.II.1.22 §II.1 source ↗

List all of the matrix equivalence classes of 1×1 matrices. Also list the similarity classes, and describe which similarity classes are contained inside of each matrix equivalence class.

Answer

There are two equivalence classes, (i) the class of rank zero matrices, of which there is one: 𝒞1={(0)}⁠, and (2) the class of rank one matrices, of which there are infinitely many: 𝒞2={(k)|k≠0}⁠.

Each 1×1 matrix is alone in its similarity class. That’s because any transformation of a one-dimensional space is multiplication by a scalar tk:V→V given by v→↦k⋅v→⁠. Thus, for any basis B=⟨β→⟩⁠, the matrix representing a transformation tk with respect to B,B is (RepB⁡(tk(β→)))=(k)⁠.

So, contained in the matrix equivalence class 𝒞1 is (obviously) the single similarity class consisting of the matrix (0)⁠. And, contained in the matrix equivalence class 𝒞2 are the infinitely many, one-member-each, similarity classes consisting of (k) for k≠0⁠.

Five.II.1.23 §II.1 source ↗

Does similarity preserve sums?

Answer

No. Here is an example that has two pairs, each of two similar matrices:

(1−112)(1003)(2/31/3−1/31/3)=(5/3−2/3−4/37/3)

and

(1−2−11)(−100−3)(−1−2−1−1)=(−5−421)

(this example is not entirely arbitrary because the center matrices on the two left sides add to the zero matrix). Note that the sums of these similar matrices are not similar

(1003)+(−100−3)=(0000)(5/3−2/3−4/37/3)+(−5−421)≠(0000)

since the zero matrix is similar only to itself.

Five.II.1.24 §II.1 source ↗

Show that if T−λI and N are similar matrices then T and N+λI are also similar.

Answer

If N=P(T−λI)P−1 then N=PTP−1−P(λI)P−1⁠. The diagonal matrix λI commutes with anything, so P(λI)P−1=PP−1(λI)=λI⁠. Thus N=PTP−1−λI and consequently N+λI=PTP−1⁠. (So not only are they similar, in fact they are similar via the same P⁠.)

Five.II.2.6 §II.2 source ↗

Repeat Example 2.5 for the matrix from Example 2.2.

Answer

Because we chose the basis vectors arbitrarily, many different answers are possible. However, here is one way to go; to diagonalize

T=(4−211)

take it as the representation of a transformation with respect to the standard basis T=Repℰ2,ℰ2⁡(t) and look for B=⟨β→1,β→2⟩ such that

RepB,B⁡(t)=(λ100λ2)

that is, such that t(β→1)=λ1 and t(β→2)=λ2⁠.

(4−211)β→1=λ1⋅β→1(4−211)β→2=λ2⋅β→2

We are looking for scalars x such that this equation

(4−211)(b1b2)=x⋅(b1b2)

has solutions b1 and b2⁠, which are not both zero. Rewrite that as a linear system

(4−x)⋅b1+−2⋅b2=01⋅b1+(1−x)⋅b2=0

If x=4 then the first equation gives that b2=0⁠, and then the second equation gives that b1=0⁠. We have disallowed the case where both b⁠’s are zero so we can assume that x≠4⁠.

⟶(−1/(4−x))ρ1+ρ2(4−x)⋅b1+−2⋅b2=0((x2−5x+6)/(4−x))⋅b2=0

Consider the bottom equation. If b2=0 then the first equation gives b1=0 or x=4⁠. The b1=b2=0 case is not allowed. The other possibility for the bottom equation is that the numerator of the fraction x2−5x+6=(x−2)(x−3) is zero. The x=2 case gives a first equation of 2b1−2b2=0⁠, and so associated with x=2 we have vectors whose first and second components are equal:

β→1=(11)(so (4−211)(11)=2⋅(11), and λ1=2).

If x=3 then the first equation is b1−2b2=0 and so the associated vectors are those whose first component is twice their second:

β→2=(21)(so (4−211)(21)=3⋅(21), and so λ2=3).

This picture

ℂwrt ℰ22→Ttℂwrt ℰ22id⁡↓id⁡↓ℂwrt B2→Dtℂwrt B2

shows how to get the diagonalization.

(2003)=(1211)−1(4−211)(1211)

Comment. This equation matches the T=PSP−1 definition under this renaming.

T=(2003)P=(1211)−1P−1=(1211)S=(4−211)

Five.II.2.7 §II.2 source ↗

Diagonalize this matrix by following the steps of Example 2.5.

(1100)

(a)

Set up the matrix-vector equation described in Lemma 2.4 and rewrite it as a linear system.

(b)

By considering solutions for that system, find two vectors to make a basis. (Consider separately the system in the x=0 and x≠0 cases. Also, recall that the zero vector cannot be a member of a basis.)

(c)

Use that basis in the similarity diagram to get the diagonal matrix as the product of three others.

Answer (3 parts)

(a)

We want a basis where the member vectors satisfy

(1100)(b1b2)=x⋅(b1b2)

so consider the resulting linear system.

b1+b2=x⋅b10=x⋅b2

(b)

If x=0 then b1+b2=0⁠, so b1=−b2⁠. That is, associated with this case are basis elements where the entries have opposite signs. (The entries must be nonzero since the zero vector cannot be a member of a basis.)

If x≠0 then b2=0 and substituting that into the top equation gives b1=x⋅b1⁠. So in this case, either b1=0 or x=1⁠. But since b2=0 in this case, b1=0 is impossible because the zero vector is not in any basis. Thus x=1 and the system becomes this.

b1+b2=1⋅b10=1⋅b2

So this case is associated with basis elements where b2=0 and b1 can have any nonzero value.

Therefore, a suitable basis is this.

B=⟨(1−1),(10)⟩

(c)

With that basis, to get the diagonalization, draw the similarity diagram.

ℂwrt ℰ22→Ttℂwrt ℰ22id⁡↓id⁡↓ℂwrt B2→Dtℂwrt B2

In this diagram,

T=(1100)

and we will calculate the matrix D using the basis B and traversing the lower-left to upper-left, to upper-right, and to lower-right path.

We first want to go from the lower-left to the upper-left. Observe that because we have chosen the basis for the top of the diagram to be the standard basis ℰ2⁠, the matrix representing that map is easy.

RepB,ℰ2⁡(id⁡)=(11−10)

The calculation for the entire path does indeed give a diagonal result (of course, the lower-left to upper-left portion appears on the right side of the product of the three matrices).

D=(11−10)−1(1100)(11−10)=(0−111)(1100)(11−10)=(0001)

Five.II.2.8 §II.2 source ↗

Follow Example 2.5 to diagonalize this matrix.

(0110)

(a)

Set up the matrix-vector equation and rewrite it as a linear system.

(b)

By considering solutions for that system, find two vectors to make a basis. (Consider separately the x=0 and x≠0 cases. Also, recall that the zero vector cannot be a member of a basis.)

(c)

With that basis use the similarity diagram to get the diagonalization as the product of three matrices.

Answer (3 parts)

(a)

The members of the basis satisfy this.

(0110)(b1b2)=x⋅(b1b2)

so this is the resulting linear system.

b2=x⋅b1b1=x⋅b2⟹−xb1+b2=0b1−xb2=0

(b)

If x=0 then b1=b2=0⁠, but the zero vector is not a member of any basis so this case cannot happen.

The other case is x≠0⁠.

−xb1+b2=0b1−xb2=0⟶(1/x)ρ1+ρ2−xb1+b2=0(−x+(1/x))b2=0

Focus on the bottom equation. Either b2=0 or −x+(1/x)=0⁠. If b2=0 then substituting into the top equation gives b1=0 since x≠0 is the assumption of this case. But the zero vector is not a member of any basis so this cannot happen.

We are left with this.

1x−x=0⟹0=1−x2x=(1−x)(1+x)x

If x=1 then the system is

−1⋅b1+b2=00=0

and so basis vectors satisfy b2=b1⁠. If x=−1 then the system is

b1+b2=00=0

and basis vectors satisfy b2=−b1⁠. Therefore, this is an example of a suitable basis.

B=⟨(11),(1−1)⟩

(c)

To get the diagonalization, draw the similarity diagram

ℂwrt ℰ22→Ttℂwrt ℰ22id⁡↓id⁡↓ℂwrt B2→Dtℂwrt B2

where T is the matrix give in the problem statement.

T=(0110)

Calculate the matrix D using the basis B by moving from the lower left to the upper left, to the upper right, and then to the lower right.

Observe that the matrix representing the lower left to upper left map is easy.

RepB,ℰ2⁡(id⁡)=(111−1)

With that, the calculation for the entire path does give a diagonal result.

D=(111−1)−1(0110)(111−1)=−12(−1−1−11)(0110)(111−1)=(100−1)

Five.II.2.9 §II.2 source ↗

Diagonalize these upper triangular matrices.

(a)

(−2102)

(b)

(5401)

Answer (2 parts)

(a)

Setting up

(−2102)(b1b2)=x⋅(b1b2)⟹(−2−x)⋅b1+b2=0(2−x)⋅b2=0

gives the two possibilities that b2=0 and x=2⁠. Following the b2=0 possibility leads to the first equation (−2−x)b1=0 with the two cases that b1=0 and that x=−2⁠. Thus, under this first possibility, we find x=−2 and the associated vectors whose second component is zero, and whose first component is free.

(−2102)(b10)=−2⋅(b10)β→1=(10)

Following the other possibility leads to a first equation of −4b1+b2=0 and so the vectors associated with this solution have a second component that is four times their first component.

(−2102)(b14b1)=2⋅(b14b1)β→2=(14)

The diagonalization is this.

(1104)(−2102)(1104)−1=(−2002)

(b)

The calculations are like those in the prior part.

(5401)(b1b2)=x⋅(b1b2)⟹(5−x)⋅b1+4⋅b2=0(1−x)⋅b2=0

The bottom equation gives the two possibilities that b2=0 and x=1⁠. Following the b2=0 possibility, and discarding the case where both b2 and b1 are zero, gives that x=5⁠, associated with vectors whose second component is zero and whose first component is free.

β→1=(10)

The x=1 possibility gives a first equation of 4b1+4b2=0 and so the associated vectors have a second component that is the negative of their first component.

β→1=(1−1)

We thus have this diagonalization.

(110−1)(5401)(110−1)−1=(5001)

Five.II.2.10 §II.2 source ↗

If we try to diagonalize the matrix of Example 2.3

N=(0010)

using the method of Example 2.5 then what goes wrong?

(a)

Draw the similarity diagram with N⁠.

(b)

Set up the matrix-vector equation described in Lemma 2.4 and rewrite it as a linear system.

(c)

By considering solutions for that system, find the trouble. (Consider separately the x=0 and x≠0 cases.)

Answer (2 parts)

(a)

In this similarity diagram

ℂwrt ℰ22→Nnℂwrt ℰ22id⁡↓id⁡↓ℂwrt B2→Dnℂwrt B2

we will try to calculate the matrix D and find a suitable basis B⁠, by traversing the lower-left to upper-left, to upper-right, and to lower-right path.

(b)

The diagram shows that where B=⟨β→1,β→2⟩ we want n(β→1)=λ1β→1 and n(β→2)=λ2β→2⁠. To determine λ1 and λ2⁠, where the matrix N represents the map with respect to the standard bases, we want to find solutions of this system.

(0010)(b1b2)=x⋅(b1b2)

The equivalent linear system

0b1+0b2=xb1b1+0b2=xb2⟹−xb1=0b1−xb2=0

By the first line, either x=0 or b1=0⁠. If x=0 then the second equation becomes b1−0=0 so b1=0 and we get this solution set.

{(01)b2|b2∈ℂ}

We can take the initial member of the basis to be β→1=e→2⁠.

If x≠0 then the first equation of the linear system gives b1=0⁠. The system’s first line becomes 0=0 and its second line becomes −xb2=0⁠. Since this case takes x≠0⁠, we get b2=0 and the set of solutions is trivial.

{(00)}

So what goes wrong is that we don’t get the basis required by Lemma 2.4.

Five.II.2.11 §II.2 source ↗

What form do the powers of a diagonal matrix have?

Answer

For any integer p⁠, we have this.

(d100⋱dn)p=(d1p00⋱dnp)

Five.II.2.12 §II.2 source ↗

Give two same-sized diagonal matrices that are not similar. Must any two different diagonal matrices come from different similarity classes?

Answer

These two are not similar

(0000)(1001)

because each is alone in its similarity class.

For the second half, these

(2003)(3002)

are similar via the matrix that changes bases from ⟨β→1,β→2⟩ to ⟨β→2,β→1⟩⁠. (Question. Are two diagonal matrices similar if and only if their diagonal entries are permutations of each others?)

Five.II.2.13 §II.2 source ↗

Give a nonsingular diagonal matrix. Can a diagonal matrix ever be singular?

Answer

Contrast these two.

(2001)(2000)

The first is nonsingular, the second is singular.

Five.II.2.14 §II.2 source ↗

Show that the inverse of a diagonal matrix is the diagonal of the inverses, if no element on that diagonal is zero. What happens when a diagonal entry is zero?

Answer

To check that the inverse of a diagonal matrix is the diagonal matrix of the inverses, just multiply.

(a1,100a2,2⋱an,n)(1/a1,1001/a2,2⋱1/an,n)

(Showing that it is a left inverse is just as easy.)

If a diagonal entry is zero then the diagonal matrix is singular; it has a zero determinant.

Five.II.2.15 §II.2 source ↗

The equation ending Example 2.5

(110−1)−1(3201)(110−1)=(3001)

is a bit jarring because for P we must take the first matrix, which is shown as an inverse, and for P−1 we take the inverse of the first matrix, so that the two −1 powers cancel and this matrix is shown without a superscript −1⁠.

(a)

Check that this nicer-appearing equation holds.

(3001)=(110−1)(3201)(110−1)−1

(b)

Is the previous item a coincidence? Or can we always switch the P and the P−1⁠?

Answer (2 parts)

(a)

The check is easy.

(110−1)(3201)=(330−1)(330−1)(110−1)−1=(3001)

(b)

It is a coincidence, in the sense that if T=PSP−1 then T need not equal P−1SP⁠. Even in the case of a diagonal matrix D⁠, the condition that D=PTP−1 does not imply that D equals P−1TP⁠. The matrices from Example 2.2 show this.

(1211)(4−211)=(605−1)(605−1)(1211)−1=(−612−611)

Five.II.2.16 §II.2 source ↗

Show that the P used to diagonalize in Example 2.5 is not unique.

Answer

The columns of the matrix are the vectors associated with the x⁠’s. The exact choice, and the order of the choice was arbitrary. We could, for instance, get a different matrix by swapping the two columns.

Five.II.2.17 §II.2 source ↗

Find a formula for the powers of this matrix. Hint: see Exercise 2.11.

(−31−42)

Answer

Diagonalizing and then taking powers of the diagonal matrix shows that

(−31−42)k=13(−11−44)+(−23)k(4−14−1).

Five.II.2.18 §II.2 source ↗

We can ask how diagonalization interacts with the matrix operations. Assume that t,s:V→V are each diagonalizable. Is ct diagonalizable for all scalars c⁠? What about t+s⁠? t∘s⁠?

Answer

Yes, ct is diagonalizable by the final theorem of this subsection.

No, t+s need not be diagonalizable. Intuitively, the problem arises when the two maps diagonalize with respect to different bases (that is, when they are not simultaneously diagonalizable). Specifically, these two are diagonalizable but their sum is not:

(1100)(−1000)

(the second is already diagonal; for the first, see Exercise 2.7). The sum is not diagonalizable because its square is the zero matrix.

The same intuition suggests that t∘s is not be diagonalizable. These two are diagonalizable but their product is not:

(1000)(0110)

(for the second, see Exercise 2.8).

Five.II.2.19 §II.2 source ↗

Show that matrices of this form are not diagonalizable.

(1c01)c≠0

Answer

If

P(1c01)P−1=(a00b)

then

P(1c01)=(a00b)P

so

(pqrs)(1c01)=(a00b)(pqrs)(pcp+qrcr+s)=(apaqbrbs)

The 1,1 entries show that a=1 and the 1,2 entries then show that pc=0⁠. Since c≠0 this means that p=0⁠. The 2,1 entries show that b=1 and the 2,2 entries then show that rc=0⁠. Since c≠0 this means that r=0⁠. But if both p and r are 0 then P is not invertible.

Five.II.2.20 §II.2 source ↗

Show that each of these is diagonalizable.

(a)

(1221)

(b)

(xyyz)x,y,z scalars

Answer (2 parts)

(a)

Using the formula for the inverse of a 2×2 matrix gives this.

(abcd)(1221)⋅1ad−bc⋅(d−b−ca)=1ad−bc(ad+2bd−2ac−bc−ab−2b2+2a2+abcd+2d2−2c2−cd−bc−2bd+2ac+ad)

Now pick scalars a,…,d so that ad−bc≠0 and 2d2−2c2=0 and 2a2−2b2=0⁠. For example, these will do.

(111−1)(1221)⋅1−2⋅(−1−1−11)=1−2(−6002)

(b)

As above,

(abcd)(xyyz)⋅1ad−bc⋅(d−b−ca)=1ad−bc(adx+bdy−acy−bcz−abx−b2y+a2y+abzcdx+d2y−c2y−cdz−bcx−bdy+acy+adz)

we are looking for scalars a,…,d so that ad−bc≠0 and −abx−b2y+a2y+abz=0 and cdx+d2y−c2y−cdz=0⁠, no matter what values x⁠, y⁠, and z have.

For starters, we assume that y≠0⁠, else the given matrix is already diagonal. We shall use that assumption because if we (arbitrarily) let a=1 then we get

−bx−b2y+y+bz=0(−y)b2+(z−x)b+y=0

and the quadratic formula gives

b=−(z−x)±(z−x)2−4(−y)(y)−2yy≠0

(note that if x⁠, y⁠, and z are real then these two b⁠’s are real as the discriminant is positive). By the same token, if we (arbitrarily) let c=1 then

dx+d2y−y−dz=0(y)d2+(x−z)d−y=0

and we get here

d=−(x−z)±(x−z)2−4(y)(−y)2yy≠0

(as above, if x,y,z∈ℝ then this discriminant is positive so a symmetric, real, 2×2 matrix is similar to a real diagonal matrix).

For a check we try x=1⁠, y=2⁠, z=1⁠.

b=0±0+16−4=∓1d=0±0+164=±1

Note that not all four choices (b,d)=(+1,+1),…,(−1,−1) satisfy ad−bc≠0⁠.

Five.II.3.22 §II.3 source ↗

This matrix has two eigenvalues λ1=3⁠, λ2=−4⁠.

(41−8−5)

Give two different diagonal form matrices with which it is similar.

Answer

We can permute the columns, by permuting the basis with which the representation is done.

(300−4)(−4003)

Five.II.3.23 §II.3 source ↗

For each, find the characteristic polynomial and the eigenvalues.

(a)

(10−94−2)

(b)

(1243)

(c)

(0370)

(d)

(0000)

(e)

(1001)

Answer (5 parts)

(a)

This

0=|10−x−94−2−x|=(10−x)(−2−x)−(−36)

simplifies to the characteristic equation x2−8x+16=0⁠. Because the equation factors into (x−4)2 there is only one eigenvalue λ1=4⁠.

(b)

0=(1−x)(3−x)−8=x2−4x−5⁠; λ1=5⁠, λ2=−1

(c)

x2−21=0⁠; λ1=21⁠, λ2=−21

(d)

x2=0⁠; λ1=0

(e)

x2−2x+1=0⁠; λ1=1

Five.II.3.24 §II.3 source ↗

For each matrix, find the characteristic equation and the eigenvalues and associated eigenvectors.

(a)

(308−1)

(b)

(32−10)

Answer (2 parts)

(a)

The characteristic equation is (3−x)(−1−x)=0⁠. Its roots, the eigenvalues, are λ1=3 and λ2=−1⁠. For the eigenvectors we consider this equation.

(3−x08−1−x)(b1b2)=(00)

For the eigenvector associated with λ1=3⁠, we consider the resulting linear system.

0⋅b1+0⋅b2=08⋅b1+−4⋅b2=0

The eigenspace is the set of vectors whose second component is twice the first component.

{(b2/2b2)|b2∈ℂ}(308−1)(b2/2b2)=3⋅(b2/2b2)

(Here, the parameter is b2 only because that is the variable that is free in the above system.) Hence, this is an eigenvector associated with the eigenvalue 3⁠.

(12)

Finding an eigenvector associated with λ2=−1 is similar. This system

4⋅b1+0⋅b2=08⋅b1+0⋅b2=0

leads to the set of vectors whose first component is zero.

{(0b2)|b2∈ℂ}(308−1)(0b2)=−1⋅(0b2)

And so this is an eigenvector associated with λ2⁠.

(01)

(b)

The characteristic equation is

0=|3−x2−1−x|=x2−3x+2=(x−2)(x−1)

and so the eigenvalues are λ1=2 and λ2=1⁠. To find eigenvectors, consider this system.

(3−x)⋅b1+2⋅b2=0−1⋅b1−x⋅b2=0

For λ1=2 we get

1⋅b1+2⋅b2=0−1⋅b1−2⋅b2=0

leading to this eigenspace and eigenvector.

{(−2b2b2)|b2∈ℂ}(−21)

For λ2=1 the system is

2⋅b1+2⋅b2=0−1⋅b1−1⋅b2=0

leading to this.

{(−b2b2)|b2∈ℂ}(−11)

Five.II.3.25 §II.3 source ↗

Find the characteristic equation, and the eigenvalues and associated eigenvectors for this matrix. Hint. The eigenvalues are complex.

(−2−152)

Answer

The characteristic equation

0=|−2−x−152−x|=x2+1

has the complex roots λ1=i and λ2=−i⁠. This system

(−2−x)⋅b1−1⋅b2=05⋅b1(2−x)⋅b2=0

For λ1=i Gauss’s Method gives this reduction.

(−2−i)⋅b1−1⋅b2=05⋅b1−(2−i)⋅b2=0⟶(−5/(−2−i))ρ1+ρ2(−2−i)⋅b1−1⋅b2=00=0

(For the calculation in the lower right get a common denominator

5−2−i−(2−i)=5−2−i−−2−i−2−i⋅(2−i)=5−(−5)−2−i

to see that it gives a 0=0 equation.) These are the resulting eigenspace and eigenvector.

{((1/(−2−i))b2b2)|b2∈ℂ}(1/(−2−i)1)

For λ2=−i the system

(−2+i)⋅b1−1⋅b2=05⋅b1−(2+i)⋅b2=0⟶(−5/(−2+i))ρ1+ρ2(−2+i)⋅b1−1⋅b2=00=0

leads to this.

{((1/(−2+i))b2b2)|b2∈ℂ}(1/(−2+i)1)

Five.II.3.26 §II.3 source ↗

Find the characteristic polynomial, the eigenvalues, and the associated eigenvectors of this matrix.

(111001001)

Answer

The characteristic equation is

0=|1−x110−x1001−x|=(1−x)2(−x)

and so the eigenvalues are λ1=1 (this is a repeated root of the equation) and λ2=0⁠. For the rest, consider this system.

(1−x)⋅b1+b2+b3=0−x⋅b2+b3=0(1−x)⋅b3=0

When x=λ1=1 then the solution set is this eigenspace.

{(b100)|b1∈ℂ}

When x=λ2=0 then the solution set is this eigenspace.

{(−b2b20)|b2∈ℂ}

So these are eigenvectors associated with λ1=1 and λ2=0⁠.

(100)(−110)

Five.II.3.27 §II.3 source ↗

For each matrix, find the characteristic equation, and the eigenvalues and associated eigenvectors.

(a)

(3−20−230005)

(b)

(0100014−178)

Answer (2 parts)

(a)

The characteristic equation is

0=|3−x−20−23−x0005−x|=x3−11x2+35x−25=(x−1)(x−5)2

and so the eigenvalues are λ1=1 and also the repeated eigenvalue λ2=5⁠. To find eigenvectors, consider this system.

(3−x)⋅b1−2⋅b2=0−2⋅b1+(3−x)⋅b2=0(5−x)⋅b3=0

For λ1=1 we get

2⋅b1−2⋅b2=0−2⋅b1+2⋅b2=04⋅b3=0

leading to this eigenspace and eigenvector.

{(b2b20)|b2∈ℂ}(110)

For λ2=5 the system is

−2⋅b1−2⋅b2=0−2⋅b1−2⋅b2=00⋅b3=0

leading to this.

{(−b2b20)+(00b3)|b2,b3∈ℂ}(−110),(001)

(b)

The characteristic equation is

0=|−x100−x14−178−x|=−x3+8x2−17x+4=−1⋅(x−4)(x2−4x+1)

and the eigenvalues are λ1=4 and (by using the quadratic equation) λ2=2+3 and λ3=2−3⁠. To find eigenvectors, consider this system.

−x⋅b1+b2=0−x⋅b2+b3=04⋅b1−17⋅b2+(8−x)⋅b3=0

Substituting x=λ1=4 gives the system

−4⋅b1+b2=0−4⋅b2+b3=04⋅b1−17⋅b2+4⋅b3=0⟶ρ1+ρ3−4⋅b1+b2=0−4⋅b2+b3=0−16⋅b2+4⋅b3=0⟶−4ρ2+ρ3−4⋅b1+b2=0−4⋅b2+b3=00=0

leading to this eigenspace and eigenvector.

V4={((1/16)⋅b3(1/4)⋅b3b3)|b2∈ℂ}(1416)

Substituting x=λ2=2+3 gives the system

(−2−3)⋅b1+b2=0(−2−3)⋅b2+b3=04⋅b1−17⋅b2+(6−3)⋅b3=0⟶(−4/(−2−3))ρ1+ρ3(−2−3)⋅b1+b2=0(−2−3)⋅b2+b3=0+(−9−43)⋅b2+(6−3)⋅b3=0

(the middle coefficient in the third equation equals the number (−4/(−2−3))−17⁠; find a common denominator of −2−3 and then rationalize the denominator by multiplying the top and bottom of the fraction by −2+3⁠)

⟶((9+43)/(−2−3))ρ2+ρ3(−2−3)⋅b1+b2=0(−2−3)⋅b2+b3=00=0

which leads to this eigenspace and eigenvector.

V2+3={((1/(2+3)2)⋅b3(1/(2+3))⋅b3b3)|b3∈ℂ}((1/(2+3)2)(1/(2+3))1)

Finally, substituting x=λ3=2−3 gives the system

(−2+3)⋅b1+b2=0(−2+3)⋅b2+b3=04⋅b1−17⋅b2+(6+3)⋅b3=0⟶(−4/(−2+3))ρ1+ρ3(−2+3)⋅b1+b2=0(−2+3)⋅b2+b3=0(−9+43)⋅b2+(6+3)⋅b3=0⟶((9−43)/(−2+3))ρ2+ρ3(−2+3)⋅b1+b2=0(−2+3)⋅b2+b3=00=0

which gives this eigenspace and eigenvector.

V2−3={((1/(−2+3)2)⋅b3(1/(2−3))⋅b3b3)|b3∈ℂ}((1/(−2+3)2)(1/(2−3))1)

Five.II.3.28 §II.3 source ↗

For each matrix, find the characteristic polynomial, and the eigenvalues and associated eigenspaces. Also find the algebraic and geometric multiplicities.

(a)

(13−4−47)

(b)

(13−3−37−3−66−2)

(c)

(23−302−3001)

Answer (3 parts)

(a)

The characteristic polynomial factors as x2−20x+75=(x−5)(x−15)⁠, so the eigenvalues are λ1=5 and λ2=15⁠. These are the associated eigenspaces.

V5={k(12)|k∈ℂ}V15={k(−21)|k∈ℂ}

For each eigenvalue, both the algebraic and geometric multiplicities are 1⁠.

(b)

The characteristic polynomial x3−6x2+32 factors into (x+2)(x−4)2⁠. The eigenvectors are λ1=−2 and λ2=4⁠. Here are the associated eigenspaces.

V−2={k(112)|k∈ℂ}V4={k1(110)+k2(10−1)|k1,k2∈ℂ}

For λ1=−2 the algebraic and geometric multiplicities are both 1⁠. For λ2=4 the algebraic and geometric multiplicities are both 2⁠.

(c)

The characteristic polynomial x3−5x2+8x−4 factors into (x−1)(x−2)2⁠. The eigenvalues are λ1=1 and λ2=2⁠. Here are the associated eigenspaces.

V1={k(−631)|k∈ℂ}V2={k(100)|k∈ℂ}

For λ1=1 the algebraic and geometric multiplicities are both 1⁠. For λ2=2 the algebraic multiplicity is 2 but the geometric multiplicity is 1⁠.

Five.II.3.29 §II.3 source ↗

Let t:𝒫2→𝒫2 be this linear map.

a0+a1x+a2x2↦(5a0+6a1+2a2)−(a1+8a2)x+(a0−2a2)x2

Find its eigenvalues and the associated eigenvectors.

Answer

With respect to the natural basis B=⟨1,x,x2⟩ the matrix representation is this.

RepB,B⁡(t)=(5620−1−810−2)

Thus the characteristic equation

0=(5−x620−1−x−810−2−x)=(5−x)(−1−x)(−2−x)−48−2⋅(−1−x)

is 0= − x3+ 2x2+ 15x− 36= − 1⋅(x+4)(x−3)2⁠. To find the associated eigenvectors, consider this system.

(5−x)⋅b1+6b2+2b3=0(−1−x)⋅b2−8b3=0b1+(−2−x)⋅b3=0

Plugging in λ1=−4 for x gives

9b1+6b2+2b3=03b2−8b3=0b1+2b3=0⟶−(1/9)ρ1+ρ3⟶(2/9)ρ2+ρ39b1+6b2+2b3=03b2−8b3=0

Here is the eigenspace and an eigenvector.

V−4={(2⋅b3(8/3)⋅b3b3)|b3∈ℂ}(28/31)

Similarly, plugging in x=λ2=3 gives

2b1+6⋅b2+2⋅b3=0−4⋅b2−8⋅b3=0b1−5⋅b3=0⟶−(1/2)ρ1+ρ3⟶−(3/4)ρ2+ρ32b1+6⋅b2+2⋅b3=0−4⋅b2−8⋅b3=0

with this eigenspace and eigenvector.

V3={(5⋅b3−2⋅b3b3)|b3∈ℂ}(5−21)

Five.II.3.30 §II.3 source ↗

Find the eigenvalues and eigenvectors of this map t:ℳ2→ℳ2⁠.

(abcd)↦(2ca+cb−2cd)

Answer

λ=1,(0001) and (2310)⁠, λ=−2,(−1010)⁠, λ=−1,(−2110)

Five.II.3.31 §II.3 source ↗

Find the eigenvalues and associated eigenvectors of the differentiation operator d/dx:𝒫3→𝒫3⁠.

Answer

Fix the natural basis B=⟨1,x,x2,x3⟩⁠. The map’s action is 1↦0⁠, x↦1⁠, x2↦2x⁠, and x3↦3x2 and its representation is easy to compute.

T=RepB,B⁡(d/dx)=(0100002000030000)B,B

We find the eigenvalues with this computation.

0=|T−xI|=|−x1000−x2000−x3000−x|=x4

Thus the map has the single eigenvalue λ=0⁠. To find the associated eigenvectors, we solve

(0100002000030000)B,B(b1b2b3b4)B=0⋅(b1b2b3b4)B⟹b2=0, b3=0, b4=0

to get this eigenspace.

{(b1000)B|b1∈ℂ}={b1+0⋅x+0⋅x2+0⋅x3|b1∈ℂ}={b1|b1∈ℂ}

Five.II.3.32 §II.3 source ↗

Prove that the eigenvalues of a triangular matrix (upper or lower triangular) are the entries on the diagonal.

Answer

The determinant of the triangular matrix T−xI is the product down the diagonal, and so it factors into the product of the terms ti,i−x⁠.

Five.II.3.33 §II.3 source ↗

This matrix has distinct eigenvalues.

(1216−10−1−2−1)

(a)

Diagonalize it.

(b)

Find a basis with respect to which this matrix has that diagonal representation.

(c)

Draw the diagram. Find the matrices P and P−1 to effect the change of basis.

Answer (3 parts)

(a)

Computing the eigenvalues gives λ1=3⁠, λ2=0⁠, and λ3=−4⁠. Thus the diagonalization of the matrix is this.

(30000000−4)

(b)

Get an eigenvector associated with each eigenvalue, and use them to make a basis.

For λ1=3 we solve this system.

(1−3216−1−30−1−2−1−3)(xyz)=(000)

Gauss’s Method suffices.

(−22106−400−1−2−40)⟶−(1/2)ρ1+ρ33ρ1+ρ2⟶(3/2)ρ1+ρ3(−221002300000)

Parametrizing gives the solution space.

S1={(−1−3/21)⋅z|z∈ℂ}

This is the λ2=0 system.

(1216−10−1−2−1)(xyz)=(000)

Gauss’s Method

(12106−100−1−2−10)⟶ρ1+ρ3−6ρ1+ρ2(12100−13−600000)

followed by parametrizing gives the solution space.

S2={(−1/13−6/131)⋅z|z∈ℂ}

The λ3=−4 system works the same way.

(1−(−4)216−1−(−4)0−1−2−1−(−4))(xyz)=(000)

Gauss’s Method

(52106300−1−230)⟶(1/5)ρ1+ρ3−(6/5)ρ1+ρ2⟶(8/3)ρ2+ρ3(521003/5−6/500000)

and parametrizing gives this.

S3={(−121)⋅z|z∈ℂ}

Thus the matrix will have that diagonal form when it represents of the same transformation but with respect to D,D for this D⁠.

D=⟨(−1−3/21),(−1/13−6/131),(−121)⟩

(c)

Call the matrix T⁠. Taking the domain and codomain spaces as ℝ3⁠, and taking the matrix to be T=Repℰ3,ℰ3⁡(t) gives this diagram.

ℂwrt ℰ33→Ttℂwrt ℰ33id⁡↓id⁡↓ℂwrt B3→Dtℂwrt B3

Reading off the diagram we have D=P−1MP where P=Repℰ3,B⁡(id⁡)⁠. To represent each vector in B with respect to E3 is easy because with respect to the standard basis each vector represents itself.

Repℰ3⁡((−1−3/21))=(−1−3/21)Repℰ3⁡((−1/13−6/131))=(−1/13−6/131)Repℰ3⁡((−121))=(−121)

Thus the matrix just concatenates the three basis vectors.

P=(−1−1/13−1−3/2−6/132111)and soP−1(−16/21−2/7−4/2113/12013/12−9/282/73/28)

Five.II.3.34 §II.3 source ↗

Find the formula for the characteristic polynomial of a 2×2 matrix.

Answer

Just expand the determinant of T−xI⁠.

|a−xcbd−x|=(a−x)(d−x)−bc=x2+(−a−d)⋅x+(ad−bc)

Five.II.3.35 §II.3 source ↗

Prove that the characteristic polynomial of a transformation is well-defined.

Answer

Any two representations of that transformation are similar, and similar matrices have the same characteristic polynomial.

Five.II.3.36 §II.3 source ↗

Prove or disprove: if all the eigenvalues of a matrix are 0 then it must be the zero matrix.

Answer

It is not true. All of the eigenvalues of this matrix are 0⁠.

(0100)

Five.II.3.37 §II.3 source ↗

(a)

Show that any non-⁠0→ vector in any nontrivial vector space can be an eigenvector. That is, given a v→≠0→ from a nontrivial V⁠, show that there is a transformation t:V→V having a scalar eigenvalue λ∈ℝ such that v→∈Vλ⁠.

(b)

What if we are given a scalar λ⁠? Can any non-⁠0→ member of any nontrivial vector space be an eigenvector associated with λ⁠?

Answer (2 parts)

(a)

Use λ=1 and the identity map.

(b)

Yes, use the transformation that multiplies all vectors by the scalar λ⁠.

Five.II.3.38 §II.3 source ↗

Suppose that t:V→V and T=RepB,B⁡(t)⁠. Prove that the eigenvectors of T associated with λ are the non-⁠0→ vectors in the kernel of the map represented (with respect to the same bases) by T−λI⁠.

Answer

If t(v→)=λ⋅v→ then v→↦0→ under the map t−λ⋅id⁡⁠.

Five.II.3.39 §II.3 source ↗

Prove that if a,…,d are all integers and a+b=c+d then

(abcd)

has integral eigenvalues, namely a+b and a−c⁠.

Answer

The characteristic equation

0=|a−xbcd−x|=(a−x)(d−x)−bc

simplifies to x2+(−a−d)⋅x+(ad−bc)⁠. Checking that the values x=a+b and x=a−c satisfy the equation (under the a+b=c+d condition) is routine.

Five.II.3.40 §II.3 source ↗

Prove that if T is nonsingular and has eigenvalues λ1,…,λn then T−1 has eigenvalues 1/λ1,…,1/λn⁠. Is the converse true?

Answer

Consider an eigenspace Vλ⁠. Any w→∈Vλ is the image w→=λ⋅v→ of some v→∈Vλ (namely, v→=(1/λ)⋅w→⁠). Thus, on Vλ (which is a nontrivial subspace) the action of t−1 is t−1(w→)=v→=(1/λ)⋅w→⁠, and so 1/λ is an eigenvalue of t−1⁠.

Five.II.3.41 §II.3 source ↗

Suppose that T is n×n and c,d are scalars.

(a)

Prove that if T has the eigenvalue λ with an associated eigenvector v→ then v→ is an eigenvector of cT+dI associated with eigenvalue cλ+d⁠.

(b)

Prove that if T is diagonalizable then so is cT+dI⁠.

Answer (2 parts)

(a)

We have (cT+dI)v→= cTv→+ dIv→= cλv→+ dv→= (cλ+d)⋅v→⁠.

(b)

Suppose that S=PTP−1 is diagonal. Then P(cT+dI)P−1= P(cT)P−1+ P(dI)P−1= cPTP−1+ dI= cS+ dI is also diagonal.

Five.II.3.42 §II.3 source ↗

Show that λ is an eigenvalue of T if and only if the map represented by T−λI is not an isomorphism.

Answer

The scalar λ is an eigenvalue if and only if the transformation t−λid⁡ is singular. A transformation is singular if and only if it is not an isomorphism (that is, a transformation is an isomorphism if and only if it is nonsingular).

Five.II.3.43 §II.3 source ↗

[Strang 80]

(a)

Show that if λ is an eigenvalue of A then λk is an eigenvalue of Ak⁠.

(b)

What is wrong with this proof generalizing that? “If λ is an eigenvalue of A and μ is an eigenvalue for B⁠, then λμ is an eigenvalue for AB⁠, for, if Ax→=λx→ and Bx→=μx→ then ABx→=Aμx→=μAx→=μλx→⁠”?

Answer (2 parts)

(a)

Where the eigenvalue λ is associated with the eigenvector x→ then Akx→= A⋯Ax→= Ak−1λx→= λAk−1x→= ⋯= λkx→⁠. (The full details require induction on k⁠.)

(b)

The eigenvector associated with λ might not be an eigenvector associated with μ⁠.

Five.II.3.44 §II.3 source ↗

Do matrix equivalent matrices have the same eigenvalues?

Answer

No. These are two same-sized, equal rank, matrices with different eigenvalues.

(1001)(1002)

Five.II.3.45 §II.3 source ↗

Show that a square matrix with real entries and an odd number of rows has at least one real eigenvalue.

Answer

The characteristic polynomial has an odd power and so has at least one real root.

Five.II.3.46 §II.3 source ↗

Diagonalize.

(−122222−3−6−6)

Answer

The characteristic polynomial x3+5x2+6x has distinct roots λ1=0⁠, λ2=−2⁠, and λ3=−3⁠. Thus the matrix can be diagonalized into this form.

(0000−2000−3)

Five.II.3.47 §II.3 source ↗

Suppose that P is a nonsingular n×n matrix. Show that the similarity transformation map tP:ℳn×n→ℳn×n sending T↦PTP−1 is an isomorphism.

Answer

We must show that it is one-to-one and onto, and that it respects the operations of matrix addition and scalar multiplication.

To show that it is one-to-one, suppose that tP(T)=tP(S)⁠, that is, suppose that PTP−1=PSP−1⁠, and note that multiplying both sides on the left by P−1 and on the right by P gives that T=S⁠. To show that it is onto, consider S∈ℳn×n and observe that S=tP(P−1SP)⁠.

The map tP preserves matrix addition since tP(T+S)= P(T+S)P−1= (PT+PS)P−1= PTP−1+ PSP−1= tP(T+S) follows from properties of matrix multiplication and addition that we have seen. Scalar multiplication is similar: tP(cT)= P(c⋅T)P−1= c⋅(PTP−1)= c⋅tP(T)⁠.

Five.II.3.48 §II.3 source ↗

[Math. Mag., Nov. 1967] Show that if A is an n square matrix and each row (column) sums to c then c is a characteristic root of A⁠. (“Characteristic root” is a synonym for eigenvalue.)

Answer

This is how the answer was given in the cited source. If the argument of the characteristic function of A is set equal to c⁠, adding the first (n−1) rows (columns) to the n⁠th row (column) yields a determinant whose n⁠th row (column) is zero. Thus c is a characteristic root of A⁠.

Five.III.1.9 §III.1 source ↗

Give the chains of range spaces and null spaces for the zero and identity transformations.

Answer

For the zero transformation, no matter what the space, the chain of range spaces is V⊃{0→}={0→}=⋯ and the chain of null spaces is {0→}⊂V=V=⋯⁠. For the identity transformation the chains are V=V=V=⋯ and {0→}={0→}=⋯⁠.

Five.III.1.10 §III.1 source ↗

For each map, give the chain of range spaces and the chain of null spaces, and the generalized range space and the generalized null space.

(a)

t0:𝒫2→𝒫2⁠, a+bx+cx2↦b+cx2

(b)

t1:ℝ2→ℝ2⁠,

(ab)↦(0a)

(c)

t2:𝒫2→𝒫2⁠, a+bx+cx2↦b+cx+ax2

(d)

t3:ℝ3→ℝ3⁠,

(abc)↦(aab)

Answer (4 parts)

(a)

Iterating t0 twice a+bx+cx2↦b+cx2↦cx2 gives

a+bx+cx2⟼t02cx2

and any higher iterate is the same map. Thus, while ℛ(t0) is the space of quadratic polynomials with no linear term {p+rx2|p,r∈ℂ}⁠, and ℛ(t02) is the space of purely-quadratic polynomials {rx2|r∈ℂ}⁠, this is where the chain stabilizes ℛ∞(t0)={rx2|r∈ℂ}⁠.

As for null spaces, 𝒩(t0) is the space of purely-linear quadratic polynomials {qx|q∈ℂ}⁠, and 𝒩(t02) is the space of linear polynomials {p+qx|p,q∈ℂ}⁠. This is the end: 𝒩∞(t0)=𝒩(t02)⁠.

(b)

The second power

(ab)⟼t1(0a)⟼t1(00)

is the zero map. Consequently, the chain of range spaces

ℝ2⊃{(0p)|p∈ℂ}⊃{0→}=⋯

and the chain of null spaces

{0→}⊂{(q0)|q∈ℂ}⊂ℝ2=⋯

each has length two. The generalized range space is the trivial subspace and the generalized null space is the entire space.

(c)

Iterates of this map cycle around

a+bx+cx2⟼t2b+cx+ax2⟼t2c+ax+bx2⟼t2a+bx+cx2⋯

and the chains of range spaces and null spaces are trivial.

𝒫2=𝒫2=⋯{0→}={0→}=⋯

Thus, obviously, generalized spaces are ℛ∞(t2)=𝒫2 and 𝒩∞(t2)={0→}⁠.

(d)

We have

(abc)↦(aab)↦(aaa)↦(aaa)↦⋯

and so the chain of range spaces

ℝ3⊃{(ppr)|p,r∈ℂ}⊃{(ppp)|p∈ℂ}=⋯

and the chain of null spaces

{0→}⊂{(00r)|r∈ℂ}⊂{(0qr)|q,r∈ℂ}=⋯

each has length two. The generalized spaces are the final ones shown above in each chain.

Five.III.1.11 §III.1 source ↗

Prove that function composition is associative (t∘t)∘t=t∘(t∘t) and so we can write t3 without specifying a grouping.

Answer

Each maps x↦t(t(t(x)))⁠.

Five.III.1.12 §III.1 source ↗

Check that a subspace must be of dimension less than or equal to the dimension of its superspace. Check that if the subspace is proper (the subspace does not equal the superspace) then the dimension is strictly less. (This is used in the proof of Lemma 1.4.)

Answer

Recall that if W is a subspace of V then we can enlarge any basis BW for W to make a basis BV for V⁠. From this the first sentence is immediate. The second sentence is also not hard: W is the span of BW and if W is a proper subspace then V is not the span of BW⁠, and so BV must have at least one vector more than does BW⁠.

Five.III.1.13 §III.1 source ↗

Prove that the generalized range space ℛ∞(t) is the entire space, and the generalized null space 𝒩∞(t) is trivial, if the transformation t is nonsingular. Is this ‘only if’ also?

Answer

It is both ‘if’ and ‘only if’. A linear map is nonsingular if and only if it preserves dimension, that is, if the dimension of its range equals the dimension of its domain. With a transformation t:V→V that means that the map is nonsingular if and only if it is onto: ℛ(t)=V (and thus ℛ(t2)=V⁠, etc).

Five.III.1.14 §III.1 source ↗

Verify the null space half of Lemma 1.4.

Answer

The null spaces form chains because because if v→∈𝒩(tj) then tj(v→)=0→ and tj+1(v→)= t(tj(v→))= t(0→)= 0→ and so v→∈𝒩(tj+1)⁠.

Now, the “further” property for null spaces follows from that fact that it holds for range spaces, along with the prior exercise. Because the dimension of ℛ(tj) plus the dimension of 𝒩(tj) equals the dimension n of the starting space V⁠, when the dimensions of the range spaces stop decreasing, so do the dimensions of the null spaces. The prior exercise shows that from this point k on, the containments in the chain are not proper — the null spaces are equal.

Five.III.1.15 §III.1 source ↗

Give an example of a transformation on a three dimensional space whose range has dimension two. What is its null space? Iterate your example until the range space and null space stabilize.

Answer

(Many examples are correct but here is one.) An example is the shift operator on triples of reals (x,y,z)↦(0,x,y)⁠. The null space is all triples that start with two zeros. The map stabilizes after three iterations.

Five.III.1.16 §III.1 source ↗

Show that the range space and null space of a linear transformation need not be disjoint. Are they ever disjoint?

Answer

The differentiation operator d/dx:𝒫1→𝒫1 has the same range space as null space. For an example of where they are disjoint — except for the zero vector — consider an identity map, or any nonsingular map.

Five.III.2.20 §III.2 source ↗

What is the index of nilpotency of the right-shift operator, here acting on the space of triples of reals?

(x,y,z)↦(0,x,y)

Answer

Three. It is at least three because ℓ2((1,1,1))=(0,0,1)≠0→⁠. It is at most three because (x, y, z)↦(0, x, y)↦(0,0, x)↦(0,0,0)⁠.

Five.III.2.21 §III.2 source ↗

For each string basis state the index of nilpotency and give the dimension of the range space and null space of each iteration of the nilpotent map.

(a)

β→1↦β→2↦0→β→3↦β→4↦0→

(b)

β→1↦β→2↦β→3↦0→β→4↦0→β→5↦0→β→6↦0→

(c)

β→1↦β→2↦β→3↦0→

Also give the canonical form of the matrix.

Answer (3 parts)

(a)

The domain has dimension four. The map’s action is that any vector in the space c1⋅β→1+ c2⋅β→2+ c3⋅β→3+ c4⋅β→4 goes to c1⋅β→2+ c2⋅0→+ c3⋅β→4+ c4⋅0→= c1⋅β→3+ c3⋅β→4⁠. The first application of the map sends two basis vectors β→2 and β→4 to zero, and therefore the null space has dimension two and the range space has dimension two. With a second application, all four basis vectors go to zero and so the null space of the second power has dimension four while the range space of the second power has dimension zero. Thus the index of nilpotency is two. This is the canonical form.

(0000100000000010)

(b)

The dimension of the domain of this map is six. For the first power the dimension of the null space is four and the dimension of the range space is two. For the second power the dimension of the null space is five and the dimension of the range space is one. Then the third iteration results in a null space of dimension six and a range space of dimension zero. The index of nilpotency is three, and this is the canonical form.

(000000100000010000000000000000000000)

(c)

The dimension of the domain is three, and the index of nilpotency is three. The first power’s null space has dimension one and its range space has dimension two. The second power’s null space has dimension two and its range space has dimension one. Finally, the third power’s null space has dimension three and its range space has dimension zero. Here is the canonical form matrix.

(000100010)

Five.III.2.22 §III.2 source ↗

Decide which of these matrices are nilpotent.

(a)

(−24−12)

(b)

(3113)

(c)

(−321−321−321)

(d)

(11430−1527)

(e)

(45−22−1933−16−1469−34−29)

Answer (7 parts)

By Lemma 1.4 the nullity has grown as large as possible by the n⁠-th iteration where n is the dimension of the domain. Thus, for the 2×2 matrices, we need only check whether the square is the zero matrix. For the 3×3 matrices, we need only check the cube.

(a)

Yes, this matrix is nilpotent because its square is the zero matrix.

(b)

No, the square is not the zero matrix.

(3113)2=(106610)

(c)

Yes, the cube is the zero matrix. In fact, the square is zero.

(d)

No, the third power is not the zero matrix.

(11430−1527)3=(2068630426826438180634)

(e)

Yes, the cube of this matrix is the zero matrix.

Another way to see that the second and fourth matrices are not nilpotent is to note that they are nonsingular.

Five.III.2.23 §III.2 source ↗

Find the canonical form of this matrix.

(0110100111000000000000000)

Answer

The table of calculations

pNp𝒩(Np)
1(0110100111000000000000000){(ru−u−vuv)|r,u,v∈ℂ}
2(0011100000000000000000000){(rs−u−vuv)|r,s,u,v∈ℂ}
2–zero matrix–ℂ5

gives these requirements of the string basis: three basis vectors map directly to zero, one more basis vector maps to zero by a second application, and the final basis vector maps to zero by a third application. Thus, the string basis has this form.

β→1↦β→2↦β→3↦0→β→4↦0→β→5↦0→

From that the canonical form is immediate.

(0000010000010000000000000)

Five.III.2.24 §III.2 source ↗

Consider the matrix from Example 2.19.

(a)

Use the action of the map on the string basis to give the canonical form.

(b)

Find the change of basis matrices that bring the matrix to canonical form.

(c)

Use the answer in the prior item to check the answer in the first item.

Answer (3 parts)

(a)

The canonical form has a 3×3 block and a 2×2 block

(0000010000010000000000010)

corresponding to the length three string and the length two string in the basis.

(b)

Assume that N is the representation of the underlying map with respect to the standard basis. Let B be the basis to which we will change. By the similarity diagram

ℂwrt ℰ22→Nnℂwrt ℰ22id⁡↓Pid⁡↓PCwrt B2→nℂwrt B2

we have that the canonical form matrix is PNP−1 where

P−1=RepB,ℰ5⁡(id⁡)=(1000001010101000011100001)

and P is the inverse of that.

P=Repℰ5,B⁡(id⁡)=(P−1)−1=(10000−111−11−1010010−11−100001)

(c)

The calculation to check this is routine.

Five.III.2.25 §III.2 source ↗

Each of these matrices is nilpotent.

(a)

(1/2−1/21/2−1/2)

(b)

(0000−110−11)

(c)

(−11−11011−11)

Put each in canonical form.

Answer (3 parts)

(a)

The calculation

pNp𝒩(Np)
1(1/2−1/21/2−1/2){(uu)|u∈ℂ}
2–zero matrix–ℂ2

shows that any map represented by the matrix must act on the string basis in this way

β→1↦β→2↦0→

because the null space after one application has dimension one and exactly one basis vector, β→2⁠, maps to zero. Therefore, this representation with respect to ⟨β→1,β→2⟩ is the canonical form.

(0010)

(b)

The calculation here is similar to the prior one.

pNp𝒩(Np)
1(0000−110−11){(uvv)|u,v∈ℂ}
2–zero matrix–ℂ3

The table shows that the string basis is of the form

β→1↦β→2↦0→β→3↦0→

because the null space after one application of the map has dimension two — β→2 and β→3 are both sent to zero — and one more iteration results in the additional vector going to zero.

(c)

The calculation

pNp𝒩(Np)
1(−11−11011−11){(u0−u)|u∈ℂ}
2(101000−10−1){(uv−u)|u,v∈ℂ}
3–zero matrix–ℂ3

shows that any map represented by this basis must act on a string basis in this way.

β→1↦β→2↦β→3↦0→

Therefore, this is the canonical form.

(000100010)

Five.III.2.26 §III.2 source ↗

Describe the effect of left or right multiplication by a matrix that is in the canonical form for nilpotent matrices.

Answer

A couple of examples

(0010)(abcd)=(00ab)(000100010)(abcdefghi)=(000abcdef)

suggest that left multiplication by a block of subdiagonal ones shifts the rows of a matrix downward. Distinct blocks

(0000100000000010)(abcdefghijklmnop)=(0000abcd0000ijkl)

act to shift down distinct parts of the matrix.

Right multiplication does an analogous thing to columns. See Exercise 2.20.

Five.III.2.27 §III.2 source ↗

Is nilpotence invariant under similarity? That is, must a matrix similar to a nilpotent matrix also be nilpotent? If so, with the same index?

Answer

Yes. Generalize the last sentence in Example 2.10. As to the index, that same last sentence shows that the index of the new matrix is less than or equal to the index of Nˆ⁠, and reversing the roles of the two matrices gives inequality in the other direction.

Another answer to this question is to show that a matrix is nilpotent if and only if any associated map is nilpotent, and with the same index. Then, because similar matrices represent the same map, the conclusion follows. This is Exercise 2.33 below.

Five.III.2.28 §III.2 source ↗

Show that the only eigenvalue of a nilpotent matrix is zero.

Answer

Observe that a canonical form nilpotent matrix has only zero eigenvalues; e.g., the determinant of this lower-triangular matrix

(−x001−x001−x)

is (−x)3⁠, the only root of which is zero. But similar matrices have the same eigenvalues and every nilpotent matrix is similar to one in canonical form.

Another way to see this is to observe that a nilpotent matrix sends all vectors to zero after some number of iterations, but that conflicts with an action on an eigenspace v→↦λv→ unless λ is zero.

Five.III.2.29 §III.2 source ↗

Is there a nilpotent transformation of index three on a two-dimensional space?

Answer

No, by Lemma 1.4 for a map on a two-dimensional space, the nullity has grown as large as possible by the second iteration.

Five.III.2.30 §III.2 source ↗

In the proof of Theorem 2.16, why isn’t the proof’s base case that the index of nilpotency is zero?

Answer

The index of nilpotency of a transformation can be zero only when the vector starting the string must be 0→⁠, that is, only when V is a trivial space.

Five.III.2.31 §III.2 source ↗

Let t:V→V be a linear transformation and suppose v→∈V is such that tk(v→)=0→ but tk−1(v→)≠0→⁠. Consider the t⁠-string ⟨v→,t(v→),…,tk−1(v→)⟩⁠.

(a)

Prove that t is a transformation on the span of the set of vectors in the string, that is, prove that t restricted to the span has a range that is a subset of the span. We say that the span is a t⁠-invariant subspace.

(b)

Prove that the restriction is nilpotent.

(c)

Prove that the t⁠-string is linearly independent and so is a basis for its span.

(d)

Represent the restriction map with respect to the t⁠-string basis.

Answer (4 parts)

(a)

Any member w→ of the span is a linear combination w→= c0⋅v→+ c1⋅t(v→)+ ⋯+ ck−1⋅tk−1(v→)⁠. But then, by the linearity of the map, t(w→)= c0⋅t(v→)+ c1⋅t2(v→)+ ⋯+ ck−2⋅tk−1(v→)+ ck−1⋅0→ is also in the span.

(b)

The operation in the prior item, when iterated k times, will result in a linear combination of zeros.

(c)

If v→=0→ then the set is empty and so is linearly independent by definition. Otherwise write c1v→+⋯+ck−1tk−1(v→)=0→ and apply tk−1 to both sides. The right side gives 0→ while the left side gives c1tk−1(v→)⁠; conclude that c1=0⁠. Continue in this way by applying tk−2 to both sides, etc.

(d)

Of course, t acts on the span by acting on this basis as a single, k⁠-long, t⁠-string.

(0000…001000…000100…00001000⋱000010)

Five.III.2.32 §III.2 source ↗

Finish the proof of Theorem 2.16.

Answer

We must check that B∪Cˆ∪{v→1,…,v→j} is linearly independent where B is a t⁠-string basis for ℛ(t)⁠, where Cˆ is a basis for 𝒩(t)⁠, and where t(v→1)=β→1,…,t(v→i)=β→i⁠. Write

0→=c1,−1v→1+c1,0β→1+c1,1t(β→1)+⋯+c1,h1th1(β→→1)+c2,−1v→2+⋯+cj,hithi(βi→)

and apply t⁠.

0→=c1,−1β→1+c1,0t(β→1)+⋯+c1,h1−1th1(β→→1)+c1,h10→+c2,−1β→2+⋯+ci,hi−1thi(βi→)+ci,hi0→

Conclude that the coefficients c1,−1, …, c1,hi−1, c2,−1, …, ci,hi−1 are all zero as B∪Cˆ is a basis. Substitute back into the first displayed equation to conclude that the remaining coefficients are zero also.

Five.III.2.33 §III.2 source ↗

Show that the terms ‘nilpotent transformation’ and ‘nilpotent matrix’, as given in Definition 2.7, fit with each other: a map is nilpotent if and only if it is represented by a nilpotent matrix. (Is it that a transformation is nilpotent if an only if there is a basis such that the map’s representation with respect to that basis is a nilpotent matrix, or that any representation is a nilpotent matrix?)

Answer

For any basis B⁠, a transformation n is nilpotent if and only if N=RepB,B⁡(n) is a nilpotent matrix. This is because only the zero matrix represents the zero map and so nj is the zero map if and only if Nj is the zero matrix.

Five.III.2.34 §III.2 source ↗

Let T be nilpotent of index four. How big can the range space of T3 be?

Answer

It can be of any size greater than or equal to one. To have a transformation that is nilpotent of index four, whose cube has range space of dimension k⁠, take a vector space, a basis for that space, and a transformation that acts on that basis in this way.

β→1↦β→2↦β→3↦β→4↦0→β→5↦β→6↦β→7↦β→8↦0→⋮β→4k−3↦β→4k−2↦β→4k−1↦β→4k↦0→⋮–𝑝𝑜𝑠𝑠𝑖𝑏𝑙𝑦 𝑜𝑡ℎ𝑒𝑟, 𝑠ℎ𝑜𝑟𝑡𝑒𝑟, 𝑠𝑡𝑟𝑖𝑛𝑔𝑠–

So the dimension of the range space of T3 can be as large as desired. The smallest that it can be is one — there must be at least one string or else the map’s index of nilpotency would not be four.

Five.III.2.35 §III.2 source ↗

Recall that similar matrices have the same eigenvalues. Show that the converse does not hold.

Answer

These two have only zero for eigenvalues

(0000)(0010)

but are not similar (they have different canonical representatives, namely, themselves).

Five.III.2.36 §III.2 source ↗

Lemma 2.1 shows that any for any linear transformation t:V→V the restriction t:ℛ∞(t)→ℛ∞(t) is one-to-one. Show that it is also onto, so it is an automorphism. Must it be the identity map?

Answer

It is onto by Lemma 1.4. It need not be the identity: consider this map t:ℝ2→ℝ2⁠.

(xy)⟼t(yx)

For that map ℛ∞(t)=ℝ2⁠, and t is not the identity.

Five.III.2.37 §III.2 source ↗

Prove that a nilpotent matrix is similar to one that is all zeros except for blocks of super-diagonal ones.

Answer

A simple reordering of the string basis will do. For instance, a map that is associated with this string basis

β→1↦β→2↦0→

is represented with respect to B=⟨β→1,β→2⟩ by this matrix

(0010)

but is represented with respect to B=⟨β→2,β→1⟩ in this way.

(0100)

Five.III.2.38 §III.2 source ↗

Prove that if a transformation has the same range space as null space. then the dimension of its domain is even.

Answer

Let t:V→V be the transformation. If rank⁡(t)=nullity⁡(t) then the equation rank⁡(t)+ nullity⁡(t)= dim⁡(V) shows that dim⁡(V) is even.

Five.III.2.39 §III.2 source ↗

Prove that if two nilpotent matrices commute then their product and sum are also nilpotent.

Answer

For the matrices to be nilpotent they must be square. For them to commute they must be the same size. Thus their product and sum are defined.

Call the matrices A and B⁠. To see that AB is nilpotent, multiply (AB)2=ABAB=AABB=A2B2⁠, and (AB)3=A3B3⁠, etc., and, as A is nilpotent, that product is eventually zero.

The sum is similar; use the Binomial Theorem.

Five.IV.1.13 §IV.1 source ↗

What are the possible minimal polynomials if a matrix has the given characteristic polynomial?

(a)

(x−3)4

(b)

(x+1)3(x−4)

(c)

(x−2)2(x−5)2

(d)

(x+3)2(x−1)(x−2)2

What is the degree of each possibility?

Answer (5 parts)

The Cayley-Hamilton Theorem Theorem 1.8 says that the minimal polynomial must contain the same linear factors as the characteristic polynomial, although possibly of lower degree but not of zero degree.

(a)

The possibilities are m1(x)=x−3⁠, m2(x)=(x−3)2⁠, m3(x)=(x−3)3⁠, and m4(x)=(x−3)4⁠. The first is a degree one polynomial, the second is degree two, the third is degree three, and the fourth is degree four.

(b)

The possibilities are m1(x)=(x+1)(x−4)⁠, m2(x)=(x+1)2(x−4)⁠, and m3(x)=(x+1)3(x−4)⁠. The first is a quadratic polynomial, that is, it has degree two. The second has degree three, and the third has degree four.

(c)

We have m1(x)=(x−2)(x−5)⁠, m2(x)=(x−2)2(x−5)⁠, m3(x)=(x−2)(x−5)2⁠, and m4(x)=(x−2)2(x−5)2⁠. They are polynomials of degree two, three, three, and four.

(d)

The possibilities are m1(x)=(x+3)(x−1)(x−2)⁠, m2(x)=(x+3)2(x−1)(x−2)⁠, m3(x)=(x+3)(x−1)(x−2)2⁠, and m4(x)=(x+3)2(x−1)(x−2)2⁠. The degree of m1 is three, the degree of m2 is four, the degree of m3 is four, and the degree of m4 is five.

Five.IV.1.14 §IV.1 source ↗

For this matrix

(0101101001011010)

find the characteristic polynomial and the minimal polynomial.

Answer

The characteristic polynomial

|T−xI|=|−x1011−x1001−x1101−x|

is x2(x2−4)⁠.

By the Caley-Hamilton Theorem, all of the distinct roots of the characteristic polynomial are also roots of the minimal polynomial so the minimal polynomial has roots 0⁠, 2⁠, and −2⁠. Thus the minimal polynomial is either x(x2−4) or x2(x2−4)⁠.

One way to find the minimal polynomial is to try the matrix in the two expressions until we find the zero matrix.

T(T2−4I)=(0101101001011010)(−20100−20220−20020−2)=(0000000000000000)

So the minimal polynomial is x(x2−4)⁠.

Five.IV.1.15 §IV.1 source ↗

Find the minimal polynomial of this matrix.

(010001100)

Answer

Its characteristic polynomial has complex roots.

|−x100−x110−x|=(1−x)⋅(x−(−12+32i))⋅(x−(−12−32i))

As the roots are distinct, the characteristic polynomial equals the minimal polynomial.

Five.IV.1.16 §IV.1 source ↗

Find the minimal polynomial of each matrix.

(a)

(300130004)

(b)

(300130003)

(c)

(300130013)

(d)

(201062002)

(e)

(221062002)

(f)

(−14000030000−4−1003−9−42−115414)

Answer (7 parts)

In each case we will use the method of Example 1.12.

(a)

Because T is triangular, T−xI is also triangular

T−xI=(3−x0013−x0004−x)

the characteristic polynomial is easy c(x)= |T− xI|= (3−x)2(4−x)= − 1⋅(x−3)2(x−4)⁠. There are only two possibilities for the minimal polynomial, m1(x)=(x−3)(x−4) and m2(x)=(x−3)2(x−4)⁠. (Note that the characteristic polynomial has a negative sign but the minimal polynomial does not since it must have a leading coefficient of one). Because m1(T) is not the zero matrix

(T−3I)(T−4I)=(000100001)(−1001−10000)=(000−100000)

the minimal polynomial is m(x)=m2(x)⁠.

(T−3I)2(T−4I)=(T−3I)⋅((T−3I)(T−4I))=(000100001)(000−100000)=(000000000)

(b)

As in the prior item, the fact that the matrix is triangular makes computation of the characteristic polynomial easy.

c(x)=|T−xI|=|3−x0013−x0003−x|=(3−x)3=−1⋅(x−3)3

There are three possibilities for the minimal polynomial m1(x)=(x−3)⁠, m2(x)=(x−3)2⁠, and m3(x)=(x−3)3⁠. We settle the question by computing m1(T)

T−3I=(000100000)

and m2(T)⁠.

(T−3I)2=(000100000)(000100000)=(000000000)

Because m2(T) is the zero matrix, m2(x) is the minimal polynomial.

(c)

Again, the matrix is triangular.

c(x)=|T−xI|=|3−x0013−x0013−x|=(3−x)3=−1⋅(x−3)3

Again, there are three possibilities for the minimal polynomial m1(x)=(x−3)⁠, m2(x)=(x−3)2⁠, and m3(x)=(x−3)3⁠. We compute m1(T)

T−3I=(000100010)

and m2(T)

(T−3I)2=(000100010)(000100010)=(000000100)

and m3(T)⁠.

(T−3I)3=(T−3I)2(T−3I)=(000000100)(000100010)=(000000000)

Therefore, the minimal polynomial is m(x)=m3(x)=(x−3)3⁠.

(d)

This case is also triangular, here upper triangular.

c(x)=|T−xI|=|2−x0106−x2002−x|=(2−x)2(6−x)=−(x−2)2(x−6)

There are two possibilities for the minimal polynomial, m1(x)=(x−2)(x−6) and m2(x)=(x−2)2(x−6)⁠. Computation shows that the minimal polynomial isn’t m1(x)⁠.

(T−2I)(T−6I)=(001042000)(−40100200−4)=(00−4000000)

It therefore must be that m(x)=m2(x)=(x−2)2(x−6)⁠. Here is a verification.

(T−2I)2(T−6I)=(T−2I)⋅((T−2I)(T−6I))=(001042000)(00−4000000)=(000000000)

(e)

The characteristic polynomial is

c(x)=|T−xI|=|2−x2106−x2002−x|=(2−x)2(6−x)=−(x−2)2(x−6)

and there are two possibilities for the minimal polynomial, m1(x)=(x−2)(x−6) and m2(x)=(x−2)2(x−6)⁠. Checking the first one

(T−2I)(T−6I)=(021042000)(−42100200−4)=(000000000)

shows that the minimal polynomial is m(x)=m1(x)=(x−2)(x−6)⁠.

(f)

The characteristic polynomial is this.

c(x)=|T−xI|=|−1−x400003−x0000−4−1−x003−9−42−x−115414−x|=(x−3)3(x+1)2

Here are the possibilities for the minimal polynomial, listed here by ascending degree: m1(x)=(x−3)(x+1)⁠, m1(x)=(x−3)2(x+1)⁠, m1(x)=(x−3)(x+1)2⁠, m1(x)=(x−3)3(x+1)⁠, m1(x)=(x−3)2(x+1)2⁠, and m1(x)=(x−3)3(x+1)2⁠. The first one doesn’t pan out

(T−3I)(T+1I)=(−44000000000−4−4003−9−4−1−115411)(04000040000−40003−9−43−115415)=(000000000000000−4−40−4−444044)

but the second one does.

(T−3I)2(T+1I)=(T−3I)((T−3I)(T+1I))=(−44000000000−4−4003−9−4−1−115411)(000000000000000−4−40−4−444044)=(0000000000000000000000000)

The minimal polynomial is m(x)=(x−3)2(x+1)⁠.

Five.IV.1.17 §IV.1 source ↗

What is the minimal polynomial of the differentiation operator d/dx on 𝒫n⁠?

Answer

We know that 𝒫n is a dimension n+1 space and that the differentiation operator is nilpotent of index n+1 (for instance, taking n=3⁠, 𝒫3= {c3x3+ c2x2+ c1x+ c0|c3, …, c0∈ℂ} and the fourth derivative of a cubic is the zero polynomial). Represent this operator using the canonical form for nilpotent transformations.

(000…01000010⋱00010)

This is an (n+1)×(n+1) matrix with an easy characteristic polynomial, c(x)=xn+1⁠. (Remark: this matrix is RepB,B⁡(d/dx) where B= ⟨xn, nxn−1, n(n−1)xn−2, …, n!⟩⁠.) To find the minimal polynomial as in Example 1.12 we consider the powers of T−0I=T⁠. But, of course, the first power of T that is the zero matrix is the power n+1⁠. So the minimal polynomial is also xn+1⁠.

Five.IV.1.18 §IV.1 source ↗

Find the minimal polynomial of matrices of this form

(λ00…01λ0001λ⋱λ000…1λ)

where the scalar λ is fixed (i.e., is not a variable).

Answer

Call the matrix T and suppose that it is n×n⁠. Because T is triangular, and so T−xI is triangular, the characteristic polynomial is c(x)=(x−λ)n⁠. To see that the minimal polynomial is the same, consider T−λI⁠.

(000…0100…0010⋱00…10)

Recognize it as the canonical form for a transformation that is nilpotent of degree n⁠; the power (T−λI)j is zero first when j is n⁠.

Five.IV.1.19 §IV.1 source ↗

What is the minimal polynomial of the transformation of 𝒫n that sends p(x) to p(x+1)⁠?

Answer

The n=3 case provides a hint. A natural basis for 𝒫3 is B=⟨1,x,x2,x3⟩⁠. The action of the transformation is

1↦1x↦x+1x2↦x2+2x+1x3↦x3+3x2+3x+1

and so the representation RepB,B⁡(t) is this upper triangular matrix.

(1111012300130001)

Because it is triangular, the fact that the characteristic polynomial is c(x)=(x−1)4 is clear. For the minimal polynomial, the candidates are m1(x)=(x−1)⁠,

T−1I=(0111002300030000)

m2(x)=(x−1)2⁠,

(T−1I)2=(0026000600000000)

m3(x)=(x−1)3⁠,

(T−1I)3=(0006000000000000)

and m4(x)=(x−1)4⁠. Because m1⁠, m2⁠, and m3 are not right, m4 must be right, as is easily verified.

In the case of a general n⁠, the representation is an upper triangular matrix with ones on the diagonal. Thus the characteristic polynomial is c(x)=(x−1)n+1⁠. One way to verify that the minimal polynomial equals the characteristic polynomial is argue something like this: say that an upper triangular matrix is 0⁠-upper triangular if there are nonzero entries on the diagonal, that it is 1⁠-upper triangular if the diagonal contains only zeroes and there are nonzero entries just above the diagonal, etc. As the above example illustrates, an induction argument will show that, where T has only nonnegative entries, Tj is j⁠-upper triangular.

Five.IV.1.20 §IV.1 source ↗

What is the minimal polynomial of the map π:ℂ3→ℂ3 projecting onto the first two coordinates?

Answer

The map twice is the same as the map once: π∘π=π⁠, that is, π2=π and so the minimal polynomial is of degree at most two since m(x)=x2−x will do. The fact that no linear polynomial will do follows from applying the maps on the left and right side of c1⋅π+c0⋅id⁡=z (where z is the zero map) to these two vectors.

(001)(100)

Thus the minimal polynomial is m⁠.

Five.IV.1.21 §IV.1 source ↗

Find a 3×3 matrix whose minimal polynomial is x2⁠.

Answer

This is one answer.

(000100000)

Five.IV.1.22 §IV.1 source ↗

What is wrong with this claimed proof of Lemma 1.9: “if c(x)=|T−xI| then c(T)=|T−TI|=0⁠” ? [Cullen]

Answer

The x must be a scalar, not a matrix.

Five.IV.1.23 §IV.1 source ↗

Verify Lemma 1.9 for 2×2 matrices by direct calculation.

Answer

The characteristic polynomial of

T=(abcd)

is (a−x)(d−x)− bc= x2− (a+d)x+ (ad−bc)⁠. Substitute

(abcd)2−(a+d)(abcd)+(ad−bc)(1001)=(a2+bcab+bdac+cdbc+d2)−(a2+adab+bdac+cdad+d2)+(ad−bc00ad−bc)

and just check each entry sum to see that the result is the zero matrix.

Five.IV.1.24 §IV.1 source ↗

Prove that the minimal polynomial of an n×n matrix has degree at most n (not n2 as a person might guess from this subsection’s opening). Verify that this maximum, n⁠, can happen.

Answer

By the Cayley-Hamilton theorem the degree of the minimal polynomial is less than or equal to the degree of the characteristic polynomial, n⁠. Example 1.6 shows that n can happen.

Five.IV.1.25 §IV.1 source ↗

Show that, on a nontrivial vector space, a linear transformation is nilpotent if and only if its only eigenvalue is zero.

Answer

Let the linear transformation be t:V→V⁠. If t is nilpotent then there is an n such that tn is the zero map, so t satisfies the polynomial p(x)=xn=(x−0)n⁠. By Lemma 1.10 the minimal polynomial of t divides p⁠, so the minimal polynomial has only zero for a root. By Cayley-Hamilton, Theorem 1.8, the characteristic polynomial has only zero for a root. Thus the only eigenvalue of t is zero.

Conversely, if a transformation t on an n⁠-dimensional space has only the single eigenvalue of zero then its characteristic polynomial is xn⁠. The Cayley-Hamilton Theorem says that a map satisfies its characteristic polynomial so tn is the zero map. Thus t is nilpotent.

Five.IV.1.26 §IV.1 source ↗

What is the minimal polynomial of a zero map or matrix? Of an identity map or matrix?

Answer

A minimal polynomial must have leading coefficient 1⁠, and so if the minimal polynomial of a map or matrix were to be a degree zero polynomial then it would be m(x)=1⁠. But the identity map or matrix equals the zero map or matrix only on a trivial vector space.

So in the nontrivial case the minimal polynomial must be of degree at least one. A zero map or matrix has minimal polynomial m(x)=x⁠, and an identity map or matrix has minimal polynomial m(x)=x−1⁠.

Five.IV.1.27 §IV.1 source ↗

What is the minimal polynomial of a diagonal matrix?

Answer

For a diagonal matrix

T=(t1,100t2,2⋱tn,n)

the characteristic polynomial is (t1,1− x)(t2,2− x)⋯(tn,n− x)⁠. Of course, some of those factors may be repeated, e.g., the matrix might have t1,1=t2,2⁠. For instance, the characteristic polynomial of

D=(300030001)

is (3−x)2(1−x)= − 1⋅(x−3)2(x−1)⁠.

To form the minimal polynomial, take the terms x−ti,i⁠, throw out repeats, and multiply them together. For instance, the minimal polynomial of D is (x−3)(x−1)⁠. To check this, note first that Theorem 1.8, the Cayley-Hamilton theorem, requires that each linear factor in the characteristic polynomial appears at least once in the minimal polynomial. One way to check the other direction — that in the case of a diagonal matrix, each linear factor need appear at most once — is to use a matrix argument. A diagonal matrix, multiplying from the left, rescales rows by the entry on the diagonal. But in a product (T−t1,1I)⋯⁠, even without any repeat factors, every row is zero in at least one of the factors.

For instance, in the product

(D−3I)(D−1I)=(D−3I)(D−1I)I=(00000000−2)(200020000)(100010001)

because the first and second rows of the first matrix D−3I are zero, the entire product will have a first row and second row that are zero. And because the third row of the middle matrix D−1I is zero, the entire product has a third row of zero.

Five.IV.1.28 §IV.1 source ↗

A projection is any transformation t such that t2=t⁠. (For instance, consider the transformation of the plane ℝ2 projecting each vector onto its first coordinate. If we project twice then we get the same result as if we project just once.) What is the minimal polynomial of a projection?

Answer

This subsection starts with the observation that the powers of a linear transformation cannot climb forever without a “repeat”, that is, that for some power n there is a linear relationship cn⋅tn+⋯+c1⋅t+c0⋅id⁡=z where z is the zero transformation. The definition of projection is that for such a map one linear relationship is quadratic, t2−t=z⁠. To finish, we need only consider whether this relationship might not be minimal, that is, are there projections for which the minimal polynomial is constant or linear?

For the minimal polynomial to be constant, the map would have to satisfy that c0⋅id⁡=z⁠, where c0=1 since the leading coefficient of a minimal polynomial is 1⁠. This is only satisfied by the zero transformation on a trivial space. This is a projection, but not an interesting one.

For the minimal polynomial of a transformation to be linear would give c1⋅t+c0⋅id⁡=z where c1=1⁠. This equation gives t=−c0⋅id⁡⁠. Coupling it with the requirement that t2=t gives t2=(−c0)2⋅id⁡=−c0⋅id⁡⁠, which gives that c0=0 and t is the zero transformation or that c0=1 and t is the identity.

Thus, except in the cases where the projection is a zero map or an identity map, the minimal polynomial is m(x)=x2−x⁠.

Five.IV.1.29 §IV.1 source ↗

The first two items of this question are review.

(a)

Prove that the composition of one-to-one maps is one-to-one.

(b)

Prove that if a linear map is not one-to-one then at least one nonzero vector from the domain maps to the zero vector in the codomain.

(c)

Verify the statement, excerpted here, that precedes Theorem 1.8.

… if a minimal polynomial m(x) for a transformation t factors as m(x)=(x−λ1)q1⋯(x−λz)qz then m(t)= (t−λ1)q1∘⋯∘(t−λz)qz is the zero map. Since m(t) sends every vector to zero, at least one of the maps t−λi sends some nonzero vectors to zero. … That is, … at least some of the λi are eigenvalues.

Answer (3 parts)

(a)

This is a property of functions in general, not just of linear functions. Suppose that f and g are one-to-one functions such that f∘g is defined. Let f∘g(x1)=f∘g(x2)⁠, so that f(g(x1))=f(g(x2))⁠. Because f is one-to-one this implies that g(x1)=g(x2)⁠. Because g is also one-to-one, this in turn implies that x1=x2⁠. Thus, in summary, f∘g(x1)=f∘g(x2) implies that x1=x2 and so f∘g is one-to-one.

(b)

If the linear map h is not one-to-one then there are unequal vectors v→1⁠, v→2 that map to the same value h(v→1)=h(v→2)⁠. Because h is linear, we have 0→= h(v→1)− h(v→2)= h(v→1−v→2) and so v→1−v→2 is a nonzero vector from the domain that h maps to the zero vector of the codomain (v→1−v→2 does not equal the zero vector of the domain because v→1 does not equal v→2⁠).

(c)

The minimal polynomial m(t) sends every vector in the domain to zero and so it is not one-to-one (except in a trivial space, which we ignore). By the first item of this question, since the composition m(t) is not one-to-one, at least one of the components t−λi is not one-to-one. By the second item, t−λi has a nontrivial null space. Because (t−λi)(v→)=0→ holds if and only if t(v→)=λi⋅v→⁠, the prior sentence gives that λi is an eigenvalue (recall that the definition of eigenvalue requires that the relationship hold for at least one nonzero v→⁠).

Five.IV.1.30 §IV.1 source ↗

True or false: for a transformation on an n dimensional space, if the minimal polynomial has degree n then the map is diagonalizable.

Answer

This is false. The natural example of a non-diagonalizable transformation works here. Consider the transformation of ℂ2 represented with respect to the standard basis by this matrix.

N=(0100)

The characteristic polynomial is c(x)=x2⁠. Thus the minimal polynomial is either m1(x)=x or m2(x)=x2⁠. The first is not right since N−0⋅I is not the zero matrix, thus in this example the minimal polynomial has degree equal to the dimension of the underlying space, and, as mentioned, we know this matrix is not diagonalizable because it is nilpotent.

Five.IV.1.31 §IV.1 source ↗

Let f(x) be a polynomial. Prove that if A and B are similar matrices then f(A) is similar to f(B)⁠.

(a)

Now show that similar matrices have the same characteristic polynomial.

(b)

Show that similar matrices have the same minimal polynomial.

(c)

Decide if these are similar.

(1323)(4−111)

Answer (4 parts)

Let A and B be similar A=PBP−1⁠. From the facts that

An=(PBP−1)n=(PBP−1)(PBP−1)⋯(PBP−1)=PB(P−1P)B(P−1P)⋯(P−1P)BP−1=PBnP−1

and c⋅A=c⋅(PBP−1)=P(c⋅B)P−1 follows the required fact that for any polynomial function f we have f(A)=Pf(B)P−1⁠. For instance, if f(x)=x2+2x+3 then

A2+2A+3I=(PBP−1)2+2⋅PBP−1+3⋅I=(PBP−1)(PBP−1)+P(2B)P−1+3⋅PP−1=P(B2+2B+3I)P−1

shows that f(A) is similar to f(B)⁠.

(a)

Taking f to be a linear polynomial we have that A−xI is similar to B−xI⁠. Similar matrices have equal determinants (since |A|= |PBP−1|= |P|⋅|B|⋅|P−1|= 1⋅|B|⋅1= |B|⁠). Thus the characteristic polynomials are equal.

(b)

As P and P−1 are invertible, f(A) is the zero matrix when, and only when, f(B) is the zero matrix.

(c)

They cannot be similar since they don’t have the same characteristic polynomial. The characteristic polynomial of the first one is x2−4x−3 while the characteristic polynomial of the second is x2−5x+5⁠.

Five.IV.1.32 §IV.1 source ↗

(a)

Show that a matrix is invertible if and only if the constant term in its minimal polynomial is not 0⁠.

(b)

Show that if a square matrix T is not invertible then there is a nonzero matrix S such that ST and TS both equal the zero matrix.

Answer (3 parts)

Suppose that m(x)= xn+ mn−1xn−1+ ⋯+ m1x+ m0 is minimal for T⁠.

(a)

For the ‘if’ argument, because Tn+⋯+m1T+m0I is the zero matrix we have that I= (Tn+⋯+m1T)/(−m0)= T⋅(Tn−1+⋯+m1I)/(−m0) and so the matrix (−1/m0)⋅(Tn−1+⋯+m1I) is the inverse of T⁠. For ‘only if’, suppose that m0=0 (we put the n=1 case aside but it is easy) so that Tn+⋯+m1T=(Tn−1+⋯+m1I)T is the zero matrix. Note that Tn−1+⋯+m1I is not the zero matrix because the degree of the minimal polynomial is n⁠. If T−1 exists then multiplying both (Tn−1+⋯+m1I)T and the zero matrix from the right by T−1 gives a contradiction.

(b)

If T is not invertible then the constant term in its minimal polynomial is zero. Thus,

Tn+⋯+m1T=(Tn−1+⋯+m1I)T=T(Tn−1+⋯+m1I)

is the zero matrix.

Five.IV.1.33 §IV.1 source ↗

(a)

Finish the proof of Lemma 1.7.

(b)

Give an example to show that the result does not hold if t is not linear.

Answer (2 parts)

(a)

For the inductive step, assume that Lemma 1.7 is true for polynomials of degree i,…,k−1 and consider a polynomial f(x) of degree k⁠. Factor f(x)=k(x−λ1)q1⋯(x−λz)qz and let k(x−λ1)q1−1⋯(x−λz)qz be cn−1xn−1+⋯+c1x+c0⁠. Substitute:

k(t−λ1)q1∘⋯∘(t−λz)qz(v→)=(t−λ1)∘(t−λ1)q1∘⋯∘(t−λz)qz(v→)=(t−λ1)(cn−1tn−1(v→)+⋯+c0v→)=f(t)(v→)

(the second equality follows from the inductive hypothesis and the third from the linearity of t⁠).

(b)

One example is to consider the squaring map s:ℝ→ℝ given by s(x)=x2⁠. It is nonlinear. The action defined by the polynomial f(t)=t2−1 changes s to f(s)=s2−1⁠, which is this map.

x⟼s2−1s∘s(x)−1=x4−1

Observe that this map differs from the map (s−1)∘(s+1)⁠; for instance, the first map takes x=5 to 624 while the second one takes x=5 to 675⁠.

Five.IV.2.15 §IV.2 source ↗

Do the check for Example 2.4.

Answer

We must check that

(3013)=N+3I=PTP−1=(1/21/2−1/41/4)(2−114)(1−212)

That calculation is easy.

Five.IV.2.16 §IV.2 source ↗

Each matrix is in Jordan form. State its characteristic polynomial and its minimal polynomial.

(a)

(3013)

(b)

(−100−1)

(c)

(20012000−1/2)

(d)

(300130013)

(e)

(3000130000300013)

(f)

(4000140000−40001−4)

(g)

(500020003)

(h)

(5000020000200003)

(i)

(5000020001200003)

Answer (8 parts)

(a)

The characteristic polynomial is c(x)=(x−3)2 and the minimal polynomial is the same.

(b)

The characteristic polynomial is c(x)=(x+1)2⁠. The minimal polynomial is m(x)=x+1⁠.

(c)

The characteristic polynomial is c(x)=(x+(1/2))(x−2)2 and the minimal polynomial is the same.

(d)

The characteristic polynomial is c(x)=(x−3)3 The minimal polynomial is the same.

(e)

The characteristic polynomial is c(x)=(x−3)4⁠. The minimal polynomial is m(x)=(x−3)2⁠.

(f)

The characteristic polynomial is c(x)=(x+4)2(x−4)2 and the minimal polynomial is the same.

(g)

The characteristic polynomial is c(x)=(x−2)2(x−3)(x−5) and the minimal polynomial is m(x)=(x−2)(x−3)(x−5)⁠.

(h)

The characteristic polynomial is c(x)=(x−2)2(x−3)(x−5) and the minimal polynomial is the same.

Five.IV.2.17 §IV.2 source ↗

Find the Jordan form from the given data.

(a)

The matrix T is 5×5 with the single eigenvalue 3⁠. The nullities of the powers are: T−3I has nullity two, (T−3I)2 has nullity three, (T−3I)3 has nullity four, and (T−3I)4 has nullity five.

(b)

The matrix S is 5×5 with two eigenvalues. For the eigenvalue 2 the nullities are: S−2I has nullity two, and (S−2I)2 has nullity four. For the eigenvalue −1 the nullities are: S+1I has nullity one.

Answer (2 parts)

(a)

The transformation t−3 is nilpotent (that is, 𝒩∞(t−3) is the entire space) and it acts on a string basis via two strings, β→1↦β→2↦β→3↦β→4↦0→ and β→5↦0→⁠. Consequently, t−3 can be represented in this canonical form.

N3=(0000010000010000010000000)

and therefore T is similar to this canonical form matrix.

J3=N3+3I=(3000013000013000013000003)

(b)

The restriction of the transformation s+1 is nilpotent on the subspace 𝒩∞(s+1)⁠, and the action on a string basis is β→1↦0→⁠. The restriction of the transformation s−2 is nilpotent on the subspace 𝒩∞(s−2)⁠, having the action on a string basis of β→2↦β→3↦0→ and β→4↦β→5↦0→⁠. Consequently the Jordan form is this.

(−1000002000012000002000012)

Five.IV.2.18 §IV.2 source ↗

Find the change of basis matrices for each example.

(a)

Example 2.12

(b)

Example 2.13

(c)

Example 2.14

Answer (4 parts)

For each, because many choices of basis are possible, many other answers are possible. Of course, the calculation to check if an answer gives that PTP−1 is in Jordan form is the arbiter of what’s correct.

(a)

Here is the arrow diagram.

ℂwrt ℰ33→Ttℂwrt ℰ33id⁡↓Pid⁡↓Pℂwrt B3→Jtℂwrt B3

The matrix to move from the lower left to the upper left is this.

P−1=(Repℰ3,B⁡(id⁡))−1=RepB,ℰ3⁡(id⁡)=(1−20101−200)

The matrix P to move from the upper right to the lower right is the inverse of P−1⁠.

(b)

We want this matrix and its inverse.

P−1=(1030140−20)

(c)

The concatenation of these bases for the generalized null spaces will do for the basis for the entire space.

B−1=⟨(−10010),(−10−101)⟩B3=⟨(11−100),(000−22),(−1−1120)⟩

The change of basis matrices are this one and its inverse.

P−1=(−1−110−10010−10−1−101100−2201020)

Five.IV.2.19 §IV.2 source ↗

Find the Jordan form and a Jordan basis for each matrix.

(a)

(−104−2510)

(b)

(5−49−7)

(c)

(400213504)

(d)

(543−10−31−21)

(e)

(973−9−7−4444)

(f)

(22−1−1−11−1−22)

(g)

(712214−1−1−215−11128)

Answer (8 parts)

The general procedure is to factor the characteristic polynomial c(x)=(x−λ1)p1(x−λ2)p2⋯ to get the eigenvalues λ1⁠, λ2⁠, etc. Then, for each λi we find a string basis for the action of the transformation t−λi when restricted to 𝒩∞(t−λi)⁠, by computing the powers of the matrix T−λiI and finding the associated null spaces, until these null spaces settle down (do not change), at which point we have the generalized null space. The dimensions of those null spaces (the nullities) tell us the action of t−λi on a string basis for the generalized null space, and so we can write the pattern of subdiagonal ones to have Nλi⁠. From this matrix, the Jordan block Jλi associated with λi is immediate Jλi=Nλi+λiI⁠. Finally, after we have done this for each eigenvalue, we put them together into the canonical form.

(a)

The characteristic polynomial of this matrix is c(x)= (−10−x)(10−x)+ 100= x2⁠, so it has only the single eigenvalue λ=0⁠.

power p(T+0⋅I)p𝒩((t−0)p)nullity
1(−104−2510){(2y/5y)|y∈ℂ}1
2(0000)ℂ22

(Thus, this transformation is nilpotent: 𝒩∞(t−0) is the entire space). From the nullities we know that t⁠’s action on a string basis is β→1↦β→2↦0→⁠. This is the canonical form matrix for the action of t−0 on 𝒩∞(t−0)=ℂ2

N0=(0010)

and this is the Jordan form of the matrix.

J0=N0+0⋅I=(0010)

Note that if a matrix is nilpotent then its canonical form equals its Jordan form.

We can find such a string basis using the techniques of the prior section.

B=⟨(10),(−10−25)⟩

We took the first basis vector so that it is in the null space of t2 but is not in the null space of t⁠. The second basis vector is the image of the first under t⁠.

(b)

The characteristic polynomial of this matrix is c(x)=(x+1)2⁠, so it is a single-eigenvalue matrix. (That is, the generalized null space of t+1 is the entire space.) We have

𝒩(t+1)={(2y/3y)|y∈ℂ}𝒩((t+1)2)=ℂ2

and so the action of t+1 on an associated string basis is β→1↦β→2↦0→⁠. Thus,

N−1=(0010)

the Jordan form of T is

J−1=N−1+−1⋅I=(−101−1)

and choosing vectors from the above null spaces gives this string basis (other choices are possible).

B=⟨(10),(69)⟩

(c)

The characteristic polynomial c(x)= (1−x)(4−x)2= − 1⋅(x−1)(x−4)2 has two roots and they are the eigenvalues λ1=1 and λ2=4⁠.

We handle the two eigenvalues separately. For λ1⁠, the calculation of the powers of T−1I yields

𝒩(t−1)={(0y0)|y∈ℂ}

and the null space of (t−1)2 is the same. Thus this set is the generalized null space 𝒩∞(t−1)⁠. The nullities show that the action of the restriction of t−1 to the generalized null space on a string basis is β→1↦0→⁠.

A similar calculation for λ2=4 gives these null spaces.

𝒩(t−4)={(0zz)|z∈ℂ}𝒩((t−4)2)={(z−yyz)|y,z∈ℂ}

(The null space of (t−4)3 is the same, as it must be because the power of the term associated with λ2=4 in the characteristic polynomial is two, and so the restriction of t−2 to the generalized null space 𝒩∞(t−2) is nilpotent of index at most two — it takes at most two applications of t−2 for the null space to settle down.) The pattern of how the nullities rise tells us that the action of t−4 on an associated string basis for 𝒩∞(t−4) is β→2↦β→3↦0→⁠.

Putting the information for the two eigenvalues together gives the Jordan form of the transformation t⁠.

(100040014)

We can take elements of the null spaces to get an appropriate basis.

B=B1⌢B4=⟨(010),(101),(055)⟩

(d)

The characteristic polynomial is c(x)= (−2−x)(4−x)2= − 1⋅(x+2)(x−4)2⁠.

For the eigenvalue λ−2 (the subscript there is −2⁠), calculation of the powers of T+2I yields this.

𝒩(t+2)={(zzz)|z∈ℂ}

The null space of (t+2)2 is the same, and so this is the generalized null space 𝒩∞(t+2)⁠. Thus the action of the restriction of t+2 to 𝒩∞(t+2) on an associated string basis is β→1↦0→⁠.

For λ2=4⁠, computing the powers of T−4I yields

𝒩(t−4)={(z−zz)|z∈ℂ}𝒩((t−4)2)={(x−zz)|x,z∈ℂ}

and so the action of t−4 on a string basis for 𝒩∞(t−4) is β→2↦β→3↦0→⁠.

Therefore the Jordan form is

(−200040014)

and a suitable basis is this.

B=B−2⌢B4=⟨(111),(0−11),(−11−1)⟩

(e)

The characteristic polynomial of this matrix is c(x)=(2−x)3=−1⋅(x−2)3⁠. This matrix has only a single eigenvalue, λ=2⁠. By finding the powers of T−2I we have

𝒩(t−2)={(−yy0)|y∈ℂ}𝒩((t−2)2)={(−y−(1/2)zyz)|y,z∈ℂ}

and

𝒩((t−2)3)=ℂ3

and so the action of t−2 on an associated string basis is β→1↦β→2↦β→3↦0→⁠. The Jordan form is this

(200120012)

and one choice of basis is this.

B=⟨(010),(7−94),(−220)⟩

(f)

The characteristic polynomial c(x)=(1−x)3=−(x−1)3 has only a single root, so the matrix has only a single eigenvalue λ=1⁠. Finding the powers of T−1I and calculating the null spaces

𝒩(t−1)={(−2y+zyz)|y,z∈ℂ}𝒩((t−1)2)=ℂ3

shows that the action of the nilpotent map t−1 on a string basis is β→1↦β→2↦0→ and β→3↦0→⁠. Therefore the Jordan form is

J=(100110001)

and an appropriate basis (a string basis associated with t−1⁠) is this.

B=⟨(010),(2−2−2),(101)⟩

(g)

The characteristic polynomial is a bit large for by-hand calculation, but just manageable c(x)= x4− 24x3+ 216x2− 864x+ 1296= (x−6)4⁠. This is a single-eigenvalue map, so the transformation t−6 is nilpotent. The null spaces

𝒩(t−6)={(−z−w−z−wzw)|z,w∈ℂ}𝒩((t−6)2)={(x−z−wzw)|x,z,w∈ℂ}

and

𝒩((t−6)3)=ℂ4

and the nullities show that the action of t−6 on a string basis is β→1↦β→2↦β→3↦0→ and β→4↦0→⁠. The Jordan form is

(6000160001600006)

and finding a suitable string basis is routine.

B=⟨(0001),(2−1−12),(33−63),(−1−110)⟩

Five.IV.2.20 §IV.2 source ↗

Find all possible Jordan forms of a transformation with characteristic polynomial (x−1)2(x+2)2⁠.

Answer

There are two eigenvalues, λ1=−2 and λ2=1⁠. The restriction of t+2 to 𝒩∞(t+2) could have either of these actions on an associated string basis.

β→1↦β→2↦0→β→1↦0→β→2↦0→

The restriction of t−1 to 𝒩∞(t−1) could have either of these actions on an associated string basis.

β→3↦β→4↦0→β→3↦0→β→4↦0→

In combination, that makes four possible Jordan forms, the two first actions, the second and first, the first and second, and the two second actions.

(−20001−20000100011)(−20000−20000100011)(−20001−20000100001)(−20000−20000100001)

Five.IV.2.21 §IV.2 source ↗

Find all possible Jordan forms of a transformation with characteristic polynomial (x−1)3(x+2)⁠.

Answer

The restriction of t+2 to 𝒩∞(t+2) can have only the action β→1↦0→⁠. The restriction of t−1 to 𝒩∞(t−1) could have any of these three actions on an associated string basis.

β→2↦β→3↦β→4↦0→β→2↦β→3↦0→β→4↦0→β→2↦0→β→3↦0→β→4↦0→

Taken together there are three possible Jordan forms, the one arising from the first action by t−1 (along with the only action from t+2⁠), the one arising from the second action, and the one arising from the third action.

(−2000010001100011)(−2000010001100001)(−2000010000100001)

Five.IV.2.22 §IV.2 source ↗

Find all possible Jordan forms of a transformation with characteristic polynomial (x−2)3(x+1) and minimal polynomial (x−2)2(x+1)⁠.

Answer

The action of t+1 on a string basis for 𝒩∞(t+1) must be β→1↦0→⁠. Because of the power of x−2 in the minimal polynomial, a string basis for t−2 has length two and so the action of t−2 on 𝒩∞(t−2) must be of this form.

β→2↦β→3↦0→β→4↦0→

Therefore there is only one Jordan form that is possible.

(−1000020001200002)

Five.IV.2.23 §IV.2 source ↗

Find all possible Jordan forms of a transformation with characteristic polynomial (x−2)4(x+1) and minimal polynomial (x−2)2(x+1)⁠.

Answer

There are two possible Jordan forms. The action of t+1 on a string basis for 𝒩∞(t+1) must be β→1↦0→⁠. There are two actions for t−2 on a string basis for 𝒩∞(t−2) that are possible with this characteristic polynomial and minimal polynomial.

β→2↦β→3↦0→β→4↦β→5↦0→β→2↦β→3↦0→β→4↦0→β→5↦0→

The resulting Jordan form matrices are these.

(−1000002000012000002000012)(−1000002000012000002000002)

Five.IV.2.24 §IV.2 source ↗

Diagonalize these.

(a)

(1100)

(b)

(0110)

Answer (2 parts)

(a)

The characteristic polynomial is c(x)=x(x−1)⁠. For λ1=0 we have

𝒩(t−0)={(−yy)|y∈ℂ}

(of course, the null space of t2 is the same). For λ2=1⁠,

𝒩(t−1)={(x0)|x∈ℂ}

(and the null space of (t−1)2 is the same). We can take this basis

B=⟨(1−1),(10)⟩

to get the diagonalization.

(11−10)−1(1100)(11−10)=(0001)

(b)

The characteristic polynomial is c(x)=x2−1=(x+1)(x−1)⁠. For λ1=−1⁠,

𝒩(t+1)={(−yy)|y∈ℂ}

and the null space of (t+1)2 is the same. For λ2=1

𝒩(t−1)={(yy)|y∈ℂ}

and the null space of (t−1)2 is the same. We can take this basis

B=⟨(1−1),(11)⟩

to get a diagonalization.

(111−1)−1(0110)(11−11)=(−1001)

Five.IV.2.25 §IV.2 source ↗

Find the Jordan matrix representing the differentiation operator on 𝒫3⁠.

Answer

The transformation d/dx:𝒫3→𝒫3 is nilpotent. Its action on the basis B=⟨x3,3x2,6x,6⟩ is x3↦3x2↦6x↦6↦0⁠. Its Jordan form is its canonical form as a nilpotent matrix.

J=(0000100001000010)

Five.IV.2.26 §IV.2 source ↗

Decide if these two are similar.

(1−14−3)(−101−1)

Answer

Yes. Each has the characteristic polynomial (x+1)2⁠. Calculations of the powers of T1+1⋅I and T2+1⋅I gives these two.

𝒩(t1+1)={(y/2y)|y∈ℂ}𝒩(t2+1)={(0y)|y∈ℂ}

(Of course, for each the null space of the square is the entire space.) The way that the nullities rise shows that each is similar to this Jordan form matrix

(−101−1)

and they are therefore similar to each other.

Five.IV.2.27 §IV.2 source ↗

Find the Jordan form of this matrix.

(0−110)

Also give a Jordan basis.

Answer

Its characteristic polynomial is c(x)=x2+1 which has complex roots x2+1=(x+i)(x−i)⁠. Because the roots are distinct, the matrix is diagonalizable and its Jordan form is that diagonal matrix.

(−i00i)

To find an associated basis we compute the null spaces.

𝒩(t+i)={(−iyy)|y∈ℂ}𝒩(t−i)={(iyy)|y∈ℂ}

For instance,

T+i⋅I=(i−11i)

and so we get a description of the null space of t+i by solving this linear system.

ix−y=0x+iy=0⟶iρ1+ρ2ix−y=00=0

(To change the relation ix=y so that the leading variable x is expressed in terms of the free variable y⁠, we can multiply both sides by −i⁠.)

As a result, one such basis is this.

B=⟨(−i1),(i1)⟩

Five.IV.2.28 §IV.2 source ↗

How many similarity classes are there for 3×3 matrices whose only eigenvalues are −3 and 4⁠?

Answer

We can count the possible classes by counting the possible canonical representatives, that is, the possible Jordan form matrices. The characteristic polynomial must be either c1(x)=(x+3)2(x−4) or c2(x)=(x+3)(x−4)2⁠. In the c1 case there are two possible actions of t+3 on a string basis for 𝒩∞(t+3)⁠.

β→1↦β→2↦0→β→1↦0→β→2↦0→

There are two associated Jordan form matrices.

(−3001−30004)(−3000−30004)

Similarly there are two Jordan form matrices that could arise out of c2⁠.

(−300040014)(−300040004)

So in total there are four possible Jordan forms.

Five.IV.2.29 §IV.2 source ↗

Prove that a matrix is diagonalizable if and only if its minimal polynomial has only linear factors.

Answer

Jordan form is unique. A diagonal matrix is in Jordan form. Thus the Jordan form of a diagonalizable matrix is its diagonalization. If the minimal polynomial has factors to some power higher than one then the Jordan form has subdiagonal 1⁠’s, and so is not diagonal.

Five.IV.2.30 §IV.2 source ↗

Give an example of a linear transformation on a vector space that has no non-trivial invariant subspaces.

Answer

One example is the transformation of ℂ that sends x to −x⁠.

Five.IV.2.31 §IV.2 source ↗

Show that a subspace is t−λ1 invariant if and only if it is t−λ2 invariant.

Answer

Apply Lemma 2.11 twice; the subspace is t−λ1⁠ invariant if and only if it is t⁠ invariant, which in turn holds if and only if it is t−λ2⁠ invariant.

Five.IV.2.32 §IV.2 source ↗

Prove or disprove: two n×n matrices are similar if and only if they have the same characteristic and minimal polynomials.

Answer

False; these two 4×4 matrices each have c(x)=(x−3)4 and m(x)=(x−3)2⁠.

(3000130000300013)(3000130000300003)

Five.IV.2.33 §IV.2 source ↗

The trace of a square matrix is the sum of its diagonal entries.

(a)

Find the formula for the characteristic polynomial of a 2×2 matrix.

(b)

Show that trace is invariant under similarity, and so we can sensibly speak of the ‘trace of a map’. (Hint: see the prior item.)

(c)

Is trace invariant under matrix equivalence?

(d)

Show that the trace of a map is the sum of its eigenvalues (counting multiplicities).

(e)

Show that the trace of a nilpotent map is zero. Does the converse hold?

Answer (5 parts)

(a)

The characteristic polynomial is this.

|a−xbcd−x|=(a−x)(d−x)−bc=ad−(a+d)x+x2−bc=x2−(a+d)x+(ad−bc)

Note that the determinant appears as the constant term.

(b)

Recall that the characteristic polynomial |T−xI| is invariant under similarity. Use the permutation expansion formula to show that the trace is the negative of the coefficient of xn−1⁠.

(c)

No, there are matrices T and S that are equivalent S=PTQ (for some nonsingular P and Q⁠) but that have different traces. An easy example is this.

PTQ=(2001)(1001)(1001)=(2001)

Even easier examples using 1×1 matrices are possible.

(d)

Put the matrix in Jordan form. By the first item, the trace is unchanged.

(e)

The first part is easy; use the third item. The converse does not hold: this matrix

(100−1)

has a trace of zero but is not nilpotent.

Five.IV.2.34 §IV.2 source ↗

To use Definition 2.10 to check whether a subspace is t⁠ invariant, we seemingly have to check all of the infinitely many vectors in a (nontrivial) subspace to see if they satisfy the condition. Prove that a subspace is t⁠ invariant if and only if its subbasis has the property that for all of its elements, t(β→) is in the subspace.

Answer

Suppose that BM is a basis for a subspace M of some vector space. Implication one way is clear; if M is t invariant then in particular, if m→∈BM then t(m→)∈M⁠. For the other implication, let BM=⟨β→1,…,β→q⟩ and note that t(m→)= t(m1β→1+⋯+mqβ→q)= m1t(β→1)+ ⋯+ mqt(β→q) is in M as any subspace is closed under linear combinations.

Five.IV.2.35 §IV.2 source ↗

Is t invariance preserved under intersection? Under union? Complementation? Sums of subspaces?

Answer

Yes, the intersection of t invariant subspaces is t⁠ invariant. Assume that M and N are t invariant. If v→∈M∩N then t(v→)∈M by the invariance of M and t(v→)∈N by the invariance of N⁠.

Of course, the union of two subspaces need not be a subspace (remember that the x⁠- and y⁠-axes are subspaces of the plane ℝ2 but the union of the two axes fails to be closed under vector addition; for instance it does not contain e→1+e→2⁠.) However, the union of invariant subsets is an invariant subset; if v→∈M∪N then v→∈M or v→∈N so t(v→)∈M or t(v→)∈N⁠.

No, the complement of an invariant subspace need not be invariant. Consider the subspace

{(x0)|x∈ℂ}

of ℂ2 under the zero transformation.

Yes, the sum of two invariant subspaces is invariant. The check is easy.

Five.IV.2.36 §IV.2 source ↗

Give a way to order the Jordan blocks if some of the eigenvalues are complex numbers. That is, suggest a reasonable ordering for the complex numbers.

Answer

One such ordering is the dictionary ordering. Order by the real component first, then by the coefficient of i⁠. For instance, 3+2i<4+1i but 4+1i<4+2i⁠.

Five.IV.2.37 §IV.2 source ↗

Let 𝒫j(ℝ) be the vector space over the reals of degree j polynomials. Show that if j≤k then 𝒫j(ℝ) is an invariant subspace of 𝒫k(ℝ) under the differentiation operator. In 𝒫7(ℝ)⁠, does any of 𝒫0(ℝ)⁠, …, 𝒫6(ℝ) have an invariant complement?

Answer

The first half is easy — the derivative of any real polynomial is a real polynomial of lower degree. The answer to the second half is ‘no’; any complement of 𝒫j(ℝ) must include a polynomial of degree j+1⁠, and the derivative of that polynomial is in 𝒫j(ℝ)⁠.

Five.IV.2.38 §IV.2 source ↗

In 𝒫n(ℝ)⁠, the vector space (over the reals) of degree n polynomials,

ℰ={p(x)∈𝒫n(ℝ)|p(−x)=p(x) for all x}

and

𝒪={p(x)∈𝒫n(ℝ)|p(−x)=−p(x) for all x}

are the even and the odd polynomials; p(x)=x2 is even while p(x)=x3 is odd. Show that they are subspaces. Are they complementary? Are they invariant under the differentiation transformation?

Answer

For the first half, show that each is a subspace and then observe that any polynomial can be uniquely written as the sum of even-powered and odd-powered terms (the zero polynomial is both). The answer to the second half is ‘no’: x2 is even while 2x is odd.

Five.IV.2.39 §IV.2 source ↗

A matrix S is the square root of another T if S2=T⁠. Show that any nonsingular matrix has a square root.

Answer

Put the matrix in Jordan form. By non-singularity, there are no zero eigenvalues on the diagonal. Ape this example:

(900190004)=(3001/630002)2

to construct a square root. Show that it holds up under similarity: if S2=T then (PSP−1)(PSP−1)=PTP−1⁠.

1 Topic: Method of Powers source ↗

Use ten iterations to estimate the largest eigenvalue of these matrices, starting from the vector with components 1 and 2⁠. Compare the answer with the one obtained by solving the characteristic equation.

(a)

(1504)

(b)

(32−10)

Answer (2 parts)

(a)

By eye, we see that the largest eigenvalue is 4⁠. Sage gives this.

sage: def eigen(M,v,num_loops=10):
....:     for p in range(num_loops):
....:         v_normalized = (1/v.norm())*v
....:         v = M*v
....:     return v 
....: 
sage: M = matrix(RDF, [[1,5], [0,4]])
sage: v = vector(RDF, [1, 2])
sage: v = eigen(M,v)
sage: (M*v).dot_product(v)/v.dot_product(v)
4.00000147259
(b)

A simple calculation shows that the largest eigenvalue is 2⁠. Sage gives this.

sage: M = matrix(RDF, [[3,2], [-1,0]])
sage: v = vector(RDF, [1, 2])
sage: v = eigen(M,v)
sage: (M*v).dot_product(v)/v.dot_product(v)
2.00097741083

2 Topic: Method of Powers source ↗

Redo the prior exercise by iterating until |v→k|−|v→k−1| has absolute value less than 0.01 At each step, normalize by dividing each vector by its length. How many iterations does it take? Are the answers significantly different?

Answer (2 parts)

(a)

Here is Sage.

sage: def eigen_by_iter(M, v, toler=0.01):
....:     dex = 0
....:     diff = 10
....:     while abs(diff)>toler:
....:         dex = dex+1
....:         v_next = M*v
....:         v_normalized = (1/v.norm())*v
....:         v_next_normalized = (1/v_next.norm())*v_next
....:         diff = (v_next_normalized-v_normalized).norm()
....:         v_prior = v_normalized
....:         v = v_next_normalized
....:     return v, v_prior, dex
....: 
sage: M = matrix(RDF, [[1,5], [0,4]])
sage: v = vector(RDF, [1, 2])
sage: v,v_prior,dex = eigen_by_iter(M,v)
sage: (M*v).norm()/v.norm()
4.00604111686
sage: dex
4
(b)

Sage takes a few more iterations on this one. This makes use of the procedure defined in the prior item.

sage: M = matrix(RDF, [[3,2], [-1,0]])
sage: v = vector(RDF, [1, 2])
sage: v,v_prior,dex = eigen_by_iter(M,v)
sage: (M*v).norm()/v.norm()
2.01585174302
sage: dex
6

3 Topic: Method of Powers source ↗

Use ten iterations to estimate the largest eigenvalue of these matrices, starting from the vector with components 1⁠, 2⁠, and 3⁠. Compare the answer with the one obtained by solving the characteristic equation.

(a)

(401−210−201)

(b)

(−122222−3−6−6)

Answer (2 parts)

(a)

The largest eigenvalue is 3⁠. Sage gives this.

sage: M = matrix(RDF, [[4,0,1], [-2,1,0], [-2,0,1]])
sage: v = vector(RDF, [1, 2, 3])
sage: v = eigen(M,v)
sage: (M*v).dot_product(v)/v.dot_product(v)
3.02362112326
(b)

The largest eigenvalue is −3⁠.

sage: M = matrix(RDF, [[-1,2,2], [2,2,2], [-3,-6,-6]])
sage: v = vector(RDF, [1, 2, 3])
sage: v = eigen(M,v)
sage: (M*v).dot_product(v)/v.dot_product(v)
-3.00941127145

4 Topic: Method of Powers source ↗

Redo the prior exercise by iterating until |v→k|−|v→k−1| has absolute value less than 0.01⁠. At each step, normalize by dividing each vector by its length. How many iterations does it take? Are the answers significantly different?

Answer (2 parts)

(a)

Sage gets this, where eigen_by_iter is defined above.

sage: M = matrix(RDF, [[4,0,1], [-2,1,0], [-2,0,1]])
sage: v = vector(RDF, [1, 2, 3])
sage: v,v_prior,dex = eigen_by_iter(M,v) 
sage: (M*v).dot_product(v)/v.dot_product(v)
3.05460392934
sage: dex
8
(b)

With this setup,

sage: M = matrix(RDF, [[-1,2,2], [2,2,2], [-3,-6,-6]])
sage: v = vector(RDF, [1, 2, 3])

Sage does not return (use <Ctrl>-c to interrupt the computation). Adding some error checking code to the routine

def eigen_by_iter(M, v, toler=0.01):
    dex = 0
    diff = 10
    while abs(diff)>toler:
        dex = dex+1
        if dex>1000:
            print "oops! probably in some loop: \nv=",v,"\nv_next=",v_next
        v_next = M*v
        if (v.norm()==0):
            print "oops! v is zero"
	    return None
        if (v_next.norm()==0):
            print "oops! v_next is zero"
            return None
        v_normalized = (1/v.norm())*v
        v_next_normalized = (1/v_next.norm())*v_next
        diff = (v_next_normalized-v_normalized).norm()
        v_prior = v_normalized
        v = v_next_normalized
    return v, v_prior, dex

gives this.

oops! probably in some loop: 
  v= (0.707106781187, -1.48029736617e-16, -0.707106781187) 
  v_next= (2.12132034356, -4.4408920985e-16, -2.12132034356) 
oops! probably in some loop: 
  v= (-0.707106781187, 1.48029736617e-16, 0.707106781187) 
  v_next= (-2.12132034356, 4.4408920985e-16, 2.12132034356)
oops! probably in some loop: 
  v= (0.707106781187, -1.48029736617e-16, -0.707106781187) 
  v_next= (2.12132034356, -4.4408920985e-16, -2.12132034356)

So it is circling.

5 Topic: Method of Powers source ↗

What happens if c1=0⁠? That is, what happens if the initial vector does not to have any component in the direction of the relevant eigenvector?

Answer

In theory, this method would produce λ2⁠. In practice, however, rounding errors in the computation introduce components in the direction of v→1⁠, and so the method will still produce λ1⁠, although it may take somewhat longer than it would have taken with a more fortunate choice of initial vector.

6 Topic: Method of Powers source ↗

How can we adapt the method of powers to find the smallest eigenvalue?

Answer

Instead of using v→k=Tv→k−1⁠, use T−1v→k=v→k−1⁠.

1 Topic: Stable Populations source ↗

For the park discussed above, what should be the initial park population in the case where the populations decline by 11% every year?

Answer

The equation

0.89I−T=(0.89000.89)−(.90.01.10.99)=(−.01−.01−.10−.10)

leads to this system.

(−.01−.01−.10−.10)(pr)=(00)

So the eigenvectors have p=−r⁠.

2 Topic: Stable Populations source ↗

What will happen to the population of the park in the event of a growth in the external population of 1% per year? Will the park growth rate lag the world or lead it? Assume that the initial park population is ten thousand, and the world population is one hundred thousand, and calculate over a ten year span.

Answer

Sage finds this.

sage: M = matrix(RDF, [[0.90,0.01], [0.10,0.99]])
sage: v = vector(RDF, [10000, 100000])
sage: for y in range(10):
....:     v[1] = v[1]*(1+.01)^y
....:     print "pop vector year",y," is",v
....:     v = M*v
....:     
pop vector year 0  is (10000.0, 100000.0)
pop vector year 1  is (10000.0, 101000.0)
pop vector year 2  is (10010.0, 103019.899)
pop vector year 3  is (10039.19899, 106111.421211)
pop vector year 4  is (10096.3933031, 110360.453787)
pop vector year 5  is (10190.3585107, 115891.187687)
pop vector year 6  is (10330.2345365, 122872.349786)
pop vector year 7  is (10525.9345807, 131525.973067)
pop vector year 8  is (10788.6008533, 142139.351965)
pop vector year 9  is (11131.1342876, 155081.09214)

So inside the park the population grows by about eleven percent while outside the park the population grows by about fifty five percent.

3 Topic: Stable Populations source ↗

The park discussed above is partially fenced so that now, every year, only 5% of the animals from inside of the park leave (still, about 1% of the animals from the outside find their way in). Under what conditions can the park maintain a stable population now?

Answer

We want that next year the population of the park is a combination of animals that stay and animals that enter from outside, pn+1=0.95pn+0.01rn⁠. And, next year’s population in the rest of the world is a combination of animals that have left the park and animals that stayed in the rest of the world rn+1=0.05pn+0.99rn⁠. The matrix equation

(pn+1rn+1)=(0.950.010.050.99)(pnrn)

means that to find eigenvalues we want to solve this.

0=|λ−0.950.010.10λ−0.99|=λ2−1.94λ−0.9405

Sage gives this.

sage: a,b,c = var('a,b,c') 
sage: qe = (x^2 - 1.94*x -.9405 == 0)
sage: print solve(qe, x)
[
x == -1/100*sqrt(18814) + 97/100,
x == 1/100*sqrt(18814) + 97/100
]
sage: n(-1/100*sqrt(18814) + 97/100)
-0.401641352540816
sage: n(1/100*sqrt(18814) + 97/100)
2.34164135254082

So the only way to have a dynamically stable population is if the park’s numbers grows by 234% every year.

4 Topic: Stable Populations source ↗

Suppose that a species of bird only lives in Canada, the United States, or in Mexico. Every year, 4% of the Canadian birds travel to the US, and 1% of them travel to Mexico. Every year, 6% of the US birds travel to Canada, and 4% go to Mexico. From Mexico, every year 10% travel to the US, and 0% go to Canada.

(a)

Give the transition matrix.

(b)

Is there a way for the three countries to have constant populations?

Answer (2 parts)

(a)

This is the recurrence.

(cn+1un+1mn+1)=(.95.060.04.90.10.01.04.90)(cnunmn)

(b)

The system

.05c−.06u=0−.04c+.10u−.10m=0−.01c−.04u+.10m=0

has infinitely many solutions.

sage: var('c,u,m')
(c, u, m)
sage: eqns = [ .05*c-.06*u  == 0,
....:         -.04*c+.10*u-.10*m == 0,
....:         -.01*c-.04*u+.10*m == 0 ]
sage: solve(eqns, c,u,m)
[[c == 30/13*r1, u == 25/13*r1, m == r1]]

1 Topic: Page Ranking source ↗

A square matrix is stochastic if the sum of the entries in each column is one. The Google matrix is computed by taking a combination G=α∗H+(1−α)∗S of two stochastic matrices. Show that G must be stochastic.

Answer

The sum of the entries in column j is ∑iαhi,j+ (1−α)si,j= ∑iαhi,j+ ∑i(1−α)si,j= α∑iαhi,j+ (1−α)∑isi,j= α⋅1+ (1−α)⋅1⁠, which is one.

2 Topic: Page Ranking source ↗

For this web of pages, the importance of each page should be equal. Verify it for α=0.85⁠.

Answer

This Sage session gives equal values.

sage: H=matrix(QQ,[[0,0,0,1], [1,0,0,0], [0,1,0,0], [0,0,1,0]])
sage: S=matrix(QQ,[[1/4,1/4,1/4,1/4], [1/4,1/4,1/4,1/4], [1/4,1/4,1/4,1/4], [1/4,1/4,1/4,1/4]])
sage: alpha=0.85     
sage: G=alpha*H+(1-alpha)*S    
sage: I=matrix(QQ,[[1,0,0,0], [0,1,0,0], [0,0,1,0], [0,0,0,1]])
sage: N=G-I
sage: 1200*N         
[-1155.00000000000  45.0000000000000  45.0000000000000  1065.00000000000]
[ 1065.00000000000 -1155.00000000000  45.0000000000000  45.0000000000000]
[ 45.0000000000000  1065.00000000000 -1155.00000000000  45.0000000000000]
[ 45.0000000000000  45.0000000000000  1065.00000000000 -1155.00000000000]
sage: M=matrix(QQ,[[-1155,45,45,1065], [1065,-1155,45,45], [45,1065,-1155,45], [45,45,1065,-1155]]) 
sage: M.echelon_form()       
[ 1  0  0 -1]
[ 0  1  0 -1]
[ 0  0  1 -1]
[ 0  0  0  0]
sage: v=vector([1,1,1,1])
sage: (v/v.norm()).n()
(0.500000000000000, 0.500000000000000, 0.500000000000000, 0.500000000000000)

3 Topic: Page Ranking source ↗

[Bryan & Leise] Give the importance ranking for this web of pages.

(a)

Use α=0.85⁠.

(b)

Use α=0.95⁠.

(c)

Observe that while p3 is linked-to from all other pages, and therefore seems important, it is not the highest ranked page. What is the highest ranked page? Explain.

Answer (4 parts)

We have this.

H=(0011/21/30001/31/201/21/31/200)

(a)

This Sage session gives the answer.

sage: H=matrix(QQ,[[0,0,1,1/2], [1/3,0,0,0], [1/3,1/2,0,1/2], [1/3,1/2,0,0]])
sage: S=matrix(QQ,[[1/4,1/4,1/4,1/4], [1/4,1/4,1/4,1/4], [1/4,1/4,1/4,1/4], [1/4,1/4,1/4,1/4]])
sage: I=matrix(QQ,[[1,0,0,0], [0,1,0,0], [0,0,1,0], [0,0,0,1]])
sage: alpha=0.85
sage: G=alpha*H+(1-alpha)*S
sage: N=G-I
sage: 1200*N
[-1155.00000000000  45.0000000000000  1065.00000000000  555.000000000000]
[ 385.000000000000 -1155.00000000000  45.0000000000000  45.0000000000000]
[ 385.000000000000  555.000000000000 -1155.00000000000  555.000000000000]
[ 385.000000000000  555.000000000000  45.0000000000000 -1155.00000000000]
sage: M=matrix(QQ,[[-1155,45,1065,555], [385,-1155,45,45], [385,555,-1155,555],
....:      [385,555,45,-1155]])
sage: M.echelon_form()
[            1             0             0 -106613/58520]
[            0             1             0        -40/57]
[            0             0             1        -57/40]
[            0             0             0             0]
sage: v=vector([106613/58520,40/57,57/40,1])
sage: (v/v.norm()).n()
(0.696483066294572, 0.268280959381099, 0.544778023143244, 0.382300367118066)
(b)

Continue the Sage session to get this.

 sage: alpha=0.95
sage: G=alpha*H+(1-alpha)*S
sage: N=G-I
sage: 1200*N
[-1185.00000000000  15.0000000000000  1155.00000000000  585.000000000000]
[ 395.000000000000 -1185.00000000000  15.0000000000000  15.0000000000000]
[ 395.000000000000  585.000000000000 -1185.00000000000  585.000000000000]
[ 395.000000000000  585.000000000000  15.0000000000000 -1185.00000000000]
sage: M=matrix(QQ,[[-1185,15,1155,585], [395,-1185,15,15], [395,585,-1185,585],
....:      [395,585,15,-1185]])
sage: M.echelon_form()
[             1              0              0 -361677/186440]
[             0              1              0         -40/59]
[             0              0              1         -59/40]
[             0              0              0              0]
sage: v=vector([361677/186440,40/59,59/40,1])
sage: (v/v.norm()).n()
(0.713196892748114, 0.249250262646952, 0.542275102671275, 0.367644137404254)
(c)

Page p3 is important, but it passes its importance on to only one page, p1⁠. So that page receives a large boost.

1 Topic: Linear Recurrences source ↗

How many months until the number of Fibonacci rabbit pairs passes a thousand? Ten thousand? A million?

Answer

We use the formula

F(n)=15[(1+52)n−(1−52)n]

As observed earlier, (1+5)/2 is larger than one while (1+5)/2 has absolute value less than one.

sage: phi = (1+5^(0.5))/2
sage: psi = (1-5^(0.5))/2
sage: phi
1.61803398874989
sage: psi
-0.618033988749895

So the value of the expression is dominated by the first term. Solving 1000=(1/5)⋅((1+5)/2)n gives this.

sage: a = ln(1000*5^(0.5))/ln(phi)
sage: a
16.0271918385296
sage: psi^(17)
-0.000280033582072583

So by the seventeenth power, the second term does not contribute enough to change the roundoff. For the ten thousand and million calculations the situation is even more extreme.

sage: b = ln(10000*5^(0.5))/ln(phi)
sage: b
20.8121638053112
sage: c = ln(1000000*5^(0.5))/ln(phi)
sage: c
30.3821077388746

The answers in these cases are 21 and 31⁠.

2 Topic: Linear Recurrences source ↗

Solve each homogeneous linear recurrence relations.

(a)

f(n)=5f(n−1)−6f(n−2)

(b)

f(n)=4f(n−2)

(c)

f(n)= 5f(n−1)− 2f(n−2)− 8f(n−3)

Answer (3 parts)

(a)

We express the relation in matrix form.

(f(n)f(n−1))=(5−610)(f(n−1)f(n−2))

The characteristic equation of the matrix

|5−λ−61−λ|=λ2−5λ+6

has roots of 2 and 3⁠. Any function of the form f(n)=c12n+c23n satisfies the recurrence.

(b)

The matrix expression of the relation is

(f(n)f(n−1))=(0410)(f(n−1)f(n−2))

and the characteristic equation

|λ2−2|=(λ−2)(λ+2)

has the two roots 2 and −2⁠. Any function of the form f(n)=c12n+c2(−2)n satisfies this recurrence.

(c)

In matrix form the relation

(f(n)f(n−1)f(n−2))=(5−2−8100010)(f(n−1)f(n−2)f(n−3))

has a characteristic equation with roots −1⁠, 2⁠, and 4⁠. Any combination of the form c1(−1)n+c22n+c34n solves the recurrence.

3 Topic: Linear Recurrences source ↗

Give a formula for the relations of the prior exercise, with these initial conditions.

(a)

f(0)=1⁠, f(1)=1

(b)

f(0)=0⁠, f(1)=1

(c)

f(0)=1⁠, f(1)=1⁠, f(2)=3⁠.

Answer (3 parts)

(a)

The solution of the homogeneous recurrence is f(n)=c12n+c23n⁠. Substituting f(0)=1 and f(1)=1 gives this linear system.

c1+c2=12c1+3c2=1

By eye we see that c1=2 and c2=−1⁠.

(b)

The solution of the homogeneous recurrence is c12n+c2(−2)n⁠. The initial conditions give this linear system.

c1+c2=02c1−2c2=1

The solution is c1=1/4⁠, c2=−1/4⁠.

(c)

The homogeneous recurrence has the solution f(n)=c1(−1)n+c22n+c34n⁠. With the initial conditions we get this linear system.

c1+c2+c3=1−c1+2c2+4c3=1c1+4c2+16c3=3

Its solution is c1=1/3⁠, c2=2/3⁠, c3=0⁠.

4 Topic: Linear Recurrences source ↗

Check that the isomorphism given between S and ℝk is a linear map.

Answer

Fix a linear homogeneous recurrence of order k⁠. Let S be the set of functions f:ℕ→ℂ satisfying the recurrence. Consider the function Φ:S→ℂk given as here.

f⟼Φ(f(0)f(1)⋮f(k−1))

This shows linearity.

Φ(a⋅f1+b⋅f2)=(af1(0)+bf2(0)⋮af1(k−1)+bf2(k−1))=a(f1(0)⋮f1(k−1))+b(f2(0)⋮f2(k−1))=a⋅Φ(f1)+b⋅Φ(f2)

5 Topic: Linear Recurrences source ↗

Show that the characteristic equation of the matrix is as stated, that is, is the polynomial associated with the relation. (Hint: expanding down the final column and using induction will work.)

Answer

We use the hint to prove this.

0=|an−1−λan−2an−3…an−k+1an−k1−λ0…0001−λ001⋮⋮⋱⋮000…1−λ|=(−1)k−1(−λk+an−1λk−1+an−2λk−2+⋯+an−k+1λ+an−k)=±(−λk+an−1λk−1+an−2λk−2+⋯+an−k+1λ+an−k)

The base step is trivial. For the inductive step, expanding down the final column gives two nonzero terms.

(−1)k−1an−k⋅1−λ⋅|an−1−λan−2an−3…an−k+11−λ0…001−λ001⋮⋮⋱000…1|

(The matrix is square so the sign in front of −λ is −1even⁠). Application of the inductive hypothesis gives the desired result.

=(−1)k−1an−k⋅1−λ⋅(−1)k−2(−λk−1+an−1λk−2+an−2λk−3+⋯+an−k+1λ0)

6 Topic: Linear Recurrences source ↗

Given a homogeneous linear recurrence relation f(n)= anf(n−1)+ ⋯+ an−kf(n−k)⁠, let r1⁠, …, rk be the roots of the associated polynomial. Prove that each function fri(n)=rkn satisfies the recurrence (without initial conditions).

Answer

This is a straightforward induction on n⁠.

7 Topic: Linear Recurrences source ↗

(This refers to the value T(64)= 18,446,744,073,709,551,615 given in the computer code below.) Transferring one disk per second, how many years would it take the priests at the Tower of Hanoi to finish the job?

Answer

Sage says that we are safe.

sage: T64 = 18446744073709551615
sage: T64_days = T64/(60*60*24)
sage: T64_days
1229782938247303441/5760
sage: T64_years = T64_days/365.25
sage: T64_years
5.84542046090626e11        
sage: age_of_universe = 13.8e9
sage: T64_years/age_of_universe
42.3581192819294

1 Topic: Coupled Oscillators source ↗

Use the formula for the cosine of a sum to give an even more general formula for simple harmonic motion.

Answer

The angle sum formula for the cosine function is cos⁡(α+β)= cos(α)cos⁡(β)− sin⁡(α)sin⁡(β)⁠. Expand Acos⁡(ωt+ϕ) to A⋅[cos(ωt)cos⁡(ϕ)− sin⁡(ωt)sin⁡(ϕ)]⁠. Then cos⁡(ϕ) and sin⁡(ϕ) do not vary with t so we get the general solution x(t)= Bcos⁡(ωt)+ Csin⁡(ωt)⁠.

2 Topic: Coupled Oscillators source ↗

Find the eigenvectors associated with the eigenvalues.

Answer

We will do the first root; the second is similar. We have this equation.

(ω02−ω2ϵ/2mϵ/2Iω02−ω2)(A1A2)=(00)

Plug in the first root ω2=ω02+ϵ/2mI⁠. The two equations are redundant so we just consider the first.

−(ϵ/2mI)⋅A1+(ϵ/2m)⋅A2=0

That’s a line through the origin in the plane, so to specify it we need only find the ratio between the first and second variables.

A2A1=(ϵ/2mI)(ϵ/2m)=mI

So this is the space of eigenvectors associated with the first eigenvalue.

{(A1m/I⋅A1)|A1∈ℂ}

3 Topic: Coupled Oscillators source ↗

Find the values of ω in the case where x(t)=A1cos⁡(ωt+ϕ) and θ(t)=A2cos⁡(ωt+(ϕ+π))⁠.

Answer

Observe that cos⁡(ωt+(ϕ+π))= − cos⁡(ωt+ϕ)⁠. Take the derivatives

dx(t)dt=−ω⋅sin⁡(ωt+ϕ)d2x(t)dt2=−ω2⋅cos⁡(ωt+ϕ)

and

dθ(t)dt=ω⋅sin⁡(ωt+ϕ)d2θ(t)dt2=ω2⋅cos⁡(ωt+ϕ)

and plug into the equations of motion (∗∗⁠).

mA1⋅(−ω2cos⁡(ωt+ϕ))+kA1⋅(cos⁡(ωt+ϕ))+ϵ2A2⋅(−cos⁡(ωt+ϕ))=0IA2⋅(ω2cos⁡(ωt+ϕ))+κA2⋅(−cos⁡(ωt+ϕ))+ϵ2A1⋅(cos⁡(ωt+ϕ))=0

Factor out cos⁡(ωt+ϕ)

(m(−ω2))⋅A1+kA1−ϵ2⋅A2=0(Iω2)⋅A2−κA2+ϵ2⋅A1=0

and divide through by m and I⁠.

(km−ω2)⋅A1−ϵ2m⋅A2=0ϵ2I⋅A1−(κI−ω2)⋅A2=0

We are assuming that k/m=ωx2 and κ/I=ωθ2 are equal, and writing ω02 for that number. Make the substitution and restate it as a matrix equation.

(ω02−ω2−ϵ/2mϵ/2Iω02−ω2)(A1A2)=(00)

We want the frequencies ω for which this system has a nontrivial solution.

|ω02−ω2−ϵ/2mϵ/2Iω2−ω02|=0

This is the same determinant as we did in the Topic body except that the second column is multiplied by −1⁠. Multiplying a row or column by a scalar multiplies the entire determinant by that scalar. So this determinant is the negative of the one in the Topic body. But we are setting it to zero so that doesn’t matter. The roots are the same as in the Topic body.

4 Topic: Coupled Oscillators source ↗

Build a Wilberforce pendulum out of a Slinky Jr and a soup can. You can drill holes in the can for bolts, either two or four of them, that you can use to adjust the moment of inertia of the can so the periods of vertical and rotational motion coincide.

Answer

See [Mewes].

1 Topic: Inner Product source ↗

Verify that ℝn along with dot product is an inner product space.

Answer

The opening paragraph of this Topic notes that dot product is linear in both inputs, and so the first condition is satisfied. Condition (2) concerns conjugation, which has no effect on reals, and so this condition is trivially satisfied. Condition (3) requires that the length (the ordinary Euclidean length in n⁠-dimensional real space) be zero if and only if the vector is 0→⁠. If the vector has components v1⁠, … vn then v12+⋯vn2=0 if and only if all the components are zero.

2 Topic: Inner Product source ↗

Show that Hermitian inner product on the two-dimensional vector space ℂ2 satisfies the three conditions in the definition of inner product.

Answer

The development in the Topic body proves (1) and (2). For (3), suppose that two members v→,w→ of ℂn have i⁠-th components vi and wi⁠, so that ⟨v→,w→⟩=v1w1+⋯+vnwn⁠. Then ⟨v→,v→⟩=v1v1+⋯+vnvn equals |v→1|2+⋯|v→n|2⁠, which clearly equals 0 if and only if each component is 0⁠.

3 Topic: Inner Product source ↗

Show that this operation is an inner product on the real vector space ℳ2×2⁠.

⟨(abcd),(efgh)⟩=ae+bf+cg+dh

Answer

Condition (1), linearity in the first argument, is easy.

⟨(a1b1c1d1)+(a2b2c2d2),(efgh)⟩=(a1+a2)e+(b1+b2)f+(c1+c2)g+(d1+d2)h=(a1e+b1f+c1g+d1h)+(a2e+b2f+c2g+d2h)=⟨(a1b1c1d1),(efgh)⟩+⟨(a2b2c2d2),(efgh)⟩

Scalar multiplication is just as routine. Condition (2) concerns conjugation, which has no effect on reals and so this condition is trivially satisfied. For (3), of course a2+b2+c2+d2≥0⁠, with equality if and only if each entry is zero.

4 Topic: Inner Product source ↗

Consider the real vector space of quadratic polynomials 𝒫2⁠. Show that this is an inner product operation: ⟨a2x2+ a1x+ a0, b2x2+ b1x+ b0⟩= a2b2+ a1b1+ a0b0⁠. Consider the resulting norm. What is ‖x2−1‖⁠?

Answer

Condition (1) is that it is linear in the first input. Verifying that this operation respects addition is routine

⟨(a2x2+a1x+a0)+(c2x2+c1x+c0),b2x2+b1x+b0⟩=(a2+c2)b2+(a1+c1)b1+(a0+c0)b0=(a2b2+a1b1+a0b0)+(c2b2+c1b1+c0b0)=⟨(a2x2+a1x+a0),b2x2+b1x+b0⟩+⟨(c2x2+c1x+c0),b2x2+b1x+b0⟩

as is verifying that it respects scalar multiplication. Because this is a real vector space, and conjugation has no effect on real number scalars, it satisfies condition (2). For condition (3), let v→=a2x2+a1x+a0⁠. Then ⟨v→,v→⟩=a22+a12+a02 is clearly non-negative, and is equal to 0 if and only if v→ is the zero polynomial.

As to the norm of x2−1⁠, we have this.

‖x2−1‖=⟨x2−1,x2−1⟩=12+02+(−1)2=2

5 Topic: Inner Product source ↗

Verify that the set of real valued functions that are continuous on the closed interval [0..1] is a vector space.

Answer

The ten conditions in the definition of a vector space are consequences of results from Calculus. For example, the first condition is closure under addition and we have that the sum of two continuous functions is a continuous function in Calculus. (Although, often a student does not see a completely precise proof until a later class.) Similarly, addition of two functions is commutative. Of course, the zero vector is the function f(x)=0⁠, which is continuous.

6 Topic: Inner Product source ↗

Show that an attempted definition of ⟨v→,w→⟩ on ℂn that is linear in both inputs loses that ⟨v→,v→⟩ gives the length.

Answer

Suppose that ⟨v→,v→⟩ is greater than or equal to 0 for all complex vectors v→⁠. Then linearity in both inputs gives 0≤⟨iv→, iv→⟩= i⟨v→, iv→⟩= i2⟨v→, v→⟩≤0⁠.

End of the book.

Exercises and answers are adapted from Linear Algebra by Jim Hefferon, which is free to read and download, and are used under the CC BY-SA 2.5 licence. They are taken from the book’s LaTeX source and set for the web, so they may differ slightly from the printed page. This adaptation is released under the same licence and is neither affiliated with nor endorsed by Jim Hefferon.