16.3  Measuring a Matrix

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 mn 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, F is R or C, all spaces are finite-dimensional, and ‖⋅‖ with no subscript is a norm in the sense of Definition 16.1.1.

16.3.1. What a matrix norm has to do

Chapter 10 §10 decided whether Am→0 by looking at one entry of one Jordan block at a time. A single number attached to A 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: Am, ∑kAk, eA. A norm that satisfies

‖AB‖≤‖A‖‖B‖

immediately gives ‖Am‖≤‖A‖m, 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 n≥1. A matrix norm on Mn(F) is a function ‖⋅‖:Mn(F)→R which is a norm on Mn(F) as a vector space (Definition 16.1.1) and which is in addition submultiplicative: for all A,B∈Mn(F),

‖AB‖≤‖A‖‖B‖.

In words: the first requirement is the ordinary one — clauses (N1), (N2), (N3) of Definition 16.1.1 — and it sees Mn(F) only as the vector space Fn2 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 ‖AB‖ far smaller than ‖A‖‖B‖, and impossible for a matrix norm to make it larger.

The condition needs a square matrix only so that AB is defined for every pair; nothing else in the definition cares about the shape.

Example (The Frobenius norm is submultiplicative)

Show that ‖AB‖F≤‖A‖F‖B‖F for all A,B∈Mn(F), so that ‖⋅‖F is a matrix norm.

Solution. The entry (AB)ij=∑kaikbkj is the inner product in Fn of the i-th row of A with the j-th column of B (up to a conjugation, which does not change moduli). By Theorem 11.1.5,

|(AB)ij|2≤(∑k|aik|2)(∑k|bkj|2).

Summing over all i and j, and noting that the two bracketed factors depend only on i and only on j respectively,

‖AB‖F2≤(∑i,k|aik|2)(∑k,j|bkj|2)=‖A‖F2‖B‖F2.

Both sides are non-negative, so taking square roots gives the claim. That ‖⋅‖F 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 A∈Mn(F) put N(A)=maxi,j|aij|. Show that N is a norm on Mn(F) but not a matrix norm for n≥2.

Solution. N is the ∞-norm of the list of n2 entries, so it is a norm (Definition 16.1.1): it vanishes only on the zero matrix, N(cA)=|c|N(A), and (N3) holds entry by entry.

Submultiplicativity fails already for n=2. Take A=B=J=(1111). Then N(J)=1, while

J2=(2222),N(J2)=2>1=N(J)N(J).

The mechanism is plain: an entry of a product is a sum of n 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 n×n the same matrix gives N(J2)=n.

Warning. Submultiplicativity is a property of the norm, not of matrices. The failure above is not caused by J being special; it is caused by N being the wrong yardstick. Exercise (C2: Repairing the entrywise maximum) shows that the single rescaling A↦nN(A) repairs it, which is a good way to see that nothing deep is being violated — only a normalization.

16.3.2. 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 A is for. It sends vectors to vectors, and if both are already measured, we can ask how much A 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 max{‖Ax‖:‖x‖=1} 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 n,m≥1, let ‖⋅‖ be a norm on Fn and also a norm on Fm, and let A∈Mm×n(F). Then the function x↦‖Ax‖ attains a maximum on the unit sphere S={x∈Fn:‖x‖=1}.

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 ‖Az‖ by a constant times ‖z‖. ② 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 z=∑j=1nzjej. By subadditivity and homogeneity,

‖Az‖≤∑j=1n|zj|‖Aej‖≤(maxj‖Aej‖)‖z‖1.

By Theorem 16.2.2, applied on Fn to the norms ‖⋅‖1 and ‖⋅‖, there is a constant c>0 with ‖z‖1≤c‖z‖ for every z. Hence ‖Az‖≤M‖z‖ with M=cmaxj‖Aej‖, a constant depending on A but not on z.

