Complex Matrices

Nearly everything we have studied in this book would remain true if the phrase real number were replaced by complex number wherever it occurs. Then we would deal with matrices with complex entries, systems of linear equations with complex coefficients (and complex solutions), determinants of complex matrices, and vector spaces with scalar multiplication by any complex number allowed. Moreover, the proofs of most theorems about (the real version of) these concepts extend easily to the complex case. It is not our intention here to give a full treatment of complex linear algebra. However, we will carry the theory far enough to give another proof of the Real Spectral Theorem (th:PrinAxes).

The set of complex numbers is denoted \(\mathbb {C}\) . We will use only the most basic properties of these numbers (mainly conjugation and absolute values), and the reader can find this material in Complex Numbers.

If \(n \ge 1\), we denote the set of all \(n\)-tuples of complex numbers by \(\mathbb {C}^n\). As with \(\RR ^n\), these \(n\)-tuples will be written either as row or column matrices and will be referred to as vectors. We define vector operations on \(\mathbb {C}^n\) as follows:

\begin{align*} [v_{1}, v_{2}, \ldots , v_{n}] + [w_{1}, w_{2}, \ldots , w_{n}] &= [v_{1} + w_{1}, v_{2} + w_{2}, \ldots , v_{n} + w_{n}] \\ u[v_{1}, v_{2}, \ldots , v_{n}] &= [uv_{1}, uv_{2}, \ldots , uv_{n}] \quad \mbox { for } u \mbox { in } \mathbb {C} \end{align*}

With these definitions, \(\mathbb {C}^n\) satisfies the axioms for a vector space (with complex scalars) given in Abstract Vector Spaces. Thus we can speak of spanning sets for \(\mathbb {C}^n\), of linearly independent subsets, and of bases. In all cases, the definitions are identical to the real case, except that the scalars are allowed to be complex numbers. In particular, the standard basis of \(\RR ^n\) remains a basis of \(\mathbb {C}^n\), called the standard basis of \(\mathbb {C}^n\).

A Generalization of the Dot Product for Complex Vectors

There is a generalization to \(\mathbb {C}^n\) of the dot product in \(\RR ^n\).

Clearly, if \(\vec {z}\) and \(\vec {w}\) actually lie in \(\RR ^n\), then \(\langle \vec {z}, \vec {w} \rangle = \vec {z} \dotp \vec {w}\) is the usual dot product.

Note that \(\langle \vec {z}, \vec {w} \rangle \) is a complex number in general, as opposed to requiring an inner product to be real as we do in Inner Product Spaces. However, if \(\vec {w} = \vec {z} = [z_{1}, z_{2}, \ldots , z_{n}]\), the definition gives \(\langle \vec {z}, \vec {z} \rangle = |z_{1}|^{2} + \ldots + |z_{n}|^{2}\) which is a nonnegative real number, equal to \(0\) if and only if \(\vec {z} = \vec {0}\). This explains the conjugation in the definition of \(\langle \vec {z}, \vec {w} \rangle \), and it gives (d) of the following theorem.

We leave (a) and (b) to the reader (Practice Problem ), and (d) has already been proved. To prove (c), write \(\vec {z} = [z_{1}, z_{2}, \ldots , z_{n}]\) and \(\vec {w} = [w_{1}, w_{2}, \ldots , w_{n}]\). Then

\begin{align*} \overline {\langle \vec {w}, \vec {z} \rangle } = (\overline {w_{1}\overline {z}_{1} + \ldots + w_{n}\overline {z}_{n}}) &= \overline {w}_{1}\overline {\overline {z}}_{1} + \ldots + \overline {w}_{n}\overline {\overline {z}}_{n} \\ &= z_{1}\overline {w}_{1} + \ldots + z_{n}\overline {w}_{n} = \langle \vec {z}, \vec {w} \rangle \end{align*}

The only properties of the norm function we will need are the following (the proofs are left to the reader):

A vector \(\vec {u}\) in \(\mathbb {C}^n\) is called a unit vector if \(\norm {\vec {u}} = 1\). Property (b) in Theorem 6 then shows that if \(\vec {z} \neq \vec {0}\) is any nonzero vector in \(\mathbb {C}^n\), then \(\vec {u} = \frac {1}{\norm {\vec {z}}}\vec {z}\) is a unit vector.

