Similarity compares two matrices of one operator,
where the same basis is used for input and output, and
the last section showed that deciding it is hard. This
section asks what happens when the two bases may be
chosen independently, as they may for any map . The
answer is strikingly simple: the rank is then the only
thing that matters. On the way we meet the column–row
factorization, which writes every matrix as a product
through a space of dimension equal to its rank, and we
collect the basic inequalities for the rank of a
composition.
The
column–row factorization
In Chapter 3 we computed the rank of a matrix by row
reduction, and the reduction gave more than a number. It
produced a basis of the column space, the pivot columns
of (Theorem 3.5.2), and
it expressed every other column in terms of them, with
coefficients read off the reduced row echelon form.
Recording both pieces of information as matrices gives a
factorization.
Theorem
4.3.1 (Column–Row Factorization)
Let have rank , with pivot columns .
Let be the matrix whose columns are the pivot
columns of , and let
consist of the first rows of the reduced row
echelon form of ,
that is, its nonzero rows. Then Moreover, the columns of are linearly
independent and the rows of are linearly
independent.
Idea. Column of is times column of , that is, a
combination of the pivot columns of with the entries of
column of the
RREF as coefficients. For a pivot column the
coefficients are a standard basis vector, and for a free
column they are exactly the coefficients that Theorem 3.5.2 (b)
provides. So the factorization is that theorem written
as a matrix identity.
Proof. Let be the RREF of , with entries . By Theorem 3.5.5,
has exactly
pivots, so its
nonzero rows are rows . Fix a column index . By Theorem 0.9.3, column
of is Case 1: is a pivot column. By Lemma 3.3.3 (a),
column of is , so the
combination is .
Case 2: is
a free column. Let be the number of pivot
columns to the left of . For , the leading
entry of row of
lies in column
, so
. The
combination is therefore , which equals by Theorem 3.5.2
(b).
In both cases column of is column of ; hence . The columns of
form a basis of
by Theorem 3.5.2 (a),
so they are independent. The rows of are the nonzero rows
of , which form a
basis of
by Theorem 3.5.3, so
they are independent. This proves the theorem.
In words: every matrix of rank is a “tall” matrix with
independent columns times a “wide” matrix with
independent rows, and the middle size is exactly . The factorization
cannot be squeezed further.
Proposition
4.3.2 (Rank as the size of the thinnest
factorization)
Let be nonzero. Then is the
smallest integer such that for some
and .
Proof. Let ; then since . By Theorem 4.3.1, a
factorization with exists. Conversely, suppose with .
By Theorem 3.5.9,
, and because is spanned by
columns (Theorem 1.6.3).
Hence .
This gives a second view of “row rank equals column
rank”. Let , and let be the matrix whose
columns form a basis of . Each column
of is for its
coordinate column (Theorem 0.9.3), so
with
. Then every row of is a row of times , a combination of the
rows of (Theorem 0.9.2 (3) and
the row form of Theorem 0.9.3), so
. Applying the same argument to gives the reverse
inequality. In the language of maps, the factorization
says that can be done in two steps, :
first a surjection onto a space of dimension , then an injection.
Example
(A column–row factorization)
Find the column–row factorization of .
Solution. Row reduce. Subtracting
times row 1 from
row 2 and times
row 1 from row 3 gives rows , , .
Subtracting row 2 from row 3 gives the RREF The pivot columns are columns and , so , and by Theorem 4.3.1
Check by rows: row 2 of is , and row 3 is . The last column of says , and indeed .
Equivalent
matrices
Go back to the change-of-basis square for a map (Theorem 4.2.3):
When , nothing ties the basis of to the basis of . So the two outer
matrices are unrelated invertible
matrices, of sizes and . Any two matrices of the same map are related
in this way, and this recurring shape deserves a
name.
Two matrices are equivalent when they describe
the same linear map, with the input basis and the output
basis each changed freely.
Definition
4.3.1 (Equivalent matrices)
Let . We say is
equivalent to if there
exist invertible matrices and such that
In words: we may multiply on the left by
any invertible matrix and on the right
by any invertible matrix, and the two
need not be related. Both sizes are forced: acts on the rows and on the columns. This is an
equivalence relation on : ; if then ; and if also ,
then , where products of
invertible matrices are invertible by Theorem 0.9.7 (part
3).
Examples.
- Row equivalence. If is obtained from by row operations,
then with
invertible (Theorem 3.4.8), so
and are equivalent with
and . Column
operations correspond to multiplying on the right, and
they are allowed too.
- Similarity. If , then
and are equivalent, with
.
- The identity and a swap. and are equivalent, since
. They
are not similar, because is similar only to
itself.
- Degenerate case. The zero matrix
is
equivalent only to itself, since . It is the
matrix of the zero map in every pair of bases.
Nonexample by minimal change. Keep
but
drop the requirement that be invertible. Then
every matrix could be sent to the zero matrix with , and the relation
would not even be symmetric: , but for any . The clause “invertible” is what makes a change of basis
rather than an arbitrary composition.
Warning. Equivalent is much
weaker than similar. The matrices and are equivalent, since
and
is invertible,
but they are not similar (see the warning in the
previous section). For square matrices, “equivalent”
allows a different change of basis on each side, and
that freedom destroys almost all the structure an
operator has.
How much structure survives? The rank does, since
multiplying by invertible matrices does not change it
(Theorem 3.5.9). The
theorem below says that nothing else survives. The
target is the simplest matrix of each rank.
For , let be the block matrix with in
positions and
everywhere else. Its first columns are , and its remaining columns are zero.
Idea. Work backwards from the goal.
says: sends the first vectors of to the first vectors of , and sends the last
vectors of
to . So ① the last
vectors of
should be a
basis of ,
and by Rank–Nullity that is exactly the right number; ②
the first vectors
of extend it to
a basis of ; ③
their images are forced to be the first vectors of , and we must show
they are independent, which follows by counting; ④
extend them to a basis of . These are the bases of
the proof of Rank–Nullity, one on each side.
Proof.
- (a) By the
Rank–Nullity Theorem (Theorem 2.3.1), . Let
be a basis of , a subspace of
by Theorem 2.2.1. By the
Basis Extension Theorem (Theorem 1.6.5) there
are such that
is a basis of ;
the new vectors may be listed first, since reordering a
basis gives a basis. Put
for .
By Theorem 2.2.4, . Since
for ,
removing these vectors does not change the span, so
spans . This list has
length , so it is a basis of by Theorem 1.6.4 (b);
in particular it is linearly independent in . By the Basis Extension
Theorem there are such that is a basis of
.
Now read off the columns of (Definition 4.1.1).
For , ,
whose -coordinate vector is
. For
, ,
whose coordinate vector is . These are the
columns of .
(b) Let ,
, and let , be the standard
bases. The -th
column of is
, the -th column of (Theorem 0.9.3), so
. By Theorem 4.1.5,
. Take , as in (a). By Theorem 4.2.3,
and both outer matrices are invertible by Theorem 4.2.1
(c). Hence is
equivalent to .
(c) If with invertible, then
by Theorem 3.5.9.
If
, then by (b) both and are equivalent to
. Since
equivalence is symmetric and transitive, is equivalent to . This proves the
theorem.
So the equivalence classes of are
exactly the sets of matrices of rank , for : there are of them,
whatever the field. Every linear map between
finite-dimensional spaces looks like “keep coordinates and delete
the rest” once the bases on both sides are chosen well.
The right pair of bases makes any single map
trivial.
The proof is also an algorithm. To find and for a concrete matrix,
compute a basis of , extend it, and
take images.
Solution. Subtracting times row 1 from row 2
gives , so , and the solutions of
satisfy ,
. Hence
, as Rank–Nullity predicts ().
Following the proof, let ,
and extend with , .
The list is independent: in
, the
second and third entries give , and then . So it is a basis
of . The images are
and , which already form a basis of since . Therefore
using Theorem 4.2.1 (c)
and Theorem 0.9.8. Check:
(the columns are , , ), and The pair is far from unique; any choice of extension
gives another.
Quick check. How many equivalence
classes does have, and which is the class of ?
Solution. Ranks in range
over , so
there are classes by Theorem 4.3.3
(c). The given matrix is nonzero and its second column
is twice the first, so its rank is , and it lies in the
class of : all
matrices of rank .
Rank
inequalities for compositions
The rank of a composition cannot exceed the rank of
either factor, and it cannot fall too far either. Both
facts come from the same observation: is applied to the subspace
. For a
subspace and a linear map , the
restriction ,
, is linear, since its values are
those of . Its
kernel is .
Theorem
4.3.4 (Rank of a Composition)
Let and be linear, with and finite-dimensional.
Then
Proof. Consider the restriction
. For , , so ; conversely every with equals . Hence
. Since is finite-dimensional (a subspace of , Theorem 1.6.7),
Rank–Nullity (Theorem 2.3.1) applied
to gives
Also , and is
finite-dimensional by Rank–Nullity for , so by
Theorem 1.6.7. This
proves the theorem.
The lower bound uses the same restriction, but now we
control its kernel instead of ignoring it.
Theorem
4.3.5 (Sylvester’s Rank Inequality)
Let and be linear, with and finite-dimensional.
Then In matrix form: for
and , .
Idea. Of the dimensions that
reach , the map
can destroy only
those lying in , and has dimension .
Rank–Nullity for the restriction turns this
sentence into an inequality.
Proof. Let , so
that as in the proof of Theorem 4.3.4,
and . By Rank–Nullity
for and then
for , where the inequality is Theorem 1.6.7.
Rearranging gives the claim. For matrices, apply it to
on and on : the matrix of in standard bases is
by Theorem 4.1.2,
and ranks of maps and matrices agree by Theorem 4.1.5.
If you worked Exercise (C1: A Lower Bound for the Rank of a Composite),
you have seen this proof already; we record the result
here as a theorem so that later arguments can cite it. A
typical use: if with and , then , that is, For instance, there are no of
rank with .
Warning. The rank of a
product depends on the order of the factors.
For and , has rank but has rank . Theorem 3.5.9 and
Sylvester’s inequality bound from both
sides; they do not determine it, and they say nothing
that would make .
Exercises
A. Check your
understanding
Exercise
(A1)
- (a) Define what it
means for to be equivalent.
- (b) State the Rank
Normal Form theorem for a matrix .
- (c) True or false:
equivalent square matrices are similar. Justify your
answer.
- (d) True or false:
for all .
Justify your answer.
- (e) How many
equivalence classes does
have?
- (f) In the proof of
the normal form for , which vectors form the end of the
basis of , and
which vectors form the start of the basis of ?
Solution.
(a) There are
invertible and with (Definition 4.3.1).
(b) If , there are
invertible
with ; and two matrices are
equivalent if and only if they have the same rank (Theorem 4.3.3).
(c) False. and both have rank , so they are
equivalent, but is similar only to
itself.
(d) False. With
and , and .
(e) The possible
ranks are ,
so there are
classes. The field plays no role.
(f) The last
vectors of
the basis of form
a basis of .
The first vectors
of the basis of
are the images of the first basis vectors of , and they form a basis
of .
B. Practice
Exercise
(B1: A column–row factorization)
Find the column–row factorization of Hence write the second and fourth columns of
as combinations
of its pivot columns.
Solution. Subtract times row 1 from row 2
and row 1 from row 3: the rows become , , . Subtract
row 2 from row 3, and then times row 2 from row 1:
The pivot columns are and , so , and by Theorem 4.3.1
Check by rows: and ; row 1 is . Reading columns and of :
and . Indeed .
Exercise
(B2: Normal form of a rank-one matrix)
Let . Find invertible and such that
, following the proof of Theorem 4.3.3.
Solution. Row 2 is times row 1, and , so . The null
space is ; setting the free variables to and gives the basis
, . Extend with :
in ,
the last two entries give and then , so is a basis
of . Its first
image is , which we extend to the
basis of (independent, since
forces ). Then Check: has columns ,
, ,
and . So . Both and are invertible by Theorem 4.2.1
(c).
Exercise
(B3: Rank-one pieces)
Let .
- (a) Use Theorem 4.3.1 to
show that if , then for some nonzero and
.
- (b) Show that if
,
then is a sum of
matrices of rank
, and that it is
not a sum of fewer than matrices of rank .
- (c) Write as a sum of two
rank-one matrices.
Solution.
(a) With ,
is a single column , and
is a single row, which we write as with . So
. Both are nonzero: the
columns of and
the rows of are
linearly independent by Theorem 4.3.1,
and a list containing a zero vector is dependent. (This
recovers one direction of Exercise (C2: Matrices of rank one).)
(b) Write with columns
of and rows of , and let , be the entries of
, . The matrix has
-entry , so by Theorem 0.9.2 (1),
Each has
rank by the direction of
Exercise (C2: Matrices of rank one),
since (a column of , whose columns are
independent) and (a row of , whose rows are
independent). If with each , then by
Exercise (C1: Rank of a sum)
applied repeatedly, .
(c) By Example (A column–row factorization),
with
and
, so
C. Going deeper
Exercise
(C1: Frobenius’s rank inequality)
Let , , and .
Write , and let be .
- (a) Show that , and deduce from
Rank–Nullity for that .
- (b) Show similarly
that .
- (c) Prove that , and deduce Frobenius’s
inequality
- (d) Deduce
Sylvester’s inequality for matrices from (c).
Hint: for (d), look for a special choice of the
middle factor.
Solution.
(a) By Theorem 0.9.2, every
vector
equals
with ; conversely every is
for some
, and then
. So . The kernel of is , and
is
finite-dimensional, so by Theorem 2.3.1,
(b) The same
argument with
in place of
gives , with kernel , so
.
(c) Every vector
lies in , so . Hence , and
by Theorem 1.6.7, . By (b) and then (a), Rearranging gives Frobenius’s
inequality.
(d) Let
and , and apply (c) with and . Then , since . This is
Theorem 4.3.5 for
matrices.
Exercise
(C2: Counting over )
- (a) Show that contains
exactly matrix of
rank , of rank , and of rank .
- (b) Deduce the sizes
of the equivalence classes of .
- (c) Compare with
similarity: how many matrices of are similar to
? Is the class
of under
equivalence a union of several similarity classes?
Solution.
(a) The only
matrix of rank is
. A matrix has
rank exactly when
its two columns span , that is (Theorem 1.6.4), when
they form a basis of . The first column
can be any of the
nonzero vectors of . The second must
lie outside the span of the first, which has elements, so there are
choices.
That gives matrices of rank . Since , the remaining matrices
have rank .
(b) By Theorem 4.3.3
(c), the classes are the sets of matrices of a fixed
rank, of sizes ,
, and .
(c) is similar only to
itself, so its similarity class has element, while its
equivalence class, the invertible matrices, has . Since similar matrices
are equivalent, each equivalence class is a union of
similarity classes, and here the class of the invertible matrices
contains the similarity class together with
others. For example is invertible but not
similar to , so
there are at least two similarity classes
inside.