Step 2: continuity. Let x,y∈Fn. By Lemma 16.1.1, applied to the vectors Ax and Ay of Fm, and then by Step 1,

|‖Ax‖−‖Ay‖|≤‖A(x−y)‖≤M‖x−y‖.

So x↦‖Ax‖ is Lipschitz with constant M, and in particular continuous.

Step 3: compactness. The set S is non-empty, since n≥1 and e1/‖e1‖∈S, the division being legal because e1≠0. 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 n,m≥1 and let a norm ‖⋅‖ be given on Fn and on Fm; we write both with the same symbol, the argument saying which is meant. The operator norm of A∈Mm×n(F) induced by those norms is

‖A‖:=max{‖Ax‖:x∈Fn, ‖x‖=1},

which exists by Lemma 16.3.1. When the vector norm is ‖⋅‖p on both sides we write ‖A‖p, and ‖A‖2 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 ‖A‖ 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 A≠0,

‖A‖=max‖x‖=1‖Ax‖=maxx≠0‖Ax‖‖x‖=max‖x‖≤1‖Ax‖.

The middle expression equals the first because ‖Ax‖/‖x‖=‖A(x/‖x‖)‖ by homogeneity, and x/‖x‖ runs over S as x runs over the non-zero vectors. The third equals the first because scaling any x with 0<‖x‖≤1 up to the sphere only increases ‖Ax‖, and x=0 contributes 0.

The smallest cases are worth seeing before any real example, because they pin the scale. For A=0 every Ax is 0, so ‖0‖=0, as (N1) demands. For m=n=1, a matrix is a scalar (a), the only unit vectors are the scalars of modulus 1, and ‖(a)‖=|a|: the operator norm of a number is its absolute value, so the definition extends the one-dimensional case rather than replacing it. And for A=I, every unit vector satisfies ‖Ix‖=‖x‖=1, so ‖I‖=1 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 ‖x‖=1 rather than over all x? What would the supremum over all x∈Fn be?

Solution. For a non-zero A it would be infinite: pick x0 with Ax0≠0; then ‖A(tx0)‖=t‖Ax0‖→∞ as t→∞. Restricting to unit vectors measures the ratio of output size to input size, which is what “how much does A magnify?” means, and by homogeneity the ratio is already known on all of Fn once it is known on the sphere.

Theorem 16.3.2 (Properties of the Operator Norm)

Let ‖⋅‖ be a norm on each of Fn, Fm, Fℓ, and let A∈Mm×n(F). Then:

  • (a) ‖Ax‖≤‖A‖‖x‖ for every x∈Fn, and ‖A‖ is the smallest constant with this property;
  • (b) ‖⋅‖ is a norm on the vector space Mm×n(F);
  • (c) ‖In‖=1;
  • (d) ‖AB‖≤‖A‖‖B‖ for every B∈Mn×ℓ(F). In particular the operator norm induced on Mn(F) by a single norm on Fn 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 x≠0 is ‖x‖ 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 x=0 both sides are 0. For x≠0, the vector u=x/‖x‖ satisfies ‖u‖=1, so ‖Au‖≤‖A‖ by Definition 16.3.2; multiplying by ‖x‖ and using homogeneity gives ‖Ax‖≤‖A‖‖x‖. If K is any constant with ‖Ax‖≤K‖x‖ for all x, then taking a unit x0 at which the maximum is attained (Lemma 16.3.1) gives ‖A‖=‖Ax0‖≤K.

  • (b) Each value is a maximum of non-negative numbers, so ‖A‖≥0. If ‖A‖=0 then Au=0 for every unit u, hence Ax=0 for every x by (a), hence A=0. For a scalar c, ‖cA‖=max‖x‖=1|c|‖Ax‖=|c|‖A‖. For subadditivity, let x be a unit vector; then ‖(A+B)x‖≤‖Ax‖+‖Bx‖≤‖A‖+‖B‖, and taking the maximum over x gives ‖A+B‖≤‖A‖+‖B‖.

  • (c) ‖Inx‖=‖x‖=1 for every unit x, so the maximum is 1.

  • (d) Let x∈Fℓ be a unit vector. By (a) applied twice, first to A and then to B,