A matrix \(A = \left [ a_{ij} \right ]\) is called a complex matrix if every entry \(a_{ij}\) is a complex number. The notion of conjugation for complex numbers extends to matrices as follows: Define the conjugate of \(A = \left [ a_{ij} \right ]\) to be the matrix

\begin{equation*} \overline {A} = \left [ \begin{array}{c} \overline {a}_{ij} \end{array}\right ] \end{equation*}

obtained from \(A\) by conjugating every entry. Then (using Appendix ??)

\begin{equation*} \overline {A + B} = \overline {A} + \overline {B} \quad \mbox { and } \quad \overline {AB} = \overline {A} \; \overline {B} \end{equation*}

holds for all (complex) matrices of appropriate size.

Transposition of complex matrices is defined just as in the real case, and the following notion is fundamental.

Observe that \(A^{H} = A^{T}\) when \(A\) is real.

The following properties of \(A^{H}\) follow easily from the rules for transposition of real matrices and extend these rules to complex matrices. Note the conjugate in property (c).

Hermitian and Unitary Matrices

If \(A\) is a real symmetric matrix, it is clear that \(A^{H} = A\). The complex matrices that satisfy this condition turn out to be the most natural generalization of the real symmetric matrices:

Hermitian matrices are easy to recognize because the entries on the main diagonal must be real, and the “reflection” of each off-diagonal entry in the main diagonal must be the conjugate of that entry.

The following theorem extends Theorem ??, and gives a very useful characterization of Hermitian matrices in terms of the standard inner product in \(\mathbb {C}^n\).

If \(A\) is Hermitian, we have \(A^T = \overline {A}\). If \(\vec {z}\) and \(\vec {w}\) are columns in \(\mathbb {C}^n\), then \(\langle \vec {z}, \vec {w} \rangle = \vec {z}^T\overline {\vec {w}}\), so

\begin{equation*} \langle A\vec {z}, \vec {w} \rangle =(A\vec {z})^T\overline {\vec {w}} = \vec {z}^TA^T\overline {\vec {w}} = \vec {z}^T\overline {A}\overline {\vec {w}} = \vec {z}^T(\overline {A\vec {w}}) = \langle \vec {z}, A\vec {w} \rangle \end{equation*}
To prove the converse, let \(\vec {e}_{j}\) denote column \(j\) of the identity matrix. If \(A = \left [ a_{ij} \right ]\), the condition gives
\begin{equation*} \overline {a}_{ij} = \langle \vec {e}_{i}, A\vec {e}_{j} \rangle = \langle A\vec {e}_{i}, \vec {e}_{j} \rangle = {a}_{ij} \end{equation*}
Hence \(\overline {A} = A^T\), so \(A\) is Hermitian.

Let \(A\) be an \(n \times n\) complex matrix. As in the real case, a complex number \(\lambda \) is called an eigenvalue of \(A\) if \(A\vec {x} = \lambda \vec {x}\) holds for some column \(\vec {x} \neq \vec {0}\) in \(\mathbb {C}^n\). In this case \(\vec {x}\) is called an eigenvector of \(A\) corresponding to \(\lambda \).

If \(A\) is an \(n \times n\) matrix, the characteristic polynomial \(c_{A}(z)\) is a polynomial of degree \(n\) and the eigenvalues of \(A\) are just the roots of \(c_{A}(z)\). In most of our examples these roots have been real numbers (in fact, the examples have been carefully chosen so this will be the case!); but it need not happen, even when the characteristic polynomial has real coefficients. For example, if \(A = \left [ \begin{array}{rr} 0 & 1 \\ -1 & 0 \end{array}\right ]\) then \(c_{A}(z) = z^{2} + 1\) has roots \(i\) and \(-i\), where \(i\) is a complex number satisfying \(i^{2} = -1\). Therefore, we have to deal with the possibility that the eigenvalues of a (real) square matrix might be complex numbers.

For a complex matrix, \(c_{A}(z)\) has complex coefficients (possibly nonreal). However, an argument like that given in Exploration ?? of The Characteristic Equation still works to show that the eigenvalues of \(A\) are the roots (possibly complex) of \(c_{A}(z)\).

