Books
Theme
Text size
Width

[ and ] go to the previous and next page.

Rank–nullity measures one matrix at a time. This section compares the ranks of matrices built from one another. Multiplying by an invertible matrix does not change the rank, and multiplying by an arbitrary matrix can only lower it:

rank(𝐴𝐵)≤min(rank𝐴,rank𝐵).

We prove this in three ways, by comparing images, by comparing kernels, and by restricting 𝐴 to the image of 𝐵, and the last proof says exactly how much rank is lost. The transpose 𝐴𝖳 has the same rank as 𝐴, although its columns are the rows of 𝐴 and live in a different space. This equality, row rank equals column rank, is the main theorem of the section. It completes the count of dimensions of the four subspaces attached to a matrix. Finally, matrices of rank one are exactly the products 𝐮𝐯𝖳, rank is subadditive, and the column–row factorization 𝐴=𝐶𝑅 writes a matrix of rank 𝑟 as a sum of 𝑟 matrices of rank one, which is the fewest possible.

Throughout, 𝐹 is a field. For 𝐴∈𝐹𝑚×𝑛 we write im𝐴=im𝑇𝐴⊆𝐹𝑚 and ker𝐴=ker𝑇𝐴⊆𝐹𝑛 for the column space and the null space (Definition 2.6.17). So rank𝐴=dimim𝐴 (Theorem 4.2.7) and dimker𝐴=𝑛−rank𝐴 (Theorem 4.2.11).

4.3.1Invertible factors

We already know that rank(𝐺𝐴)=rank𝐴 for invertible 𝐺, because rref(𝐺𝐴)=rref(𝐴) (Proposition 3.5.2(c)), and that rank(𝐴𝐻)=rank𝐴 for invertible 𝐻, from the uniqueness of the rank normal form (Exercise 3.5.41(c)). The second proof was roundabout, because row reduction cannot see column operations. Now that rank is a dimension, both statements have short proofs. They are also the reason the rank of a linear map can be read off any of its matrices: changing bases multiplies [𝑇]𝒞︀ℬ︀ on the left and on the right by invertible matrices (Section 4.5), and that cannot change the rank.

Intuition (Outputs and inputs). In 𝐺𝐴 the matrix 𝐺 acts after 𝐴, on the outputs 𝐴𝐱. An invertible 𝐺 is injective, so it never sends a nonzero output to 𝟎: 𝐺𝐴𝐱=𝟎 only when 𝐴𝐱=𝟎, and ker(𝐺𝐴)=ker𝐴. But 𝐺 moves the outputs, so im(𝐺𝐴), the set of vectors 𝐺𝐲 with 𝐲∈im𝐴, is in general a different subspace of 𝐹𝑚. It has the same dimension, since an injective map preserves dimension.

In 𝐴𝐻 the matrix 𝐻 acts before 𝐴, on the inputs. An invertible 𝐻 is surjective, so 𝐴𝐻 still sees every input, and im(𝐴𝐻)=im𝐴. But the kernel moves: 𝐴𝐻𝐱=𝟎 exactly when 𝐻𝐱∈ker𝐴. So each of 𝐺𝐴 and 𝐴𝐻 keeps one of the two subspaces and moves the other by an invertible map.

Proposition 4.3.1 (Invertible factors preserve rank). Let 𝐴∈𝐹𝑚×𝑛, and let 𝐺∈𝐹𝑚×𝑚 and 𝐻∈𝐹𝑛×𝑛 be invertible.
  1. (a)ker(𝐺𝐴)=ker𝐴, and rank(𝐺𝐴)=rank𝐴.
  2. (b)im(𝐴𝐻)=im𝐴, and rank(𝐴𝐻)=rank𝐴.
In particular, rank(𝐺𝐴𝐻)=rank𝐴.

The first proof is the lecture’s. It compares images, using the restriction identity im(𝑇𝑆)=im(𝑇↾im𝑆).

Proof (by restriction, comparing images).routine

Uses. Theorem 2.4.2 (The matrix of a composite), Proposition 2.5.4 (Left inverses give uniqueness), Proposition 2.5.5 (Right inverses give existence), Proposition 4.2.3 (Properties of restriction), Proposition 4.2.6 (Injective maps preserve dimension), Proposition 2.6.15 (Kernels and images of composites), Theorem 4.2.7 (The image of 𝑇𝐴)

Techniques. Restrict to a subspace (Technique 4.2.5)

Hint 1. Write 𝑇𝐺𝐴=𝑇𝐺∘𝑇𝐴. By Proposition 4.2.3(b), im(𝐺𝐴) is the image of 𝑇𝐺 restricted to im𝐴. What do you know about 𝑇𝐺?
Hint 2. 𝑇𝐺 is injective, so its restriction to im𝐴 is an injective linear map onto im(𝐺𝐴). For 𝐴𝐻, what is im𝑇𝐻?
Sketch. im(𝐺𝐴)=im(𝑇𝐺↾im𝐴), and 𝑇𝐺 is injective because 𝑇𝐺−1 is a left inverse, so the restriction carries im𝐴 injectively onto im(𝐺𝐴) and the dimensions agree. On the other side, im(𝐴𝐻)=im(𝑇𝐴↾im𝐻)=im(𝑇𝐴↾𝐹𝑛)=im𝐴, because 𝑇𝐻 is surjective. The kernel statement in (a) is Proposition 2.6.15(a).

By Theorem 2.4.2, 𝑇𝐺𝐴=𝑇𝐺∘𝑇𝐴 and 𝑇𝐴𝐻=𝑇𝐴∘𝑇𝐻. Since 𝐺−1𝐺=𝐼𝑚 and 𝐻𝐻−1=𝐼𝑛, the map 𝑇𝐺 is injective (Proposition 2.5.4) and 𝑇𝐻 is surjective (Proposition 2.5.5).

(a) By Proposition 4.2.3(b), im(𝐺𝐴)=im(𝑇𝐺↾im𝐴). The restriction 𝑇𝐺↾im𝐴 is linear and injective (Proposition 4.2.3(a), (c)), so Proposition 4.2.6(a), applied to it, gives dimim(𝐺𝐴)=dimim𝐴. That is, rank(𝐺𝐴)=rank𝐴 (Theorem 4.2.7). Since 𝑇𝐺 is injective, ker(𝐺𝐴)=ker𝐴 by Proposition 2.6.15(a).

(b) By Proposition 4.2.3(b), im(𝐴𝐻)=im(𝑇𝐴↾im𝐻). Since im𝐻=𝐹𝑛, this restriction is 𝑇𝐴 itself, so im(𝐴𝐻)=im𝐴, and in particular rank(𝐴𝐻)=rank𝐴.

Finally, rank(𝐺𝐴𝐻)=rank(𝐺(𝐴𝐻))=rank(𝐴𝐻)=rank𝐴 by (a) and (b). □

The lecture remarks that the same statements follow from kernels and rank–nullity. Here is that proof.

Proof (by kernels and rank–nullity).routine

Uses. Proposition 2.6.15 (Kernels and images of composites), Theorem 4.2.11 (Rank–nullity), Proposition 4.2.3 (Properties of restriction), Proposition 4.2.6 (Injective maps preserve dimension), Proposition 4.1.21 (Dimensions of subspaces)

Techniques. Dimension counting (Technique 4.2.24)

Hint 1. 𝐺𝐴 and 𝐴𝐻 have 𝑛 columns, as 𝐴 does. By rank–nullity it suffices to compare dimker(𝐺𝐴) and dimker(𝐴𝐻) with dimker𝐴.
Hint 2. 𝐺 can be cancelled from 𝐺𝐴𝐱=𝟎. For 𝐴𝐻, show that 𝐱∈ker(𝐴𝐻) if and only if 𝐻𝐱∈ker𝐴, and use 𝑇𝐻 to move one kernel onto the other.
Sketch. If 𝐺𝐴𝐱=𝟎, then 𝐴𝐱=𝐺−1𝐺𝐴𝐱=𝟎, so ker(𝐺𝐴)=ker𝐴, and rank–nullity gives rank(𝐺𝐴)=𝑛−dimker𝐴=rank𝐴. Next, 𝐱∈ker(𝐴𝐻) if and only if 𝐻𝐱∈ker𝐴, so 𝑇𝐻 restricts to an injective linear map from ker(𝐴𝐻) onto ker𝐴. The kernels have the same dimension, so rank(𝐴𝐻)=rank𝐴 by rank–nullity, and then im(𝐴𝐻)⊆im𝐴 with equal dimensions forces im(𝐴𝐻)=im𝐴.

(a) By Proposition 2.6.15(a), ker𝐴⊆ker(𝐺𝐴). Conversely, if 𝐺𝐴𝐱=𝟎, then 𝐴𝐱=𝐺−1(𝐺𝐴𝐱)=𝟎. So ker(𝐺𝐴)=ker𝐴. The matrices 𝐺𝐴 and 𝐴 both have 𝑛 columns, so by Theorem 4.2.11(a),

rank(𝐺𝐴)=𝑛−dimker(𝐺𝐴)=𝑛−dimker𝐴=rank𝐴.

(b) For 𝐱∈𝐹𝑛, we have 𝐴𝐻𝐱=𝟎 if and only if 𝐻𝐱∈ker𝐴. So 𝑇𝐻 maps ker(𝐴𝐻) into ker𝐴. Since 𝐻−1𝐻=𝐼𝑛, 𝑇𝐻 is injective, so its restriction to ker(𝐴𝐻) is linear and injective (Proposition 4.2.3). It maps onto ker𝐴: if 𝐲∈ker𝐴, then 𝐱=𝐻−1𝐲 satisfies 𝐴𝐻𝐱=𝐴𝐲=𝟎 and 𝐻𝐱=𝐲. By Proposition 4.2.6(a), dimker(𝐴𝐻)=dimker𝐴, and since 𝐴𝐻 has 𝑛 columns, Theorem 4.2.11(a) gives

rank(𝐴𝐻)=𝑛−dimker(𝐴𝐻)=𝑛−dimker𝐴=rank𝐴.

Finally, im(𝐴𝐻)⊆im𝐴 by Proposition 2.6.15(b), and these subspaces of 𝐹𝑚 have the same dimension, so they are equal (Proposition 4.1.21(b)).

As in the first proof, rank(𝐺𝐴𝐻)=rank𝐴 follows from (a) and (b). □

The two proofs split the work in opposite ways. In the first, the left factor needs an argument (an injective map preserves dimension) and the right factor is free (a surjective map on the inputs changes nothing). In the second, the left factor is free (the kernel does not change) and the right factor needs the argument. The same split reappears in Section 4.3.2, where an arbitrary product is compared with its factors.

By Exercise 3.5.41(d), two matrices 𝐴,𝐵∈𝐹𝑚×𝑛 satisfy 𝐵=𝐺𝐴𝐻 for some invertible 𝐺 and 𝐻 if and only if they have the same rank. So the rank is the only information about 𝐴 that row and column operations together preserve. In particular, similar matrices 𝐵=𝑃𝐴𝑃−1 have the same rank, a fact we return to in Section 4.5.

Exercise 4.3.2.routine Let 𝐴=(1224), 𝐺=(10−21), and 𝐻=(1−201), over ℚ. Compute 𝐺𝐴 and 𝐴𝐻, and find the kernel and the image of each of 𝐴, 𝐺𝐴, and 𝐴𝐻. Which subspaces agree, and which do not?

Solution.

Uses. Proposition 4.3.1 (Invertible factors preserve rank), Proposition 2.6.18 (The column space is spanned by the columns)

Hint 1. 𝐺 adds −2 times row 1 to row 2, and 𝐻 adds −2 times column 1 to column 2.
Hint 2. All three matrices have rank 1, so each kernel and each image is a line. Find a spanning vector for each.
Sketch. 𝐺𝐴=(1200) and 𝐴𝐻=(1020). The kernels are span((−2,1)), span((−2,1)), span(𝐞2), and the images are span((1,2)), span(𝐞1), span((1,2)), for 𝐴, 𝐺𝐴, 𝐴𝐻 respectively. So 𝐺𝐴 keeps the kernel of 𝐴 and moves the image, and 𝐴𝐻 keeps the image and moves the kernel.

Multiplying out,

𝐺𝐴=(122−44−8)=(1200),𝐴𝐻=(1−2+22−4+4)=(1020).

The columns of 𝐴 are (1,2) and 2(1,2), so im𝐴=span((1,2)) (Proposition 2.6.18), and 𝐴𝐱=𝟎 means 𝑥1+2𝑥2=0, so ker𝐴=span((−2,1)). For 𝐺𝐴, the columns are 𝐞1 and 2𝐞1, so im(𝐺𝐴)=span(𝐞1), and 𝐺𝐴𝐱=𝟎 again means 𝑥1+2𝑥2=0, so ker(𝐺𝐴)=span((−2,1)). For 𝐴𝐻, the columns are (1,2) and 𝟎, so im(𝐴𝐻)=span((1,2)), and 𝐴𝐻𝐱=𝑥1(1,2) vanishes exactly when 𝑥1=0, so ker(𝐴𝐻)=span(𝐞2).

So ker(𝐺𝐴)=ker𝐴 but im(𝐺𝐴)≠im𝐴, and im(𝐴𝐻)=im𝐴 but ker(𝐴𝐻)≠ker𝐴, as Proposition 4.3.1 and Intuition (Outputs and inputs) predict. All three matrices have rank 1. Note that ker(𝐴𝐻)=𝐻−1(ker𝐴): indeed 𝐻−1=(1201) and 𝐻−1(−2,1)=(0,1).

4.3.2The rank of a product

Intuition (A pipeline is no wider than its narrowest stage). For 𝐱∈𝐹𝑝, 𝐴𝐵𝐱=𝐴(𝐵𝐱): first 𝐵, then 𝐴. Every output of 𝐴𝐵 is an output of 𝐴, so im(𝐴𝐵)⊆im𝐴; the product cannot reach a direction that 𝐴 cannot reach. On the other hand, 𝐴 only ever receives the vectors 𝐵𝐱, which fill a space of dimension rank𝐵, and a linear map does not raise dimension. So

rank(𝐴𝐵)≤rank𝐴andrank(𝐴𝐵)≤rank𝐵.

The first bound is about where the outputs land and is immediate. The second is about how many independent directions enter 𝐴, and it needs a dimension argument.

Theorem 4.3.3 (Rank of a product).
  1. (a)Let 𝐴∈𝐹𝑚×𝑛 and 𝐵∈𝐹𝑛×𝑝. Then im(𝐴𝐵)⊆im𝐴 and ker𝐵⊆ker(𝐴𝐵), and

    rank(𝐴𝐵)≤min(rank𝐴,rank𝐵).
  2. (b)More generally, let 𝑈 and 𝑉 be finite-dimensional, and let 𝑆∈L(𝑈,𝑉) and 𝑇∈L(𝑉,𝑊). Then

    rank(𝑇𝑆)≤min(rank𝑇,rank𝑆).

Part (a) is part (b) for 𝑆=𝑇𝐵 and 𝑇=𝑇𝐴, since 𝑇𝐴𝐵=𝑇𝐴∘𝑇𝐵 (Theorem 2.4.2) and rank𝐴=rank𝑇𝐴 (Definition 4.2.14); the two containments in (a) are Proposition 2.6.15. So we prove (b). In both proofs below, the bound by rank𝑇 comes from im(𝑇𝑆)⊆im𝑇. They differ in how they prove the bound by rank𝑆.

Proof (comparing images).routine

Uses. Proposition 2.6.15 (Kernels and images of composites), Theorem 4.2.11 (Rank–nullity), Proposition 4.1.21 (Dimensions of subspaces), Proposition 4.2.3 (Properties of restriction), Proposition 2.1.16 (Linear maps preserve spans and linear relations), Corollary 4.1.11 (Spanning lists are long)

Techniques. Restrict to a subspace (Technique 4.2.5)

Hint 1. im(𝑇𝑆)⊆im𝑇, which gives one inequality at once.
Hint 2. For the other, im(𝑇𝑆)=im(𝑇↾im𝑆) consists of the values of 𝑇 on im𝑆. Take a basis of im𝑆 and push it forward with 𝑇.
Hint 3. A space spanned by 𝑟 vectors has dimension at most 𝑟.
Sketch. im(𝑇𝑆)⊆im𝑇 gives rank(𝑇𝑆)≤rank𝑇. If 𝐬1,…,𝐬𝑟 is a basis of im𝑆, then im(𝑇𝑆)=im(𝑇↾im𝑆) is spanned by 𝑇(𝐬1),…,𝑇(𝐬𝑟), so its dimension is at most 𝑟=rank𝑆.

Since 𝑈 and 𝑉 are finite-dimensional, im𝑆, im𝑇, and im(𝑇𝑆) are finite-dimensional (Theorem 4.2.11(b)), so all three ranks are defined.

The bound by rank𝑇. By Proposition 2.6.15(b), im(𝑇𝑆)⊆im𝑇, so rank(𝑇𝑆)≤rank𝑇 by Proposition 4.1.21(a).

The bound by rank𝑆. Let 𝑟=rank𝑆 and let 𝐬1,…,𝐬𝑟 be a basis of im𝑆 (the empty list if 𝑟=0). By Proposition 4.2.3(b) and Proposition 2.1.16(a),

im(𝑇𝑆)=im(𝑇↾im𝑆)=span(𝑇(𝐬1),…,𝑇(𝐬𝑟)).

So im(𝑇𝑆) is spanned by a list of length 𝑟, and dimim(𝑇𝑆)≤𝑟 by Corollary 4.1.11(a). That is, rank(𝑇𝑆)≤rank𝑆. □

The second proof moves the comparison to the domain 𝑈, where 𝑆 and 𝑇𝑆 have a common domain and rank–nullity applies to both.

Proof (comparing kernels, with rank–nullity).routine

Uses. Proposition 2.6.15 (Kernels and images of composites), Proposition 4.1.21 (Dimensions of subspaces), Theorem 4.2.11 (Rank–nullity)

Techniques. Dimension counting (Technique 4.2.24)

Hint 1. The bound by rank𝑇 is the image comparison im(𝑇𝑆)⊆im𝑇. For the bound by rank𝑆, compare kernels instead: which of ker𝑆 and ker(𝑇𝑆) contains the other?
Hint 2. ker𝑆⊆ker(𝑇𝑆), both subspaces of 𝑈. Apply rank–nullity to 𝑆 and to 𝑇𝑆.
Sketch. rank(𝑇𝑆)≤rank𝑇 because im(𝑇𝑆)⊆im𝑇. Since ker𝑆⊆ker(𝑇𝑆), rank–nullity on 𝑈 gives rank(𝑇𝑆)=dim𝑈−dimker(𝑇𝑆)≤dim𝑈−dimker𝑆=rank𝑆.

As in the first proof, im(𝑇𝑆)⊆im𝑇 (Proposition 2.6.15(b)), so rank(𝑇𝑆)≤rank𝑇 (Proposition 4.1.21(a)).

By Proposition 2.6.15(a), ker𝑆⊆ker(𝑇𝑆). Both are subspaces of the finite-dimensional space 𝑈, so dimker𝑆≤dimker(𝑇𝑆) (Proposition 4.1.21(a)). By Theorem 4.2.11(b), applied to 𝑆 and to 𝑇𝑆, which are both defined on 𝑈,

rank(𝑇𝑆)=dim𝑈−dimker(𝑇𝑆)≤dim𝑈−dimker𝑆=rank𝑆,