‖ABx‖≤‖A‖‖Bx‖≤‖A‖‖B‖‖x‖=‖A‖‖B‖.

Taking the maximum over unit x gives ‖AB‖≤‖A‖‖B‖. 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.

16.3.3. The three formulas

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.

Theorem 16.3.3 (Computing the 1-, 2- and Infinity-Norms)

Let A=(aij)∈Mm×n(F) with m,n≥1. Then

  • (a) ‖A‖1=max1≤j≤n ∑i=1m|aij|, the largest absolute column sum;
  • (b) ‖A‖∞=max1≤i≤m ∑j=1n|aij|, the largest absolute row sum;
  • (c) ‖A‖2=σ1(A), the largest singular value of A (Definition 13.8.1).

Idea. Each part is an upper bound and a witness. ① For ‖⋅‖1 the upper bound is subadditivity of the modulus with the column sums pulled out, and the witness is a standard basis vector ek 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 ‖⋅‖2 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 C=maxj∑i|aij| and R=maxi∑j|aij|.

(a). Let ‖x‖1=1. Then

‖Ax‖1=∑i|∑jaijxj|≤∑j|xj|∑i|aij|≤C∑j|xj|=C,

where the first inequality is subadditivity of the modulus together with an exchange of two finite sums. So ‖A‖1≤C. For the reverse, choose k with ∑i|aik|=C. The vector ek has ‖ek‖1=1 and Aek is the k-th column of A, so ‖Aek‖1=C. Hence ‖A‖1≥C.

(b). Let ‖x‖∞=1. For each i,

|∑jaijxj|≤∑j|aij||xj|≤∑j|aij|≤R,

so ‖Ax‖∞≤R and ‖A‖∞≤R. For the reverse, choose k with ∑j|akj|=R and define x by xj=akj―/|akj| when akj≠0 and xj=1 otherwise. Every xj has modulus 1, so ‖x‖∞=1; and akjxj=|akj| for every j, the choice being legal exactly because we checked akj≠0 before dividing. Hence the k-th entry of Ax is ∑j|akj|=R, giving ‖Ax‖∞≥R.

(c). Let A=UΣV∗ be a singular value decomposition (Theorem 13.8.1), with U,V unitary and Σ carrying σ1≥⋯≥σp≥0 on its diagonal, p=min(m,n). Let ‖x‖2=1 and put y=V∗x. A unitary matrix is an isometry for ‖⋅‖2 (Theorem 11.7.3, (a) and (f)), so ‖y‖2=1 and ‖Ax‖2=‖UΣy‖2=‖Σy‖2. The vector Σy has entries σiyi for i≤p and zeros below, so

‖Ax‖22=∑i=1pσi2|yi|2≤σ12∑i=1n|yi|2=σ12.

Hence ‖A‖2≤σ1. For the reverse, take x=v1, the first column of V; then y=V∗v1=e1, so ‖Av1‖2=‖Σe1‖2=σ1, and ‖v1‖2=1 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 σ1=max{‖Ax‖:‖x‖=1} and said that Chapter 16 would give that quantity a name, the operator norm ‖A‖2, “and this computation is why ‖A‖2=σ1”. 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 ‖A‖F 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 A=(3045)∈M2(R), compute ‖A‖1, ‖A‖∞, ‖A‖F and ‖A‖2.

Solution. The column sums of absolute values are 3+4=7 and 0+5=5, so ‖A‖1=7 by Theorem 16.3.3 (a). The row sums are 3 and 9, so ‖A‖∞=9 by (b). The Frobenius norm is 9+16+25=52.

For the spectral norm, form

A⊤A=(25202025).