It is at this point that the advantage of working with complex numbers becomes apparent. The real numbers are incomplete in the sense that the characteristic polynomial of a real matrix may fail to have all its roots real. However, this difficulty does not occur for the complex numbers. The so-called fundamental theorem of algebra ensures that every polynomial of positive degree with complex coefficients has a complex root. Hence every square complex matrix \(A\) has a (complex) eigenvalue. Indeed (see ??), \(c_{A}(z)\) factors completely as follows:

\begin{equation*} c_{A}(z) = (z -\lambda _{1})(z -\lambda _{2}) \cdots (z -\lambda _{n}) \end{equation*}

where \(\lambda _{1}, \lambda _{2}, \ldots , \lambda _{n}\) are the eigenvalues of \(A\) (with possible repetitions due to multiple roots).

The next result extends Theorem ??, which asserts that eigenvectors of a symmetric real matrix corresponding to distinct eigenvalues are orthogonal. In the complex context, two \(n\)-tuples \(\vec {z}\) and \(\vec {w}\) in \(\mathbb {C}^n\) are said to be orthogonal if \(\langle \vec {z}, \vec {w} \rangle = 0\).

Let \(\lambda \) and \(\mu \) be eigenvalues of \(A\) with (nonzero) eigenvectors \(\vec {z}\) and \(\vec {w}\). Then \(A\vec {z} = \lambda \vec {z}\) and \(A\vec {w} = \mu \vec {w}\), so Theorem 16 gives

\begin{equation} \label {eigenvalEq} \lambda \langle \vec {z}, \vec {w} \rangle = \langle \lambda \vec {z}, \vec {w} \rangle = \langle A\vec {z}, \vec {w} \rangle = \langle \vec {z}, A\vec {w} \rangle = \langle \vec {z}, \mu \vec {w} \rangle = \overline {\mu } \langle \vec {z}, \vec {w} \rangle \end{equation}
If \(\mu = \lambda \) and \(\vec {w} = \vec {z}\), this becomes \(\lambda \langle \vec {z}, \vec {z} \rangle = \overline {\lambda } \langle \vec {z}, \vec {z} \rangle \). Because \(\langle \vec {z}, \vec {z} \rangle = \norm {\vec {z}}^{2} \neq 0\), this implies \(\lambda = \overline {\lambda }\). Thus \(\lambda \) is real, proving (1). Similarly, \(\mu \) is real, so equation (1) gives \(\lambda \langle \vec {z}, \vec {w} \rangle = \mu \langle \vec {z}, \vec {w} \rangle \). If \(\lambda \neq \mu \), this implies \(\langle \vec {z}, \vec {w} \rangle = 0\), proving (2).

Symmetric real matrices are Hermitian, and so the result follows immediately from Theorem 20.

The Real Spectral Theorem (??) asserts that every real symmetric matrix \(A\) is orthogonally diagonalizable—that is \(U^{T}AU\) is diagonal where \(U\) is an orthogonal matrix \((U^{-1} = U^{T})\). The next theorem identifies the complex analogs of these orthogonal real matrices.

If \(A = \left [ \begin{array}{cccc} | & | & & | \\ \vec {c}_{1} & \vec {c}_{2} & \cdots & \vec {c}_{n} \\ | & | & & | \end{array}\right ]\) is a complex matrix with \(j\)th column \(\vec {c}_{j}\), then \(A^T\overline {A} = \left [ \langle \vec {c}_{i}, \vec {c}_{j}\rangle \right ]\), as in Theorem ??. Now (a) \(\Leftrightarrow \) (b) follows, and (a) \(\Leftrightarrow \) (c) is proved in the same way.

Thus a real matrix is unitary if and only if it is orthogonal.

Given a real symmetric matrix \(A\), we saw in Orthogonal Matrices and Symmetric Matrices a procedure for finding an orthogonal matrix \(Q\) such that \(Q^{T}AQ\) is diagonal (see Example ??). The following example illustrates Theorem 20 and shows that the technique works for complex matrices.

Unitary Diagonalization

An \(n \times n\) complex matrix \(A\) is called unitarily diagonalizable if \(U^{H}AU\) is diagonal for some unitary matrix \(U\). As Example 29 suggests, we are going to prove that every Hermitian matrix is unitarily diagonalizable. However, with only a little extra effort, we can get a very important theorem that has this result as an easy consequence.

A complex matrix is called upper triangular if every entry below the main diagonal is zero. We owe the following theorem to Issai Schur.