as desired. □

Both inequalities can be strict, even at the same time.

Example 4.3.4 (A product of rank zero). Let 𝑁=(0100), so 𝑁𝐞1=𝟎 and 𝑁𝐞2=𝐞1. Then rank𝑁=1, but 𝑁2=0 has rank 0<min(rank𝑁,rank𝑁). In 𝑁𝑁, the left factor kills exactly what the right factor produces: im𝑁=span(𝐞1)=ker𝑁. The order of the factors matters too: with 𝑀=(1000), 𝑀𝑁=𝑁 has rank 1 but 𝑁𝑀=0, so rank(𝐴𝐵) and rank(𝐵𝐴) can differ even for square matrices.

What a product loses

The example shows where the rank goes. 𝐴𝐵 applies 𝐴 to the vectors of im𝐵, and 𝐴 kills those that lie in ker𝐴. Restricting 𝐴 to im𝐵 makes this exact. This is the lecture’s restriction view of a composite, and it gives the third proof of the bound by rank𝐵.

Proposition 4.3.5 (Rank of a product, exactly). Let 𝑆∈L(𝑈,𝑉) and 𝑇∈L(𝑉,𝑊), with 𝑈 finite-dimensional.

  1. (a)rank(𝑇𝑆)=rank𝑆−dim(im𝑆∩ker𝑇).
  2. (b)rank(𝑇𝑆)=rank𝑆 if and only if im𝑆∩ker𝑇={𝟎}. If 𝑉 is also finite-dimensional, then rank(𝑇𝑆)=rank𝑇 if and only if im𝑆+ker𝑇=𝑉.

For matrices 𝐴∈𝐹𝑚×𝑛 and 𝐵∈𝐹𝑛×𝑝, (a) reads

rank(𝐴𝐵)=rank𝐵−dim(im𝐵∩ker𝐴).
(4.3.1)

Proof (by restriction).standard

Uses. Theorem 4.2.11 (Rank–nullity), Proposition 4.2.3 (Properties of restriction), Exercise 4.2.44, Proposition 2.6.15 (Kernels and images of composites), Proposition 4.1.21 (Dimensions of subspaces), Exercise 2.6.39, Theorem 2.4.2 (The matrix of a composite)

Techniques. Restrict to a subspace (Technique 4.2.5), Dimension counting (Technique 4.2.24)

Hint 1. im(𝑇𝑆)=im(𝑇↾im𝑆). What is the kernel of 𝑇↾im𝑆?
Hint 2. Apply rank–nullity to 𝑇↾im𝑆, whose domain is the finite-dimensional space im𝑆.
Hint 3. For the last statement, im(𝑇𝑆)⊆im𝑇; when do two finite-dimensional nested subspaces have the same dimension? Then use Exercise 2.6.39(b).
Sketch. The restriction 𝑅=𝑇↾im𝑆 has domain im𝑆, of dimension rank𝑆, kernel im𝑆∩ker𝑇, and image im(𝑇𝑆). Rank–nullity for 𝑅 is (a), and (a) gives the first half of (b). For the second half, im(𝑇𝑆)⊆im𝑇, so the ranks agree exactly when the images agree, which by Exercise 2.6.39(b) means im𝑆+ker𝑇=𝑉.

(a) Since 𝑈 is finite-dimensional, so is im𝑆, and dimim𝑆=rank𝑆 (Theorem 4.2.11(b)). Let 𝑅=𝑇↾im𝑆:im𝑆→𝑊. By Proposition 4.2.3, 𝑅 is linear, ker𝑅=im𝑆∩ker𝑇, and im𝑅=im(𝑇𝑆). Rank–nullity for 𝑅 (Theorem 4.2.11(b)) gives

rank𝑆=dimim𝑆=dim(im𝑆∩ker𝑇)+dimim(𝑇𝑆),

which is (a). (This is Exercise 4.2.44(a) for the subspace im𝑆.)

(b) By (a), rank(𝑇𝑆)=rank𝑆 if and only if dim(im𝑆∩ker𝑇)=0, that is, im𝑆∩ker𝑇={𝟎}. Now let 𝑉 be finite-dimensional, so that im𝑇 is finite-dimensional. Since im(𝑇𝑆)⊆im𝑇 (Proposition 2.6.15(b)), Proposition 4.1.21(b) shows that rank(𝑇𝑆)=rank𝑇 if and only if im(𝑇𝑆)=im𝑇, and by Exercise 2.6.39(b) this holds if and only if im𝑆+ker𝑇=𝑉.

The matrix statement is (a) for 𝑆=𝑇𝐵 and 𝑇=𝑇𝐴 (Theorem 2.4.2). □

Intuition (Three proofs of rank(𝑇𝑆)≤rank𝑆). The image proof works in the middle space 𝑉: it pushes a basis of im𝑆 forward with 𝑇 and observes that at most rank𝑆 vectors come out. The kernel proof works in the domain 𝑈: whatever 𝑆 kills, 𝑇𝑆 kills too, and rank–nullity converts a larger kernel into a smaller image. The restriction proof is the sharpest, because it computes the loss instead of bounding it: 𝑇𝑆 loses exactly the directions that 𝑆 produces and 𝑇 destroys, im𝑆∩ker𝑇. It answers the question of when the inequality is an equality, and it also gives lower bounds for rank(𝑇𝑆), which are the subject of Exercise 4.3.30 and Exercise 4.3.31. When all you need is the inequality, the image proof is the one to remember.

Example 4.3.6 (Powers of a square matrix). Let 𝐴∈𝐹𝑛×𝑛. By (4.3.1) with 𝐵=𝐴,

rank(𝐴2)=rank𝐴−dim(im𝐴∩ker𝐴).

So rank(𝐴2)=rank𝐴 if and only if im𝐴∩ker𝐴={𝟎}. For the projection 𝑃=(1000), im𝑃=span(𝐞1) and ker𝑃=span(𝐞2) meet only in 𝟎, and indeed 𝑃2=𝑃. For 𝑁 of Example 4.3.4, im𝑁=ker𝑁, and rank(𝑁2)=1−1=0.

With 𝐵=𝐴𝑘 instead, the formula gives rank(𝐴𝑘+1)=rank(𝐴𝑘)−dim(im𝐴𝑘∩ker𝐴). The subspaces im𝐴𝑘 shrink as 𝑘 grows, so the drops rank(𝐴𝑘)−rank(𝐴𝑘+1) can only decrease. This is Exercise 4.2.52 in the language of images, and by Exercise 4.2.51 the ranks stop decreasing by 𝑘=𝑛.

Exercise 4.3.7.routine Let 𝐴∈𝐹𝑚×𝑛 and 𝐵∈𝐹𝑛×𝑚 with 𝑚<𝑛.

  1. (a)Show that 𝐵𝐴∈𝐹𝑛×𝑛 is not invertible.
  2. (b)Give an example with 𝐴𝐵=𝐼𝑚.
  3. (c)Show that if 𝐴𝐵=𝐼𝑚, then rank(𝐵𝐴)=𝑚.

Solution.

Uses. Theorem 4.3.3 (Rank of a product), Proposition 3.5.2 (First properties of rank), Corollary 3.5.25 (The invertible matrix theorem), Theorem 2.4.34 (Block multiplication)

Hint 1. rank(𝐵𝐴)≤rank𝐴, and 𝐴 has only 𝑚 rows.
Hint 2. For (b), try block matrices built from 𝐼𝑚 and zeros.
Hint 3. For (c), sandwich 𝐼𝑚=(𝐴𝐵)(𝐴𝐵)=𝐴(𝐵𝐴)𝐵.
Sketch. (a) rank(𝐵𝐴)≤rank𝐴≤𝑚<𝑛, so 𝐵𝐴 is not invertible. (b) 𝐴=(𝐼𝑚0), 𝐵=(𝐼𝑚0). (c) 𝑚=rank(𝐼𝑚)=rank(𝐴(𝐵𝐴)𝐵)≤rank(𝐵𝐴)≤𝑚.

(a) By Theorem 4.3.3 and Proposition 3.5.2(a), rank(𝐵𝐴)≤rank𝐴≤𝑚<𝑛. An invertible 𝑛×𝑛 matrix has rank 𝑛 (Corollary 3.5.25), so 𝐵𝐴 is not invertible.

(b) Let 𝐴=(𝐼𝑚0)∈𝐹𝑚×𝑛 and 𝐵=(𝐼𝑚0)∈𝐹𝑛×𝑚, where the zero blocks have 𝑛−𝑚 columns and 𝑛−𝑚 rows. By Theorem 2.4.34, 𝐴𝐵=𝐼𝑚𝐼𝑚+0⋅0=𝐼𝑚. (Here 𝐵𝐴=(𝐼𝑚000), of rank 𝑚.)

(c) If 𝐴𝐵=𝐼𝑚, then 𝐼𝑚=(𝐴𝐵)(𝐴𝐵)=𝐴(𝐵𝐴)𝐵. By Theorem 4.3.3, applied twice,

𝑚=rank(𝐼𝑚)=rank(𝐴(𝐵𝐴)𝐵)≤rank(𝐵𝐴)≤𝑚,

so rank(𝐵𝐴)=𝑚.

4.3.3Row rank equals column rank

Intuition (Why this is surprising). Consider

𝐴=(1214240612−12)∈ℚ3×4.

Its four columns are vectors in ℚ3. They satisfy 𝐚2=2𝐚1 and 𝐚4=3𝐚1+𝐚3, and 𝐚1,𝐚3 are independent, so exactly two of the columns are independent. Its three rows are vectors in ℚ4. They satisfy one relation, row1−row2+row3=𝟎, and any two of them are independent, so again exactly two of the rows are independent.

Nothing in the definitions suggests that these two numbers should agree. The columns satisfy two independent relations and the rows only one; the columns live in 𝐹3 and the rows in 𝐹4. For a 3×100 matrix, one count is about 100 vectors in 𝐹3 and the other is about 3 vectors in 𝐹100. Yet the two counts always agree. Equivalently, the numbers of independent relations always differ by exactly the difference between the numbers of columns and rows: here 2−1=4−3.

The columns of 𝐴𝖳 are the rows of 𝐴, written as columns (Definition 2.4.21). So im𝐴𝖳⊆𝐹𝑛 is the span of the rows of 𝐴 (Proposition 2.6.18), and rank𝐴𝖳 is its dimension.

Theorem 4.3.8 (Row rank equals column rank). For every 𝐴∈𝐹𝑚×𝑛,

rank𝐴𝖳=rank𝐴.

That is, the span of the rows of 𝐴, viewed as vectors in 𝐹𝑛, has the same dimension as the span of its columns in 𝐹𝑚.

The first proof reduces 𝐴 to its rank normal form, where the statement is obvious.

Proof (via the rank normal form).routine

Uses. Theorem 3.5.27 (Rank normal form), Proposition 2.4.23 (Transpose of a product), Proposition 2.5.17 (Algebra of inverses), Proposition 4.3.1 (Invertible factors preserve rank), Proposition 2.6.18 (The column space is spanned by the columns), Theorem 4.2.7 (The image of 𝑇𝐴)

Techniques. Transpose to trade rows for columns (Technique 3.2.27)

Hint 1. The statement is obvious for 𝑁𝑘=(𝐼𝑘000). Reduce to it.
Hint 2. Transpose 𝐺𝐴𝐻=𝑁𝑘. Which factors of the result are invertible, and what is the rank of 𝑁𝑘𝖳?
Sketch. With 𝐺𝐴𝐻=𝑁𝑘 and 𝑘=rank𝐴, transposing gives 𝐻𝖳𝐴𝖳𝐺𝖳=𝑁𝑘𝖳, with 𝐻𝖳 and 𝐺𝖳 invertible. By Proposition 4.3.1, rank𝐴𝖳=rank𝑁𝑘𝖳, and 𝑁𝑘𝖳 has columns 𝐞1,…,𝐞𝑘 followed by zero columns, so its rank is 𝑘.

Let 𝑘=rank𝐴. By Theorem 3.5.27 there are invertible 𝐺∈𝐹𝑚×𝑚 and 𝐻∈𝐹𝑛×𝑛 with 𝐺𝐴𝐻=𝑁𝑘. Transposing with Proposition 2.4.23,

𝐻𝖳𝐴𝖳𝐺𝖳=(𝐺𝐴𝐻)𝖳=𝑁𝑘𝖳.

The matrices 𝐻𝖳 and 𝐺𝖳 are invertible (Proposition 2.5.17(d)), so rank𝐴𝖳=rank(𝐻𝖳𝐴𝖳𝐺𝖳)=rank𝑁𝑘𝖳 by Proposition 4.3.1. The matrix 𝑁𝑘𝖳∈𝐹𝑛×𝑚 has 1 in positions (1,1),…,(𝑘,𝑘) and 0 elsewhere, so its columns are 𝐞1,…,𝐞𝑘∈𝐹𝑛 followed by 𝑚−𝑘 zero columns. Hence im𝑁𝑘𝖳=span(𝐞1,…,𝐞𝑘) (Proposition 2.6.18), which has the basis 𝐞1,…,𝐞𝑘. So rank𝑁𝑘𝖳=𝑘 (Theorem 4.2.7), and rank𝐴𝖳=𝑘=rank𝐴. □

The second proof uses only row operations. They change the column space of 𝐴 but keep its row space, and the reduced form exposes a basis of the row space.

Proof (row operations keep the row space, and RREF exposes it).standard

Uses. Theorem 3.3.10 (𝐺𝐴=rref(𝐴)), Proposition 2.4.23 (Transpose of a product), Proposition 2.5.17 (Algebra of inverses), Proposition 4.3.1 (Invertible factors preserve rank), Proposition 3.3.8 (Properties of the output), Lemma 3.3.16 (Reading the columns of a matrix in RREF), Proposition 2.6.18 (The column space is spanned by the columns), Theorem 1.5.5 (Independence is decided at zero), Theorem 4.2.7 (The image of 𝑇𝐴)

Techniques. Reduce to RREF (Technique 3.3.22), Transpose to trade rows for columns (Technique 3.2.27)

Hint 1. Transpose 𝐺𝐴=𝐴𝑟, and recall from Proposition 4.3.1(b) that an invertible factor on the right does not change the image.
Hint 2. So 𝐴 and 𝐴𝑟 have the same row space. Which rows of 𝐴𝑟 are nonzero, and why are they linearly independent?
Hint 3. Row 𝑠 of 𝐴𝑟 has a 1 in the pivot column 𝑝𝑠, and every other row has a 0 there.
Sketch. From 𝐺𝐴=𝐴𝑟 we get 𝐴𝑟𝖳=𝐴𝖳𝐺𝖳 with 𝐺𝖳 invertible, so im𝐴𝑟𝖳=im𝐴𝖳 by Proposition 4.3.1(b): 𝐴 and 𝐴𝑟 have the same row space. The nonzero rows of 𝐴𝑟 are its first 𝑘=rank𝐴 rows. They span its row space, and they are independent because row 𝑠 is the only one with a nonzero entry in column 𝑝𝑠. So the row space has dimension 𝑘.

Let 𝑘=rank𝐴, and let 𝐺∈𝐹𝑚×𝑚 be invertible with 𝐺𝐴=𝐴𝑟 (Theorem 3.3.10). Transposing, 𝐴𝑟𝖳=𝐴𝖳𝐺𝖳 (Proposition 2.4.23), and 𝐺𝖳 is invertible (Proposition 2.5.17(d)). By Proposition 4.3.1(b),

im𝐴𝑟𝖳=im𝐴𝖳.

That is, 𝐴 and 𝐴𝑟 have the same row space.

By Proposition 3.3.8(a), rows 𝑘+1,…,𝑚 of 𝐴𝑟 are zero. Let 𝛒1,…,𝛒𝑘∈𝐹𝑛 be the transposes of rows 1,…,𝑘 of 𝐴𝑟. These are the first 𝑘 columns of 𝐴𝑟𝖳, and its other columns are zero, so im𝐴𝑟𝖳=span(𝛒1,…,𝛒𝑘) (Proposition 2.6.18). We show that 𝛒1,…,𝛒𝑘 is linearly independent. Let 𝑝1<⋯<𝑝𝑘 be the pivot columns of 𝐴𝑟. Column 𝑝𝑡 of 𝐴𝑟 is 𝐞𝑡 (Lemma 3.3.16(b)), so entry 𝑝𝑡 of 𝛒𝑠 is 1 if 𝑠=𝑡 and 0 if 𝑠≠𝑡. Hence if 𝑐1𝛒1+⋯+𝑐𝑘𝛒𝑘=𝟎, then entry 𝑝𝑡 of the left side is 𝑐𝑡, so 𝑐𝑡=0 for every 𝑡 (Theorem 1.5.5).

So 𝛒1,…,𝛒𝑘 is a basis of im𝐴𝑟𝖳=im𝐴𝖳, and rank𝐴𝖳=dimim𝐴𝖳=𝑘=rank𝐴 (Theorem 4.2.7). □

Both proofs reduce 𝐴 to a form in which the theorem can be seen, but they teach different things. The first uses row and column operations together and lands on 𝑁𝑘, whose transpose is again of the same shape. The second uses only row operations, and it rests on a contrast. Row operations change the column space but keep the relations among the columns (Corollary 3.2.17), hence the number of independent columns. They keep the row space itself. So reduction preserves both numbers, and in 𝐴𝑟 both are visibly the number of pivots. The second proof also computes something: the nonzero rows of rref(𝐴) form a basis of the row space of 𝐴. A third proof, from the column–row factorization, is in Section 4.3.6, and for real matrices Exercise 4.3.33 gives a fourth. The transpose of a matrix is the matrix of the dual map, and the coordinate-free form of the theorem is that a linear map and its dual have the same rank (Section 13.1). That point of view is beyond the course.

Remark 4.3.9 (The two halves of the product inequality). With the transpose available, the two bounds in Theorem 4.3.3 are mirror images of each other (Technique 3.2.27). If we only knew the easy bound rank(𝑋𝑌)≤rank𝑋, which comes from im(𝑋𝑌)⊆im𝑋, then

rank(𝐴𝐵)=rank((𝐴𝐵)𝖳)=rank(𝐵𝖳𝐴𝖳)≤rank𝐵𝖳=rank𝐵.

So the bound that needed a dimension argument follows from the bound about images, once rows and columns are known to have the same rank.

The theorem adds a condition about rows to the Second Fundamental Theorem.

Corollary 4.3.10 (Independent rows). Let 𝐴∈𝐹𝑚×𝑛. The following are equivalent.
  1. (i)The rows of 𝐴 are linearly independent.
  2. (ii)rank𝐴=𝑚.
  3. (iii)For every 𝐲∈𝐹𝑚, the equation 𝐴𝐱=𝐲 has a solution.
In particular, a square matrix is invertible if and only if its rows are linearly independent.

Proof.routine

Uses. Theorem 4.3.8 (Row rank equals column rank), Theorem 3.5.12 (First Fundamental Theorem, completed), Theorem 3.5.13 (Second Fundamental Theorem, completed), Corollary 3.5.25 (The invertible matrix theorem)

Techniques. Transpose to trade rows for columns (Technique 3.2.27)

