Chapter 8: The spectral theorems
8.4Isometries and unitary operators
A rotation of the plane moves almost every vector but changes no length and no angle. Section 6.3 met the matrices that behave like this, the unitary matrices, through a condition on their columns, . This section starts from the other end, from linear maps that preserve lengths, and shows that preserving lengths forces the rest: inner products are preserved, , and the columns of every matrix in orthonormal bases are orthonormal. For operators on a finite-dimensional space these maps are the unitary operators of Table 8.1.1, the analogues of the complex numbers on the unit circle, and over the spectral theorem describes them completely. Over the description is the classical one, by rotations in perpendicular planes and reflections. We derive it in the second half of the section, together with the distance-preserving maps of that are not linear.
Throughout, and are inner product spaces over , and inner products are linear in the first slot.
8.4.1Isometries
Definition 8.4.1 (Isometry). Let and be inner product spaces over . A linear map is an isometry if
Since is linear,
- (a)Every unitary matrix gives an isometry of (Exercise 6.3.37). In particular, the rotations and the reflections of Example 6.3.35 are isometries of .
(b)More generally, let have orthonormal columns. Then (Theorem 6.3.29), so for every ,
For instance, maps isometrically onto a plane in ; it sends to , and both have length .
- (c)If is a subspace of , the inclusion , , is an isometry.
- (d)Let be an orthonormal basis of . By Theorem 6.3.22(b), the coordinate map is an isometry from to with the dot product, and so is its inverse
.
Non-example 8.4.3 (Maps that are not isometries).
- The shear sends to , which has length .
- An orthogonal projection with is not an isometry. For , the vector is nonzero and orthogonal to , so by Theorem 6.3.3.
- The columns of are orthogonal but have length . The matrix is times an orthogonal matrix, so it multiplies every length by . It preserves right angles without being an isometry (compare Exercise 8.4.34).
Intuition (Why lengths control angles). Over , the law of cosines
computes the inner product of and from the three side lengths of the triangle with vertices , , . A linear map sends this triangle to the triangle with vertices , , , whose third side is
The main theorem lists the equivalent forms of the definition. Parts of it appeared as exercises in Section 6.6 (Exercise 6.1.51 and Exercise 6.6.40).
Theorem 8.4.4 (Characterizations of isometries). Let and be finite-dimensional inner product spaces with , and let . The following are equivalent.
- (I) is an isometry.
- (II) preserves inner products:
for all . - (III).
- (IV) maps every orthonormal list in to an orthonormal list in .
- (V) maps some orthonormal basis of to an orthonormal list in .
- (VI)For some orthonormal bases of and of , the columns of are orthonormal, that is, .
- (VII)For all orthonormal bases of and of , the columns of are orthonormal.
Proof (by polarization).
Uses. Proposition 6.1.32 (Polarization identities), Definition 6.6.8 (Adjoint), Proposition 6.1.37 (Inner products distinguish vectors), Theorem 6.4.20 (Orthonormal bases exist), Proposition 6.3.18 (Coordinates in an orthonormal basis), Corollary 6.3.6 (Pythagoras for orthogonal lists), Theorem 6.6.15 (Existence of adjoints), Theorem 2.4.2 (The matrix of a composite), Theorem 2.3.20 (Linear maps correspond to matrices), Theorem 6.3.29 (Orthonormal columns)
Techniques. Proving TFAE (Technique 2.7.1), Polarization (Technique 6.1.33), Move an operator across the inner product (Technique 6.6.20), Translate between maps and matrices (Technique 2.3.27)
We prove (I) (II) (III) (IV) (V) (I), and then (III) (VII) (VI) (III).
(I) (II). Let . For every scalar , linearity gives
If , Proposition 6.1.32(a) with gives
If , the same computation with the four values in Proposition 6.1.32(b) gives
(II) (III). Let . By the definition of the adjoint (Definition 6.6.8) and (II), for every
Conjugating both sides,
(III) (IV). Let be orthonormal in . Then
so is orthonormal.
(IV) (V). has an orthonormal basis (Theorem 6.4.20(a)), and by (IV) its image is orthonormal.
(V) (I). Let be an orthonormal basis of with orthonormal, and let . By Proposition 6.3.18, with , so
(III) (VII). Let be orthonormal bases of and , and . By Theorem 6.6.15(b), , and by Theorem 2.4.2,
By Theorem 6.3.29, says exactly that the columns of are orthonormal.
(VII) (VI). Orthonormal bases of and exist (Theorem 6.4.20(a)).
(VI) (III). The computation in (III) (VII), read backwards, gives
The polarization step is the only place where an identity for single vectors becomes one for pairs. There is a second route from lengths to that hides the pairs inside a result about self-adjoint operators.
Proof (through the quadratic form of ).
Uses. Proposition 6.6.18 (Properties of the adjoint), Definition 6.6.8 (Adjoint), Proposition 8.1.4 (Self-adjoint operators are determined by their diagonal values)
Techniques. Move an operator across the inner product (Technique 6.6.20)
We prove the cycle (I) (III) (II) (I); the remaining implications are as in the first proof.
Let . By Proposition 6.6.18,
If (I) holds, this vanishes for every , and by Proposition 8.1.4; that is (III). Conversely, if (III) holds, then for all ,
which is (II), and taking in (II) gives (I). □
The first proof recovers inner products from lengths; the second recovers the operator from its values . They rest on the same fact (the proof of Proposition 8.1.4 is a polarization argument), but the second shows the shape of the condition. An isometry is a map whose “defect” has zero quadratic form, and Section 8.3 studies maps whose is merely compared with .
Corollary 8.4.6 (Properties of isometries). Let be an isometry.
- (a) is injective, and it preserves inner products, distances, and orthogonality.
- (b)If is also an isometry, so is .
- (c)If and are finite-dimensional, then , and is the orthogonal projection .
- (d)If , then is bijective, , and is an isometry.
Uses. Remark 8.4.5 (Without finite dimension), Theorem 8.4.4 (Characterizations of isometries), Proposition 4.2.6 (Injective maps preserve dimension), Corollary 6.6.29 (Orthogonal projections are the self-adjoint idempotents), Proposition 6.6.18 (Properties of the adjoint), Theorem 4.2.11 (Rank–nullity)
Techniques. Test injectivity at (Technique 2.6.12)
(a) If , then , so ; hence is injective. Inner products are preserved by Remark 8.4.5, distances because
(b)
(c) Since is injective,
so is the orthogonal projection onto (Corollary 6.6.29). Clearly , and gives . So .
(d) By (a) and rank–nullity (Theorem 4.2.11),
For example, the matrix of Example 8.4.2(b) has
and is the projection onto the plane , whose unit normal is . Indeed with . So is a left inverse of but not a right inverse (compare Exercise 6.6.41).
The next corollary makes precise the remark in Section 6.3 that an -dimensional inner product space “looks like” with the dot product.
Uses. Theorem 6.4.20 (Orthonormal bases exist), Theorem 2.2.1 (Linear maps are determined by a basis), Theorem 8.4.4 (Characterizations of isometries), Corollary 8.4.6 (Properties of isometries), Example 8.4.2 (Isometries)
Techniques. Construct a map by choosing it on a basis (Technique 2.2.4)
So everything about lengths and angles in an -dimensional inner product space can be computed in with the dot product, once an orthonormal basis is fixed. What does depend on the space is how particular operators sit inside it.
Exercise 8.4.8. Let be finite-dimensional and .
- (a)Prove or disprove: if for every vector of some orthonormal basis of , then is an isometry.
- (b)Prove or disprove: if for every unit vector , then is an isometry.
Uses. Theorem 8.4.4 (Characterizations of isometries), Definition 8.4.1 (Isometry)
(a) False. Let have matrix . Then , a unit vector, but while . The images of the basis have length but are not orthogonal to each other, and (V) of Theorem 8.4.4 needs both.
(b) True. Let . Then is a unit vector, so by hypothesis and homogeneity of the norm,
For both sides are . So is an isometry.
8.4.2Unitary operators
An isometry is injective. In finite dimensions it is therefore invertible, and its inverse is again an isometry; in infinite dimensions it need not be (Example 8.4.12). The invertible ones get a name.
In Table 8.1.1 the unitary row was described by . The two descriptions agree in every dimension.
Proposition 8.4.10 (Unitary means ). Let be an inner product space and . Then is unitary if and only if has an adjoint and
In that case , and is unitary as well.
Uses. Remark 8.4.5 (Without finite dimension), Definition 6.6.8 (Adjoint), Proposition 6.6.9 (Adjoints are unique and linear), Definition 2.5.12 (Invertible)
Techniques. Find the adjoint by isolating the input (Technique 6.6.11)
Suppose is unitary. By Remark 8.4.5, preserves inner products, so for all ,
So is an adjoint of (Definition 6.6.8), and by Proposition 6.6.9, . Hence .
Conversely, suppose . Then is invertible with (Definition 2.5.12), and for every ,
So is unitary. Finally, is invertible, and
Theorem 8.4.11 (Characterizations of unitary operators). Let be a finite-dimensional inner product space and . The following are equivalent.
- (I) is an isometry.
- (II) is unitary.
- (III).
- (IV).
- (V) is an isometry.
- (VI) maps some orthonormal basis of to an orthonormal basis of .
- (VII) maps every orthonormal basis of to an orthonormal basis of .
- (VIII) is a unitary matrix for some orthonormal basis of , or equivalently for every one.
Uses. Theorem 8.4.4 (Characterizations of isometries), Corollary 8.4.6 (Properties of isometries), Proposition 8.4.10 (Unitary means ), Remark 4.2.20 (One-sided inverses of operators), Proposition 6.6.18 (Properties of the adjoint), Corollary 6.3.14 (Counting orthogonal vectors), Definition 6.3.33 (Unitary and orthogonal matrices)
Techniques. Proving TFAE (Technique 2.7.1)
We prove (I) (II) (IV) (III) (I), then (IV) (V), and finally (I) (VI) (VII) and (III) (VIII).
(I) (II). By Corollary 8.4.6(d) with , is bijective, with ; so is an invertible isometry.
(II) (IV). Proposition 8.4.10.
(IV) (III). is a right inverse of the operator on the finite-dimensional space , hence also a left inverse (Remark 4.2.20).
(III) (I). (III) (I) of Theorem 8.4.4.
(IV) (V). By Proposition 6.6.18, . So by (I) (III) of Theorem 8.4.4 applied to , is an isometry if and only if .
(I) (VI) (VII). Let . An orthonormal list of length in is a basis (Corollary 6.3.14). So (VI) says the same as (V) of Theorem 8.4.4, and (VII) says the same as its (IV) restricted to lists of length . Thus (I) (VII) by (I) (IV) there, (VII) (VI) is clear, and (VI) (I) is (V) (I) there.
(III) (VIII). With , (VI) and (VII) of Theorem 8.4.4 say that
So for operators on finite-dimensional spaces “isometry”, “unitary”, and “unitary matrix in an orthonormal basis” mean the same thing, and we use them interchangeably. For a real matrix , the operator on with the dot product is orthogonal exactly when is an orthogonal matrix.
Beyond 21-242 Infinite-dimensional spaces are not covered in 21-242.
Example 8.4.12 (The shift). On , define
padding the shorter polynomial with zero coefficients. This is the dot product of coefficient vectors, so it is an inner product, and is an orthonormal list. Let be multiplication by , which shifts the coefficients one step to the right. Then for every , so is an isometry, but it is not surjective: every has constant term , so . So is an isometry of an inner product space into itself that is not unitary, and (I) (II) of Theorem 8.4.11 fails without finite dimension. Exercise 8.4.40 computes its adjoint, which is a left inverse but not a right inverse.
Isometries have two more invariants that we can read off right away.
- (a)If is an isometry and is an eigenvalue of , then .
- (b)If is finite-dimensional and is unitary, then . In particular when .
Proof (by norms and by ).
Uses. Definition 7.2.6 (Eigenvalue, eigenvector), Definition 5.4.18 (Determinant of an operator), Theorem 8.4.11 (Characterizations of unitary operators), Theorem 5.4.14 (The transpose), Theorem 5.4.10 (Multiplicativity), Definition 5.3.6 (Determinant)
(a) Let with . Then
(b) Let be an orthonormal basis of and , which is a unitary matrix by Theorem 8.4.11; by Definition 5.4.18. The determinant is a polynomial in the entries with integer coefficients (Definition 5.3.6), and complex conjugation respects sums and products, so , where conjugates every entry. With Theorem 5.4.14,
So . Over , is a real number of absolute value . □
Proof (through the eigenvalues, over ).
Uses. Theorem 7.5.22 (Trace and determinant from eigenvalues), Theorem 8.4.11 (Characterizations of unitary operators), Definition 5.4.18 (Determinant of an operator)
For unitary operators, (a) is also Corollary 8.1.20(c), which was proved through normality. The norm argument here needs neither normality nor an adjoint, so it applies to every isometry, including the shift of Example 8.4.12 (which turns out to have no eigenvalues at all). Neither condition is sufficient for unitarity: the shear has determinant and the single eigenvalue , and it is not an isometry.
Invariant subspaces and the spectral form
For a normal operator, the orthogonal complement of an invariant subspace spanned by eigenvectors is invariant; this is what made the induction in the second proof of Theorem 8.2.5 work. For a unitary operator the same conclusion holds for every invariant subspace, with a proof that needs no eigenvectors. (By the spectral theorem it holds for normal operators on complex spaces too, but that is a consequence of the theorem, not a step towards it.)
Proof (directly).
Uses. Corollary 4.2.18 (Injective if and only if surjective), Corollary 8.4.6 (Properties of isometries), Theorem 8.4.11 (Characterizations of unitary operators), Definition 4.2.1 (Restriction)
Since is injective, so is , and is finite-dimensional, so is surjective (Corollary 4.2.18). That is, . Now let and . Write with ; since preserves inner products (Corollary 8.4.6(a)),
So , and is -invariant. The restrictions of to and to preserve norms, so they are isometries of finite-dimensional spaces into themselves, hence unitary by Theorem 8.4.11. □
Proof (through the adjoint).
Uses. Proposition 8.4.10 (Unitary means ), Exercise 6.6.38, Corollary 4.2.18 (Injective if and only if surjective), Theorem 8.4.11 (Characterizations of unitary operators)
Techniques. Move an operator across the inner product (Technique 6.6.20)
The first proof uses only that preserves inner products and maps onto itself. The second places the lemma in the pattern of Section 6.6: invariance passes from on to on , and for unitary the adjoint is the inverse, which has the same invariant subspaces. The lemma fails for isometries that are not onto: the shift of Example 8.4.12 maps (the polynomials with constant term ) into itself, but is not -invariant.
Over every operator has an eigenvector, so the lemma splits off one eigenline at a time. Alternatively, unitary operators are normal, and the complex spectral theorem applies. Either way:
Theorem 8.4.15 (Spectral form of unitary operators). Let be a finite-dimensional complex inner product space and . The following are equivalent.
- (I) is unitary.
- (II) is normal and every eigenvalue of has absolute value .
- (III) has an orthonormal basis consisting of eigenvectors of , and every eigenvalue of has absolute value .
In matrix form: is unitary if and only if for some unitary and some
Proof (from the complex spectral theorem).
Uses. Proposition 8.4.10 (Unitary means ), Corollary 8.4.13 (Eigenvalues and determinant), Theorem 8.2.5 (Complex spectral theorem), Theorem 8.4.11 (Characterizations of unitary operators), Theorem 7.4.2 (Diagonalizable means an eigenbasis), Exercise 6.3.37, Example 6.3.35 (Unitary and orthogonal matrices)
Techniques. Proving TFAE (Technique 2.7.1)
(I) (II). By Proposition 8.4.10, , so is normal, and its eigenvalues have absolute value by Corollary 8.4.13(a).
(II) (III). Theorem 8.2.5 gives an orthonormal eigenbasis, and the eigenvalues have absolute value by hypothesis.
(III) (I). Let be an orthonormal basis with and . Then
so maps an orthonormal basis to an orthonormal basis, and is unitary by (VI) (II) of Theorem 8.4.11.
Matrices. If is unitary, apply (I) (III) to on with the dot product. Let have the orthonormal eigenvectors as columns; it is unitary, and by Theorem 7.4.2,
Proof ((I) (III) by induction, without the spectral theorem).
Uses. Theorem 7.3.27 (Existence of eigenvalues over ), Corollary 8.4.13 (Eigenvalues and determinant), Lemma 8.4.14 (Invariant complements), Corollary 6.4.25 (Dimension of the orthogonal complement), Theorem 6.4.24 (Orthogonal decomposition)
Techniques. Induction on dimension via an invariant subspace (Technique 7.5.12)
The first proof inherits everything from normal operators. The second shows why the unitary case is easier to prove than the normal one: the invariance of the complement comes from Lemma 8.4.14, which needs neither eigenvectors nor adjoint identities. Over the second proof breaks at the first step, because a rotation of has no eigenvector. Repairing it with two-dimensional invariant subspaces gives the normal form of Section 8.4.4.
(a)The rotation
with has no real eigenvalue, but as a complex matrix With , which is unitary,
. The same works for every , which reflects the fact that rotations of the plane commute with each other. (b)Let be the cyclic shift, , with indices read modulo as in Definition 6.3.40. It is a permutation matrix, hence orthogonal, and . For the Fourier vector
, so . Since is unitary (Proposition 6.3.41), with and
. The eigenvalues of are all the th roots of unity. For , is the rotation of by about the line through (Example 8.4.28).
Uses. Theorem 8.4.15 (Spectral form of unitary operators), Theorem 2.2.1 (Linear maps are determined by a basis), Proposition 6.3.34 (Characterizations of unitary matrices)
An operator is determined by its values on the basis (Theorem 2.2.1), so there is exactly one such . Let ; its columns are orthonormal and . The columns of are eigenvectors for and , so
It has an orthonormal eigenbasis with eigenvalues and of absolute value , so it is unitary by (III) (I) of Theorem 8.4.15. Next,
the conjugate transpose of , as it must be. As a check, , since .
Unitary operators can move any orthonormal list to any other of the same length. This is the precise form of the statement that an inner product space “looks the same” in all directions.
Proposition 8.4.19 (Moving orthonormal lists). Let be finite-dimensional.
- (a)If and are orthonormal lists in , there is a unitary with for .
- (b)If is a subspace of and is an isometry, there is a unitary with .
Uses. Theorem 6.4.20 (Orthonormal bases exist), Theorem 2.2.1 (Linear maps are determined by a basis), Theorem 8.4.11 (Characterizations of unitary operators), Theorem 8.4.4 (Characterizations of isometries)
Techniques. Construct a map by choosing it on a basis (Technique 2.2.4), Extend a basis of a subspace (Technique 4.1.18)
(a) By Theorem 6.4.20(b), the lists extend to orthonormal bases and of , both of length . By Theorem 2.2.1 there is with for all . It maps an orthonormal basis to an orthonormal basis, so it is unitary by Theorem 8.4.11.
(b) Let be an orthonormal basis of . By (IV) of Theorem 8.4.4, is orthonormal in , so (a) gives a unitary with for all . The linear maps and agree on a basis of , so they are equal by the uniqueness in Theorem 2.2.1. □
In particular, any unit vector can be moved to any other by a unitary operator. When a problem is unchanged by unitary operators (it involves only lengths and inner products), we may therefore assume that a given unit vector is . Section 8.4.3 writes down an explicit operator that does this.
8.4.3Householder reflections
The simplest unitary operators other than scalars are reflections across hyperplanes.
Example 8.4.20 (Householder reflections). Let be a unit vector, and let
Then
which keeps the component of in and negates the component along . This is the operator of Exercise 8.1.30, so is Hermitian and unitary, and . Directly:
Extending to an orthonormal basis of eigenvectors, with eigenvalues , gives (Theorem 7.5.22).
For instance, for ,
the orthogonal matrix of Example 6.3.35(d). So that matrix is a reflection, across the plane .
The point of Householder reflections is that they can move a given vector to any other vector of the same length.
Proposition 8.4.21 (A reflection exchanging two vectors). Let with , , and . Let
Proof (by computation).
Uses. Example 8.4.20 (Householder reflections), Proposition 6.1.29 (Expanding a norm), Proposition 6.1.4 (Sesquilinearity)
Proof (by decomposing along the mirror).
Uses. Example 8.4.20 (Householder reflections), Proposition 6.1.4 (Sesquilinearity)
Techniques. Orthogonal decomposition (Technique 6.2.4)
First,
since and . So , and fixes it, while is negated (Example 8.4.20). Hence
and in the same way
The second proof is the picture in Figure 8.4.1: the mirror is the perpendicular bisector of the segment from to , which passes through because and are equally long. The condition cannot be dropped over : since is Hermitian, is always real (Proposition 6.6.30), so no Householder reflection maps to a with . Multiplying by a suitable scalar of absolute value first repairs this (Exercise).
Take , where (and if ). Then
Uses. Proposition 8.4.21 (A reflection exchanging two vectors), Example 8.4.20 (Householder reflections)
Let , so . Both targets have length and real dot products with , so Proposition 8.4.21 applies with , , and
For : and , so
Indeed
The matrix is symmetric by inspection. Its columns , , have squared lengths and pairwise dot products , , , all ; so is orthogonal.
For : and , so
and
(In floating-point computation one chooses the sign of the target that makes large, here the second, to avoid subtracting nearly equal numbers.)
8.4.4Real orthogonal operators Beyond 21-242
The normal form of real orthogonal operators, orientation, and rigid motions are not covered in lecture; see LADW Chapters 5 and 6 and LADR 7D.
Over an orthogonal operator need not have an eigenvector: a rotation of the plane has none. So the spectral form of Theorem 8.4.15 has no direct real analogue. Instead, we will show that every orthogonal operator is assembled from rotations of perpendicular planes and from on perpendicular lines. We first list the orthogonal operators of the plane.
Lemma 8.4.23 (Orthogonal matrices). Let be orthogonal. Then there is exactly one such that is one of
and , . The matrix is the reflection across the line spanned by
Uses. Definition 6.3.33 (Unitary and orthogonal matrices), Corollary 6.3.14 (Counting orthogonal vectors), Proposition 6.3.18 (Coordinates in an orthonormal basis), Definition 5.1.14 (The determinant)
Techniques. Take an inner product to isolate a coefficient (Technique 6.3.12)
The columns and of are orthonormal (Definition 6.3.33). Since , there is exactly one with
The first coefficient is , and , so . For we get and for we get ; the determinants are and . Since is determined by the first column and the sign by the determinant, the description is unique.
Write . By the angle-subtraction formulas,
□
So an orthogonal operator of the plane is a rotation if its determinant is and a reflection across a line if it is . Rotations by and are , the only orthogonal matrices that are scalars. In higher dimensions we need one more lemma, the block version of “a triangular unitary matrix is diagonal” (Exercise 6.3.57(c)).
Uses. Proposition 6.3.34 (Characterizations of unitary matrices), Definition 6.3.33 (Unitary and orthogonal matrices), Definition 2.4.33 (Block matrix)
Techniques. A sum of nonnegative terms (Technique 6.1.9), Induction on matrix size (Technique 5.5.17)
We induct on . If , then . Let , let be , and write
with and block upper triangular with diagonal blocks . The first columns of are orthonormal and have zero entries below , so adding their squared lengths gives , the sum over all entries of . By Proposition 6.3.34, the first rows of are orthonormal too, so the sum of over all entries of these rows is also . These rows consist of the entries of and of , so
A sum of nonnegative terms vanishes only if each term does, so . Now and
so and are unitary. By the induction hypothesis applied to , its blocks above the diagonal vanish and are unitary. □
Theorem 8.4.25 (Normal form of orthogonal operators). Let be a real inner product space with , and let be orthogonal. Then has an orthonormal basis in which is block diagonal with diagonal blocks
where , , and
Proof (from the real Schur form).
Uses. Theorem 7.5.27 (Real block upper-triangular form), Theorem 8.4.11 (Characterizations of unitary operators), Lemma 8.4.24 (Block-triangular unitary matrices), Lemma 8.4.23 (Orthogonal matrices), Theorem 4.5.17 (Change of basis)
Techniques. Prove it for triangular matrices first (Technique 7.5.23)
By Theorem 7.5.27(b), has an orthonormal basis such that is block upper triangular with diagonal blocks of size or , where no diagonal block has a real eigenvalue. By Theorem 8.4.11, is orthogonal, so by Lemma 8.4.24 it is block diagonal and each diagonal block is orthogonal.
A orthogonal matrix has , so it is or . A diagonal block is orthogonal without real eigenvalues. By Lemma 8.4.23, it is or for some ; the reflection has the eigenvalue , and , have real eigenvalues, so with
with . Finally, reordering (keeping the two vectors of each rotation block together) permutes the diagonal blocks, so we can list the rotation blocks first, in increasing order of angle, then the blocks , then the blocks .
For the matrix form, apply this to on with the dot product, and let have the vectors of as columns. Then is orthogonal and
Proof (by induction on dimension).
Uses. Theorem 7.3.31 (Small invariant subspaces), Lemma 8.4.14 (Invariant complements), Lemma 8.4.23 (Orthogonal matrices), Corollary 8.4.13 (Eigenvalues and determinant), Theorem 6.4.24 (Orthogonal decomposition), Remark 7.2.5 (Invariant subspaces and block matrices)
Techniques. Induction on dimension via an invariant subspace (Technique 7.5.12)
We induct on , proving the statement without the ordering of the blocks; reordering the basis as in the first proof then arranges them. By Theorem 7.3.31, has an invariant subspace with . By Lemma 8.4.14, is -invariant and , are orthogonal.
The piece . If , then is spanned by an eigenvector, whose eigenvalue is real of absolute value (Corollary 8.4.13(a)), so , a block or in a unit basis vector of . If , let be an orthonormal basis of ; the matrix of in it is orthogonal, so by Lemma 8.4.23 it is or . If it is , the unit eigenvectors of the lemma (taken as combinations of ) form an orthonormal basis of in which the matrix is . If it is with , it is . Otherwise it is a rotation block, and replacing by if necessary makes , as in the first proof.
The piece . If we are done. Otherwise , and the induction hypothesis gives an orthonormal basis of in which has the required block form.
By Theorem 6.4.24, , so the two bases together form an orthonormal basis of . Both summands are invariant, so in this basis is block diagonal, with the blocks of followed by those of (Remark 7.2.5). □
The first proof spends its effort in the real Schur form and needs only bookkeeping afterwards; the second needs only a small invariant subspace and does the splitting by hand, using that invariant complements come for free for orthogonal operators. A third proof complexifies and pairs conjugate eigenvectors, as in Proposition 7.3.34; Exercise carries it out.
The numbers in the normal form are determined by . The complex eigenvalues of , counted with multiplicity, are
8.4.5Orientation and rotations of Beyond 21-242
By Corollary 8.4.13(b), every orthogonal operator has determinant or . The sign records whether preserves orientation: whether it maps a right-handed orthonormal basis to a right-handed one, or reverses the handedness as a mirror does. In the plane, Lemma 8.4.23 says the two kinds are the rotations and the reflections, and the “F” in the picture below shows the difference.
The determinant is a continuous function of the entries that takes only the values on , so no path of orthogonal matrices leads from a rotation to a reflection. Conversely, by Theorem 8.4.25 every rotation is with an even number of blocks , which pair up into blocks ; letting all the angles shrink to moves continuously to inside . This is the precise sense in which the rotations are the orthogonal operators that can be “performed physically”.
In three dimensions a rotation has more structure than the word suggests: it always turns about a fixed axis.
Theorem 8.4.27 (Euler’s rotation theorem). Let be a real inner product space with , and let be a rotation. Then has an axis: a unit vector with . Moreover, for every orthonormal basis of starting with ,
for some , and .
Proof (by a determinant identity).
Uses. Theorem 8.4.11 (Characterizations of unitary operators), Theorem 5.4.10 (Multiplicativity), Theorem 5.4.14 (The transpose), Definition 5.1.1 (Determinant function), Theorem 7.2.9 (Eigenvalues and invertibility), Lemma 8.4.14 (Invariant complements), Lemma 8.4.23 (Orthogonal matrices), Proposition 5.3.17 (Block-diagonal matrices)
An axis. Let for some orthonormal basis ; it is orthogonal with . Using , Theorem 5.4.10, and Theorem 5.4.14,
Since
The matrix. Let be an orthonormal basis. The line is -invariant, so by Lemma 8.4.14 the plane
Proof (by counting eigenvalues).
Uses. Theorem 7.1.32 (Fundamental Theorem of Algebra), Proposition B.1.4 (Roots of real polynomials come in conjugate pairs), Corollary 7.1.43 (Odd degree gives a real root), Theorem 7.5.22 (Trace and determinant from eigenvalues), Corollary 8.4.13 (Eigenvalues and determinant), Corollary 7.3.13 (Eigenvalues of an operator)
Techniques. Complexify, then take real and imaginary parts (Technique 7.3.35)
We show that is an eigenvalue; the rest is as in the first proof. Let be the matrix of in an orthonormal basis, regarded as a complex unitary matrix, and let be its eigenvalues counted with multiplicity, the roots of (Theorem 7.1.32). Each satisfies (Corollary 8.4.13(a)), and (Theorem 7.5.22). Since is a real polynomial of odd degree, it has a real root (Corollary 7.1.43), say , and its nonreal roots come in conjugate pairs (Proposition B.1.4).
If , then , so and . If all three are real, each is , and their product is , so they are not all . In both cases is a root of , hence an eigenvalue of (Corollary 7.3.13). □
Proof (from the normal form).
Uses. Theorem 8.4.25 (Normal form of orthogonal operators), Proposition 5.3.17 (Block-diagonal matrices)
The first proof is the shortest and shows what matters: in any odd dimension ,
The trace formula identifies a rotation of quickly. Its axis is a kernel vector of , and the angle satisfies . The cosine determines only up to sign; replacing by replaces by , so the sense of the rotation depends on an orientation of the plane .
Example 8.4.28 (Axis and angle). Let
Its columns have squared lengths and pairwise dot products , , , all , so is orthogonal. Expanding along the first row,
Squaring, , the cyclic shift of Example 8.4.16(b) for : it turns about the same axis through , and its trace is .
Since is a group, the composite of two rotations of , about any two axes, is again a rotation about some axis, which is not obvious from the pictures (Exercise).
8.4.6Rigid motions Beyond 21-242
We now drop linearity and ask for maps that preserve distances.
Definition 8.4.29 (Rigid motion). Let be a real inner product space. A function , not assumed linear, is a rigid motion if
Translations and orthogonal operators are rigid motions, and so are their composites . These are affine maps (Definition 4.7.35) with orthogonal linear part. The theorem of this subsection is that there are no others. In Problem 6.8 the linearity of was proved by polarization; here we give a different proof that uses distances only, through the observation that in an inner product space a point is determined by its distances to the two ends of a segment that it lies on.
Lemma 8.4.30 (Points between two points). Let be an inner product space, let , and let . If
then .
Uses. Theorem 6.2.11 (Triangle inequality), Proposition 6.1.27 (Basic properties of the norm)
If , the first hypothesis gives , so
By Theorem 6.2.11(b), one of is a nonnegative real multiple of the other. If with , then , so
Theorem 8.4.31 (Distance-preserving maps are affine). Let and be real inner product spaces, and let satisfy
In particular, if is finite-dimensional, every rigid motion of has the form with orthogonal, and is bijective. The rigid motions of are the maps with and .
Proof (by betweenness).
Uses. Lemma 8.4.30 (Points between two points), Theorem 2.1.3 (Checking linearity), Definition 4.7.35 (Affine map), Theorem 8.4.11 (Characterizations of unitary operators), Proposition 2.3.17 (Every linear map is some )
Techniques. Choose a base point (Technique 4.7.4)
Step 1: segments go to segments. Let , , and . Then and
By Lemma 8.4.30 in ,
Step 2: lines go to lines. We show that
For , write
Step 3: is linear. Let
For , Step 1 with gives
So is linear (Theorem 2.1.3), and
Finite dimension. If is finite-dimensional, is orthogonal by (I) (II) of Theorem 8.4.11, so it is bijective, and so is , with inverse . For , is for a matrix (Proposition 2.3.17), which is orthogonal by Theorem 8.4.11. □
The proof used the inner product only through Lemma 8.4.30, that is, through the equality case of the triangle inequality: the “midpoint” of and , defined by distances alone, is unique. This is the idea behind the Mazur–Ulam theorem, which says that every surjective distance-preserving map between real normed spaces is affine. In a normed space such as with
In homogeneous coordinates (Section 4.7.6), the rigid motion of is the matrix
and by Theorem 4.7.42 composites and inverses correspond to products and inverses:
So the rigid motions of form a group, the Euclidean group. A rigid motion is proper if ; these are the motions of a solid body in space, and in the plane they are the translations and the rotations about a point (Exercise 4.7.62).
Exercises
Exercise 8.4.32. Decide which of the following linear maps are isometries, and which of the isometries are unitary.
- (a),
. - (b),
. - (c),
, where with . - (d),
, with as in (c). - (e), , with the Frobenius inner product.
Uses. Theorem 8.4.4 (Characterizations of isometries), Theorem 8.4.11 (Characterizations of unitary operators), Definition 8.4.1 (Isometry), Proposition 6.1.16 (Frobenius inner product)
(a) The matrix of is , with columns of squared length and and dot product . So is an isometry by Theorem 8.4.4 (with the standard bases). It maps into , so it is not an operator, and not unitary; its image is a plane.
(b) The matrix is . Its columns and have squared length , and their dot product is
(c) and , so
(d) For , , while has . So is not an isometry.
(e) is linear, and by Proposition 6.1.16,
Exercise 8.4.33. Prove or disprove each statement. Here and are finite-dimensional.
- (a)If is an isometry, then so is .
- (b)If is complex and every eigenvalue of has absolute value , then is unitary.
- (c)If is complex, is normal, and for every eigenvector of , then is unitary.
- (d)Statement (c) with “complex” replaced by “real”.
Uses. Corollary 8.4.6 (Properties of isometries), Theorem 8.4.15 (Spectral form of unitary operators), Theorem 8.2.5 (Complex spectral theorem), Example 8.1.9 (Normal operators)
(a) False. Let , , an isometry. Its adjoint is , since
(b) False. The shear on has the single eigenvalue but maps to , of length .
(c) True. By Theorem 8.2.5, has an orthonormal basis with . By hypothesis , and every eigenvalue of is one of the . So is unitary by (III) (I) of Theorem 8.4.15.
(d) False. Let have matrix . It is normal, since
Uses. Theorem 8.4.4 (Characterizations of isometries), Theorem 6.4.20 (Orthonormal bases exist)
Techniques. Plug in clever test vectors (Technique 6.1.40)
If , then . Otherwise let be an orthonormal basis of (Theorem 6.4.20). For , , so by hypothesis. Also
so by hypothesis and the orthogonality of and ,
So for one and all . If , then vanishes on a basis and . If , let . Then are pairwise orthogonal unit vectors, so is an isometry by (V) (I) of Theorem 8.4.4, and . Conversely, preserves orthogonality since
- (a)Show that an orthogonal matrix with integer entries is a signed permutation matrix: each row and each column has exactly one nonzero entry, which is or .
- (b)Count the orthogonal matrices with integer entries, and count those with determinant (for ).
Uses. Definition 6.3.33 (Unitary and orthogonal matrices), Theorem 6.3.29 (Orthonormal columns), Definition 5.2.29 (Permutation matrix), Theorem 5.4.10 (Multiplicativity)
Techniques. A sum of nonnegative terms (Technique 6.1.9)
(a) Let be the th column. It is a unit vector (Theorem 6.3.29), so with every . A nonzero integer has square at least , so exactly one entry of is nonzero, and it is ; say it is in row . If with , then , contradicting orthogonality. So is injective, hence a permutation of , and every row also contains exactly one nonzero entry.
(b) Conversely, every signed permutation matrix has orthonormal columns (they are for distinct ), so it is orthogonal with integer entries. Such a matrix is determined by the permutation ( choices) and the signs of its nonzero entries ( choices), so there are of them. Let . The map negates the first row, so it maps signed permutation matrices to signed permutation matrices, it is its own inverse, and (Theorem 5.4.10). So it exchanges those of determinant and those of determinant (the only possible values), and exactly half, , have determinant .
Exercise 8.4.36. Let
- (a)Show that , and that the sum equals if and only if every row of sums to .
- (b)Show that
, with equality if and only if every . Give an example of equality for .
Uses. Theorem 6.2.1 (Cauchy–Schwarz inequality), Theorem 8.4.11 (Characterizations of unitary operators), Definition 6.3.33 (Unitary and orthogonal matrices), Proposition 6.3.34 (Characterizations of unitary matrices)
Techniques. Apply Cauchy–Schwarz to a well-chosen inner product (Technique 6.2.21)
Let , so , and since is orthogonal.
(a) The th entry of is the th row sum, so
If every row sums to , then and the sum is . Conversely, if the sum is , then by Theorem 6.2.11(b) (the case of equality ), one of and is a nonnegative multiple of the other. They have the same nonzero length, so : every row sums to .
(b) Let be the th row of , a unit vector by Proposition 6.3.34, and let
and summing over the rows gives
Uses. Theorem 8.4.15 (Spectral form of unitary operators), Theorem 7.5.22 (Trace and determinant from eigenvalues), Theorem B.1.11 (Triangle inequality in )
By Theorem 8.4.15, with unitary and
If with , then . Conversely, if , equality holds in the triangle inequality, so by its equality case all nonzero are positive real multiples of one another. They all have absolute value , so they are all equal to one with . Then and
If is real orthogonal with , it is unitary as a complex matrix, and , so with , that is, .
Exercise 8.4.38. Let be a finite-dimensional complex inner product space and unitary.
- (a)Show that there is a self-adjoint , with all eigenvalues in , such that for (Definition 8.2.23).
- (b)Show that has a unitary square root: a unitary with .
- (c)Suppose and has distinct eigenvalues. Show that has exactly square roots , and that all of them are unitary.
Uses. Theorem 8.4.15 (Spectral form of unitary operators), Theorem 8.2.18 (Spectral decomposition), Definition 8.2.23 (Functions of an operator), Proposition 8.2.24 (Properties of ), Theorem 2.2.1 (Linear maps are determined by a basis), Definition 7.2.13 (Eigenspace)
Techniques. Apply a function to the eigenvalues (Technique 8.2.25)
(a) By Theorem 8.4.15 and Theorem 8.2.18,
(b) Let
(c) Let be an orthonormal eigenbasis of with and distinct . Suppose . Then , so
Exercise 8.4.39. A student argues: “Let be an orthogonal operator on a finite-dimensional real inner product space. Then , so is normal, and by the spectral theorem has an orthonormal basis of eigenvectors of .”
- (a)Find the error, and give a counterexample to the conclusion.
- (b)Show that an orthogonal has an orthonormal eigenbasis if and only if , and describe such operators geometrically.
Uses. Theorem 8.2.12 (Real spectral theorem), Theorem 8.2.5 (Complex spectral theorem), Proposition 8.4.10 (Unitary means ), Corollary 8.4.13 (Eigenvalues and determinant), Proposition 8.2.2 (Necessary conditions)
(a) The normal spectral theorem (Theorem 8.2.5) is a theorem about complex spaces. Over the spectral theorem (Theorem 8.2.12) requires to be self-adjoint, and an orthogonal operator need not be. The rotation is orthogonal, but its characteristic polynomial has no real root, so it has no eigenvectors at all.
(b) By Theorem 8.2.12 (with Proposition 8.2.2 for the forward direction), has an orthonormal eigenbasis if and only if . Since (Proposition 8.4.10), this says , that is, . In that case the eigenvalues are real of absolute value (Corollary 8.4.13(a)), so with the two eigenspaces orthogonal. So fixes the subspace and negates its orthogonal complement ; it is the reflection through , that is, (compare Exercise 8.1.30).
Exercise 8.4.40. Let be multiplication by on , with the inner product of Example 8.4.12, and define by
- (a)Show that is the adjoint of .
- (b)Show that and .
- (c)Show that has no eigenvalues, and that the only eigenvalue of is .
Uses. Example 8.4.12 (The shift), Definition 6.6.8 (Adjoint), Definition 7.2.6 (Eigenvalue, eigenvector), Definition 6.4.11 (Orthogonal projection)
(a) Let and (finite sums, padded with zeros). Then
So is the adjoint of (Definition 6.6.8).
(b)
(c) If , then , while is or has degree . So never holds with , and has no eigenvalue. For : the constant satisfies , so is an eigenvalue. Suppose with . If , then is or has degree , so cannot have degree ; hence , , and , which forces to be constant, contradicting . So is a nonzero constant, , and .
Exercise 8.4.41. Show that
is a rotation of , and find its axis and its angle.
Uses. Definition 8.4.26 (Rotation), Theorem 8.4.27 (Euler’s rotation theorem), Theorem 6.3.29 (Orthonormal columns)
The columns of are , , , each of squared length , with dot products , , and . So the columns of are orthonormal and is orthogonal. Expanding along the first row,
so and is a rotation (Definition 8.4.26).
For the axis,
Exercise 8.4.42. Let .
- (a)Show that if , then is an eigenvalue of .
- (b)Show that if and is odd, then is an eigenvalue of .
- (c)Give a with no real eigenvalue, and a with that does not have the eigenvalue , or show that none exists.
Uses. Theorem 5.4.10 (Multiplicativity), Theorem 5.4.14 (The transpose), Definition 5.1.1 (Determinant function), Theorem 7.2.9 (Eigenvalues and invertibility), Lemma 8.4.23 (Orthogonal matrices), Proposition 7.3.15 (Block-triangular matrices)
(a) Using , Theorem 5.4.10, and Theorem 5.4.14,
So
(b) Similarly,
(c) Let
Exercise 8.4.43. For , let
- (a)Show that . In words: reflecting across a line and then across a line through rotates the plane by twice the angle from to .
- (b)Deduce that every rotation is a product of two reflections, and that the first of them (the one applied first) can be any reflection.
Uses. Definition 2.4.1 (Matrix product)
(a) Let . Multiplying out,
So, since rotations compose by adding angles (, by the addition formulas),
The mirror of is the line at angle (Lemma 8.4.23, or directly), and that of is at angle . The angle from the second mirror to the first is , and the composite rotates by , twice that.
(b) Given and any reflection , let . By (a), . Since every reflection of the plane through is some , the reflection applied first can be chosen freely.
Exercise 8.4.44. Let and . A reflection of means a Householder matrix .
- (a)Show that if , then and .
- (b)Show that is a product of exactly reflections. (This is the Cartan–Dieudonné theorem for , in its sharp form.)
- (c)Deduce that every rotation of other than is a product of two reflections, and never of an odd number of reflections.
Uses. Example 8.4.20 (Householder reflections), Proposition 8.4.21 (A reflection exchanging two vectors), Corollary 6.4.25 (Dimension of the orthogonal complement), Theorem 4.2.11 (Rank–nullity), Theorem 5.4.10 (Multiplicativity), Theorem 8.4.27 (Euler’s rotation theorem)
Techniques. Induction on dimension via an invariant subspace (Technique 7.5.12), Parity via an invariant (Technique 5.2.19)
(a) Let
(b) By (a), no product of fewer than reflections equals , so it suffices to show that is a product of at most reflections. We use strong induction on . If , then , the empty product. Let , and let , of dimension . Pick , so . Since and all entries are real, Proposition 8.4.21 gives a reflection ,
because and preserves inner products. So , , and . Hence fixes every vector of , which has dimension since . So
(c) Let , . By Theorem 8.4.27, in an orthonormal basis starting with an axis has matrix with , that is, . Then is invertible (its determinant is
Exercise 8.4.45. On , let
- (a)Show that
for all , but that is not affine. - (b)Which step of the proof of Theorem 8.4.31 fails for ? Find all points with
.
Uses. Theorem 8.4.31 (Distance-preserving maps are affine), Lemma 8.4.30 (Points between two points), Theorem 4.7.37 (Affine maps preserve affine combinations), Corollary B.1.12 (Two terms, and the reverse triangle inequality)
(a) For ,
(b) Step 1 of the proof, which rests on Lemma 8.4.30, fails: in an inner product space the point is determined by its distances to and , but not for . Let . The conditions are
Both force and , so ; then , and both conditions hold exactly when . So the solutions are the points with , a whole segment; the midpoint is one of them, and is another. The max norm has no strict equality case in the triangle inequality:
Uses. Theorem 8.4.31 (Distance-preserving maps are affine), Proposition 4.7.20 (Affine independence through differences), Theorem 2.2.1 (Linear maps are determined by a basis)
Techniques. Choose a base point (Technique 4.7.4)
By Theorem 8.4.31, and with orthogonal. For each ,
By Proposition 4.7.20, the vectors
In , three points are affinely independent exactly when they are not collinear. Applying the first part with shows that a rigid motion fixing three non-collinear points is the identity. (Two points do not suffice: the reflection across the line through them fixes both.)
More practice
The practice bank has 15 exercises on this section’s topics.