We use induction on \(n\), mirroring the form of the proof of ??. If \(n = 1\), \(A\) is already upper triangular. If \(n > 1\), assume the theorem is valid for \((n - 1) \times (n - 1)\) complex matrices. Let \(\lambda _{1}\) be an eigenvalue of \(A\), and let \(\vec {y}_{1}\) be an eigenvector with \(\norm {\vec {y}_{1}} = 1\). Next, the (complex analog of the) Gram-Schmidt process provides \(\vec {y}_{2}, \ldots , \vec {y}_{n}\) such that \(\{\vec {y}_{1}, \vec {y}_{2}, \ldots , \vec {y}_{n}\}\) is an orthonormal basis of \(\mathbb {C}^n\). If \(U_{1} = \left [ \begin{array}{cccc} | & | & & | \\ \vec {y}_{1} & \vec {y}_{2} & \cdots & \vec {y}_{n} \\ | & | & & | \end{array}\right ]\) is the matrix with these vectors as its columns, then

\begin{equation*} U_{1}^HAU_{1} = \left [ \begin{array}{cc} \lambda _{1} & X_{1} \\ 0 & A_{1} \end{array}\right ] \end{equation*}
in block form. Now apply induction to find a unitary \((n - 1) \times (n - 1)\) matrix \(W_{1}\) such that \(W_{1}^HA_{1}W_{1} = T_{1}\) is upper triangular. Then \(U_{2} = \left [ \begin{array}{cc} 1 & 0 \\ 0 & W_{1} \end{array}\right ]\) is a unitary \(n \times n\) matrix. Hence \(U = U_{1}U_{2}\) is unitary (using Theorem 25), and
\begin{align*} U^HAU &= U_{2}^H(U_{1}^HAU_{1})U_{2} \\ &= \left [ \begin{array}{cc} 1 & 0 \\ 0 & W_{1}^H \end{array}\right ] \left [ \begin{array}{cc} \lambda _{1} & X_{1} \\ 0 & A_{1} \end{array}\right ] \left [ \begin{array}{cc} 1 & 0 \\ 0 & W_{1} \end{array}\right ] = \left [ \begin{array}{cc} \lambda _{1} & X_{1}W_{1} \\ 0 & T_{1} \end{array}\right ] \end{align*}

is upper triangular. Finally, \(A\) and \(U^{H}AU = T\) have the same eigenvalues by (the complex version of) Theorem ??, and they are the diagonal entries of \(T\) because \(T\) is upper triangular.

The fact that similar matrices have the same traces and determinants gives the following consequence of Schur’s theorem.

Schur’s theorem asserts that every complex matrix can be “unitarily triangularized.” However, we cannot substitute “unitarily diagonalized” here. In fact, if \(A = \left [ \begin{array}{cc} 1 & 1 \\ 0 & 1 \end{array}\right ]\), there is no invertible complex matrix \(U\) at all such that \(U^{-1}AU\) is diagonal. However, the situation is much better for Hermitian matrices.

By Schur’s theorem, let \(U^{H}AU = T\) be upper triangular where \(U\) is unitary. Since \(A\) is Hermitian, this gives

\begin{equation*} T^H = (U^HAU)^H = U^HA^HU^{HH} = U^HAU = T \end{equation*}
This means that \(T\) is both upper and lower triangular. Hence \(T\) is actually diagonal.

The Real Spectral Theorem asserts that a real matrix \(A\) is symmetric if and only if it is orthogonally diagonalizable (that is, \(P^{T}AP\) is diagonal for some real orthogonal matrix \(P\)). Theorem 35 is the complex analog of half of this result. However, the converse is false for complex matrices: There exist unitarily diagonalizable matrices that are not Hermitian.

There is a very simple way to characterize those complex matrices that are unitarily diagonalizable. To this end, an \(n \times n\) complex matrix \(N\) is called normal if \(NN^{H} = N^{H}N\). It is clear that every Hermitian or unitary matrix is normal, as is the matrix \(\left [ \begin{array}{rr} 0 & 1 \\ -1 & 0 \end{array}\right ]\) in Example 37. In fact we have the following result.

Assume first that \(U^{H}AU = D\), where \(U\) is unitary and \(D\) is diagonal. Then \(DD^{H} = D^{H}D\) as is easily verified. Because \(DD^{H} = U^{H}(AA^{H})U\) and \(D^{H}D = U^{H}(A^{H}A)U\), it follows by cancellation that \(AA^{H} = A^{H}A\).