Hint. The rows of 𝐴 are the columns of 𝐴𝖳∈𝐹𝑛×𝑚. Apply the First Fundamental Theorem to 𝐴𝖳.
Sketch. The rows of 𝐴 are independent iff the 𝑚 columns of 𝐴𝖳 are, iff rank𝐴𝖳=𝑚 (Theorem 3.5.12 for 𝐴𝖳), iff rank𝐴=𝑚, iff 𝑇𝐴 is surjective (Theorem 3.5.13).
The rows of 𝐴, written as columns, are the 𝑚 columns of 𝐴𝖳∈𝐹𝑛×𝑚. By (III) ⇔ (VI) of Theorem 3.5.12 for 𝐴𝖳, they are linearly independent if and only if rank𝐴𝖳=𝑚, that is, rank𝐴=𝑚 by Theorem 4.3.8. This proves (i) ⇔ (ii), and (ii) ⇔ (iii) is (VI) ⇔ (II) of Theorem 3.5.13. For 𝑚=𝑛, condition (ii) says that 𝐴 is invertible (Corollary 3.5.25). □

Read as a statement about equations, (i) ⇔ (iii) is natural in one direction. If some combination of the rows vanishes, say 𝑧1(row1)+⋯+𝑧𝑚(row𝑚)=𝟎 with not all 𝑧𝑖=0, then every solvable right side satisfies 𝑧1𝑦1+⋯+𝑧𝑚𝑦𝑚=0, so not every 𝐲 is reachable. The corollary says that this is the only obstruction.

Exercise 4.3.11.routine Over 𝔽5, let 𝐴=(120321100423).

  1. (a)Find rank𝐴, and a basis of the row space of 𝐴 from rref(𝐴).
  2. (b)Find a basis of the row space consisting of rows of 𝐴, and express the remaining row in terms of it.

Solution.

Uses. Theorem 4.3.8 (Row rank equals column rank), Proposition 4.3.1 (Invertible factors preserve rank), Corollary 4.1.11 (Spanning lists are long)

Hint 1. Subtract twice row 1 from row 2, then scale by 2−1=3.
Hint 2. For (b), the row space has dimension rank𝐴. Are rows 1 and 2 multiples of each other?
Sketch. rref(𝐴)=(104401320000), so rank𝐴=2 and (1,0,4,4),(0,1,3,2) is a basis of the row space. Rows 1 and 2 are not proportional, so they are independent, and two independent vectors in a 2-dimensional space form a basis. Row 3 is row 1 plus twice row 2.

(a) Adding −2 times row 1 to row 2 gives (0,1−4,1,−6)=(0,2,1,4) in 𝔽5, and scaling by 2−1=3 gives (0,6,3,12)=(0,1,3,2). Adding −4 times this row to row 3 gives (0,0,2−12,3−8)=(0,0,0,0), and adding −2 times it to row 1 gives (1,0,−6,−1)=(1,0,4,4). So

rref(𝐴)=(104401320000),

rank𝐴=2, and by the second proof of Theorem 4.3.8 the nonzero rows (1,0,4,4) and (0,1,3,2) form a basis of the row space of 𝐴.

(b) The row space has dimension 2. Rows 1 and 2 of 𝐴 are not proportional (a multiple 𝑐(1,2,0,3) with third entry 1 does not exist), so they are linearly independent, and two independent vectors in a 2-dimensional space form a basis. Row 3 must be a combination of them, say row 3=𝑎 (row 1) +𝑏 (row 2). The first entries give 0=𝑎+2𝑏 and the second entries 4=2𝑎+𝑏; so 𝑎=−2𝑏=3𝑏 and 4=7𝑏=2𝑏 in 𝔽5, giving 𝑏=2 and 𝑎=1. Indeed (1,2,0,3)+2(2,1,1,0)=(5,4,2,3)=(0,4,2,3) in 𝔽5.

4.3.4The four fundamental subspaces

A matrix 𝐴∈𝐹𝑚×𝑛 has a kernel and an image, and so does its transpose. Together they are the four subspaces attached to 𝐴.

Definition 4.3.12 (The four fundamental subspaces). Let 𝐴∈𝐹𝑚×𝑛. The four fundamental subspaces of 𝐴 are

  • the column space im𝐴⊆𝐹𝑚, the span of the columns of 𝐴;
  • the null space ker𝐴⊆𝐹𝑛;
  • the row space im𝐴𝖳⊆𝐹𝑛, the span of the rows of 𝐴 written as columns;
  • the left null space ker𝐴𝖳⊆𝐹𝑚.

The column rank of 𝐴 is dimim𝐴, and the row rank is dimim𝐴𝖳.

So Theorem 4.3.8 says that the row rank equals the column rank, and both equal rank𝐴. Table 4.3.1 summarizes the dimensions and bases found below. The name “left null space” comes from writing the condition on the left: since (𝐴𝖳𝐲)𝖳=𝐲𝖳𝐴 (Proposition 2.4.23), a vector 𝐲∈𝐹𝑚 lies in ker𝐴𝖳 exactly when 𝐲𝖳𝐴=𝟎𝖳, that is,

𝑦1(row1of𝐴)+⋯+𝑦𝑚(row𝑚of𝐴)=𝟎.

So the null space records the linear relations among the columns of 𝐴, and the left null space records the linear relations among the rows, or equivalently among the equations of the system 𝐴𝐱=𝐲.

Proposition 4.3.13 (Dimensions of the four subspaces). Let 𝐴∈𝐹𝑚×𝑛 have rank 𝑟, and let 𝐺∈𝐹𝑚×𝑚 be invertible with 𝐺𝐴=𝐴𝑟.
  1. (a)dimim𝐴=𝑟, dimker𝐴=𝑛−𝑟, dimim𝐴𝖳=𝑟, and dimker𝐴𝖳=𝑚−𝑟.
  2. (b)The pivot columns of 𝐴 form a basis of im𝐴; the basic solutions form a basis of ker𝐴; the nonzero rows of 𝐴𝑟, written as columns, form a basis of im𝐴𝖳; and the last 𝑚−𝑟 rows of 𝐺, written as columns, form a basis of ker𝐴𝖳.

Proof.routine

Uses. Theorem 4.2.7 (The image of 𝑇𝐴), Theorem 4.2.11 (Rank–nullity), Theorem 4.3.8 (Row rank equals column rank), Remark 3.4.11 (A basis of the kernel), Exercise 3.3.57, Corollary 4.1.11 (Spanning lists are long)

Hint 1. Rank–nullity for 𝐴 gives two of the dimensions. Rank–nullity for 𝐴𝖳∈𝐹𝑛×𝑚, together with Theorem 4.3.8, gives the other two.
Hint 2. For the basis of ker𝐴𝖳, Exercise 3.3.57 shows that the last 𝑚−𝑟 rows of 𝐺 span the relations among the rows of 𝐴. How many of them are there?
Sketch. (a) dimim𝐴=𝑟 and dimker𝐴=𝑛−𝑟 are Theorem 4.2.7 and rank–nullity. For 𝐴𝖳∈𝐹𝑛×𝑚, rank𝐴𝖳=𝑟 and rank–nullity gives dimker𝐴𝖳=𝑚−𝑟. (b) The first two bases are Theorem 4.2.7 and Remark 3.4.11, and the third is the second proof of Theorem 4.3.8. The last 𝑚−𝑟 rows of 𝐺, transposed, span ker𝐴𝖳 by Exercise 3.3.57, and 𝑚−𝑟 vectors spanning a space of dimension 𝑚−𝑟 form a basis.

(a) By Theorem 4.2.7 and Theorem 4.2.11(a), dimim𝐴=𝑟 and dimker𝐴=𝑛−𝑟. The transpose 𝐴𝖳∈𝐹𝑛×𝑚 has 𝑚 columns and rank 𝑟 (Theorem 4.3.8), so dimim𝐴𝖳=𝑟 and dimker𝐴𝖳=𝑚−𝑟, by the same two results applied to 𝐴𝖳.

(b) The pivot columns of 𝐴 form a basis of im𝐴 by Theorem 4.2.7, and the basic solutions form a basis of ker𝐴 by Remark 3.4.11. The nonzero rows of 𝐴𝑟, transposed, form a basis of im𝐴𝖳, as shown in the second proof of Theorem 4.3.8. Finally, 𝐲∈ker𝐴𝖳 if and only if 𝐲𝖳𝐴=𝟎𝖳, that is, 𝐲𝖳 is a relation among the rows of 𝐴. By Exercise 3.3.57, these relations are exactly the linear combinations of the last 𝑚−𝑟 rows of 𝐺. So the transposes of these 𝑚−𝑟 rows span ker𝐴𝖳, a space of dimension 𝑚−𝑟 by (a), and they form a basis by Corollary 4.1.11(b). □

SubspaceDimensionBasis, from 𝐺𝐴=𝐴𝑟
im𝐴⊆𝐹𝑚𝑟the pivot columns of 𝐴
ker𝐴⊆𝐹𝑛𝑛−𝑟the basic solutions
im𝐴𝖳⊆𝐹𝑛𝑟the nonzero rows of 𝐴𝑟
ker𝐴𝖳⊆𝐹𝑚𝑚−𝑟the last 𝑚−𝑟 rows of 𝐺
Table 4.3.1. The four fundamental subspaces of 𝐴∈𝐹𝑚×𝑛 of rank 𝑟 (Proposition 4.3.13). Rows are written as columns. The two subspaces of 𝐹𝑛 have dimensions adding to 𝑛, and the two subspaces of 𝐹𝑚 have dimensions adding to 𝑚.

Example 4.3.14 (The four subspaces of a 3×4 matrix). Take the matrix of Intuition (Why this is surprising),

𝐴=(1214240612−12),𝐴𝑟=(120300110000).

Here 𝑚=3, 𝑛=4, 𝑟=2, and the pivot columns are 1 and 3.

  • Column space. Columns 1 and 3 of 𝐴: im𝐴=span((1,2,1),(1,0,−1)), of dimension 2.
  • Null space. The free columns are 2 and 4, with basic solutions (−2,1,0,0) and (−3,0,−1,1), so dimker𝐴=2.
  • Row space. The nonzero rows of 𝐴𝑟: im𝐴𝖳=span((1,2,0,3),(0,0,1,1)), of dimension 2.
  • Left null space. The relation row1−row2+row3=𝟎 gives (1,−1,1)∈ker𝐴𝖳, and dimker𝐴𝖳=3−2=1, so ker𝐴𝖳=span((1,−1,1)).

The dimensions are 2+2=4=𝑛 in 𝐹𝑛 and 2+1=3=𝑚 in 𝐹𝑚.

In this example something more is true. Every vector of the row space has dot product 0 with every vector of the null space: for instance (1,2,0,3)⋅(−3,0,−1,1)=−3+3=0. Likewise (1,2,1)⋅(1,−1,1)=0=(1,0,−1)⋅(1,−1,1). This always happens. Writing 𝐮⋅𝐱=𝑢1𝑥1+⋯+𝑢𝑛𝑥𝑛, entry 𝑖 of 𝐴𝐱 is (row 𝑖) ⋅𝐱, so 𝐱∈ker𝐴 exactly when every vector of the row space has dot product 0 with 𝐱 (Exercise 3.3.50(a)); applied to 𝐴𝖳, the same holds for the column space and the left null space. Over ℝ, the row space and the null space are then orthogonal complements in ℝ𝑛: they meet only in 𝟎, and together they span ℝ𝑛. So ℝ𝑛 splits into two perpendicular pieces, of dimensions 𝑟 and 𝑛−𝑟, and ℝ𝑚 into two perpendicular pieces of dimensions 𝑟 and 𝑚−𝑟. This picture is behind least squares (Section 6.5), and it is proved in Section 6.6 (Theorem 6.6.21). Over some other fields, such as 𝔽2 and ℂ, the dimensions still add up, but the two pieces can overlap.

Non-example 4.3.15 (Row space and null space can overlap). Over 𝔽2, let 𝐴=(11). Its row space is span((1,1)), and (1,1) lies in ker𝐴, since 1+1=0. So im𝐴𝖳=ker𝐴: the two subspaces of 𝔽22 coincide, although their dimensions 1+1=2 add up correctly.

Over ℂ, let 𝐴=(1𝑖). Then (1,𝑖)∈im𝐴𝖳, and 𝐴(1,𝑖)=1+𝑖2=0, so (1,𝑖)∈ker𝐴 as well. Over ℂ the cure is to replace 𝐴𝖳 by the conjugate transpose 𝐴∗ of Definition 6.3.26: the column space of 𝐴∗ meets ker𝐴 only in 𝟎 (Section 6.6). Over 𝔽2 no such repair exists: for every symmetric bilinear pairing on 𝔽2𝑛 with 𝑛≥2, some nonzero vector pairs to 0 with itself (Section 12.4).

4.3.5Rank-one matrices and sums

The simplest nonzero matrices have rank one. In

(2−14−2)=(12)(2−1),

every column is a multiple of (1,2) and every row is a multiple of (2,−1), and the matrix is the product of these two vectors.

Proposition 4.3.16 (Rank-one matrices). Let 𝐴∈𝐹𝑚×𝑛. Then rank𝐴=1 if and only if 𝐴=𝐮𝐯𝖳 for some nonzero 𝐮∈𝐹𝑚 and 𝐯∈𝐹𝑛. In that case im𝐴=span(𝐮) and im𝐴𝖳=span(𝐯), and if also 𝐴=𝐮′𝐯′𝖳, then 𝐮′=𝑐𝐮 and 𝐯′=𝑐−1𝐯 for some nonzero 𝑐∈𝐹.

Proof.routine

Uses. Proposition 2.4.8 (Four views of the product), Proposition 2.6.18 (The column space is spanned by the columns), Proposition 1.2.35 (Zero products), Theorem 4.2.7 (The image of 𝑇𝐴), Proposition 2.4.23 (Transpose of a product)

Hint 1. Column 𝑗 of 𝐮𝐯𝖳 is 𝑣𝑗𝐮. What does this say about im(𝐮𝐯𝖳)?
Hint 2. Conversely, if rank𝐴=1, then im𝐴=span(𝐮) for some 𝐮≠𝟎. Write each column of 𝐴 as a multiple of 𝐮.
Hint 3. For uniqueness, 𝐮′ lies in im𝐴=span(𝐮). Then compare 𝐮𝐯𝖳 with 𝐮′𝐯′𝖳 in a row where 𝐮 is nonzero.
Sketch. Column 𝑗 of 𝐮𝐯𝖳 is 𝑣𝑗𝐮, so im(𝐮𝐯𝖳)=span(𝐮) when 𝐯≠𝟎, a line. If rank𝐴=1, then im𝐴=span(𝐮), so column 𝑗 is 𝑣𝑗𝐮 for some 𝑣𝑗, and 𝐴=𝐮𝐯𝖳. Transposing gives the row space. If 𝐮𝐯𝖳=𝐮′𝐯′𝖳, then span(𝐮′)=span(𝐮), so 𝐮′=𝑐𝐮, and row 𝑖 with 𝑢𝑖≠0 gives 𝐯=𝑐𝐯′.

If. Let 𝐮≠𝟎 and 𝐯≠𝟎. By the column view (Proposition 2.4.8(b)), column 𝑗 of 𝐮𝐯𝖳 is 𝑣𝑗𝐮. So every column lies in span(𝐮), and im(𝐮𝐯𝖳)⊆span(𝐮) by Proposition 2.6.18. Some 𝑣𝑗 is nonzero, and then 𝑣𝑗𝐮≠𝟎 (Proposition 1.2.35), so im(𝐮𝐯𝖳) is a nonzero subspace of the line span(𝐮), hence equal to it. So rank(𝐮𝐯𝖳)=1 (Theorem 4.2.7) and im𝐴=span(𝐮).

Only if. Let rank𝐴=1. Then im𝐴 is one-dimensional, so im𝐴=span(𝐮) for any nonzero 𝐮∈im𝐴. Each column 𝐚𝑗 lies in im𝐴, so 𝐚𝑗=𝑣𝑗𝐮 for some 𝑣𝑗∈𝐹. With 𝐯=(𝑣1,…,𝑣𝑛), column 𝑗 of 𝐮𝐯𝖳 is 𝑣𝑗𝐮=𝐚𝑗, so 𝐴=𝐮𝐯𝖳. And 𝐯≠𝟎, since 𝐴≠0.

The row space and uniqueness. By Proposition 2.4.23, 𝐴𝖳=𝐯𝐮𝖳, so im𝐴𝖳=span(𝐯) by the first paragraph. If also 𝐴=𝐮′𝐯′𝖳 with 𝐮′,𝐯′ nonzero, then span(𝐮′)=im𝐴=span(𝐮), so 𝐮′=𝑐𝐮 with 𝑐≠0. Then 𝐮(𝐯−𝑐𝐯′)𝖳=𝐮𝐯𝖳−𝐮′𝐯′𝖳=0. Choose 𝑖 with 𝑢𝑖≠0. By the row view (Proposition 2.4.8(c)), row 𝑖 of the left side is 𝑢𝑖(𝐯−𝑐𝐯′)𝖳, which is therefore zero, so 𝐯=𝑐𝐯′, that is, 𝐯′=𝑐−1𝐯. □

The rank normal form gives the same factorization in one line: if 𝐺𝐴𝐻=𝑁1=𝐞1𝐞1𝖳, then 𝐴=(𝐺−1𝐞1)(𝐞1𝖳𝐻−1), a column of 𝐺−1 times a row of 𝐻−1. The form 𝐮𝐯𝖳 is worth reaching for whenever rank one appears.

Technique 4.3.17 (Write a rank-one matrix as 𝐮𝐯𝖳). A matrix of rank one is 𝐮𝐯𝖳 with 𝐮,𝐯≠𝟎 (Proposition 4.3.16). In this form products collapse to scalars, because 𝐯𝖳𝐰 is a 1×1 matrix and can be moved out of the way:

(𝐮𝐯𝖳)𝐱=(𝐯𝖳𝐱)𝐮,(𝐮𝐯𝖳)(𝐰𝐳𝖳)=(𝐯𝖳𝐰)𝐮𝐳𝖳,(𝐮𝐯𝖳)2=(𝐯𝖳𝐮)𝐮𝐯𝖳.

Reach for it when a hypothesis says that a matrix has rank one, when a matrix differs from a known one by a rank-one matrix (as 𝐼𝑛+𝐮𝐯𝖳 or 𝐴+𝐮𝐯𝖳 do), and when a matrix is to be split into simple pieces. The canonical example is Proposition 4.3.16, and Example 4.3.18 shows the collapse at work. Later instances are the matrix determinant lemma for det(𝐴+𝐮𝐯𝖳) (Theorem 5.5.27) and the singular value decomposition, which writes every matrix as a sum of rank-one pieces in the best possible way (Section 9.1).

Example 4.3.18 (Powers of a rank-one matrix). Let 𝐴=𝐮𝐯𝖳∈𝐹𝑛×𝑛 have rank one. Its diagonal entries are 𝑢𝑖𝑣𝑖, so 𝐯𝖳𝐮=𝑢1𝑣1+⋯+𝑢𝑛𝑣𝑛=tr𝐴. By Technique 4.3.17,

𝐴2=𝐮(𝐯𝖳𝐮)𝐯𝖳=(tr𝐴)𝐴,

and by induction 𝐴𝑘=(tr𝐴)𝑘−1𝐴 for 𝑘≥1. So a rank-one matrix is nilpotent if and only if tr𝐴=0, in which case already 𝐴2=0; and 𝐴2=𝐴 if and only if tr𝐴=1, since 𝐴≠0. For instance (2−14−2)=(12)(2−1) has trace 0, and its square is 0; while (1100)=𝐞1(11) has trace 1 and equals its square.

