Page 324 - Matrix Analysis & Applied Linear Algebra
P. 324

320              Chapter 5                    Norms, Inner Products, and Orthogonality

                   5.6 UNITARY AND ORTHOGONAL MATRICES



                                    The purpose of this section is to examine square matrices whose columns (or
                                    rows) are orthonormal. The standard inner product and the euclidean 2-norm
                                    are the only ones used in this section, so distinguishing subscripts are omitted.




                                                  Unitary and Orthogonal Matrices

                                       •   A unitary matrix is defined to be a complex matrix U n×n whose
                                                                                            n
                                           columns (or rows) constitute an orthonormal basis for C .
                                       •   An orthogonal matrix is defined to be a real matrix P n×n whose
                                                                                            n
                                           columns (or rows) constitute an orthonormal basis for   .



                                        Unitary and orthogonal matrices have some nice features, one of which is
                                    the fact that they are easy to invert. To see why, notice that the columns of

                                    U n×n = u 1 | u 2 |· · · |u n  are an orthonormal set if and only if


                                                           1  when i = j,                   −1
                                                                                                  ∗
                                         [U U] = u u j =                 ⇐⇒ U U = I ⇐⇒ U      = U .
                                                                               ∗
                                                    ∗
                                           ∗
                                               ij   i      0  when i  = j,
                                                        ∗
                                    Notice that because U U = I ⇐⇒ UU = I, the columns of U are orthonor-
                                                                       ∗
                                    mal if and only if the rows of U are orthonormal, and this is why the definitions
                                    of unitary and orthogonal matrices can be stated either in terms of orthonormal
                                    columns or orthonormal rows.
                                        Another nice feature is that multiplication by a unitary matrix doesn’t
                                    change the length of a vector. Only the direction can be altered because
                                                       2                       2         n
                                                            ∗
                                                                       ∗
                                                               ∗
                                                  	Ux	 = x U Ux = x x = 	x	       ∀ x ∈C .         (5.6.1)
                                        Conversely, if (5.6.1) holds, then U must be unitary. To see this, set x = e i
                                    in (5.6.1) to observe u u i = 1 for each i, and then set x = e j + e k for j  = k
                                                        ∗
                                                        i
                                                                     ∗
                                                  ∗
                                    to obtain 0 = u u k + u u j = 2 Re (u u k ) . By setting x = e j +ie k in (5.6.1)
                                                         ∗
                                                  j      k           j
                                    it also follows that 0 = 2 Im (u u k ) , so u u k = 0 for each j  = k, and thus
                                                                          ∗
                                                                ∗
                                                                j         j
                                    (5.6.1) guarantees that U is unitary.
                                        In the case of orthogonal matrices, everything is real so that ( ) ∗  can be
                                                 T
                                    replaced by ( ) . Below is a summary of these observations.
   319   320   321   322   323   324   325   326   327   328   329