This matrix is 25I+20(0110), whose eigenvalues are 25±20, namely 45 and 5. So the singular values are σ1=35 and σ2=5, and ‖A‖2=35≈6.708 by (c).

Two checks come for free. The product σ1σ2=15 equals |detA|=15, and σ12+σ22=50=‖A‖F2. The four numbers 7, 9, 52≈7.071 and 35≈6.708 are all different, which is the point: “the size of A” is not a well-posed phrase until the norm is named.

Example (The oblique projection, finally measured)

Chapter 13 §12 studied a projection P of F2 onto a line, with matrix (1x00) in a suitable orthonormal basis, and found that the largest value of ‖Pv‖ over unit v is 1+|x|2, which over R is 1/sin⁡θ for the angle θ between image and kernel. Confirm this with Theorem 16.3.3 (c) for x=−1, the matrix P′=(1−100) of Example (Two Oblique Projections, Measured).

Solution. (P′)⊤P′=(1−1−11), with eigenvalues 2 and 0; so σ1=2 and ‖P′‖2=2. This agrees with 1+|−1|2=2, and with 1/sin⁡θ for θ=π/4, the angle between the two lines. An orthogonal projection, by contrast, has ‖PU‖2=1 whenever U≠{0}, since it fixes every unit vector of U 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 A∈Mm×n(F) be non-zero, of rank r. Then

‖A‖2≤‖A‖F≤r‖A‖2,

and the left inequality is an equality if and only if r=1.

Idea. Write both norms in terms of the singular values: ‖A‖2=σ1 and ‖A‖F2=σ12+⋯+σr2. The claim is then the statement that one term of a sum of r non-negative numbers is at most the sum, which is at most r times the largest — and that the first is an equality exactly when there is only one term.

Proof. By Lemma 11.7.5 and A=UΣV∗ (Theorem 13.8.1), ‖A‖F=‖Σ‖F, so ‖A‖F2=∑i=1pσi2. Exactly r of the σi are non-zero (Theorem 13.8.1), and each satisfies σi≤σ1. Hence

σ12≤∑i=1pσi2≤rσ12,

which is the displayed chain after taking square roots and using ‖A‖2=σ1 (Theorem 16.3.3 (c)). The left inequality is an equality exactly when σi=0 for all i≥2, that is, when r=1.

So the two norms never differ by a factor worse than r≤min(m,n), and they agree exactly on the rank-one matrices. What they do not agree on is more interesting.

16.3.4. 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 ‖In‖=1 by Theorem 16.3.2 (c), whereas

‖In‖F=1+1+⋯+1⏟n=n.

So for n≥2 the Frobenius norm is not induced by any norm on Fn 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 ‖I‖=1 is forced: the identity stretches nothing, so a norm that measures stretching must give it 1. A norm that assigns the identity the value n is measuring something else — the amount of material in the matrix, which for In grows with n.

Warning. Two different objects are called ‖A‖∞ in the wild. In this book ‖A‖∞ is the operator norm induced by ‖⋅‖∞ on vectors, that is, the largest row sum. The entrywise maximum N(A)=maxi,j|aij| of Example (The entrywise maximum is not submultiplicative) is a different function, and it is not even submultiplicative. Whenever a source writes ‖A‖∞, check which one is meant by evaluating it at I2: the row-sum norm gives 1, the entrywise maximum gives 1 as well — so evaluate instead at J, where the row-sum norm gives 2 and the entrywise maximum gives 1.

Quick check. The Frobenius norm satisfies ‖Ax‖2≤‖A‖F‖x‖2 for every x, by Proposition 16.3.4 and Theorem 16.3.2 (a). Does that make it an operator norm?

Solution. No. Being compatible with a vector norm — satisfying ‖Ax‖≤‖A‖‖x‖ — is weaker than being induced by it. The induced norm is the smallest constant that works (Theorem 16.3.2 (a)), and for A=In with n≥2 the smallest constant is 1 while ‖In‖F=n. 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.