Exercise 4.3.19.routine Let 𝐴=(311622−3−1−1).

  1. (a)Over ℚ, show that rank𝐴=1 by writing 𝐴=𝐮𝐯𝖳, and compute 𝐴10.
  2. (b)Over 𝔽2, compute 𝐴2.

Solution.

Uses. Proposition 4.3.16 (Rank-one matrices), Example 4.3.18 (Powers of a rank-one matrix)

Techniques. Write a rank-one matrix as 𝐮𝐯𝖳 (Technique 4.3.17)

Hint 1. Every row is a multiple of (3,1,1).
Hint 2. 𝐴𝑘=(tr𝐴)𝑘−1𝐴. Over 𝔽2, reduce the trace modulo 2.
Sketch. (a) 𝐴=(1,2,−1)(3,1,1)𝖳, so rank𝐴=1, tr𝐴=3+2−1=4, and 𝐴10=49𝐴=262144𝐴. (b) Over 𝔽2, tr𝐴=0, so 𝐴2=0.

(a) The rows are 1, 2, and −1 times (3,1,1), so 𝐴=𝐮𝐯𝖳 with 𝐮=(1,2,−1) and 𝐯=(3,1,1), both nonzero, and rank𝐴=1 by Proposition 4.3.16. Its trace is 3+2−1=4, so by Example 4.3.18, 𝐴10=49𝐴=262144𝐴.

(b) Over 𝔽2, 𝐴=(111000111)=(1,0,1)(1,1,1)𝖳, of rank one, with tr𝐴=1+0+1=0. So 𝐴2=(tr𝐴)𝐴=0. Directly: (1,1,1)⋅(1,0,1)=2=0 in 𝔽2.

For sums, the image of 𝐴+𝐵 lies inside im𝐴+im𝐵, which gives an upper bound but no exact formula.

Proposition 4.3.20 (Rank is subadditive). For 𝐴,𝐵∈𝐹𝑚×𝑛,

rank(𝐴+𝐵)≤rank𝐴+rank𝐵.

Proof (comparing images).routine

Uses. Proposition 1.4.24 (The sum of two subspaces), Proposition 4.1.21 (Dimensions of subspaces), Corollary 4.1.11 (Spanning lists are long)

Hint 1. (𝐴+𝐵)𝐱=𝐴𝐱+𝐵𝐱. Which subspace contains im(𝐴+𝐵)?
Hint 2. Concatenate a basis of im𝐴 and a basis of im𝐵.
Sketch. Every (𝐴+𝐵)𝐱=𝐴𝐱+𝐵𝐱 lies in im𝐴+im𝐵, which is spanned by a basis of im𝐴 followed by a basis of im𝐵, a list of length rank𝐴+rank𝐵.

For every 𝐱∈𝐹𝑛, (𝐴+𝐵)𝐱=𝐴𝐱+𝐵𝐱∈im𝐴+im𝐵, so im(𝐴+𝐵)⊆im𝐴+im𝐵, a subspace of 𝐹𝑚 (Proposition 1.4.24). Let 𝐚1,…,𝐚𝑘 and 𝐛1,…,𝐛𝑙 be bases of im𝐴 and im𝐵. Every vector of im𝐴+im𝐵 is a sum of a combination of the 𝐚𝑖 and a combination of the 𝐛𝑗, so the list 𝐚1,…,𝐚𝑘,𝐛1,…,𝐛𝑙 spans im𝐴+im𝐵, and dim(im𝐴+im𝐵)≤𝑘+𝑙 (Corollary 4.1.11(a)). Hence

rank(𝐴+𝐵)≤dim(im𝐴+im𝐵)≤𝑘+𝑙=rank𝐴+rank𝐵

by Proposition 4.1.21(a). □

Proof (comparing kernels).routine

Uses. Proposition 4.2.26 (Subspaces that must meet), Theorem 4.2.11 (Rank–nullity), Proposition 4.1.21 (Dimensions of subspaces)

Techniques. Dimension counting (Technique 4.2.24)

Hint 1. ker𝐴∩ker𝐵⊆ker(𝐴+𝐵).
Hint 2. How large must ker𝐴∩ker𝐵 be? Use Proposition 4.2.26 in 𝐹𝑛.
Sketch. ker𝐴∩ker𝐵⊆ker(𝐴+𝐵), and by Proposition 4.2.26, dim(ker𝐴∩ker𝐵)≥(𝑛−rank𝐴)+(𝑛−rank𝐵)−𝑛. So dimker(𝐴+𝐵)≥𝑛−rank𝐴−rank𝐵, and rank–nullity gives the claim.

If 𝐴𝐱=𝐵𝐱=𝟎, then (𝐴+𝐵)𝐱=𝟎, so ker𝐴∩ker𝐵⊆ker(𝐴+𝐵). By Proposition 4.1.21(a), Proposition 4.2.26 in 𝐹𝑛, and Theorem 4.2.11(a),

dimker(𝐴+𝐵)≥dim(ker𝐴∩ker𝐵)≥(𝑛−rank𝐴)+(𝑛−rank𝐵)−𝑛=𝑛−rank𝐴−rank𝐵.

By Theorem 4.2.11(a) again, rank(𝐴+𝐵)=𝑛−dimker(𝐴+𝐵)≤rank𝐴+rank𝐵. □

As for products, the image proof works in the target and the kernel proof in the domain. On the kernel side, the inequality is the statement that two large subspaces of 𝐹𝑛 must overlap (Proposition 4.2.26).

By induction on 𝑘, a sum of 𝑘 matrices of rank one has rank at most 𝑘. The next subsection shows that this bound is attained: every matrix of rank 𝑟 is a sum of 𝑟 matrices of rank one.

4.3.6The column–row factorization

The pivot columns of 𝐴 form a basis of im𝐴 (Theorem 4.2.7), so every column of 𝐴 is a combination of them, and the coefficients can be read off rref(𝐴) (Lemma 3.3.16(d), Corollary 3.2.17). Recording all these combinations at once is a matrix product.

Theorem 4.3.21 (Column–row factorization). Let 𝐴∈𝐹𝑚×𝑛 have rank 𝑟≥1 and pivot columns 𝑝1<⋯<𝑝𝑟. Let 𝐶∈𝐹𝑚×𝑟 be the matrix whose columns are the pivot columns 𝐚𝑝1,…,𝐚𝑝𝑟 of 𝐴, and let 𝑅∈𝐹𝑟×𝑛 consist of the first 𝑟 rows of rref(𝐴), which are its nonzero rows. Then

𝐴=𝐶𝑅.

Moreover, the columns of 𝐶 form a basis of im𝐴, and the rows of 𝑅, written as columns, form a basis of im𝐴𝖳. In particular rank𝐶=rank𝑅=𝑟.

Proof.routine

Uses. Theorem 3.3.10 (𝐺𝐴=rref(𝐴)), Proposition 3.3.8 (Properties of the output), Theorem 2.4.34 (Block multiplication), Lemma 3.3.16 (Reading the columns of a matrix in RREF), Definition 2.4.1 (Matrix product), Theorem 4.2.7 (The image of 𝑇𝐴), Theorem 4.3.8 (Row rank equals column rank), Theorem 3.3.12 (Matrices in RREF)

Hint 1. Write 𝐴=𝐺−1𝐴𝑟, where 𝐴𝑟=(𝑅0) in block form.
Hint 2. Split 𝐺−1 into its first 𝑟 columns and the rest, and multiply in blocks. Then identify column 𝑠 of 𝐺−1 by looking at column 𝑝𝑠 of 𝐴=𝐺−1𝐴𝑟.
Sketch. 𝐺𝐴=𝐴𝑟=(𝑅0). With 𝐺−1=(𝐶′𝐶″) split after 𝑟 columns, block multiplication gives 𝐴=𝐶′𝑅. Column 𝑠 of 𝐶′ is 𝐺−1𝐞𝑠, and since column 𝑝𝑠 of 𝐴𝑟 is 𝐞𝑠, column 𝑝𝑠 of 𝐴 is also 𝐺−1𝐞𝑠; so 𝐶′=𝐶. The basis statements are Theorem 4.2.7 and the second proof of Theorem 4.3.8.

Let 𝐺∈𝐹𝑚×𝑚 be invertible with 𝐺𝐴=𝐴𝑟 (Theorem 3.3.10). By Proposition 3.3.8(a), rows 𝑟+1,…,𝑚 of 𝐴𝑟 are zero, so in block form 𝐴𝑟=(𝑅0), with a zero block of 𝑚−𝑟 rows. Split 𝐺−1=(𝐶′𝐶″), where 𝐶′∈𝐹𝑚×𝑟 consists of the first 𝑟 columns of 𝐺−1. By Theorem 2.4.34,

𝐴=𝐺−1𝐴𝑟=(𝐶′𝐶″)(𝑅0)=𝐶′𝑅+𝐶″0=𝐶′𝑅.

Column 𝑠 of 𝐶′ is 𝐺−1𝐞𝑠. Column 𝑝𝑠 of 𝐴𝑟 is 𝐞𝑠 (Lemma 3.3.16(b)), so column 𝑝𝑠 of 𝐴=𝐺−1𝐴𝑟 is 𝐺−1𝐞𝑠 (Definition 2.4.1). Hence column 𝑠 of 𝐶′ is 𝐚𝑝𝑠, that is, 𝐶′=𝐶, and 𝐴=𝐶𝑅.

The columns of 𝐶 form a basis of im𝐴 by Theorem 4.2.7, so rank𝐶=𝑟. The rows of 𝑅 are the nonzero rows of 𝐴𝑟, and the second proof of Theorem 4.3.8 showed that, written as columns, they form a basis of im𝐴𝖳. Finally, 𝑅 satisfies (R1)–(R4) with 𝑟 pivots, so it is in RREF and rank𝑅=𝑟 (Theorem 3.3.12). □

Read column by column (Proposition 2.4.8(b)), 𝐴=𝐶𝑅 says that column 𝑗 of 𝐴 is the combination of the pivot columns of 𝐴 whose coefficients are column 𝑗 of rref(𝐴):

𝐚𝑗=𝑅1𝑗𝐚𝑝1+⋯+𝑅𝑟𝑗𝐚𝑝𝑟.

This is Corollary 3.2.17 once more: the free columns of rref(𝐴) record how the free columns of 𝐴 are built from its pivot columns.

Example 4.3.22 (A column–row factorization). For the matrix 𝐴 of Example 4.3.14, the pivot columns are 1 and 3, so 𝐴=𝐶𝑅 with

𝐶=(11201−1),𝑅=(12030011).

Column 4 of 𝑅 is (3,1), and indeed 𝐚4=3𝐚1+𝐚3; column 2 of 𝑅 is (2,0), and 𝐚2=2𝐚1. By the column–row view of the product (Proposition 2.4.8(d)), 𝐴 is a sum of two matrices of rank one:

𝐴=(121)(1203)+(10−1)(0011)=(120324061203)+(0011000000−1−1).
(1214240612−12)⇝(120300110000)
The reduction of Example 4.3.22. At the end, the pivot columns listed in the summary pick out the columns of 𝐴 that form 𝐶, and the nonzero rows of the final matrix form 𝑅. Use “Edit” to try your own matrix: the number of pivots is both the number of columns of 𝐶 and the number of rows of 𝑅.

The column–row factorization is one of many factorizations of 𝐴 through 𝐹𝑟, and the next proposition says that all of them carry the same information.

Proposition 4.3.23 (Factorizations through 𝐹𝑘). Let 𝐴=𝑋𝑌 with 𝑋∈𝐹𝑚×𝑘 and 𝑌∈𝐹𝑘×𝑛.
  1. (a)rank𝐴≤𝑘.
  2. (b)If 𝑘=rank𝐴, then the columns of 𝑋 form a basis of im𝐴, and the rows of 𝑌, written as columns, form a basis of im𝐴𝖳.

Proof.routine

Uses. Theorem 4.3.3 (Rank of a product), Proposition 3.5.2 (First properties of rank), Proposition 2.6.18 (The column space is spanned by the columns), Proposition 4.1.21 (Dimensions of subspaces), Corollary 4.1.11 (Spanning lists are long), Proposition 2.4.23 (Transpose of a product), Theorem 4.3.8 (Row rank equals column rank)

Techniques. Dimension counting (Technique 4.2.24), Transpose to trade rows for columns (Technique 3.2.27)

Hint 1. 𝑋 has only 𝑘 columns.
Hint 2. For (b), im𝐴⊆im𝑋, and im𝑋 is spanned by 𝑘 vectors. Compare dimensions. For 𝑌, transpose.
Sketch. (a) rank𝐴≤rank𝑋≤𝑘. (b) im𝐴⊆im𝑋, and im𝑋 is spanned by the 𝑘 columns of 𝑋, so dimim𝑋≤𝑘=dimim𝐴. Hence im𝐴=im𝑋, and 𝑘 vectors spanning a 𝑘-dimensional space form a basis. Apply the same argument to 𝐴𝖳=𝑌𝖳𝑋𝖳, of rank 𝑘 by Theorem 4.3.8.

(a) By Theorem 4.3.3 and Proposition 3.5.2(a), rank𝐴=rank(𝑋𝑌)≤rank𝑋≤𝑘.

(b) Let 𝑘=rank𝐴. By Theorem 4.3.3, im𝐴⊆im𝑋, and im𝑋 is spanned by the 𝑘 columns of 𝑋 (Proposition 2.6.18), so dimim𝑋≤𝑘=dimim𝐴 (Corollary 4.1.11(a)). Also dimim𝐴≤dimim𝑋 (Proposition 4.1.21(a)), so the dimensions agree and im𝐴=im𝑋 by Proposition 4.1.21(b). This is a space of dimension 𝑘 spanned by the 𝑘 columns of 𝑋; so these columns form a basis of im𝐴 (Corollary 4.1.11(b)). Next, 𝐴𝖳=𝑌𝖳𝑋𝖳 (Proposition 2.4.23), with 𝑌𝖳∈𝐹𝑛×𝑘, and rank𝐴𝖳=𝑘 (Theorem 4.3.8). The argument just given, applied to 𝐴𝖳, shows that the columns of 𝑌𝖳, which are the rows of 𝑌 written as columns, form a basis of im𝐴𝖳. □

So in any factorization 𝐴=𝑋𝑌 with inner size rank𝐴, the left factor is a basis of the column space and the right factor a basis of the row space. Two such factorizations differ by an invertible 𝑟×𝑟 matrix (Exercise 4.3.44), and the column–row factorization is the one in which the right factor is in RREF. Every matrix of rank 𝑟 is therefore a sum of 𝑟 rank-one matrices, and Proposition 4.3.20 says that fewer will not do.

Corollary 4.3.24 (Rank as a number of rank-one pieces). Let 𝐴∈𝐹𝑚×𝑛. Then rank𝐴 is the smallest 𝑘≥0 such that 𝐴 is a sum of 𝑘 matrices of rank one (the empty sum being 0).

Proof.routine

Uses. Theorem 4.3.21 (Column–row factorization), Proposition 2.4.8 (Four views of the product), Proposition 4.3.16 (Rank-one matrices), Proposition 4.3.20 (Rank is subadditive)

Techniques. Write a rank-one matrix as 𝐮𝐯𝖳 (Technique 4.3.17)

Hint 1. For a decomposition with rank𝐴 terms, write 𝐴=𝐶𝑅 and use the column–row view of the product.
Hint 2. For the lower bound, use Proposition 4.3.20 and induction on the number of terms.
Sketch. If 𝑟=rank𝐴≥1, then 𝐴=𝐶𝑅=∑𝑠=1𝑟𝐜𝑠𝑅𝑠,⋅, where 𝐜𝑠 is column 𝑠 of 𝐶 and 𝑅𝑠,⋅ is row 𝑠 of 𝑅; both are nonzero, so each term has rank one. Conversely, a sum of 𝑘 rank-one matrices has rank at most 𝑘 by subadditivity.

Let 𝑟=rank𝐴. If 𝑟=0, then 𝐴=0 is the empty sum. If 𝑟≥1, write 𝐴=𝐶𝑅 as in Theorem 4.3.21. By the column–row view (Proposition 2.4.8(d)),

𝐴=𝐜1𝑅1,⋅+⋯+𝐜𝑟𝑅𝑟,⋅,

where 𝐜𝑠 is column 𝑠 of 𝐶 and 𝑅𝑠,⋅ is row 𝑠 of 𝑅. Each 𝐜𝑠 is nonzero, as a vector of a basis, and each row of 𝑅 is nonzero, as it contains a pivot. So each term has rank one by Proposition 4.3.16, and 𝐴 is a sum of 𝑟 matrices of rank one.

Conversely, if 𝐴=𝐴1+⋯+𝐴𝑘 with rank𝐴𝑖=1 for each 𝑖, then by Proposition 4.3.20 and induction on 𝑘, rank𝐴≤rank𝐴1+⋯+rank𝐴𝑘=𝑘. So no sum of fewer than 𝑟 rank-one matrices equals 𝐴. □

The lower bound can also be seen from Proposition 4.3.23(a): a sum ∑𝑖=1𝑘𝐮𝑖𝐯𝑖𝖳 is the product of the 𝑚×𝑘 matrix with columns 𝐮𝑖 and the 𝑘×𝑛 matrix with rows 𝐯𝑖𝖳, so its rank is at most 𝑘. Storing 𝐶 and 𝑅 takes 𝑟(𝑚+𝑛) numbers instead of the 𝑚𝑛 entries of 𝐴, which is far fewer when 𝑟 is small. Low-rank approximation (Section 16.3) exploits this to compress data.

The factorization also gives a third proof of Theorem 4.3.8, the one in LADR. It never row reduces, and its idea is to prove an inequality for every matrix and then apply it to the transpose.

Proof of Theorem 4.3.8 (via a column–row factorization).routine

Uses. Theorem 4.2.7 (The image of 𝑇𝐴), Proposition 2.4.8 (Four views of the product), Proposition 2.4.23 (Transpose of a product), Theorem 4.3.3 (Rank of a product), Proposition 3.5.2 (First properties of rank), Proposition 2.4.22 (Transpose is a linear involution)

Techniques. Transpose to trade rows for columns (Technique 3.2.27)

Hint 1. It suffices to prove rank𝐴𝖳≤rank𝐴 for every matrix 𝐴. Why?
Hint 2. Write 𝐴=𝐶𝑅 with 𝐶 consisting of a basis of im𝐴, so that 𝑅 has only 𝑟=rank𝐴 rows. What does 𝐴𝖳=𝑅𝖳𝐶𝖳 say about rank𝐴𝖳?
Sketch. Let 𝑟=rank𝐴 and let 𝐶 have as columns a basis of im𝐴. Every column of 𝐴 is a combination of these, so 𝐴=𝐶𝑅 with 𝑅∈𝐹𝑟×𝑛. Then 𝐴𝖳=𝑅𝖳𝐶𝖳, and rank𝐴𝖳≤rank𝑅𝖳≤𝑟, since 𝑅𝖳 has 𝑟 columns. So rank𝐴𝖳≤rank𝐴 for every 𝐴; applied to 𝐴𝖳, it gives the reverse inequality.

