Chapter 12 §11 gave the Frobenius norm its own symbol
and explained the subscript: a bare already
meant the norm of a vector, and a matrix, it
warned, would shortly be asked to carry both kinds of
measurement. This is that moment. A matrix is a list of
numbers, and it
is also a map; measuring the list is easy, measuring the
map is what we usually want, and the two answers are not
the same number. This section builds the second
measurement, computes it in the three cases anyone
computes, and says exactly where the first one
fails.
Throughout, is or , all spaces
are finite-dimensional, and with no
subscript is a norm in the sense of Definition 16.1.1.
What a
matrix norm has to do
Chapter 10 §10 decided whether by looking
at one entry of one Jordan block at a time. A single
number attached to would be a far better
instrument — but only if it interacts with matrix
multiplication, because the questions we want to answer
are questions about products: , , . A norm that
satisfies
immediately gives , and with that one line a
geometric series is in reach. A norm that does not
satisfy it gives nothing at all about powers. So the
property is not a technical nicety; it is the entire
reason for the definition.
A matrix norm is a norm on matrices that never
lets a product be bigger than the product of the
sizes.
Definition
16.3.1 (Matrix Norm)
Let . A
matrix norm on is a function
which is a norm on as a vector space
(Definition 16.1.1) and which
is in addition submultiplicative: for
all ,
In words: the first requirement is the ordinary one —
clauses (N1), (N2), (N3) of Definition 16.1.1 — and it
sees only as
the vector space in disguise. The
second requirement is the one that notices the
multiplication. Note that submultiplicativity is an
inequality, not an equality: it is easy for a norm to
make
far smaller than , and
impossible for a matrix norm to make it larger.
The condition needs a square matrix only so that
is defined for
every pair; nothing else in the definition cares about
the shape.
Example
(The Frobenius norm is submultiplicative)
Show that for all , so
that
is a matrix norm.
Solution. The entry is the inner product in of the -th row of with the -th column of (up to a conjugation,
which does not change moduli). By Theorem 11.1.5,
Summing over all and , and noting that the
two bracketed factors depend only on and only on respectively,
Both sides are non-negative, so taking square roots
gives the claim. That is a norm
was settled in Chapter 11: it is the norm induced by the
Frobenius inner product (Definition 11.1.3),
so Theorem 11.1.2 and
Corollary 11.1.6
apply.
Now the minimal change that breaks it. Keep the idea
“read the matrix as a list of numbers” but measure the
list with the maximum instead of the Euclidean
length.
Example
(The entrywise maximum is not submultiplicative)
For put . Show that is a norm on but not a matrix
norm for .
Solution. is the -norm of the list
of entries, so
it is a norm (Definition 16.1.1): it
vanishes only on the zero matrix, , and (N3) holds entry by entry.
Submultiplicativity fails already for . Take . Then , while
The mechanism is plain: an entry of a product is a
sum of terms,
each as large as the largest entry, and a maximum taken
entry by entry has no way of charging for the length of
that sum. For the same matrix gives .
Warning.
Submultiplicativity is a property of the norm,
not of matrices. The failure above is not
caused by being
special; it is caused by being the wrong
yardstick. Exercise (C2: Repairing the entrywise maximum)
shows that the single rescaling repairs
it, which is a good way to see that nothing deep is
being violated — only a normalization.
The norm a
vector norm induces
The definition above is a filter, not a construction:
it says which norms are usable, not where to find one.
The construction comes from remembering what is for. It sends
vectors to vectors, and if both are already measured, we
can ask how much
magnifies.
This is exactly the question Chapter 13 §08 asked of
the singular value decomposition, and the question
Chapter 13 §12 asked of an oblique projection, in each
case computing and then promising that the number would
get a name here.
Before naming it, we must know that the maximum
exists. “Take the largest value of a function over
infinitely many vectors” is not something one may write
without a reason, and the reason is compactness, already
packaged for us in Corollary 16.2.5.
Lemma
16.3.1 (The Largest Stretch Is Attained)
Let ,
let be
a norm on and
also a norm on ,
and let . Then the function
attains a maximum on the unit sphere .
Idea. Two steps, and the second is
already packaged for us. ① The function is Lipschitz,
hence continuous: Lemma 16.1.1 converts a
difference of norms into the norm of a difference, and a
crude expansion in the standard basis bounds by a constant
times . ②
The unit sphere is closed and bounded, hence compact,
and a continuous real function on it attains its bounds
— which is exactly what Corollary 16.2.5 says.
Proof. Step 1: a Lipschitz
bound. Write . By subadditivity and
homogeneity,
By Theorem 16.2.2, applied on
to the norms
and
, there
is a constant with for every . Hence with , a constant depending on
but not on .
Step 2: continuity. Let . By Lemma 16.1.1, applied to the
vectors and
of , and then by Step
1,
So is Lipschitz with constant , and in particular
continuous.
Step 3: compactness. The set is non-empty, since
and ,
the division being legal because . By Corollary 16.2.5 the unit
sphere of any norm on a finite-dimensional space is
compact, and a continuous real-valued function on it
attains a maximum. That corollary is where this
chapter’s imported facts (A3), the compactness of closed
bounded sets, and (A4), the extreme value theorem, are
cashed in; here we only apply it, to the continuous
function of Step 2. This proves the lemma.
The operator norm of a matrix is the worst
stretching it can do to a unit vector.
Definition
16.3.2 (Operator Norm)
Let
and let a norm be given on
and on ; we write both with
the same symbol, the argument saying which is meant. The
operator norm of
induced by those norms is
which exists by Lemma 16.3.1. When the
vector norm is on both
sides we write , and is also
called the spectral norm.
Clause by clause: the maximum is over
unit vectors only, so a matrix that
stretches one long vector enormously is not thereby
given a large norm — what counts is the factor.
The norm on the source and the norm on the target may be
different, and then depends on
both; the notation hides that, which is why the vector
norm must always be named. And the word
induced is doing real work: not every
matrix norm arises this way, as the end of this section
shows.
Three equivalent descriptions are worth having at
once. For ,
The middle expression equals the first because by homogeneity, and
runs
over as runs over the non-zero
vectors. The third equals the first because scaling any
with up
to the sphere only increases , and contributes .
The smallest cases are worth seeing before any real
example, because they pin the scale. For every is , so , as (N1)
demands. For , a matrix is a scalar , the only unit
vectors are the scalars of modulus , and : the operator norm of a number is its
absolute value, so the definition extends the
one-dimensional case rather than replacing it. And for
, every unit
vector satisfies , so for
every choice of vector norm — a fact used so
often below that it is worth fixing now.
Quick check. Why is the maximum in
Definition 16.3.2 taken over
rather than over all ? What would the
supremum over all be?
Solution. For a non-zero it would be infinite:
pick with
; then
as . Restricting
to unit vectors measures the ratio of output
size to input size, which is what “how much does magnify?” means, and
by homogeneity the ratio is already known on all of
once it is
known on the sphere.
Theorem
16.3.2 (Properties of the Operator Norm)
Let
be a norm on each of , , , and let .
Then:
- (a) for every
, and
is the
smallest constant with this property;
- (b) is a norm on
the vector space ;
- (c) ;
- (d) for every . In particular the operator norm
induced on
by a single norm on is a matrix norm in
the sense of Definition 16.3.1.
Idea. Everything is the definition
plus homogeneity. The one move worth naming is in (a): a
general
is times
a unit vector, so the inequality on the sphere upgrades
to an inequality everywhere; and once (a) is available,
(d) is two applications of it.
Proof.
(a) For both sides are
. For , the vector
satisfies , so by Definition 16.3.2;
multiplying by and using
homogeneity gives . If is any constant with
for all , then taking a unit
at which the
maximum is attained (Lemma 16.3.1) gives .
(b) Each value is
a maximum of non-negative numbers, so . If then for every unit
, hence for every by (a), hence . For a scalar
, . For subadditivity, let be a unit vector; then
, and taking the maximum over gives .
(c) for every unit , so the maximum is
.
(d) Let be a unit
vector. By (a) applied twice, first to and then to ,
Taking the maximum over unit gives . This proves the
theorem.
Part (c) is small and will do a large amount of work:
it is the single property that separates the operator
norms from the rest of the matrix norms.
A definition by maximization is not yet a
computation. For the three standard vector norms the
maximization can be carried out completely, and in two
of the three cases the answer is read straight off the
entries.
Idea. Each part is an upper bound
and a witness. ① For the upper
bound is subadditivity of the modulus with the column
sums pulled out, and the witness is a standard basis
vector
selecting the heaviest column. ② For the
upper bound is the same estimate read along a row, and
the witness is a vector of unimodular entries chosen to
rotate every term of the heaviest row into the positive
reals, so that no cancellation occurs. ③ For the
witness is the first right singular vector, and the
upper bound is the singular value decomposition: unitary
factors are invisible to the Euclidean norm, so only the
diagonal matrix is left, and a
diagonal matrix stretches most along its largest entry,
because it scales each coordinate separately.
Proof. Write and .
(a). Let . Then
where the first inequality is subadditivity of the
modulus together with an exchange of two finite sums. So
.
For the reverse, choose with . The vector has and
is the -th column of , so .
Hence .
(b). Let .
For each ,
so and .
For the reverse, choose with and define by when and otherwise. Every
has modulus
, so ;
and for every , the choice being legal
exactly because we checked before
dividing. Hence the -th entry of is , giving .
(c). Let be a singular value
decomposition (Theorem 13.8.1), with
unitary and
carrying
on its diagonal, . Let and put
. A
unitary matrix is an isometry for (Theorem 11.7.3, (a)
and (f)), so and .
The vector has entries
for
and zeros
below, so
Hence . For the reverse, take , the first
column of ; then
, so , and because
the columns of a unitary matrix are orthonormal. This
proves the theorem.
Part (c) is the promise this section was written to
keep. Chapter 13 §08 computed and said that Chapter 16 would
give that quantity a name, the operator norm , “and this
computation is why ”.
The name is Definition 16.3.2; the
reason has now been separated from the name, so that the
equality is a theorem about two independently defined
numbers rather than a definition dressed up as one.
Chapter 13 §12 measured an oblique projection by the
same maximum and deferred the word twice; Example (The oblique projection, finally measured)
below supplies it. And Chapter 12 §11, when it wrote
with an
obligatory subscript, explained that the bare symbol was
reserved because a matrix would soon be asked to be
measured both as a list and as a map: parts (a)–(c) are
that second measurement, and the rest of this section is
the difference between the two.
Example
(One matrix, four numbers)
For , compute , , and .
Solution. The column sums of
absolute values are and , so by Theorem 16.3.3 (a). The row
sums are and
, so by
(b). The Frobenius norm is .
For the spectral norm, form
This matrix is , whose eigenvalues are
, namely
and . So the singular values
are and , and
by (c).
Two checks come for free. The product
equals , and . The four
numbers , , and
are all different, which is the point:
“the size of ” is
not a well-posed phrase until the norm is named.
Example
(The oblique projection, finally measured)
Chapter 13 §12 studied a projection of onto a line, with
matrix in a
suitable orthonormal basis, and found that the largest
value of
over unit is
, which over is for the
angle
between image and kernel. Confirm this with Theorem 16.3.3 (c) for , the matrix of Example (Two Oblique Projections, Measured).
Solution. , with eigenvalues and ; so and
. This agrees with , and with for , the angle
between the two lines. An orthogonal projection, by
contrast, has whenever , since it fixes every unit vector of
and shortens the
rest.
The two matrix norms now in play compare in a way
worth recording, because the Frobenius norm is far
cheaper to compute and the spectral norm is usually the
one wanted.
Proposition
16.3.4 (Frobenius Against Spectral)
Let be non-zero, of rank . Then
and the left inequality is an equality if and only if
.
Idea. Write both norms in terms of
the singular values:
and . The claim is
then the statement that one term of a sum of non-negative numbers is
at most the sum, which is at most times the largest — and
that the first is an equality exactly when there is only
one term.
Proof. By Lemma 11.7.5 and
(Theorem 13.8.1), , so . Exactly of the are non-zero (Theorem 13.8.1), and
each satisfies . Hence
which is the displayed chain after taking square
roots and using (Theorem 16.3.3 (c)). The
left inequality is an equality exactly when for all
, that is,
when .
So the two norms never differ by a factor worse than
, and they agree exactly on the
rank-one matrices. What they do not agree on is
more interesting.
What
submultiplicativity does not buy
The Frobenius norm is a matrix norm: Example (The Frobenius norm is submultiplicative)
proved it submultiplicative. It is tempting to conclude
that it is therefore the operator norm induced by some
vector norm — perhaps a cleverly chosen one. It is not,
and the obstruction is a single line.
Warning. Submultiplicative
does not mean induced. Every operator norm
satisfies by Theorem 16.3.2 (c),
whereas
So for
the Frobenius norm is not induced by
any norm on
whatsoever — not by a bad choice, but by no choice at
all. It is a matrix norm and not an operator norm, and
the two phrases must be kept apart.
The identity is the cheapest witness, and it is
decisive because is forced: the identity stretches nothing, so
a norm that measures stretching must give it . A norm that assigns
the identity the value is measuring
something else — the amount of material in the matrix,
which for
grows with .
Warning. Two different
objects are called in the
wild. In this book is the
operator norm induced by on
vectors, that is, the largest row sum. The entrywise
maximum of Example (The entrywise maximum is not submultiplicative)
is a different function, and it is not even
submultiplicative. Whenever a source writes , check
which one is meant by evaluating it at : the row-sum norm
gives , the
entrywise maximum gives as well — so evaluate
instead at ,
where the row-sum norm gives and the entrywise
maximum gives .
Solution. No. Being
compatible with a vector norm — satisfying
— is weaker than being
induced by it. The induced norm is the
smallest constant that works (Theorem 16.3.2 (a)), and for
with
the
smallest constant is while .
Compatibility says the number is big enough; being
induced says it is exactly right.
Where this leaves us: we have a supply of matrix
norms (Definition 16.3.1), among
them a distinguished family, the induced ones (Definition 16.3.2), and
inside that family three that are computable (Theorem 16.3.3). The next
section asks which quantity all of these norms are
simultaneously bounding from above, and finds that the
answer is the spectral radius.
Exercises
A. Check your
understanding
Exercise
(A1)
- (a) Define a
matrix norm on , and define the
operator norm induced by a vector
norm.
- (b) State the three
formulas of Theorem 16.3.3.
- (c) Which analysis
facts make Definition 16.3.2
legitimate, and where is each used?
- (d) Determine
whether the following statement is correct, and justify
your answer: every matrix norm on is the operator
norm induced by some norm on .
- (e) Determine
whether the following statement is correct, and justify
your answer: for every and every operator norm, .
Solution.
(a) A matrix norm
is a norm on the vector space that also
satisfies for all (Definition 16.3.1). Given
norms on and
, the operator
norm of is (Definition 16.3.2).
(b) is the
largest absolute column sum, the
largest absolute row sum, and .
(c) Two: the
compactness of closed bounded sets, fact (A3) of the
chapter introduction, used for the unit sphere ; and the extreme value
theorem, fact (A4), used for .
Both enter through Corollary 16.2.5, quoted in
Step 3 of Lemma 16.3.1. The continuity
needed for the second is proved, not imported, and rests
on Theorem 16.2.2 and Lemma 16.1.1.
(d) Incorrect.
The Frobenius norm is a matrix norm, and for , whereas every operator norm has by Theorem 16.3.2 (c).
(e) Incorrect;
only holds, by
Theorem 16.3.2 (d). Take
, for
which but , so .
B. Practice
Exercise
(B1: Four numbers for one matrix)
Let . Compute , , and , and verify
that .
Solution. Absolute column sums:
and , so . Absolute
row sums:
and , so
. Also .
For the spectral norm,
with trace
and determinant , so its eigenvalues solve , giving and . Hence , and by Theorem 16.3.3 (c). Finally
, as
Proposition 16.3.4 requires;
the inequality is strict because .
Exercise
(B2: A rank-one matrix)
Let
and be
non-zero and put .
- (a) Prove that .
- (b) Verify the
formula for and in by computing
directly.
Solution.
(a) For we have , so
. By Theorem 11.1.5,
, with equality when
. Hence for unit the quantity is at most
and attains that value, so .
(b) Here and , of
trace and
determinant . So the eigenvalues are and , giving . And , in agreement. Note also , as Proposition 16.3.4 predicts
for a rank-one matrix.
Exercise
(B3: The star trick)
Let . Prove that . Hint: the eigenvalues of
.
Solution. By Definition 13.8.1 the
eigenvalues of are , where the are the singular
values of padded
with zeros if . The matrix is positive
semidefinite, hence Hermitian, so by Corollary 12.4.3
it is
with unitary and
. Then , whose eigenvalues are , so the
singular values of are the and its
largest is . By Theorem 16.3.3 (c), applied
to and
then to ,
This proves the claim. Notice that it is an equality,
whereas Theorem 16.3.2 (d) alone
would give only .
C. Going deeper
Exercise
(C1: The spectral norm between the other two)
Let .
Hint for (a): apply Theorem 11.1.5 to
the splitting inside one
row.
Solution.
- (a) Write and
for the largest column and row sums (Theorem 16.3.3). Let . For each
, Theorem 11.1.5
applied to the vectors with entries
and gives
Summing over
and exchanging the two finite sums,
Taking the maximum over unit gives , as
claimed.
- (b) There and ,
so the bound is , while . The bound overestimates by about 18 per
cent. It is sharp in general: for all three
numbers are .
Exercise
(C2: Repairing the entrywise maximum)
For let , as in Example (The entrywise maximum is not submultiplicative).
- (a) Prove that ,
and deduce that is a matrix norm.
- (b) Deduce that
is not an
operator norm for any .
- (c) Explain why (a)
and (b) together show that submultiplicativity and are
independent requirements.
Solution.
- (a) For all ,
so . Put . It is
a norm, being a positive multiple of the norm , and
(b) for , while every operator norm gives the identity
the value (Theorem 16.3.2
(c)).
(c) itself is a norm with
that is
not submultiplicative; is submultiplicative
with .
So neither property implies the other, and an operator
norm is precisely a matrix norm that also passes the
identity test — a necessary condition which, as already
showed, is still not sufficient on its own for large.
Exercise
(C3: When is a bound an equality?)
Let be non-zero.
- (a) Prove that
if and only if .
- (b) Give an example
with
where , the worst case allowed by
Proposition 16.3.4.
- (c) Is there a
non-zero with
? Justify your answer.
Solution.
(a) This is the
equality case of Proposition 16.3.4, proved
there: and , so the two agree exactly when for every
, which by
Theorem 13.8.1 says
.
(b) Take . Its singular
values are , so
and
. Any matrix with two equal non-zero
singular values does the same, for instance any unitary
matrix.
(c) No. always, since every term of the sum
is non-negative. The Frobenius norm counts all the
singular values and the spectral norm only the largest,
so the Frobenius norm can never be the smaller of the
two.