16.3.5. Exercises

A. Check your understanding

Exercise (A1)
  • (a) Define a matrix norm on Mn(F), 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 Mn(R) is the operator norm induced by some norm on Rn.
  • (e) Determine whether the following statement is correct, and justify your answer: for every A∈Mn(F) and every operator norm, ‖A2‖=‖A‖2.

Solution.

  • (a) A matrix norm is a norm on the vector space Mn(F) that also satisfies ‖AB‖≤‖A‖‖B‖ for all A,B (Definition 16.3.1). Given norms on Fn and Fm, the operator norm of A∈Mm×n(F) is ‖A‖=max{‖Ax‖:‖x‖=1} (Definition 16.3.2).

  • (b) ‖A‖1 is the largest absolute column sum, ‖A‖∞ the largest absolute row sum, and ‖A‖2=σ1(A).

  • (c) Two: the compactness of closed bounded sets, fact (A3) of the chapter introduction, used for the unit sphere S; and the extreme value theorem, fact (A4), used for x↦‖Ax‖. 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 ‖In‖F=n≠1 for n≥2, whereas every operator norm has ‖In‖=1 by Theorem 16.3.2 (c).

  • (e) Incorrect; only ≤ holds, by Theorem 16.3.2 (d). Take A=(0100), for which ‖A‖2=1 but A2=0, so ‖A2‖2=0.

B. Practice

Exercise (B1: Four numbers for one matrix)

Let A=(2−122)∈M2(R). Compute ‖A‖1, ‖A‖∞, ‖A‖F and ‖A‖2, and verify that ‖A‖2≤‖A‖F.

Solution. Absolute column sums: 2+2=4 and 1+2=3, so ‖A‖1=4. Absolute row sums: 2+1=3 and 2+2=4, so ‖A‖∞=4. Also ‖A‖F=4+1+4+4=13.

For the spectral norm,

A⊤A=(8225),

with trace 13 and determinant 40−4=36, so its eigenvalues solve λ2−13λ+36=0, giving λ=9 and λ=4. Hence σ1=3, σ2=2 and ‖A‖2=3 by Theorem 16.3.3 (c). Finally 3≤13, as Proposition 16.3.4 requires; the inequality is strict because rankA=2.

Exercise (B2: A rank-one matrix)

Let u∈Fm and v∈Fn be non-zero and put A=uv∗.

  • (a) Prove that ‖A‖2=‖u‖2‖v‖2.
  • (b) Verify the formula for u=(1,2) and v=(2,−1) in R2 by computing σ1(A) directly.

Solution.

  • (a) For x∈Fn we have Ax=(v∗x)u, so ‖Ax‖2=|v∗x|‖u‖2. By Theorem 11.1.5, |v∗x|=|⟨x,v⟩|≤‖x‖2‖v‖2, with equality when x=v/‖v‖2. Hence for unit x the quantity ‖Ax‖2 is at most ‖u‖2‖v‖2 and attains that value, so ‖A‖2=‖u‖2‖v‖2.

  • (b) Here A=(2−14−2) and A⊤A=(20−10−105), of trace 25 and determinant 100−100=0. So the eigenvalues are 25 and 0, giving σ1=5. And ‖u‖2‖v‖2=5⋅5=5, in agreement. Note also ‖A‖F=5=‖A‖2, as Proposition 16.3.4 predicts for a rank-one matrix.

Exercise (B3: The star trick)

Let A∈Mm×n(F). Prove that ‖A∗A‖2=‖A‖22. Hint: the eigenvalues of A∗A.

Solution. By Definition 13.8.1 the eigenvalues of A∗A are σ12≥⋯≥σn2≥0, where the σi are the singular values of A padded with zeros if n>m. The matrix A∗A is positive semidefinite, hence Hermitian, so by Corollary 12.4.3 it is WDW∗ with W unitary and D=diag(σ12,…,σn2). Then (A∗A)∗(A∗A)=WD2W∗, whose eigenvalues are σi4, so the singular values of A∗A are the σi2 and its largest is σ12. By Theorem 16.3.3 (c), applied to A∗A and then to A,