We first show that rank𝐴𝖳≤rank𝐴 for every 𝐴∈𝐹𝑚×𝑛. Let 𝑟=rank𝐴. If 𝑟=0, then 𝐴=0, so 𝐴𝖳=0 and both ranks are 0. Otherwise let 𝐜1,…,𝐜𝑟 be a basis of im𝐴 (for instance the pivot columns, Theorem 4.2.7), and let 𝐶∈𝐹𝑚×𝑟 have these columns. Each column 𝐚𝑗 of 𝐴 lies in im𝐴, so 𝐚𝑗=𝑥1𝑗𝐜1+⋯+𝑥𝑟𝑗𝐜𝑟 for some scalars 𝑥𝑠𝑗. Let 𝑅=(𝑥𝑠𝑗)∈𝐹𝑟×𝑛. By the column view (Proposition 2.4.8(b)), column 𝑗 of 𝐶𝑅 is ∑𝑠𝑥𝑠𝑗𝐜𝑠=𝐚𝑗, so 𝐴=𝐶𝑅. Transposing, 𝐴𝖳=𝑅𝖳𝐶𝖳 (Proposition 2.4.23), and by Theorem 4.3.3 and Proposition 3.5.2(a),

rank𝐴𝖳≤rank𝑅𝖳≤𝑟=rank𝐴,

since 𝑅𝖳∈𝐹𝑛×𝑟 has only 𝑟 columns.

Applying this to 𝐴𝖳, whose transpose is 𝐴 (Proposition 2.4.22), gives rank𝐴≤rank𝐴𝖳. Hence rank𝐴𝖳=rank𝐴. □

Intuition (Comparing the proofs that row rank equals column rank). The normal-form proof reduces 𝐴 to the simplest matrix of its rank and observes that transposing does not disturb it; the work was done in proving Theorem 3.5.27 and Proposition 4.3.1. The row-operation proof is the computational one: it finds a basis of the row space, and it explains the theorem through what row operations preserve. The factorization proof is the shortest, given Theorem 4.3.3, and it shows the mechanism most directly: every row of 𝐴=𝐶𝑅 is a combination of the 𝑟 rows of 𝑅, so the rows cannot span more than 𝑟 dimensions. Its last step, proving an inequality for all matrices and then applying it to the transpose to get the reverse inequality, is a pattern worth remembering.

Exercise 4.3.25.routine Find the column–row factorization of

𝐴=(1𝑖02𝑖−112𝑖1+𝑖−1+𝑖12+2𝑖)∈ℂ3×4,

and write 𝐴 as a sum of rank𝐴 matrices of rank one.

Solution.

Uses. Theorem 4.3.21 (Column–row factorization), Corollary 4.3.24 (Rank as a number of rank-one pieces)

Hint 1. Clear the first column with row 1. Row 3 is the sum of rows 1 and 2.
Hint 2. The pivot columns are 1 and 3.
Sketch. rref(𝐴)=(1𝑖0200100000), with pivot columns 1,3. So 𝐶=(10𝑖11+𝑖1), 𝑅=(1𝑖020010), and 𝐴=(1,𝑖,1+𝑖)(1𝑖02)+(0,1,1)(0010).

Adding −𝑖 times row 1 to row 2 gives (𝑖−𝑖,−1−𝑖2,1,2𝑖−2𝑖)=(0,0,1,0). Row 3 is the sum of rows 1 and 2 of 𝐴, so adding −1 times each of them to it gives 𝟎. Hence

rref(𝐴)=(1𝑖0200100000),

with pivot columns 1 and 3, and rank𝐴=2. By Theorem 4.3.21,

𝐴=𝐶𝑅,𝐶=(10𝑖11+𝑖1),𝑅=(1𝑖020010).

Check: column 2 of 𝐶𝑅 is 𝑖(1,𝑖,1+𝑖)=(𝑖,−1,−1+𝑖) and column 4 is 2(1,𝑖,1+𝑖)=(2,2𝑖,2+2𝑖), the columns of 𝐴. By Corollary 4.3.24 and its proof,

𝐴=(1𝑖1+𝑖)(1𝑖02)+(011)(0010).

Exercises

Exercise 4.3.26.routine Over ℚ, let

𝐴=(101011112),𝐵1=(1011−10),𝐵2=(100100).
  1. (a)Find ker𝐴, and for 𝑖=1,2 find im𝐵𝑖∩ker𝐴. Use (4.3.1) to find rank(𝐴𝐵𝑖), and check by computing 𝐴𝐵𝑖.
  2. (b)Is there a matrix 𝐵∈ℚ3×2 of rank 2 with 𝐴𝐵=0?

Solution.

Uses. Proposition 4.3.5 (Rank of a product, exactly), Theorem 4.2.11 (Rank–nullity)

Hint 1. Row 3 of 𝐴 is the sum of rows 1 and 2, so rank𝐴=2 and ker𝐴 is a line.
Hint 2. For (b), how small can rank(𝐴𝐵) be when rank𝐵=2 and dimker𝐴=1?
Sketch. ker𝐴=span((1,1,−1)). im𝐵1 contains (1,1,−1), so im𝐵1∩ker𝐴=ker𝐴 and rank(𝐴𝐵1)=2−1=1; im𝐵2=span(𝐞1,𝐞2) misses it, so rank(𝐴𝐵2)=2. (b) No: rank(𝐴𝐵)=2−dim(im𝐵∩ker𝐴)≥2−1=1.

(a) 𝐴𝐱=𝟎 means 𝑥1+𝑥3=0, 𝑥2+𝑥3=0, and 𝑥1+𝑥2+2𝑥3=0, the third being the sum of the first two. So ker𝐴=span((1,1,−1)) and rank𝐴=2.

The first column of 𝐵1 is (1,1,−1), so ker𝐴⊆im𝐵1 and im𝐵1∩ker𝐴=ker𝐴, of dimension 1. The columns of 𝐵1 are independent, so rank𝐵1=2, and (4.3.1) gives rank(𝐴𝐵1)=2−1=1. Indeed

𝐴𝐵1=(1−101−111+1−21)=(000101),

which has rank 1.

im𝐵2=span(𝐞1,𝐞2) consists of the vectors with third entry 0, and the only such multiple of (1,1,−1) is 𝟎. So im𝐵2∩ker𝐴={𝟎} and rank(𝐴𝐵2)=2. Indeed 𝐴𝐵2=(100111), the first two columns of 𝐴, which are independent.

(b) No. If rank𝐵=2, then by (4.3.1), rank(𝐴𝐵)=2−dim(im𝐵∩ker𝐴)≥2−dimker𝐴=1, so 𝐴𝐵≠0.

Exercise 4.3.27.routine Let 𝐴=(110000111111).

  1. (a)Over 𝔽2, find bases of the four fundamental subspaces of 𝐴, and check the dimensions against Proposition 4.3.13. Show that im𝐴𝖳=ker𝐴.
  2. (b)Over ℚ, do the same, and show that im𝐴𝖳∩ker𝐴={𝟎}.

Solution.

Uses. Proposition 4.3.13 (Dimensions of the four subspaces), Definition 4.3.12 (The four fundamental subspaces)

Hint 1. Over either field, row 3 is the sum of rows 1 and 2, and the first two rows are already in RREF.
Hint 2. Over 𝔽2, compare the basic solutions with the rows. Over ℚ, apply 𝐴 to a general vector 𝑎(1,1,0,0)+𝑏(0,0,1,1) of the row space.
Sketch. Over both fields rref(𝐴)=(110000110000), 𝑟=2, pivot columns 1,3. Over 𝔽2: im𝐴=span((1,0,1),(0,1,1)), ker𝐴=span((1,1,0,0),(0,0,1,1))=im𝐴𝖳, ker𝐴𝖳=span((1,1,1)). Over ℚ: ker𝐴=span((−1,1,0,0),(0,0,−1,1)), ker𝐴𝖳=span((1,1,−1)), and 𝐴(𝑎,𝑎,𝑏,𝑏)=(2𝑎,2𝑏,2𝑎+2𝑏) vanishes only for 𝑎=𝑏=0.

Over any field, row 3 of 𝐴 is the sum of rows 1 and 2, so subtracting rows 1 and 2 from row 3 gives

rref(𝐴)=(110000110000),

with 𝑟=2, pivot columns 1,3, and free columns 2,4. Here 𝑚=3 and 𝑛=4, so Proposition 4.3.13 predicts dimensions 2, 2, 2, 1.

(a) Over 𝔽2:

  • im𝐴 has the basis (1,0,1),(0,1,1), columns 1 and 3 of 𝐴.
  • The basic solutions are (−1,1,0,0)=(1,1,0,0) and (0,0,−1,1)=(0,0,1,1), a basis of ker𝐴.
  • im𝐴𝖳 has the basis (1,1,0,0),(0,0,1,1), the nonzero rows of rref(𝐴).
  • Row 1 plus row 2 plus row 3 is (2,2,2,2)=𝟎, so (1,1,1)∈ker𝐴𝖳, which has dimension 3−2=1; so ker𝐴𝖳=span((1,1,1)).

The bases of ker𝐴 and im𝐴𝖳 are the same list, so im𝐴𝖳=ker𝐴. Each row is “orthogonal” to itself: (1,1,0,0)⋅(1,1,0,0)=2=0 in 𝔽2.

(b) Over ℚ, im𝐴=span((1,0,1),(0,1,1)) and im𝐴𝖳=span((1,1,0,0),(0,0,1,1)) as before. The basic solutions are now (−1,1,0,0) and (0,0,−1,1), and row 1 plus row 2 minus row 3 is 𝟎, so ker𝐴𝖳=span((1,1,−1)). A vector of the row space is 𝐰=(𝑎,𝑎,𝑏,𝑏), and 𝐴𝐰=(2𝑎,2𝑏,2𝑎+2𝑏). This is 𝟎 only if 𝑎=𝑏=0, since 2≠0 in ℚ. So im𝐴𝖳∩ker𝐴={𝟎}. The two computations differ only in whether 1+1=0.

Exercise 4.3.28.routine For each item, find a matrix with the given properties, or show that none exists.

  1. (a)𝐴∈ℝ3×3 with im𝐴=span((1,0,1)) and ker𝐴=span((1,1,0),(0,1,1)).
  2. (b)𝐴∈𝐹2×4 with ker𝐴=span(𝐞1).
  3. (c)𝐴∈𝐹3×3 whose row space contains (1,1,0) and whose null space contains (1,0,0).
  4. (d)𝐴∈𝐹2×2 with im𝐴=ker𝐴.
  5. (e)𝐴∈𝐹3×3 with im𝐴=ker𝐴.

Solution.

Uses. Proposition 4.3.16 (Rank-one matrices), Proposition 4.3.13 (Dimensions of the four subspaces), Theorem 4.2.11 (Rank–nullity), Exercise 3.3.50

Hint 1. In (a), 𝐴 has rank 1, so 𝐴=𝐮𝐯𝖳. What must 𝐮 be, and which condition on 𝐯 makes ker𝐴 correct?
Hint 2. In (b) and (e), count dimensions. In (c), every row has dot product 0 with every vector of the null space (Exercise 3.3.50(a)).
Sketch. (a) 𝐴=(1,0,1)(1,−1,1)𝖳. (b) None: dimker𝐴≥4−2=2. (c) None: (1,1,0)⋅(1,0,0)=1≠0. (d) 𝐴=(0100). (e) None: 𝑟=3−𝑟 is impossible.

(a) Such an 𝐴 has rank 1, so 𝐴=𝐮𝐯𝖳, and we may take 𝐮=(1,0,1) (Proposition 4.3.16). Then 𝐴𝐱=(𝐯𝖳𝐱)𝐮, so ker𝐴={𝐱|𝐯𝖳𝐱=0}. We need 𝐯𝖳(1,1,0)=0 and 𝐯𝖳(0,1,1)=0, that is, 𝑣1+𝑣2=0=𝑣2+𝑣3; take 𝐯=(1,−1,1). Then

𝐴=(101)(1−11)=(1−110001−11).

Its kernel contains (1,1,0) and (0,1,1), which are independent, and has dimension 3−1=2 (Theorem 4.2.11), so it is their span; and im𝐴=span((1,0,1)).

(b) None. By Theorem 4.2.11, dimker𝐴=4−rank𝐴≥4−2=2, while dimspan(𝐞1)=1.

(c) None, over every field. By Exercise 3.3.50(a), every vector of the row space has dot product 0 with every vector of the null space, but (1,1,0)⋅(1,0,0)=1≠0.

(d) 𝐴=(0100): its image is spanned by its nonzero column 𝐞1, and 𝐴(𝑥1,𝑥2)=(𝑥2,0) vanishes exactly when 𝑥2=0, so ker𝐴=span(𝐞1) too.

(e) None. If im𝐴=ker𝐴, then rank𝐴=dimker𝐴=3−rank𝐴 by Theorem 4.2.11, so 2rank𝐴=3, which no integer satisfies. In general, a square matrix with im𝐴=ker𝐴 must have even size.

Exercise 4.3.29.routine Let 𝐴,𝐵∈𝐹𝑚×𝑛. Prove or disprove.

  1. (a)If 𝐴 and 𝐵 have the same four fundamental subspaces, then 𝐴=𝐵.
  2. (b)If 𝐴 and 𝐵 have the same row space, then ker𝐴=ker𝐵.
  3. (c)If 𝐴 and 𝐵 have the same column space, then ker𝐴𝖳=ker𝐵𝖳.
  4. (d)If ker𝐴=ker𝐵, then im𝐴=im𝐵.

Solution.

Hint 1. For (a), think of invertible matrices.
Hint 2. For (b), Exercise 3.3.50(a) describes ker𝐴 in terms of the row space alone. Then transpose for (c).
Hint 3. For (d), find two 2×2 matrices of rank 1 with the same kernel.
Sketch. (a) False: 𝐼2 and (1101) are both invertible. (b) True: ker𝐴 is the set of 𝐱 with 𝐰⋅𝐱=0 for all 𝐰 in the row space. (c) True: apply (b) to 𝐴𝖳 and 𝐵𝖳, whose row spaces are im𝐴 and im𝐵. (d) False: (1000) and (0010).

(a) False. 𝐼2 and 𝐵=(1101) are both invertible, so both have column space and row space 𝐹2 and null space and left null space {𝟎}; but 𝐼2≠𝐵.

(b) True. By Exercise 3.3.50(a), ker𝐴 is the set of 𝐱∈𝐹𝑛 with 𝐰⋅𝐱=0 for every 𝐰 in the row space of 𝐴. This description involves only the row space, so equal row spaces give equal null spaces.

(c) True. The row space of 𝐴𝖳 is im(𝐴𝖳)𝖳=im𝐴 (Proposition 2.4.22), and likewise for 𝐵. So if im𝐴=im𝐵, then 𝐴𝖳 and 𝐵𝖳 have the same row space, and (b) gives ker𝐴𝖳=ker𝐵𝖳.

(d) False. Let 𝐴=(1000) and 𝐵=(0010). Both send (𝑥1,𝑥2) to a vector that is zero exactly when 𝑥1=0, so ker𝐴=ker𝐵=span(𝐞2), but im𝐴=span(𝐞1) and im𝐵=span(𝐞2). Parts (b)–(d) show that the row space determines the null space and the column space determines the left null space, but that the null space does not determine the column space.

Exercise 4.3.30.standard (Sylvester’s rank inequality.) Let 𝐴∈𝐹𝑚×𝑛 and 𝐵∈𝐹𝑛×𝑝.

  1. (a)Show that rank(𝐴𝐵)≥rank𝐴+rank𝐵−𝑛.
  2. (b)Deduce that if 𝐴𝐵=0, then rank𝐴+rank𝐵≤𝑛.
  3. (c)Show that for 𝐴1,…,𝐴𝑘∈𝐹𝑛×𝑛,

    rank(𝐴1𝐴2⋯𝐴𝑘)≥rank𝐴1+⋯+rank𝐴𝑘−(𝑘−1)𝑛.
  4. (d)Deduce that if 𝐴∈𝐹𝑛×𝑛 has rank 𝑛−1, then rank(𝐴𝑘)≥𝑛−𝑘 for every 𝑘≥1, and show by an example that equality can hold for every 𝑘≤𝑛.

Solution.

Uses. Proposition 4.3.5 (Rank of a product, exactly), Theorem 4.2.11 (Rank–nullity), Proposition 4.1.21 (Dimensions of subspaces)

Techniques. Dimension counting (Technique 4.2.24)

Hint 1. By (4.3.1), rank(𝐴𝐵)=rank𝐵−dim(im𝐵∩ker𝐴). Bound the subtracted term.
Hint 2. For (c), induct on 𝑘, applying (a) to 𝐴1⋯𝐴𝑘−1 and 𝐴𝑘.
Hint 3. For the example in (d), take the shift 𝑁 with 𝑁𝐞1=𝟎 and 𝑁𝐞𝑗=𝐞𝑗−1.
Sketch. (a) dim(im𝐵∩ker𝐴)≤dimker𝐴=𝑛−rank𝐴. (b) Put rank(𝐴𝐵)=0. (c) Induction: rank(𝑃𝐴𝑘)≥rank𝑃+rank𝐴𝑘−𝑛 with 𝑃=𝐴1⋯𝐴𝑘−1. (d) With all 𝐴𝑖=𝐴, 𝑘(𝑛−1)−(𝑘−1)𝑛=𝑛−𝑘; the shift has rank𝑁𝑘=𝑛−𝑘.

(a) By (4.3.1), rank(𝐴𝐵)=rank𝐵−dim(im𝐵∩ker𝐴). Since im𝐵∩ker𝐴⊆ker𝐴, its dimension is at most dimker𝐴=𝑛−rank𝐴 (Proposition 4.1.21(a), Theorem 4.2.11(a)). Hence

rank(𝐴𝐵)≥rank𝐵−(𝑛−rank𝐴)=rank𝐴+rank𝐵−𝑛.

(b) If 𝐴𝐵=0, then 0=rank(𝐴𝐵)≥rank𝐴+rank𝐵−𝑛.

(c) By induction on 𝑘. For 𝑘=1 there is nothing to prove. For 𝑘≥2, let 𝑃=𝐴1⋯𝐴𝑘−1. By (a), rank(𝑃𝐴𝑘)≥rank𝑃+rank𝐴𝑘−𝑛, and by the induction hypothesis,

rank𝑃≥rank𝐴1+⋯+rank𝐴𝑘−1−(𝑘−2)𝑛.

Substituting the second inequality into the first gives the claim for 𝑘.

(d) With 𝐴1=⋯=𝐴𝑘=𝐴, (c) gives rank(𝐴𝑘)≥𝑘(𝑛−1)−(𝑘−1)𝑛=𝑛−𝑘. For the shift 𝑁∈𝐹𝑛×𝑛 with 𝑁𝐞1=𝟎 and 𝑁𝐞𝑗=𝐞𝑗−1 for 𝑗≥2, the power 𝑁𝑘 sends 𝐞𝑗 to 𝐞𝑗−𝑘 for 𝑗>𝑘 and to 𝟎 for 𝑗≤𝑘, so im𝑁𝑘=span(𝐞1,…,𝐞𝑛−𝑘) and rank𝑁𝑘=𝑛−𝑘 for 𝑘≤𝑛. In particular rank𝑁=𝑛−1.

Sylvester’s inequality is the lower bound that the exact formula hides: the product loses at most dimker𝐴 dimensions, because the lost directions lie in ker𝐴.