Conversely, assume \(A\) is normal—that is, \(AA^{H} = A^{H}A\). By Schur’s theorem, let \(U^{H}AU = T\), where \(T\) is upper triangular and \(U\) is unitary. Then \(T\) is normal too:

\begin{equation*} TT^H = U^H(AA^H)U = U^H(A^HA)U = T^HT \end{equation*}
Hence it suffices to show that a normal \(n \times n\) upper triangular matrix \(T\) must be diagonal. We induct on \(n\); it is clear if \(n = 1\). If \(n > 1\) and \(T = \left [ t_{ij} \right ]\), then equating \((1, 1)\)-entries in \(TT^{H}\) and \(T^{H}T\) gives
\begin{equation*} |t_{11}|^2 + |t_{12}|^2 + \ldots + |t_{1n}|^2 = |t_{11}|^2 \end{equation*}
This implies \(t_{12} = t_{13} = \ldots = t_{1n} = 0\), so \(T = \left [ \begin{array}{cc} t_{11} & 0 \\ 0 & T_{1} \end{array}\right ]\) in block form. Hence \(T = \left [ \begin{array}{cc} \overline {t}_{11} & 0 \\ 0 & T_{1}^H \end{array}\right ]\) so \(TT^{H} = T^{H}T\) implies \(T_{1}T_{1}^H = T_{1}T_{1}^H\). Thus \(T_{1}\) is diagonal by induction, and the proof is complete.

We conclude this section by using Schur’s theorem (Theorem 32) to prove a famous theorem about matrices. Recall that the characteristic polynomial of a square matrix \(A\) is defined by \(c_{A}(z) = \mbox {det}(zI - A)\), and that the eigenvalues of \(A\) are just the roots of \(c_{A}(z)\).

If \(p(z)\) is any polynomial with complex coefficients, then \(p(P^{-1}AP) = P^{-1}p(A)P\) for any invertible complex matrix \(P\). Hence, by Schur’s theorem, we may assume that \(A\) is upper triangular. Then the eigenvalues \(\lambda _{1}, \lambda _{2}, \ldots , \lambda _{n}\) of \(A\) appear along the main diagonal, so

\begin{equation*} c_{A}(z) = (z - \lambda _{1})(z - \lambda _{2})(z - \lambda _{3}) \cdots (z -\lambda _{n}) \end{equation*}
Thus
\begin{equation*} c_{A}(A) = (A - \lambda _{1}I)(A - \lambda _{2}I)(A - \lambda _{3}I) \cdots (A - \lambda _{n}I) \end{equation*}
Note that each matrix \(A - \lambda _{i}I\) is upper triangular. Now observe:
(a)
\(A - \lambda _{1}I\) has zero first column because column 1 of \(A\) is \((\lambda _{1}, 0, 0, \ldots , 0)^{T}\).
(b)
Then \((A - \lambda _{1}I)(A - \lambda _{2}I)\) has the first two columns zero because the second column of \((A - \lambda _{2}I)\) is \((b, 0, 0, \ldots , 0)^{T}\) for some constant \(b\).
(c)
Next \((A - \lambda _{1}I)(A - \lambda _{2}I)(A - \lambda _{3}I)\) has the first three columns zero because column 3 of \((A -\lambda _{3}I)\) is \((c, d, 0, \ldots , 0)^{T}\) for some constants \(c\) and \(d\).

Continuing in this way we see that \((A - \lambda _{1}I)(A - \lambda _{2}I)(A - \lambda _{3}I) \cdots (A - \lambda _{n}I)\) has all \(n\) columns zero; that is, \(c_{A}(A) = 0\).

Practice Problems