‖A∗A‖2=σ12=‖A‖22.

This proves the claim. Notice that it is an equality, whereas Theorem 16.3.2 (d) alone would give only ‖A∗A‖2≤‖A∗‖2‖A‖2.

C. Going deeper

Exercise (C1: The spectral norm between the other two)

Let A∈Mm×n(F).

Hint for (a): apply Theorem 11.1.5 to the splitting |aij|=|aij|1/2⋅|aij|1/2|xj| inside one row.

Solution.

  • (a) Write C=‖A‖1 and R=‖A‖∞ for the largest column and row sums (Theorem 16.3.3). Let ‖x‖2=1. For each i, Theorem 11.1.5 applied to the vectors with entries |aij|1/2 and |aij|1/2|xj| gives

    |∑jaijxj|2≤(∑j|aij|)(∑j|aij||xj|2)≤R∑j|aij||xj|2.

    Summing over i and exchanging the two finite sums,

    ‖Ax‖22≤R∑j|xj|2∑i|aij|≤RC∑j|xj|2=RC.

    Taking the maximum over unit x gives ‖A‖22≤RC, as claimed.

  • (b) There ‖A‖1=7 and ‖A‖∞=9, so the bound is 63≈7.937, while ‖A‖2=35≈6.708. The bound overestimates by about 18 per cent. It is sharp in general: for A=In all three numbers are 1.
Exercise (C2: Repairing the entrywise maximum)

For A∈Mn(F) let N(A)=maxi,j|aij|, as in Example (The entrywise maximum is not submultiplicative).

  • (a) Prove that N(AB)≤nN(A)N(B), and deduce that A↦nN(A) is a matrix norm.
  • (b) Deduce that nN(⋅) is not an operator norm for any n≥2.
  • (c) Explain why (a) and (b) together show that submultiplicativity and ‖I‖=1 are independent requirements.

Solution.

  • (a) For all i,j,

    |(AB)ij|=|∑kaikbkj|≤∑k|aik||bkj|≤nN(A)N(B),

    so N(AB)≤nN(A)N(B). Put ‖A‖=nN(A). It is a norm, being a positive multiple of the norm N, and

    ‖AB‖=nN(AB)≤n2N(A)N(B)=‖A‖‖B‖.

  • (b) ‖In‖=nN(In)=n≠1 for n≥2, while every operator norm gives the identity the value 1 (Theorem 16.3.2 (c)).

  • (c) N itself is a norm with N(In)=1 that is not submultiplicative; nN is submultiplicative with nN(In)=n. 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 ‖⋅‖F already showed, is still not sufficient on its own for n large.

Exercise (C3: When is a bound an equality?)

Let A∈Mm×n(F) be non-zero.

  • (a) Prove that ‖A‖2=‖A‖F if and only if rankA=1.
  • (b) Give an example with m=n=2 where ‖A‖F=2‖A‖2, the worst case allowed by Proposition 16.3.4.
  • (c) Is there a non-zero A with ‖A‖F<‖A‖2? Justify your answer.

Solution.

  • (a) This is the equality case of Proposition 16.3.4, proved there: ‖A‖F2=∑iσi2 and ‖A‖22=σ12, so the two agree exactly when σi=0 for every i≥2, which by Theorem 13.8.1 says rankA=1.

  • (b) Take A=I2. Its singular values are 1,1, so ‖A‖2=1 and ‖A‖F=2. Any matrix with two equal non-zero singular values does the same, for instance any 2×2 unitary matrix.

  • (c) No. ‖A‖2=σ1≤(∑iσi2)1/2=‖A‖F 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.

Found a typo or a mistake? Report it, or read this page’s source

The text of this book is licensed under CC BY-NC-SA 4.0.