Exercise 4.3.31.hard (Frobenius’s rank inequality.) Let 𝐴∈𝐹𝑚×𝑛, 𝐵∈𝐹𝑛×𝑝, and 𝐶∈𝐹𝑝×𝑞.

  1. (a)Show that

    rank(𝐴𝐵)+rank(𝐵𝐶)≤rank𝐵+rank(𝐴𝐵𝐶).
  2. (b)Deduce Sylvester’s inequality (Exercise 4.3.30(a)).
  3. (c)Let 𝐴∈𝐹𝑛×𝑛 and 𝑑𝑘=rank(𝐴𝑘)−rank(𝐴𝑘+1) for 𝑘≥0. Deduce from (a) that 𝑑0≥𝑑1≥𝑑2≥⋯ (compare Example 4.3.6).

Solution.

Uses. Proposition 4.3.5 (Rank of a product, exactly), Proposition 2.6.15 (Kernels and images of composites), Proposition 4.1.21 (Dimensions of subspaces)

Techniques. Dimension counting (Technique 4.2.24)

Hint 1. Rewrite the inequality as rank(𝐵𝐶)−rank(𝐴(𝐵𝐶))≤rank𝐵−rank(𝐴𝐵): what 𝐴 destroys of im(𝐵𝐶) is at most what it destroys of im𝐵.
Hint 2. Express both sides with (4.3.1). How are im(𝐵𝐶) and im𝐵 related?
Hint 3. For (b), take 𝐵=𝐼𝑛. For (c), take the three matrices 𝐴, 𝐴𝑘, 𝐴.
Sketch. By (4.3.1), rank𝐵−rank(𝐴𝐵)=dim(im𝐵∩ker𝐴) and rank(𝐵𝐶)−rank(𝐴𝐵𝐶)=dim(im(𝐵𝐶)∩ker𝐴). Since im(𝐵𝐶)⊆im𝐵, the second is at most the first. (b) 𝐵=𝐼𝑛. (c) With 𝐴,𝐴𝑘,𝐴: 2rank(𝐴𝑘+1)≤rank(𝐴𝑘)+rank(𝐴𝑘+2), which rearranges to 𝑑𝑘+1≤𝑑𝑘.

(a) By (4.3.1), applied to 𝐴 and 𝐵 and to 𝐴 and 𝐵𝐶,

rank𝐵−rank(𝐴𝐵)=dim(im𝐵∩ker𝐴),rank(𝐵𝐶)−rank(𝐴𝐵𝐶)=dim(im(𝐵𝐶)∩ker𝐴).

By Proposition 2.6.15(b), im(𝐵𝐶)⊆im𝐵, so im(𝐵𝐶)∩ker𝐴⊆im𝐵∩ker𝐴, and Proposition 4.1.21(a) gives

rank(𝐵𝐶)−rank(𝐴𝐵𝐶)≤rank𝐵−rank(𝐴𝐵).

Rearranging gives (a). In words: of the directions in im(𝐵𝐶), which is part of im𝐵, the matrix 𝐴 can destroy no more than it destroys of all of im𝐵.

(b) Let 𝐴∈𝐹𝑚×𝑛 and 𝐶∈𝐹𝑛×𝑞, and take 𝐵=𝐼𝑛 in (a): rank𝐴+rank𝐶≤rank𝐼𝑛+rank(𝐴𝐶)=𝑛+rank(𝐴𝐶).

(c) Apply (a) to the three 𝑛×𝑛 matrices 𝐴, 𝐴𝑘, 𝐴:

rank(𝐴𝑘+1)+rank(𝐴𝑘+1)≤rank(𝐴𝑘)+rank(𝐴𝑘+2).

Rearranged, rank(𝐴𝑘+1)−rank(𝐴𝑘+2)≤rank(𝐴𝑘)−rank(𝐴𝑘+1), that is, 𝑑𝑘+1≤𝑑𝑘. (For 𝑘=0, 𝐴0=𝐼𝑛.) By rank–nullity, 𝑑𝑘=dimker𝐴𝑘+1−dimker𝐴𝑘, so this is Exercise 4.2.52(b) once more.

Exercise 4.3.32.hard Let 0≤𝑎≤min(𝑚,𝑛) and 0≤𝑏≤min(𝑛,𝑝) be integers. Show that there are 𝐴∈𝐹𝑚×𝑛 of rank 𝑎 and 𝐵∈𝐹𝑛×𝑝 of rank 𝑏 with rank(𝐴𝐵)=𝑟 if and only if

max(0,𝑎+𝑏−𝑛)≤𝑟≤min(𝑎,𝑏).

For 𝑚=𝑛=𝑝=5, 𝑎=3, and 𝑏=4, list the possible values of rank(𝐴𝐵) and give a pair 𝐴, 𝐵 for each.

Solution.

Uses. Theorem 4.3.3 (Rank of a product), Exercise 4.3.30, Lemma 2.4.24 (Matrix units), Proposition 2.6.18 (The column space is spanned by the columns), Corollary 1.5.17 (Sublists), Theorem 4.2.7 (The image of 𝑇𝐴)

Techniques. Finding counterexamples and minimal examples (Technique 1.4.18)

Hint 1. The bounds are Theorem 4.3.3 and Sylvester’s inequality. For the converse, use matrices made of 0s and a few 1s.
Hint 2. Let 𝐴 have 1 in positions (1,1),…,(𝑎,𝑎) and 0 elsewhere. Let 𝐵 have, in its columns 1,…,𝑏, distinct standard basis vectors 𝐞𝑡1,…,𝐞𝑡𝑏. What is 𝐴𝐵?
Hint 3. 𝐴𝐵 keeps exactly the columns 𝐞𝑡𝑘 with 𝑡𝑘≤𝑎. How many of the 𝑡𝑘 can you place in {1,…,𝑎}?
Sketch. Necessity: Theorem 4.3.3 and Exercise 4.3.30(a). Sufficiency: 𝐴=∑𝑖≤𝑎𝐞𝑖𝐞𝑖𝖳 and 𝐵=∑𝑘≤𝑏𝐞𝑡𝑘𝐞𝑘𝖳 with distinct 𝑡𝑘; then 𝐴𝐵=∑𝑘:𝑡𝑘≤𝑎𝐞𝑡𝑘𝐞𝑘𝖳 has rank |{𝑘|𝑡𝑘≤𝑎}|. Take 𝑡𝑘=𝑘 for 𝑘≤𝑟 and 𝑡𝑘=𝑎+𝑘−𝑟 for 𝑟<𝑘≤𝑏; this needs 𝑟≤𝑎 and 𝑎+𝑏−𝑟≤𝑛. For (5;3,4), 𝑟∈{2,3}.

Necessity. By Theorem 4.3.3, rank(𝐴𝐵)≤min(𝑎,𝑏), and by Exercise 4.3.30(a), rank(𝐴𝐵)≥𝑎+𝑏−𝑛; ranks are nonnegative.

Sufficiency. Let 𝑟 satisfy the bounds. A matrix whose nonzero columns are distinct standard basis vectors has rank equal to the number of those columns, since its image is spanned by them (Proposition 2.6.18), they are independent as part of the standard basis (Corollary 1.5.17), and rank is the dimension of the image (Theorem 4.2.7). Let

𝐴=∑𝑖=1𝑎𝐞𝑖𝐞𝑖𝖳∈𝐹𝑚×𝑛,

whose columns are 𝐞1,…,𝐞𝑎∈𝐹𝑚 followed by zeros, so rank𝐴=𝑎. Choose distinct indices 𝑡1,…,𝑡𝑏∈{1,…,𝑛}, specified below, and let

𝐵=∑𝑘=1𝑏𝐞𝑡𝑘𝐞𝑘𝖳∈𝐹𝑛×𝑝,

whose columns are 𝐞𝑡1,…,𝐞𝑡𝑏∈𝐹𝑛 followed by zeros, so rank𝐵=𝑏. By Lemma 2.4.24(c), 𝐞𝑖𝐞𝑖𝖳𝐞𝑡𝑘𝐞𝑘𝖳 is 𝐞𝑖𝐞𝑘𝖳 if 𝑡𝑘=𝑖 and 0 otherwise, so

𝐴𝐵=∑𝑘:𝑡𝑘≤𝑎𝐞𝑡𝑘𝐞𝑘𝖳,

whose nonzero columns are the distinct vectors 𝐞𝑡𝑘 with 𝑡𝑘≤𝑎. So rank(𝐴𝐵) is the number of 𝑘 with 𝑡𝑘≤𝑎.

Now take 𝑡𝑘=𝑘 for 1≤𝑘≤𝑟 and 𝑡𝑘=𝑎+(𝑘−𝑟) for 𝑟<𝑘≤𝑏. These are distinct, since 𝑟≤𝑎. They lie in {1,…,𝑛}, since the largest is 𝑎+𝑏−𝑟≤𝑛 when 𝑟≥𝑎+𝑏−𝑛. Exactly 𝑡1,…,𝑡𝑟 are at most 𝑎, as 𝑟≤𝑏 and 𝑡𝑘>𝑎 for 𝑘>𝑟. So rank(𝐴𝐵)=𝑟.

The example. With 𝑛=5, 𝑎=3, 𝑏=4: max(0,3+4−5)=2 and min(3,4)=3, so 𝑟∈{2,3}. Take 𝐴=(𝐼3000)∈𝐹5×5. For 𝑟=3, take 𝑡=(1,2,3,4), that is, 𝐵=(𝐼4000); then 𝐴𝐵=𝐴 has rank 3. For 𝑟=2, take 𝑡=(1,2,4,5), that is, 𝐵 has columns 𝐞1,𝐞2,𝐞4,𝐞5,𝟎; then 𝐴𝐵 has columns 𝐞1,𝐞2,𝟎,𝟎,𝟎, of rank 2.

Exercise 4.3.33.standard Let 𝐴∈𝐹𝑚×𝑛.

  1. (a)Show that, over any field, rank(𝐴𝖳𝐴)=rank𝐴−dim(im𝐴∩ker𝐴𝖳).
  2. (b)Let 𝐹=ℝ. Show that im𝐴∩ker𝐴𝖳={𝟎}, and deduce that rank(𝐴𝖳𝐴)=rank𝐴.
  3. (c)Let 𝐹=ℝ. Without using Theorem 4.3.8, deduce from (b) that rank𝐴𝖳=rank𝐴 and rank(𝐴𝐴𝖳)=rank𝐴.

Solution.

Uses. Proposition 4.3.5 (Rank of a product, exactly), Theorem 4.3.3 (Rank of a product), Proposition 2.4.23 (Transpose of a product), Proposition 2.4.22 (Transpose is a linear involution)

Techniques. A sum of nonnegative terms (Technique 6.1.9)

Hint 1. (a) is (4.3.1) for the product 𝐴𝖳𝐴.
Hint 2. For (b), if 𝐲=𝐴𝐱 and 𝐴𝖳𝐲=𝟎, compute 𝐲𝖳𝐲 in two ways.
Hint 3. For (c), rank𝐴=rank(𝐴𝖳𝐴)≤rank𝐴𝖳. Now apply the same to 𝐴𝖳.
Sketch. (a) is (4.3.1). (b) If 𝐲=𝐴𝐱∈ker𝐴𝖳, then ∑𝑖𝑦𝑖2=𝐲𝖳𝐲=𝐱𝖳𝐴𝖳𝐲=0, so 𝐲=𝟎. (c) rank𝐴=rank(𝐴𝖳𝐴)≤rank𝐴𝖳 by Theorem 4.3.3, and the same for 𝐴𝖳 gives the reverse; then rank(𝐴𝐴𝖳)=rank((𝐴𝖳)𝖳𝐴𝖳)=rank𝐴𝖳=rank𝐴.

(a) This is (4.3.1) with 𝐴𝖳 in the role of 𝐴 and 𝐴 in the role of 𝐵.

(b) Let 𝐲∈im𝐴∩ker𝐴𝖳, say 𝐲=𝐴𝐱 with 𝐴𝖳𝐲=𝟎. By Proposition 2.4.23, 𝐲𝖳=𝐱𝖳𝐴𝖳, so

𝑦12+⋯+𝑦𝑚2=𝐲𝖳𝐲=𝐱𝖳(𝐴𝖳𝐲)=𝐱𝖳𝟎=0.

A sum of squares of real numbers is 0 only if every term is 0 (Technique 6.1.9), so 𝐲=𝟎. By (a), rank(𝐴𝖳𝐴)=rank𝐴.

(c) By (b) and Theorem 4.3.3, rank𝐴=rank(𝐴𝖳𝐴)≤rank𝐴𝖳. Applying this to the real matrix 𝐴𝖳, whose transpose is 𝐴 (Proposition 2.4.22), gives rank𝐴𝖳≤rank𝐴. So rank𝐴𝖳=rank𝐴. Finally, (b) for 𝐴𝖳 gives rank(𝐴𝐴𝖳)=rank((𝐴𝖳)𝖳𝐴𝖳)=rank𝐴𝖳=rank𝐴.

Over ℂ the same argument works with the conjugate transpose 𝐴∗ in place of 𝐴𝖳, since 𝑦1𝑦1+⋯+𝑦𝑚𝑦𝑚=|𝑦1|2+⋯+|𝑦𝑚|2; this is Proposition 6.5.9. Exercise 4.3.34 shows what goes wrong with 𝐴𝖳.

Exercise 4.3.34.standard Prove or disprove.

  1. (a)Over 𝔽2, rank(𝐴𝖳𝐴)=rank𝐴 for every matrix 𝐴.
  2. (b)Over ℂ, rank(𝐴𝖳𝐴)=rank𝐴 for every matrix 𝐴.
  3. (c)Over ℚ, rank(𝐴𝖳𝐴)=rank𝐴 for every matrix 𝐴.
  4. (d)Over 𝔽2, if the columns of 𝐴 are linearly independent, then 𝐴𝖳𝐴 is invertible.

Which step of Exercise 4.3.33(b) fails in the false cases?

Solution.

Hint 1. Look for a nonzero vector 𝐲 with 𝑦12+⋯+𝑦𝑚2=0, and use it as a one-column matrix.
Hint 2. For (c), a rational matrix is a real matrix, and its rank does not depend on which of the two fields we compute in.
Sketch. (a) False: 𝐴=(1,1)𝖳 gives 𝐴𝖳𝐴=1+1=0. (b) False: 𝐴=(1,𝑖)𝖳 gives 1+𝑖2=0. (c) True, by Exercise 4.3.33 and Exercise 3.3.54. (d) False: 𝐴=(101001) gives 𝐴𝖳𝐴=(0001). The failing step is “a sum of squares is 0 only if each term is 0”.

(a) False. Let 𝐴=(11)∈𝔽22×1, of rank 1. Then 𝐴𝖳𝐴=(1⋅1+1⋅1)=(0), of rank 0.

(b) False. Let 𝐴=(1𝑖)∈ℂ2×1, of rank 1. Then 𝐴𝖳𝐴=(1+𝑖2)=(0).

(c) True. A matrix with rational entries is also a real matrix, and by Exercise 3.3.54(a) its reduced row echelon form, hence its rank, is the same over ℚ and over ℝ. The same holds for 𝐴𝖳𝐴, which also has rational entries. So rank(𝐴𝖳𝐴)=rank𝐴 over ℚ follows from Exercise 4.3.33(b) over ℝ.

(d) False. Over 𝔽2, let 𝐴=(101001), whose two columns are independent. Then

𝐴𝖳𝐴=(1+1001)=(0001),

which has rank 1 and is not invertible (Corollary 3.5.25).

In each false case, the vector 𝐲=𝐴𝐞1, that is, (1,1), (1,𝑖), or (1,1,0), is nonzero, satisfies 𝐲𝖳𝐲=0, and lies in im𝐴∩ker𝐴𝖳, which is exactly what Exercise 4.3.33(a) says must happen when the rank drops. The proof of Exercise 4.3.33(b) used that a sum of squares vanishes only when every term does, which is a property of ℝ (and of its subfields), not of fields in general.

Exercise 4.3.35.standard Find the error in the following argument. Is its conclusion true? Give a matrix for which the false step fails, and repair the argument.

Claim. For every 𝐴∈𝐹𝑚×𝑛, rank𝐴𝖳=rank𝐴.

Proof. Let 𝐴𝑟=rref(𝐴), with 𝑟 pivots. Row operations do not change the column space, so im𝐴=im𝐴𝑟. The columns of 𝐴𝑟 lie in span(𝐞1,…,𝐞𝑟) and include 𝐞1,…,𝐞𝑟, so im𝐴=span(𝐞1,…,𝐞𝑟), and rank𝐴=𝑟. Row operations do not change the row space either, and the 𝑟 nonzero rows of 𝐴𝑟 are linearly independent, so rank𝐴𝖳=𝑟 as well.

Solution.

Uses. Non-example 4.2.8 (The columns of 𝐴𝑟 are the wrong vectors), Proposition 4.3.1 (Invertible factors preserve rank), Theorem 4.2.7 (The image of 𝑇𝐴)

Techniques. Finding counterexamples and minimal examples (Technique 1.4.18)

Hint 1. Which operations preserve the column space: row operations or column operations?
Hint 2. Try 𝐴=(1224).
Hint 3. The conclusion rank𝐴=𝑟 is still right. Why do im𝐴 and im𝐴𝑟 have the same dimension even though they are different?
Sketch. The false step is im𝐴=im𝐴𝑟: row operations preserve the row space and the relations among the columns, not the column space. For 𝐴=(1224), im𝐴=span((1,2)) but im𝐴𝑟=span(𝐞1). The claim is true. Repair: 𝐴=𝐺−1𝐴𝑟, so rank𝐴=rank𝐴𝑟=𝑟 by Proposition 4.3.1(a) (or directly by Theorem 4.2.7); the rest of the argument is the second proof of Theorem 4.3.8.

The error. The step “row operations do not change the column space” is false. Row operations preserve the row space and the linear relations among the columns, but they move the columns themselves. For 𝐴=(1224) we have 𝐴𝑟=(1200), so im𝐴=span((1,2)) while im𝐴𝑟=span(𝐞1) (Non-example 4.2.8). In particular the intermediate conclusion im𝐴=span(𝐞1,…,𝐞𝑟) is false for this 𝐴.

The conclusion. The claim is true; it is Theorem 4.3.8. Only the first half of the argument needs repair. With 𝐺 invertible and 𝐺𝐴=𝐴𝑟, Proposition 4.3.1(a) gives rank𝐴=rank(𝐺𝐴)=rank𝐴𝑟=𝑟, although im𝐴≠im𝐴𝑟 in general. (This is also the definition of rank𝐴 as the number of pivots, together with Theorem 4.2.7.) The second half is correct as stated: 𝐴 and 𝐴𝑟 have the same row space, and the nonzero rows of 𝐴𝑟 are independent. It is the second proof of Theorem 4.3.8.

The two halves are not symmetric: row operations keep the row space as a set, but keep the column space only up to an invertible map, so only its dimension survives.

Exercise 4.3.36.standard Prove or disprove. All matrices are over a field 𝐹.

  1. (a)rank(𝐴𝐵)=rank(𝐵𝐴) for all 𝐴,𝐵∈𝐹𝑛×𝑛.
  2. (b)If 𝐴∈𝐹𝑚×𝑛 and 𝐵∈𝐹𝑛×𝑝 both have rank 𝑛, then rank(𝐴𝐵)=𝑛.
  3. (c)If 𝐴,𝐵∈𝐹𝑚×𝑛 and im𝐴∩im𝐵={𝟎}, then rank(𝐴+𝐵)=rank𝐴+rank𝐵.
  4. (d)If 𝐴∈𝐹𝑛×𝑛 and rank(𝐴2)=rank𝐴, then rank(𝐴𝑘)=rank𝐴 for every 𝑘≥1.