In each case, compute the norm of the complex vector.
(a)
\((1, 1 - i, -2, i)\)
\[\answer {\sqrt {6}}\]
(b)
\((1 - i, 1 + i, 1, -1)\)
(c)
\((2 + i, 1 - i, 2, 0, -i)\)
\[\answer {\sqrt {13}}\]
(d)
\((-2, -i, 1 + i, 1 - i, 2i)\)
In each case, determine whether the two vectors are orthogonal.
(a)
\(\begin{bmatrix} 4\\ -3i\\ 2 + i\end{bmatrix}\), \(\begin{bmatrix} i\\ 2\\ 2 - 4i\end{bmatrix}\) YesNo
(b)
\(\begin{bmatrix} i\\ -i\\ 2 + i \end{bmatrix}\), \(\begin{bmatrix} i\\ i\\ 2 -i \end{bmatrix}\) YesNo
(c)
\(\begin{bmatrix} 1\\ 1\\ i\\ i\end{bmatrix}\), \(\begin{bmatrix} 1\\ i\\ -i\\ 1\end{bmatrix}\) YesNo
(d)
\(\begin{bmatrix} 4 + 4i\\ 2 + i\\ 2i\end{bmatrix}\), \(\begin{bmatrix} -1 + i\\ 2\\ 3 - 2i\end{bmatrix}\) YesNo
A subset \(U\) of \(\mathbb {C}^n\) is called a complex subspace of \(\mathbb {C}^n\) if it contains \(0\) and if, given \(\vec {v}\) and \(\vec {w}\) in \(U\), both \(\vec {v} + \vec {w}\) and \(z\vec {v}\) lie in \(U\) (\(z\) any complex number). In each case, determine whether \(U\) is a complex subspace of \(\mathbb {C}^3\).
(a)
\(U = \{(w, \overline {w}, 0) \mid w \mbox { in } \mathbb {C}\}\)
Not a subspace. For example, \(i(0, 0, 1) = (0, 0, i)\) is not in \(U\).
(b)
\(U = \{(w, 2w, a) \mid w \mbox { in } \mathbb {C}, a \mbox { in } \RR \}\)
(c)
\(U = \RR ^3\) SubspaceNot a Subspace
(d)
\(U = \{(v + w, v - 2w, v) \mid v, w \mbox { in } \mathbb {C}\}\)
In each case, find a basis over \(\mathbb {C}\), and determine the dimension of the complex subspace \(U\) of \(\mathbb {C}^3\) (see the previous exercise).
(a)
\(U = \{(w, v + w, v - iw) \mid v, w \mbox { in } \mathbb {C}\}\)
Basis \(\{(i, 0, 2), (1, 0, -1)\}\); dimension \(2\)
(b)
\(U = \{(iv + w, 0, 2v - w) \mid v, w \mbox { in } \mathbb {C}\}\)
(c)
\(U = \{(u, v, w) \mid iu - 3v + (1 - i)w = 0;\ \\ u, v, w \mbox { in } \mathbb {C}\}\)
Basis \(\{(1, 0, -2i), (0, 1, 1 - i)\}\); dimension \(2\)
(d)
\(U = \{(u, v, w) \mid 2u + (1 + i)v - iw = 0;\ \\ u, v, w \mbox { in } \mathbb {C}\}\)
In each case, determine whether the given matrix is Hermitian, unitary, or normal.
(a)
\(\left [ \begin{array}{rr} 1 & -i \\ i & i \end{array}\right ]\)
Normal only
(b)
\(\left [ \begin{array}{rr} 2 & 3 \\ -3 & 2 \end{array}\right ]\)
(c)
\(\left [ \begin{array}{rr} 1 & i \\ -i & 2 \end{array}\right ]\)
Hermitian (and normal), not unitary
(d)
\(\left [ \begin{array}{rr} 1 & -i \\ i & -1 \end{array}\right ]\)
(e)
\(\frac {1}{\sqrt {2}} \left [ \begin{array}{rr} 1 & -1 \\ 1 & 1 \end{array}\right ]\)
(f)
\(\left [ \begin{array}{cc} 1 & 1 + i \\ 1 + i & i \end{array}\right ]\)
None of these adjectives apply.
(g)
\(\left [ \begin{array}{cc} 1 + i & 1 \\ -i & -1 + i \end{array}\right ]\)
(h)
\(\frac {1}{\sqrt {2}|z|}\left [ \begin{array}{rr} z & z \\ \overline {z} & -\overline {z} \end{array}\right ]\), \(z \neq 0\)
Unitary (and normal); Hermitian if and only if \(z\) is real
Show that a matrix \(N\) is normal if and only if \(\overline {N}N^T = N^T\overline {N}\).
Let \(A = \left [ \begin{array}{cc} z & \overline {v} \\ v & w \end{array}\right ]\) where \(v\), \(w\), and \(z\) are complex numbers. Characterize in terms of \(v\), \(w\), and \(z\) when \(A\) is
(a)
Hermitian
(b)
unitary
(c)
normal.
In each case, find a unitary matrix \(U\) such that \(U^{H}AU\) is diagonal.
(a)
\(A = \left [ \begin{array}{rr} 1 & i\\ -i & 1 \end{array}\right ]\)
(b)
\(A = \left [ \begin{array}{cc} 4 & 3 - i \\ 3 + i & 1 \end{array}\right ]\)
\(U = \frac {1}{\sqrt {14}}\left [ \begin{array}{cc} -2 & 3 - i \\ 3 + i & 2 \end{array}\right ]\), \(U^HAU = \left [ \begin{array}{rr} -1 & 0 \\ 0 & 6 \end{array}\right ]\)
(c)
\(A = \left [ \begin{array}{rr} a & b\\ -b & a \end{array}\right ]\); \(a\), \(b\), real
(d)
\(A = \left [ \begin{array}{cc} 2 & 1 + i\\ 1 - i & 3 \end{array}\right ]\)
\(U = \frac {1}{\sqrt {3}}\left [ \begin{array}{cc} 1 + i & 1 \\ -1 & 1 - i \end{array}\right ]\), \(U^HAU = \left [ \begin{array}{rr} 1 & 0 \\ 0 & 4 \end{array}\right ]\)
(e)
\(A = \left [ \begin{array}{ccc} 1 & 0 & 1 + i\\ 0 & 2 & 0 \\ 1 - i & 0 & 0 \end{array}\right ]\)
(f)
\(A = \left [ \begin{array}{ccc} 1 & 0 & 0\\ 0 & 1 & 1 + i\\ 0 & 1 - i & 2 \end{array}\right ]\)
\(U = \frac {1}{\sqrt {3}}\left [ \begin{array}{ccc} \sqrt {3} & 0 & 0 \\ 0 & 1 + i & 1 \\ 0 & -1 & 1 - i \end{array}\right ]\), \(U^HAU = \left [ \begin{array}{rrr} 1 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 3 \end{array}\right ]\)
Show that \(\langle A \vec {x}, \vec {y} \rangle = \langle \vec {x}, A^{H}\vec {y}\rangle \) holds for all \(n \times n\) matrices \(A\) and for all \(n\)-tuples \(\vec {x}\) and \(\vec {y}\) in \(\mathbb {C}^n\).
(a)
Prove (a) and (b) of Theorem 3.
(b)
Prove Theorem 6.
\(\norm { \lambda Z }^2 = \langle \lambda Z, \lambda Z \rangle = \lambda \overline {\lambda } \langle Z, Z \rangle = |\lambda |^2 \norm { Z }^2\)
(c)
Prove Theorem 12.
a.
Show that \(A\) is Hermitian if and only if \(\overline {A} = A^T\).
b.
Show that the diagonal entries of any Hermitian matrix are real.
If the \((k, k)\)-entry of \(A\) is \(a_{kk}\), then the \((k, k)\)-entry of \(\overline {A}\) is \(\overline {a}_{kk}\) so the \((k, k)\)-entry of \((\overline {A})^T = A^{H}\) is \(\overline {a}_{kk}\). This equals \(a\), so \(a_{kk}\) is real.
(a)
Show that every complex matrix \(Z\) can be written uniquely in the form \(Z = A + iB\), where \(A\) and \(B\) are real matrices.
(b)
If \(Z = A + iB\) as in (a), show that \(Z\) is Hermitian if and only if \(A\) is symmetric, and \(B\) is skew-symmetric (that is, \(B^{T} = -B\)).
If \(Z\) is any complex \(n \times n\) matrix, show that \(ZZ^{H}\) and \(Z + Z^{H}\) are Hermitian.
A complex matrix \(B\) is called skew-Hermitian if \(B^{H} = -B\).
a.
Show that \(Z - Z^{H}\) is skew-Hermitian for any square complex matrix \(Z\).
b.
If \(B\) is skew-Hermitian, show that \(B^{2}\) and \(iB\) are Hermitian.
Show that \((B^2)^H = B^HB^H = (-B)(-B) = B^2\); \((iB)^H = \overline {i}B^H = (-i)(-B) = iB\).
c.
If \(B\) is skew-Hermitian, show that the eigenvalues of \(B\) are pure imaginary (\(i \lambda \) for real \(\lambda \)).
d.
Show that every \(n \times n\) complex matrix \(Z\) can be written uniquely as \(Z = A + B\), where \(A\) is Hermitian and \(B\) is skew-Hermitian.
If \(Z = A + B\), as given, first show that \(Z^{H} = A - B\), and hence that \(A = \frac {1}{2}(Z + Z^{H})\) and \(B = \frac {1}{2}(Z - Z^{H})\).
Let \(U\) be a unitary matrix. Show that:
(a)
\(\norm { U\vec {x} } = \norm { \vec {x} }\) for all columns \(\vec {x}\) in \(\mathbb {C}^n\).
(b)
\(|\lambda | = 1\) for every eigenvalue \(\lambda \) of \(U\).
(a)
If \(Z\) is an invertible complex matrix, show that \(Z^{H}\) is invertible and that \((Z^{H})^{-1} = (Z^{-1})^{H}\).
(b)
Show that the inverse of a unitary matrix is again unitary.
If \(U\) is unitary, \((U^{-1})^{-1} = (U^{H})^{-1} = (U^{-1})^{H}\), so \(U^{-1}\) is unitary.
(c)
If \(U\) is unitary, show that \(U^{H}\) is unitary.
Let \(Z\) be an \(m \times n\) matrix such that \(Z^{H}Z = I_{n}\) (for example, \(Z\) is a unit column in \(\mathbb {C}^n\)).
(a)
Show that \(V = ZZ^{H}\) is Hermitian and satisfies
\(V^{2} = V\).
(b)
Show that \(U = I - 2ZZ^{H}\) is both unitary and Hermitian (so \(U^{-1} = U^{H} = U\)).
(a)
If \(N\) is normal, show that \(zN\) is also normal for all complex numbers \(z\).
(b)
Show that (a) fails if normal is replaced by Hermitian.
\(H = \left [ \begin{array}{rr} 1 & i \\ -i & 0 \end{array}\right ]\) is Hermitian but \(iH = \left [ \begin{array}{rr} i & -1 \\ 1 & 0 \end{array}\right ]\) is not.
Show that a real \(2 \times 2\) normal matrix is either symmetric or has the form \(\left [ \begin{array}{rr} a & b \\ -b & a \end{array}\right ]\).
If \(A\) is Hermitian, show that all the coefficients of \(c_{A}(z)\) are real numbers.
(a)
If \(A = \left [ \begin{array}{rr} 1 & 1 \\ 0 & 1 \end{array}\right ]\), show that \(U^{-1}AU\) is not diagonal for any invertible complex matrix \(U\).
(b)
If \(A = \left [ \begin{array}{rr} 0 & 1 \\ -1 & 0 \end{array}\right ]\), show that \(U^{-1}AU\) is not upper triangular for any real invertible matrix \(U\).
Let \(U = \left [ \begin{array}{rr} a & b \\ c & d \end{array}\right ]\) be real and invertible, and assume that \(U^{-1}AU = \left [ \begin{array}{rr} \lambda & \mu \\ 0 & v \end{array}\right ]\). Then \(AU = U\left [ \begin{array}{rr} \lambda & \mu \\ 0 & v \end{array}\right ]\), and first column entries are \(c = a\lambda \) and \(-a = c\lambda \). Hence \(\lambda \) is real (\(c\) and \(a\) are both real and are not both \(0\)), and \((1 + \lambda ^{2})a = 0\). Thus \(a = 0\), \(c = a\lambda = 0\), a contradiction.
If \(A\) is any \(n \times n\) matrix, show that \(U^{H}AU\) is lower triangular for some unitary matrix \(U\).
If \(A\) is a \(3 \times 3\) matrix, show that \(A^{2} = 0\) if and only if there exists a unitary matrix \(U\) such that \(U^{H}AU\) has the form \(\left [ \begin{array}{rrr} 0 & 0 & u \\ 0 & 0 & v \\ 0 & 0 & 0 \end{array}\right ]\) or the form \(\left [ \begin{array}{rrr} 0 & u & v \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{array}\right ]\).
If \(A^{2} = A\), show that rank \(A = \mbox {tr}A\). [Hint: Use Schur’s theorem.]

Text Source

This section was adapted from Section 8.6 of Keith Nicholson’s Linear Algebra with Applications. (CC-BY-NC-SA)

W. Keith Nicholson, Linear Algebra with Applications, Lyryx 2018, Open Edition, pp. 445–456.