Solution.

Uses. Proposition 4.3.5 (Rank of a product, exactly), Theorem 4.2.11 (Rank–nullity), Proposition 2.6.15 (Kernels and images of composites), Example 4.3.4 (A product of rank zero)

Techniques. Finding counterexamples and minimal examples (Technique 1.4.18)

Hint 1. For (a), see Example 4.3.4.
Hint 2. For (b), what are ker𝐴 and im𝐵 when both ranks are 𝑛?
Hint 3. For (c), the columns may avoid each other while the rows do not. Try 2×2 matrices with a single nonzero entry each, in the same column.
Hint 4. For (d), use (4.3.1) with 𝐴𝑘 and 𝐴: the loss is dim(im𝐴𝑘∩ker𝐴).
Sketch. (a) False: 𝑀=(1000), 𝑁=(0100) give 𝑀𝑁=𝑁, 𝑁𝑀=0. (b) True: ker𝐴={𝟎}, so rank(𝐴𝐵)=rank𝐵=𝑛. (c) False: (1000)+(0010)=(1010) has rank 1. (d) True: im𝐴∩ker𝐴={𝟎}, and im𝐴𝑘⊆im𝐴, so rank(𝐴𝑘+1)=rank(𝐴𝑘)−dim(im𝐴𝑘∩ker𝐴)=rank(𝐴𝑘).

(a) False. With 𝑀=(1000) and 𝑁=(0100), we have 𝑀𝑁=𝑁, of rank 1, and 𝑁𝑀=0 (Example 4.3.4).

(b) True. 𝐴 has 𝑛 columns and rank 𝑛, so dimker𝐴=0 (Theorem 4.2.11). By (4.3.1), rank(𝐴𝐵)=rank𝐵−dim(im𝐵∩ker𝐴)=𝑛−0=𝑛.

(c) False. Let 𝐴=(1000) and 𝐵=(0010). Then im𝐴=span(𝐞1) and im𝐵=span(𝐞2) meet only in 𝟎, but 𝐴+𝐵=(1010) has rank 1<2. Here the row spaces coincide. (If both im𝐴∩im𝐵={𝟎} and im𝐴𝖳∩im𝐵𝖳={𝟎}, then the equality does hold, but that needs more work.)

(d) True. By (4.3.1) with 𝐵=𝐴, the hypothesis says im𝐴∩ker𝐴={𝟎}. For 𝑘≥1, (4.3.1) applied to 𝐴 and 𝐴𝑘 gives

rank(𝐴𝑘+1)=rank(𝐴𝑘)−dim(im𝐴𝑘∩ker𝐴).

Since 𝐴𝑘=𝐴∘𝐴𝑘−1 as maps, im𝐴𝑘⊆im𝐴 (Proposition 2.6.15(b)), so im𝐴𝑘∩ker𝐴⊆im𝐴∩ker𝐴={𝟎}, and rank(𝐴𝑘+1)=rank(𝐴𝑘). By induction, rank(𝐴𝑘)=rank𝐴 for all 𝑘≥1.

Exercise 4.3.37.standard Let 𝐮,𝐯∈𝐹𝑛, and let 𝑠=𝐯𝖳𝐮.

  1. (a)Show that if 1+𝑠≠0, then 𝐼𝑛+𝐮𝐯𝖳 is invertible, with

    (𝐼𝑛+𝐮𝐯𝖳)−1=𝐼𝑛−11+𝑠𝐮𝐯𝖳.
  2. (b)Show that if 1+𝑠=0, then 𝐼𝑛+𝐮𝐯𝖳 is not invertible, and that its rank is 𝑛−1.
  3. (c)Let 𝐽∈𝐹𝑛×𝑛 have every entry equal to 1. For which 𝑐∈𝐹 is 𝐼𝑛+𝑐𝐽 invertible, and what is its inverse? Is 𝐼4+𝐽 invertible over ℚ? Over 𝔽5?

Solution.

Uses. Theorem 3.5.23 (One-sided inverses of square matrices), Proposition 4.3.20 (Rank is subadditive), Proposition 4.3.16 (Rank-one matrices), Corollary 3.5.25 (The invertible matrix theorem)

Techniques. Write a rank-one matrix as 𝐮𝐯𝖳 (Technique 4.3.17), Guess the inverse, then check both products (Technique 2.5.16)

Hint 1. Multiply out, using 𝐮𝐯𝖳𝐮𝐯𝖳=𝑠𝐮𝐯𝖳.
Hint 2. For (b), apply 𝐼𝑛+𝐮𝐯𝖳 to 𝐮. For the rank, compare 𝐼𝑛 with (𝐼𝑛+𝐮𝐯𝖳)−𝐮𝐯𝖳.
Hint 3. 𝐽=𝟏𝟏𝖳, where 𝟏 is the vector of ones.
Sketch. (a) (𝐼+𝐮𝐯𝖳)(𝐼−(1+𝑠)−1𝐮𝐯𝖳)=𝐼+𝐮𝐯𝖳−(1+𝑠)−1(1+𝑠)𝐮𝐯𝖳=𝐼. (b) (𝐼+𝐮𝐯𝖳)𝐮=(1+𝑠)𝐮=𝟎 with 𝐮≠𝟎; and 𝑛=rank𝐼≤rank(𝐼+𝐮𝐯𝖳)+1. (c) 𝑐𝐽=(𝑐𝟏)𝟏𝖳, 𝑠=𝑐𝑛: invertible iff 1+𝑐𝑛≠0, inverse 𝐼−𝑐(1+𝑐𝑛)−1𝐽. 𝐼4+𝐽: invertible over ℚ (inverse 𝐼−𝐽/5), not over 𝔽5.

(a) Let 𝑁=𝐮𝐯𝖳. Then 𝑁2=𝐮(𝐯𝖳𝐮)𝐯𝖳=𝑠𝑁 (Technique 4.3.17), so

(𝐼𝑛+𝑁)(𝐼𝑛−(1+𝑠)−1𝑁)=𝐼𝑛+𝑁−(1+𝑠)−1(𝑁+𝑠𝑁)=𝐼𝑛+𝑁−𝑁=𝐼𝑛.

A one-sided inverse of a square matrix is two-sided (Theorem 3.5.23), so 𝐼𝑛+𝑁 is invertible with the stated inverse.

(b) Suppose 1+𝑠=0. Then 𝐮≠𝟎, since otherwise 𝑠=0. Now (𝐼𝑛+𝑁)𝐮=𝐮+𝐮(𝐯𝖳𝐮)=(1+𝑠)𝐮=𝟎, so 𝐼𝑛+𝑁 has a nonzero vector in its kernel and is not invertible (Corollary 3.5.25). Its rank is at most 𝑛−1. Also 𝐯≠𝟎, since 𝐯𝖳𝐮=−1, so 𝑁 has rank one (Proposition 4.3.16), and by Proposition 4.3.20,

𝑛=rank𝐼𝑛=rank((𝐼𝑛+𝑁)+(−𝑁))≤rank(𝐼𝑛+𝑁)+1.

So rank(𝐼𝑛+𝑁)=𝑛−1.

(c) 𝐽=𝟏𝟏𝖳, so 𝑐𝐽=𝐮𝐯𝖳 with 𝐮=𝑐𝟏 and 𝐯=𝟏, and 𝑠=𝟏𝖳(𝑐𝟏)=𝑐𝑛, where 𝑛 means 𝑛⋅1∈𝐹. By (a) and (b), 𝐼𝑛+𝑐𝐽 is invertible if and only if 1+𝑐𝑛≠0, and then

(𝐼𝑛+𝑐𝐽)−1=𝐼𝑛−𝑐1+𝑐𝑛𝐽.

For 𝐼4+𝐽, 𝑐=1 and 1+𝑐𝑛=5. Over ℚ, 5≠0, so (𝐼4+𝐽)−1=𝐼4−𝐽/5. Over 𝔽5, 5=0, so 𝐼4+𝐽 is not invertible, and indeed (𝐼4+𝐽)𝟏=𝟏+4𝟏=5𝟏=𝟎.

Exercise 4.3.38.standard Let 𝐴,𝐵∈𝐹𝑛×𝑛 and suppose that rank(𝐴𝐵−𝐵𝐴)=1. Show that (𝐴𝐵−𝐵𝐴)2=0. Give an example in which rank(𝐴𝐵−𝐵𝐴)=1.

Solution.

Uses. Proposition 4.3.16 (Rank-one matrices), Example 4.3.18 (Powers of a rank-one matrix), Proposition 2.4.28 (Cyclicity of the trace), Proposition 2.4.27 (The trace is linear)

Techniques. Write a rank-one matrix as 𝐮𝐯𝖳 (Technique 4.3.17), Cyclicity of trace (Technique 2.4.29)

Hint 1. A rank-one matrix 𝑋 satisfies 𝑋2=(tr𝑋)𝑋.
Hint 2. What is the trace of a commutator?
Sketch. 𝑋=𝐴𝐵−𝐵𝐴 has rank one, so 𝑋2=(tr𝑋)𝑋 by Example 4.3.18, and tr𝑋=tr(𝐴𝐵)−tr(𝐵𝐴)=0. Example: 𝐴=(1000), 𝐵=(0100), 𝐴𝐵−𝐵𝐴=𝐵.

Let 𝑋=𝐴𝐵−𝐵𝐴. Since rank𝑋=1, Example 4.3.18 gives 𝑋2=(tr𝑋)𝑋. By linearity and cyclicity of the trace (Proposition 2.4.27, Proposition 2.4.28), tr𝑋=tr(𝐴𝐵)−tr(𝐵𝐴)=0. Hence 𝑋2=0.

For an example, let 𝐴=(1000) and 𝐵=(0100). Then 𝐴𝐵=𝐵 and 𝐵𝐴=0, so 𝐴𝐵−𝐵𝐴=𝐵, which has rank 1, and indeed 𝐵2=0.

The hypothesis cannot be weakened to rank(𝐴𝐵−𝐵𝐴)≤2: with 𝐴=𝐞1𝐞2𝖳 and 𝐵=𝐞2𝐞1𝖳 in 𝐹2×2, 𝐴𝐵−𝐵𝐴=(100−1), which has rank 2 and square 𝐼2 over every field.

Exercise 4.3.39.routine
  1. (a)Show that rank(𝐴+𝐵)≥|rank𝐴−rank𝐵| for 𝐴,𝐵∈𝐹𝑚×𝑛.
  2. (b)Show that changing one entry of a matrix changes its rank by at most 1.
  3. (c)Let 𝑀𝑡∈ℝ5×5 have every diagonal entry equal to 1 and every other entry equal to 𝑡. Show that rank𝑀𝑡≥4 for every 𝑡≠1.

Solution.

Uses. Proposition 4.3.20 (Rank is subadditive), Proposition 4.3.16 (Rank-one matrices), Lemma 2.4.24 (Matrix units), Proposition 4.3.1 (Invertible factors preserve rank)

Techniques. Write a rank-one matrix as 𝐮𝐯𝖳 (Technique 4.3.17)

Hint 1. Write 𝐴=(𝐴+𝐵)+(−𝐵) and apply subadditivity.
Hint 2. Changing entry (𝑖,𝑗) by 𝑐 adds 𝑐𝐞𝑖𝐞𝑗𝖳.
Hint 3. 𝑀𝑡=(1−𝑡)𝐼5+𝑡𝐽, where 𝐽 is the all-ones matrix, of rank one.
Sketch. (a) rank𝐴≤rank(𝐴+𝐵)+rank(−𝐵) and symmetrically. (b) The change adds 𝑐𝐞𝑖𝐞𝑗𝖳, of rank at most 1; apply (a). (c) 𝑀𝑡=(1−𝑡)𝐼5+𝑡𝐽; for 𝑡≠1, (1−𝑡)𝐼5 has rank 5 and 𝑡𝐽 rank at most 1, so rank𝑀𝑡≥4.

(a) −𝐵 has the same image as 𝐵, so rank(−𝐵)=rank𝐵. By Proposition 4.3.20, rank𝐴=rank((𝐴+𝐵)+(−𝐵))≤rank(𝐴+𝐵)+rank𝐵, so rank(𝐴+𝐵)≥rank𝐴−rank𝐵. Exchanging 𝐴 and 𝐵 gives rank(𝐴+𝐵)≥rank𝐵−rank𝐴.

(b) Changing entry (𝑖,𝑗) of 𝐴 from 𝑎𝑖𝑗 to 𝑎𝑖𝑗+𝑐 replaces 𝐴 by 𝐴′=𝐴+𝑐𝐞𝑖𝐞𝑗𝖳 (Lemma 2.4.24). The matrix 𝑐𝐞𝑖𝐞𝑗𝖳 has rank 0 (if 𝑐=0) or 1 (Proposition 4.3.16). By (a) and Proposition 4.3.20, rank𝐴−1≤rank𝐴′≤rank𝐴+1.

(c) Let 𝐽 be the all-ones matrix, 𝐽=𝟏𝟏𝖳, of rank 1. Then 𝑀𝑡=(1−𝑡)𝐼5+𝑡𝐽. For 𝑡≠1, (1−𝑡)𝐼5 is invertible, of rank 5, and rank(𝑡𝐽)≤1, so by (a), rank𝑀𝑡≥5−1=4. (For 𝑡=1, 𝑀1=𝐽 has rank 1. By Exercise 4.3.37(c), applied to 𝐼5+(𝑡/(1−𝑡))𝐽, the rank is exactly 4 only for 𝑡=−1/4.)

Exercise 4.3.40.hard Let 𝑛≥3, and let 𝐴∈𝐹𝑛×𝑛 have entries 𝑎𝑖𝑗=(𝑖−𝑗)2, read as elements of 𝐹.

  1. (a)Show that rank𝐴≤3 over every field, by writing 𝐴 as a sum of three matrices of rank at most one.
  2. (b)Show that rank𝐴=3 if 1+1≠0 in 𝐹.
  3. (c)Find rank𝐴 over 𝔽2.

Solution.

Uses. Proposition 4.3.20 (Rank is subadditive), Proposition 4.3.16 (Rank-one matrices), Proposition 4.3.1 (Invertible factors preserve rank), Theorem 4.3.8 (Row rank equals column rank), Theorem 3.5.13 (Second Fundamental Theorem, completed), Proposition 2.6.15 (Kernels and images of composites), Theorem 3.4.24 (Polynomial interpolation), Corollary 3.5.16 (Lists of vectors in 𝐹𝑚), Theorem 2.4.34 (Block multiplication)

Techniques. Write a rank-one matrix as 𝐮𝐯𝖳 (Technique 4.3.17)

Hint 1. (𝑖−𝑗)2=𝑖2−2𝑖𝑗+𝑗2. Let 𝟏, 𝐰, 𝐪 be the vectors with entries 1, 𝑖, 𝑖2.
Hint 2. Write 𝐴=𝑋𝑀𝑋𝖳 with 𝑋=(𝟏𝐰𝐪)∈𝐹𝑛×3 and a 3×3 matrix 𝑀. When is 𝑀 invertible, and when are the columns of 𝑋 independent?
Hint 3. Over 𝔽2, −2=0 and 𝑖2=𝑖.
Sketch. (a) 𝐴=𝐪𝟏𝖳−2𝐰𝐰𝖳+𝟏𝐪𝖳. (b) 𝐴=𝑋𝑀𝑋𝖳 with 𝑀=(0010−20100), invertible when 2≠0. Rows 1,2,3 of 𝑋 form the Vandermonde matrix for the distinct nodes 1,2,3, so rank𝑋=3, and 𝑇𝑋𝖳 is onto 𝐹3. Hence im𝐴=im(𝑋𝑀)=im𝑋, of dimension 3. (c) Over 𝔽2, 𝑎𝑖𝑗=𝑖+𝑗, so 𝐴=𝐰𝟏𝖳+𝟏𝐰𝖳 has rank 2.

Let 𝟏,𝐰,𝐪∈𝐹𝑛 have entries 1, 𝑖, and 𝑖2 in position 𝑖.

(a) Since (𝑖−𝑗)2=𝑖2⋅1−2𝑖𝑗+1⋅𝑗2, entry by entry

𝐴=𝐪𝟏𝖳−2𝐰𝐰𝖳+𝟏𝐪𝖳.

Each term has rank at most one (Proposition 4.3.16), so rank𝐴≤3 by Proposition 4.3.20.

(b) Let 𝑋=(𝟏𝐰𝐪)∈𝐹𝑛×3 and 𝑀=(0010−20100). By block multiplication (Theorem 2.4.34), 𝑋𝑀𝑋𝖳=𝐪𝟏𝖳−2𝐰𝐰𝖳+𝟏𝐪𝖳=𝐴. Suppose 2≠0 in 𝐹.

𝑀 is invertible. 𝑀2=(100040001), which is invertible since 4≠0; so 𝑀 has the inverse 𝑀(𝑀2)−1.

rank𝑋=3. The nodes 1,2,3 of 𝐹 are distinct, since their differences are 1 and 2, both nonzero. The first three rows of 𝑋 form the matrix (111124139), which is invertible by Theorem 3.4.24, so its columns are independent. A relation among the columns of 𝑋 restricts to a relation among the columns of this 3×3 matrix, so the columns of 𝑋 are independent, and rank𝑋=3 (Corollary 3.5.16(a)).

Conclusion. By Theorem 4.3.8, rank𝑋𝖳=3, so 𝑇𝑋𝖳:𝐹𝑛→𝐹3 is surjective. Hence im𝐴=im(𝑋𝑀𝑋𝖳)=im(𝑋𝑀) (Proposition 2.6.15(b)), and rank(𝑋𝑀)=rank𝑋=3 (Proposition 4.3.1). So rank𝐴=3.

(c) Over 𝔽2, −2=0, and 𝑖2 has the same parity as 𝑖, so 𝐪=𝐰 and 𝑎𝑖𝑗=𝑖+𝑗. So 𝐴=𝐰𝟏𝖳+𝟏𝐰𝖳 has rank at most 2. Rows 1 and 2 are (0,1,0,1,…) and (1,0,1,0,…), which are independent, so rank𝐴≥2 by Theorem 4.3.8. Hence rank𝐴=2 over 𝔽2.

Exercise 4.3.41.standard
  1. (a)Let 𝐴∈𝐹𝑚×𝑛. Show that rank𝐴≤1 if and only if 𝑎𝑖𝑗𝑎𝑘𝑙=𝑎𝑖𝑙𝑎𝑘𝑗 for all 𝑖,𝑘∈{1,…,𝑚} and 𝑗,𝑙∈{1,…,𝑛}.
  2. (b)Find all 𝑡∈ℝ for which (1𝑡2𝑡42𝑡) has rank 1.

Solution.

Hint 1. For “only if”, write 𝐴=𝐮𝐯𝖳, so 𝑎𝑖𝑗=𝑢𝑖𝑣𝑗.
Hint 2. For “if”, assume 𝐴≠0 and fix an entry 𝑎𝑘𝑙≠0. Solve the condition for 𝑎𝑖𝑗.
Sketch. If 𝐴=𝐮𝐯𝖳, both sides are 𝑢𝑖𝑢𝑘𝑣𝑗𝑣𝑙. Conversely, if 𝑎𝑘𝑙≠0, the condition gives 𝑎𝑖𝑗=𝑎𝑖𝑙(𝑎𝑘𝑗/𝑎𝑘𝑙), so 𝐴=𝐮𝐯𝖳 with 𝐮 column 𝑙 and 𝑣𝑗=𝑎𝑘𝑗/𝑎𝑘𝑙. (b) 4−𝑡2=0 and 2𝑡2−8=0, and the remaining condition holds for all 𝑡: 𝑡=±2.

(a) Only if. If 𝐴=0, both sides are 0. If rank𝐴=1, then 𝐴=𝐮𝐯𝖳 (Proposition 4.3.16), so 𝑎𝑖𝑗=𝑢𝑖𝑣𝑗, and both sides equal 𝑢𝑖𝑢𝑘𝑣𝑗𝑣𝑙.

If. Suppose the condition holds. If 𝐴=0, then rank𝐴=0. Otherwise fix 𝑘,𝑙 with 𝑎𝑘𝑙≠0. For all 𝑖,𝑗, the condition gives

𝑎𝑖𝑗=𝑎𝑖𝑙⋅𝑎𝑘𝑗𝑎𝑘𝑙.

Let 𝐮=(𝑎1𝑙,…,𝑎𝑚𝑙), column 𝑙 of 𝐴, and 𝐯=(𝑎𝑘1,…,𝑎𝑘𝑛)/𝑎𝑘𝑙. Then 𝑎𝑖𝑗=𝑢𝑖𝑣𝑗 for all 𝑖,𝑗, so 𝐴=𝐮𝐯𝖳, and rank𝐴≤1 by Proposition 4.3.16 (in fact =1, since 𝑢𝑘=𝑎𝑘𝑙≠0 and 𝑣𝑙=1).

(b) The matrix is nonzero, so it has rank 1 exactly when all the conditions in (a) hold. With rows 1,2 and columns (1,2), (1,3), (2,3) they read 1⋅4=𝑡⋅𝑡, 1⋅2𝑡=2⋅𝑡, and 𝑡⋅2𝑡=2⋅4 (the other choices of indices repeat these or are trivial). The second always holds, and the first and third both say 𝑡2=4. So 𝑡=2 or 𝑡=−2, giving (122244) and (1−22−24−4).

Exercise 4.3.42.hard Let 𝐴∈𝐹𝑛×𝑛. Show that

rank𝐴+rank(𝐼𝑛−𝐴)=𝑛

if and only if 𝐴2=𝐴.

Solution.

Uses. Proposition 4.3.20 (Rank is subadditive), Exercise 4.3.30, Theorem 4.2.11 (Rank–nullity), Proposition 4.1.21 (Dimensions of subspaces)

Techniques. Dimension counting (Technique 4.2.24)

Hint 1. One inequality always holds: 𝐼𝑛=𝐴+(𝐼𝑛−𝐴).
Hint 2. If 𝐴2=𝐴, then 𝐴(𝐼𝑛−𝐴)=0. Which inequality bounds rank𝐴+rank(𝐼𝑛−𝐴) from above?
Hint 3. Conversely, show ker𝐴⊆im(𝐼𝑛−𝐴), compare dimensions, and conclude 𝐴(𝐼𝑛−𝐴)=0.
Sketch. Always 𝑛=rank𝐼𝑛≤rank𝐴+rank(𝐼𝑛−𝐴). If 𝐴2=𝐴, then 𝐴(𝐼𝑛−𝐴)=0, and Sylvester gives rank𝐴+rank(𝐼𝑛−𝐴)≤𝑛. Conversely, 𝐴𝐱=𝟎 implies 𝐱=(𝐼𝑛−𝐴)𝐱, so ker𝐴⊆im(𝐼𝑛−𝐴); the hypothesis says dimker𝐴=𝑛−rank𝐴=rank(𝐼𝑛−𝐴), so ker𝐴=im(𝐼𝑛−𝐴), that is, 𝐴(𝐼𝑛−𝐴)=0.

By Proposition 4.3.20, 𝑛=rank𝐼𝑛=rank(𝐴+(𝐼𝑛−𝐴))≤rank𝐴+rank(𝐼𝑛−𝐴) for every 𝐴.

If. Suppose 𝐴2=𝐴. Then 𝐴(𝐼𝑛−𝐴)=𝐴−𝐴2=0, so by Exercise 4.3.30(b), rank𝐴+rank(𝐼𝑛−𝐴)≤𝑛. With the inequality above, equality holds.

Only if. Suppose rank𝐴+rank(𝐼𝑛−𝐴)=𝑛. If 𝐴𝐱=𝟎, then 𝐱=𝐱−𝐴𝐱=(𝐼𝑛−𝐴)𝐱, so ker𝐴⊆im(𝐼𝑛−𝐴). By Theorem 4.2.11 and the hypothesis, dimker𝐴=𝑛−rank𝐴=rank(𝐼𝑛−𝐴)=dimim(𝐼𝑛−𝐴). By Proposition 4.1.21(b), ker𝐴=im(𝐼𝑛−𝐴). So 𝐴(𝐼𝑛−𝐴)𝐱=𝟎 for every 𝐱, that is, 𝐴−𝐴2=0.

So a matrix is idempotent exactly when the ranks of 𝐴 and 𝐼𝑛−𝐴 are as small as subadditivity allows. Geometrically, 𝐹𝑛 then splits into im𝐴, on which 𝐴 acts as the identity, and ker𝐴=im(𝐼𝑛−𝐴); this is the decomposition of Section 4.4.

Exercise 4.3.43.hard Let 𝐴∈𝐹𝑛×𝑛 satisfy 𝐴2=𝐴, and let 𝑟=rank𝐴≥1.

  1. (a)Let 𝐴=𝐶𝑅 be a factorization with 𝐶∈𝐹𝑛×𝑟 and 𝑅∈𝐹𝑟×𝑛, as in Theorem 4.3.21. Show that 𝑅𝐶=𝐼𝑟.
  2. (b)Deduce that tr𝐴=𝑟⋅1, the sum of 𝑟 copies of 1 in 𝐹.
  3. (c)Deduce that over ℚ, ℝ, or ℂ, the rank of an idempotent matrix equals its trace. Show by an example that over 𝔽2 the rank of an idempotent is not determined by its trace.

Solution.

Uses. Theorem 4.3.21 (Column–row factorization), Proposition 3.5.9 (Left inverses), Proposition 3.5.8 (Right inverses from the pivots), Proposition 2.4.28 (Cyclicity of the trace), Proposition 4.3.23 (Factorizations through 𝐹𝑘)

Techniques. Cyclicity of trace (Technique 2.4.29)

Hint 1. 𝐶 has independent columns and 𝑅 has independent rows. What one-sided inverses do they have?
Hint 2. Multiply 𝐶𝑅𝐶𝑅=𝐶𝑅 on the left by a left inverse of 𝐶 and on the right by a right inverse of 𝑅.
Hint 3. tr(𝐶𝑅)=tr(𝑅𝐶).
Sketch. rank𝐶=𝑟 = number of columns, so 𝐿𝐶=𝐼𝑟 for some 𝐿; rank𝑅=𝑟 = number of rows, so 𝑅𝑍=𝐼𝑟 for some 𝑍. Multiplying 𝐶(𝑅𝐶)𝑅=𝐶𝑅 by 𝐿 on the left and 𝑍 on the right gives 𝑅𝐶=𝐼𝑟. Then tr𝐴=tr(𝐶𝑅)=tr(𝑅𝐶)=tr𝐼𝑟=𝑟⋅1. Over 𝔽2, 𝐼2 and 0 both have trace 0.

(a) By Theorem 4.3.21, rank𝐶=𝑟, the number of columns of 𝐶, and rank𝑅=𝑟, the number of rows of 𝑅. By Proposition 3.5.9 and Proposition 3.5.8, there are 𝐿∈𝐹𝑟×𝑛 with 𝐿𝐶=𝐼𝑟 and 𝑍∈𝐹𝑛×𝑟 with 𝑅𝑍=𝐼𝑟. Since 𝐴2=𝐴, 𝐶(𝑅𝐶)𝑅=𝐶𝑅. Multiplying on the left by 𝐿 and on the right by 𝑍,

𝑅𝐶=(𝐿𝐶)(𝑅𝐶)(𝑅𝑍)=𝐿(𝐶𝑅𝐶𝑅)𝑍=𝐿(𝐶𝑅)𝑍=(𝐿𝐶)(𝑅𝑍)=𝐼𝑟.

(The same argument works for any factorization 𝐴=𝑋𝑌 with inner size 𝑟, by Proposition 4.3.23.)

(b) By Proposition 2.4.28, tr𝐴=tr(𝐶𝑅)=tr(𝑅𝐶)=tr𝐼𝑟=𝑟⋅1.

(c) In ℚ, ℝ, and ℂ, 𝑟⋅1 is the integer 𝑟, so tr𝐴=rank𝐴 by (b); the idempotent 𝐴=0 has trace and rank 0. Over 𝔽2, the idempotents 𝐼2 and 0 have trace 1+1=0 and 0, but ranks 2 and 0. In general, (b) determines rank𝐴 only modulo the characteristic of 𝐹. For example, the rank of the idempotent 𝐴=(1100) is 1=tr𝐴, which (b) predicts over any field.

Exercise 4.3.44.hard Let 𝐴∈𝐹𝑚×𝑛 have rank 𝑟≥1, and suppose 𝐴=𝐶𝑅=𝐶′𝑅′ with 𝐶,𝐶′∈𝐹𝑚×𝑟 and 𝑅,𝑅′∈𝐹𝑟×𝑛.

  1. (a)Show that there is exactly one invertible 𝑃∈𝐹𝑟×𝑟 with 𝐶′=𝐶𝑃, and that then 𝑅′=𝑃−1𝑅.
  2. (b)Deduce that the factorization of Theorem 4.3.21 is the only factorization 𝐴=𝐶′𝑅′ with inner size 𝑟 in which 𝑅′ is in RREF.

Solution.

Uses. Proposition 4.3.23 (Factorizations through 𝐹𝑘), Proposition 2.4.8 (Four views of the product), Proposition 3.5.9 (Left inverses), Proposition 3.5.8 (Right inverses from the pivots), Corollary 3.5.21 (Invertibility and row equivalence), Corollary 3.3.18 (RREF and row equivalence)

Techniques. Uniqueness via two candidates (Technique 1.1.8)

Hint 1. By Proposition 4.3.23(b), the columns of 𝐶 and of 𝐶′ are both bases of im𝐴. Express one basis in terms of the other.
Hint 2. 𝐶 has a left inverse, so 𝐶𝑋=𝐶𝑌 implies 𝑋=𝑌.
Hint 3. For (b), 𝑅′=𝑃−1𝑅 says that 𝑅′ and 𝑅 are row equivalent.
Sketch. The columns of 𝐶 form a basis of im𝐴 and the columns of 𝐶′ lie in it, so 𝐶′=𝐶𝑃 for a unique 𝑃; likewise 𝐶=𝐶′𝑄, and 𝐶𝑃𝑄=𝐶 gives 𝑃𝑄=𝐼𝑟 by cancelling 𝐶, so 𝑃 is invertible. Then 𝐶𝑅=𝐶𝑃𝑅′ gives 𝑅=𝑃𝑅′. (b) If 𝑅′ is in RREF, then 𝑅′=𝑃−1𝑅 is row equivalent to the RREF matrix 𝑅, so 𝑅′=𝑅, and then 𝑃=𝐼𝑟 and 𝐶′=𝐶.

By Proposition 4.3.23(b), the columns of 𝐶 form a basis of im𝐴, and so do the columns of 𝐶′. Since rank𝐶=𝑟 is the number of columns of 𝐶, 𝐶 has a left inverse 𝐿 (Proposition 3.5.9), so 𝐶𝑋=𝐶𝑌 implies 𝑋=𝐿𝐶𝑋=𝐿𝐶𝑌=𝑌. That is, 𝐶 can be cancelled on the left.

(a) Each column of 𝐶′ lies in im𝐴, so it is a combination of the columns of 𝐶; collecting the coefficients into the columns of 𝑃∈𝐹𝑟×𝑟 gives 𝐶′=𝐶𝑃 (Proposition 2.4.8(b)). Symmetrically, 𝐶=𝐶′𝑄 for some 𝑄. Then 𝐶(𝑃𝑄)=𝐶′𝑄=𝐶=𝐶𝐼𝑟, so 𝑃𝑄=𝐼𝑟 by cancellation. Likewise 𝐶′(𝑄𝑃)=𝐶′ and 𝐶′ can be cancelled, so 𝑄𝑃=𝐼𝑟. So 𝑃 is invertible. It is unique: if also 𝐶′=𝐶𝑃′, then 𝐶𝑃=𝐶𝑃′ gives 𝑃=𝑃′. Finally, 𝐶𝑅=𝐴=𝐶′𝑅′=𝐶𝑃𝑅′, so 𝑅=𝑃𝑅′ by cancellation, and 𝑅′=𝑃−1𝑅.

(b) Let 𝑅 be the first 𝑟 rows of rref(𝐴) and suppose 𝑅′ is in RREF. By (a), 𝑅′=𝑃−1𝑅 with 𝑃−1 invertible, so 𝑅′∼𝑅 (Corollary 3.5.21); both are in RREF, so 𝑅′=𝑅 by Corollary 3.3.18(c). Then 𝑅=𝑃𝑅, and since 𝑅 has a right inverse 𝑍 (Proposition 3.5.8), 𝐼𝑟=𝑅𝑍=𝑃𝑅𝑍=𝑃. So 𝐶′=𝐶. (Problem 3.17 proved the uniqueness of 𝐶 for this 𝑅 by the same cancellation.)

Exercise 4.3.45.challenge Let 𝐴∈𝐹𝑛×𝑛.

  1. (a)Suppose rank(𝐴2)=rank𝐴=𝑟≥1, and let 𝐴=𝐶𝑅 be a factorization with inner size 𝑟. Show that 𝑅𝐶∈𝐹𝑟×𝑟 is invertible.
  2. (b)Under the hypothesis of (a), show that 𝑋=𝐶(𝑅𝐶)−2𝑅 satisfies

    𝐴𝑋𝐴=𝐴,𝑋𝐴𝑋=𝑋,𝐴𝑋=𝑋𝐴.
  3. (c)Conversely, show that if some 𝑋∈𝐹𝑛×𝑛 satisfies the three equations in (b), then rank(𝐴2)=rank𝐴.

Solution.

Uses. Proposition 4.3.23 (Factorizations through 𝐹𝑘), Proposition 3.5.9 (Left inverses), Proposition 3.5.8 (Right inverses from the pivots), Theorem 4.3.3 (Rank of a product), Proposition 4.3.1 (Invertible factors preserve rank), Corollary 3.5.25 (The invertible matrix theorem)

Hint 1. 𝐴2=𝐶(𝑅𝐶)𝑅. Use one-sided inverses of 𝐶 and 𝑅 to show rank(𝑅𝐶)=rank(𝐴2).
Hint 2. For (b), compute 𝐴𝑋 and 𝑋𝐴 first; both equal 𝐶(𝑅𝐶)−1𝑅.
Hint 3. For (c), 𝐴=𝐴𝑋𝐴=𝐴2𝑋.
Sketch. (a) With 𝐿𝐶=𝐼𝑟 and 𝑅𝑍=𝐼𝑟, 𝑅𝐶=𝐿(𝐶𝑅𝐶𝑅)𝑍=𝐿𝐴2𝑍, so rank(𝑅𝐶)≥rank(𝐴2)=𝑟 and 𝑅𝐶 is invertible. (b) 𝐴𝑋=𝐶𝑅𝐶(𝑅𝐶)−2𝑅=𝐶(𝑅𝐶)−1𝑅 and 𝑋𝐴=𝐶(𝑅𝐶)−2𝑅𝐶𝑅=𝐶(𝑅𝐶)−1𝑅; then 𝐴𝑋𝐴=𝐶(𝑅𝐶)−1𝑅𝐶𝑅=𝐶𝑅 and 𝑋𝐴𝑋=𝐶(𝑅𝐶)−1𝑅𝐶(𝑅𝐶)−2𝑅=𝑋. (c) rank𝐴=rank(𝐴2𝑋)≤rank(𝐴2)≤rank𝐴.

By Proposition 4.3.23(b), 𝐶 has 𝑟 independent columns and 𝑅 has 𝑟 independent rows, so rank𝐶=rank𝑅=𝑟, and there are 𝐿 and 𝑍 with 𝐿𝐶=𝐼𝑟 and 𝑅𝑍=𝐼𝑟 (Proposition 3.5.9, Proposition 3.5.8).

(a) 𝐴2=𝐶𝑅𝐶𝑅, so 𝐿𝐴2𝑍=(𝐿𝐶)(𝑅𝐶)(𝑅𝑍)=𝑅𝐶. By Theorem 4.3.3, rank(𝑅𝐶)=rank(𝐿𝐴2𝑍)≤rank(𝐴2), and also rank(𝐴2)=rank(𝐶(𝑅𝐶)𝑅)≤rank(𝑅𝐶). So rank(𝑅𝐶)=rank(𝐴2)=𝑟, and the 𝑟×𝑟 matrix 𝑅𝐶 is invertible (Corollary 3.5.25).

(b) Let 𝐾=𝑅𝐶, invertible by (a), and 𝑋=𝐶𝐾−2𝑅. Then

𝐴𝑋=𝐶𝑅𝐶𝐾−2𝑅=𝐶𝐾𝐾−2𝑅=𝐶𝐾−1𝑅,𝑋𝐴=𝐶𝐾−2𝑅𝐶𝑅=𝐶𝐾−2𝐾𝑅=𝐶𝐾−1𝑅,

so 𝐴𝑋=𝑋𝐴. Next, 𝐴𝑋𝐴=𝐶𝐾−1𝑅𝐶𝑅=𝐶𝐾−1𝐾𝑅=𝐶𝑅=𝐴, and

𝑋𝐴𝑋=𝐶𝐾−1𝑅𝐶𝐾−2𝑅=𝐶𝐾−1𝐾𝐾−2𝑅=𝐶𝐾−2𝑅=𝑋.

(c) If 𝐴𝑋𝐴=𝐴 and 𝐴𝑋=𝑋𝐴, then 𝐴=𝐴(𝑋𝐴)=𝐴(𝐴𝑋)=𝐴2𝑋. By Theorem 4.3.3, rank𝐴=rank(𝐴2𝑋)≤rank(𝐴2)≤rank𝐴. So rank(𝐴2)=rank𝐴.

The matrix 𝑋 is called the group inverse of 𝐴. If 𝐴 is invertible, then 𝑟=𝑛, 𝐶 and 𝑅 are invertible, and 𝑋=𝐶(𝑅𝐶)−2𝑅=𝐶𝐶−1𝑅−1𝐶−1𝑅−1𝑅=𝐴−1. For the projection 𝐴=(1000)=𝐞1𝐞1𝖳, 𝑅𝐶=(1) and 𝑋=𝐴. For 𝑁=(0100), rank(𝑁2)=0<1, and no such 𝑋 exists.

More practice

The practice bank has 31 exercises on this section’s topics.