Tutorials › Linear Algebra › Unitary Diagonalization Criteria for Complex Matrices

Normal Operators and Complex Inner Product Spaces · Tutorial 871 of 1000

Unitary Diagonalization Criteria for Complex Matrices

Use normality to decide when a complex matrix has an orthonormal eigenbasis, and distinguish that condition from ordinary diagonalizability.

Advanced 10 min read

What You'll Learn

  • Define normal and unitarily diagonalizable complex matrices.
  • State the unitary diagonalization criterion for finite-dimensional complex matrices.
  • Prove both directions of the criterion.
  • Use the criterion to distinguish ordinary diagonalizability from unitary diagonalizability.
  • Check a proposed diagonalization by multiplying the matrices back.
  • Avoid confusing transpose with conjugate transpose when testing normality.

Normality and Unitary Diagonalization

Diagonalization expresses a matrix using its eigenvalues and eigenvectors. Unitary diagonalization is a more restrictive form: it requires the eigenvectors to be chosen as an orthonormal basis. That extra requirement matters in complex inner product spaces, where orthonormal coordinates preserve lengths and inner products. The exact test is normality, a condition stated directly in terms of a matrix and its adjoint.

Throughout, let \(\mathbf{A}\in\mathbb{C}^{n\times n}\), with \(n\geq 1\). The adjoint of a matrix is its conjugate transpose, denoted \(\mathbf{A}^{*}\). A square matrix \(\mathbf{U}\) is unitary if \(\mathbf{U}^{*}\mathbf{U}=\mathbf{U}\mathbf{U}^{*}=\mathbf{I}_n\). Its columns then form an orthonormal basis of \(\mathbb{C}^n\).

Definition: A matrix \(\mathbf{A}\in\mathbb{C}^{n\times n}\) is normal if \(\mathbf{A}^{*}\mathbf{A}=\mathbf{A}\mathbf{A}^{*}\). It is unitarily diagonalizable if there are a unitary matrix \(\mathbf{U}\) and a diagonal matrix \(\mathbf{D}\) such that \(\mathbf{A}=\mathbf{U}\mathbf{D}\mathbf{U}^{*}\). Equivalently, \(\mathbf{U}^{*}\mathbf{A}\mathbf{U}=\mathbf{D}\).

The columns of \(\mathbf{U}\) are the eigenvectors used in this representation. Indeed, from \(\mathbf{A}=\mathbf{U}\mathbf{D}\mathbf{U}^{*}\) and \(\mathbf{U}^{*}\mathbf{U}=\mathbf{I}_n\), we get \(\mathbf{A}\mathbf{U}=\mathbf{U}\mathbf{D}\). If the \(j\)th diagonal entry of \(\mathbf{D}\) is \(\lambda_j\), the \(j\)th column \(\mathbf{u}_j\) of \(\mathbf{U}\) satisfies \(\mathbf{A}\mathbf{u}_j=\lambda_j\mathbf{u}_j\). The additional condition is that these eigenvectors are orthonormal, not merely linearly independent.

Standard result (Complex Spectral Theorem): For a finite-dimensional complex inner product space, a linear operator is normal if and only if it has an orthonormal basis of eigenvectors. Equivalently, a complex square matrix is normal if and only if its columns can be represented by a unitary matrix of eigenvectors. The proof of this theorem was established in the earlier material on the Complex Spectral Theorem.

The Unitary Diagonalization Criterion

Theorem (Unitary diagonalization criterion): A matrix \(\mathbf{A}\in\mathbb{C}^{n\times n}\) is unitarily diagonalizable if and only if it is normal.

Proof. First suppose that \(\mathbf{A}\) is unitarily diagonalizable, so \(\mathbf{A}=\mathbf{U}\mathbf{D}\mathbf{U}^{*}\) for a unitary \(\mathbf{U}\) and a diagonal \(\mathbf{D}\). By the adjoint laws, \(\mathbf{A}^{*}=\mathbf{U}\mathbf{D}^{*}\mathbf{U}^{*}\). Therefore

$$ \begin{aligned} \mathbf{A}^{*}\mathbf{A} &=\mathbf{U}\mathbf{D}^{*}\mathbf{U}^{*}\mathbf{U}\mathbf{D}\mathbf{U}^{*} =\mathbf{U}\mathbf{D}^{*}\mathbf{D}\mathbf{U}^{*},\\ \mathbf{A}\mathbf{A}^{*} &=\mathbf{U}\mathbf{D}\mathbf{U}^{*}\mathbf{U}\mathbf{D}^{*}\mathbf{U}^{*} =\mathbf{U}\mathbf{D}\mathbf{D}^{*}\mathbf{U}^{*}. \end{aligned} $$

Every diagonal matrix commutes with its adjoint: if its diagonal entries are \(d_1,\dots,d_n\), then both \(\mathbf{D}^{*}\mathbf{D}\) and \(\mathbf{D}\mathbf{D}^{*}\) have diagonal entries \(|d_1|^2,\dots,|d_n|^2\). Hence \(\mathbf{A}^{*}\mathbf{A}=\mathbf{A}\mathbf{A}^{*}\), so \(\mathbf{A}\) is normal.

Conversely, suppose that \(\mathbf{A}\) is normal. By the Complex Spectral Theorem, \(\mathbb{C}^n\) has an orthonormal basis \(\mathbf{u}_1,\dots,\mathbf{u}_n\) consisting of eigenvectors of \(\mathbf{A}\). Let \(\mathbf{U}\) be the matrix with these vectors as columns, and let \(\mathbf{D}\) be diagonal with their corresponding eigenvalues on the diagonal. Orthonormality gives \(\mathbf{U}^{*}\mathbf{U}=\mathbf{I}_n\), and because the \(n\) vectors form a basis, \(\mathbf{U}\) is unitary. The eigenvector equations give \(\mathbf{A}\mathbf{U}=\mathbf{U}\mathbf{D}\). Multiplying on the right by \(\mathbf{U}^{*}\) yields \(\mathbf{A}=\mathbf{U}\mathbf{D}\mathbf{U}^{*}\). Thus \(\mathbf{A}\) is unitarily diagonalizable. Both directions hold, proving the equivalence. \(\square\)

The first direction is a direct computation. The second uses the Complex Spectral Theorem, which is why the field matters: the theorem guarantees an orthonormal eigenbasis for a normal matrix over \(\mathbb{C}\). The diagonal entries may repeat; the argument requires an orthonormal eigenbasis, not distinct eigenvalues.

Worked Applications

Worked Example: Constructing a Unitary Diagonalization

Verify that the matrix below is unitarily diagonalizable, and compute it from a unitary matrix and a diagonal matrix:

$$ \mathbf{U}=\frac{1}{\sqrt{2}} \begin{bmatrix}1&1\\1&-1\end{bmatrix}, \qquad \mathbf{D}= \begin{bmatrix}1&0\\0&i\end{bmatrix}. $$

The entries of \(\mathbf{U}\) are real, so \(\mathbf{U}^{*}=\mathbf{U}^{\top}=\mathbf{U}\). Direct multiplication gives

$$ \mathbf{U}^{*}\mathbf{U} =\frac12\begin{bmatrix}1&1\\1&-1\end{bmatrix} \begin{bmatrix}1&1\\1&-1\end{bmatrix} =\begin{bmatrix}1&0\\0&1\end{bmatrix}. $$

Thus \(\mathbf{U}\) is unitary. Multiplying the factors gives

$$ \begin{aligned} \mathbf{U}\mathbf{D}\mathbf{U}^{*} &=\frac12 \begin{bmatrix}1&i\\1&-i\end{bmatrix} \begin{bmatrix}1&1\\1&-1\end{bmatrix}\\ &=\frac12 \begin{bmatrix}1+i&1-i\\1-i&1+i\end{bmatrix} =\mathbf{A}. \end{aligned} $$

Since this is a unitary diagonalization, the criterion says \(\mathbf{A}\) is normal. As a check, \(\mathbf{D}^{*}=\begin{bmatrix}1&0\\0&-i\end{bmatrix}\), so \(\mathbf{D}^{*}\mathbf{D}=\mathbf{D}\mathbf{D}^{*}=\mathbf{I}_2\). The theorem’s computation then gives \(\mathbf{A}^{*}\mathbf{A}=\mathbf{A}\mathbf{A}^{*}=\mathbf{I}_2\). The columns of \(\mathbf{U}\) are eigenvectors: \(\mathbf{A}\mathbf{u}_1=\mathbf{u}_1\) and \(\mathbf{A}\mathbf{u}_2=i\mathbf{u}_2\).

Worked Example: Diagonalizable but Not Unitarily Diagonalizable

Consider the real matrix, regarded as a complex matrix,

$$ \mathbf{A}=\begin{bmatrix}1&1\\0&2\end{bmatrix}. $$

Its columns are not the issue here; find eigenvectors by solving the eigenvalue equations. The vectors \(\mathbf{v}_1=\begin{bmatrix}1\\0\end{bmatrix}\) and \(\mathbf{v}_2=\begin{bmatrix}1\\1\end{bmatrix}\) satisfy \(\mathbf{A}\mathbf{v}_1=\mathbf{v}_1\) and \(\mathbf{A}\mathbf{v}_2=2\mathbf{v}_2\). They are linearly independent, so they give an eigenbasis. With \(\mathbf{P}=\begin{bmatrix}1&1\\0&1\end{bmatrix}\) and \(\mathbf{D}=\begin{bmatrix}1&0\\0&2\end{bmatrix}\), this is ordinary diagonalization:

$$ \mathbf{P}^{-1}=\begin{bmatrix}1&-1\\0&1\end{bmatrix}, \qquad \mathbf{P}\mathbf{D}\mathbf{P}^{-1} =\begin{bmatrix}1&2\\0&2\end{bmatrix} \begin{bmatrix}1&-1\\0&1\end{bmatrix} =\begin{bmatrix}1&1\\0&2\end{bmatrix} =\mathbf{A}. $$

But \(\mathbf{P}\) is not unitary: its columns have dot product \(1\), not \(0\). More decisively, compute the normality test. Because \(\mathbf{A}\) is real, \(\mathbf{A}^{*}=\mathbf{A}^{\top}\), and

$$ \mathbf{A}^{*}\mathbf{A} =\begin{bmatrix}1&0\\1&2\end{bmatrix} \begin{bmatrix}1&1\\0&2\end{bmatrix} =\begin{bmatrix}1&1\\1&5\end{bmatrix}, \qquad \mathbf{A}\mathbf{A}^{*} =\begin{bmatrix}1&1\\0&2\end{bmatrix} \begin{bmatrix}1&0\\1&2\end{bmatrix} =\begin{bmatrix}2&2\\2&4\end{bmatrix}. $$

These products differ, so \(\mathbf{A}\) is not normal and, by the criterion, cannot be unitarily diagonalized. This example shows why ordinary diagonalizability is not enough: a basis of eigenvectors need not be orthonormal.

Worked Example: A Real Matrix Diagonalized over the Complex Numbers

Consider the rotation matrix \(\mathbf{R}=\begin{bmatrix}0&-1\\1&0\end{bmatrix}\). Test normality and find a unitary diagonalization over \(\mathbb{C}\).

The matrix is real, so \(\mathbf{R}^{*}=\mathbf{R}^{\top}=\begin{bmatrix}0&1\\-1&0\end{bmatrix}\). Multiplication in either order gives

$$ \mathbf{R}^{*}\mathbf{R} =\begin{bmatrix}1&0\\0&1\end{bmatrix}, \qquad \mathbf{R}\mathbf{R}^{*} =\begin{bmatrix}1&0\\0&1\end{bmatrix}. $$

It is normal. For \(\lambda=i\), the vector \(\begin{bmatrix}1\\-i\end{bmatrix}\) is an eigenvector because \(\mathbf{R}\begin{bmatrix}1\\-i\end{bmatrix}=\begin{bmatrix}i\\1\end{bmatrix}=i\begin{bmatrix}1\\-i\end{bmatrix}\). For \(\lambda=-i\), \(\begin{bmatrix}1\\i\end{bmatrix}\) is an eigenvector because \(\mathbf{R}\begin{bmatrix}1\\i\end{bmatrix}=\begin{bmatrix}-i\\1\end{bmatrix}=-i\begin{bmatrix}1\\i\end{bmatrix}\). Normalize both vectors and place them in the columns of \(\mathbf{U}\):

$$ \mathbf{U}=\frac{1}{\sqrt{2}} \begin{bmatrix}1&1\\-i&i\end{bmatrix}, \qquad \mathbf{D}=\begin{bmatrix}i&0\\0&-i\end{bmatrix}. $$

The columns have norm \(1\), and their inner product is \(\tfrac12(1+i^2)=0\), so \(\mathbf{U}^{*}\mathbf{U}=\mathbf{I}_2\). Multiplying back verifies the diagonalization:

$$ \mathbf{U}\mathbf{D}\mathbf{U}^{*} =\frac12 \begin{bmatrix}i&-i\\1&1\end{bmatrix} \begin{bmatrix}1&i\\1&-i\end{bmatrix} =\begin{bmatrix}0&-1\\1&0\end{bmatrix} =\mathbf{R}. $$

The eigenvalues are not real, but they are valid complex eigenvalues. Over \(\mathbb{R}\), this matrix has no real eigenvectors; over \(\mathbb{C}\), the displayed orthonormal eigenvectors give a unitary diagonalization.

Common Pitfall: Using the Transpose in Place of the Adjoint

For a complex matrix, the adjoint is the conjugate transpose, not merely the transpose. A tempting but wrong calculation is to take a complex symmetric matrix, observe \(\mathbf{A}^{\top}=\mathbf{A}\), and conclude that \(\mathbf{A}^{\top}\mathbf{A}=\mathbf{A}\mathbf{A}^{\top}\) makes it normal. Normality is defined using \(\mathbf{A}^{*}\), so that calculation does not test the required condition.

For a concrete correction, let \(\mathbf{B}=\begin{bmatrix}1&i\\i&0\end{bmatrix}\). It is symmetric, so the transpose-based comparison gives equal products automatically. But its adjoint is \(\mathbf{B}^{*}=\begin{bmatrix}1&-i\\-i&0\end{bmatrix}\), and the correct products are

$$ \mathbf{B}^{*}\mathbf{B} =\begin{bmatrix}2&i\\-i&1\end{bmatrix}, \qquad \mathbf{B}\mathbf{B}^{*} =\begin{bmatrix}2&-i\\i&1\end{bmatrix}. $$

They are unequal, so \(\mathbf{B}\) is not normal and cannot be unitarily diagonalized. The conjugation of each entry is essential to the criterion.

How to Use the Criterion

For a complex square matrix, the criterion separates two questions. Ordinary diagonalizability asks whether there is some eigenbasis, equivalently a representation \(\mathbf{A}=\mathbf{P}\mathbf{D}\mathbf{P}^{-1}\) with \(\mathbf{P}\) invertible. Unitary diagonalizability asks whether that eigenbasis can be chosen orthonormal. Normality answers the second question exactly. It does not say that every diagonalizable matrix is normal, as the second example demonstrates.

In practice, when the goal is to decide unitary diagonalizability, compute \(\mathbf{A}^{*}\mathbf{A}\) and \(\mathbf{A}\mathbf{A}^{*}\) and compare them. Equality proves normality and, by the theorem, guarantees a unitary diagonalization. Inequality rules one out. The test concerns the matrix as a whole; checking a few eigenvectors or finding a diagonalization with an arbitrary invertible matrix does not replace it.

PropertyRequired representation or testWhat it guarantees
Normal\(\mathbf{A}^{*}\mathbf{A}=\mathbf{A}\mathbf{A}^{*}\)Unitary diagonalizability over \(\mathbb{C}\)
Ordinarily diagonalizable\(\mathbf{A}=\mathbf{P}\mathbf{D}\mathbf{P}^{-1}\), with \(\mathbf{P}\) invertibleAn eigenbasis, not necessarily orthonormal
Unitarily diagonalizable\(\mathbf{A}=\mathbf{U}\mathbf{D}\mathbf{U}^{*}\), with \(\mathbf{U}\) unitaryNormality and an orthonormal eigenbasis

Check Your Understanding

Use the normality test and the distinction between ordinary and unitary diagonalization.

  1. State the unitary diagonalization criterion for a matrix in \(\mathbb{C}^{n\times n}\).
  2. Why does \(\mathbf{A}=\mathbf{U}\mathbf{D}\mathbf{U}^{*}\) imply \(\mathbf{A}^{*}\mathbf{A}=\mathbf{A}\mathbf{A}^{*}\)?
  3. Can the matrix \(\begin{bmatrix}3&1\\0&4\end{bmatrix}\) be unitarily diagonalized? Explain which test decides.
  4. What is the difference between the matrices \(\mathbf{P}\) and \(\mathbf{U}\) in ordinary and unitary diagonalizations?
  5. For a complex matrix, which operation must be used to form the adjoint, and why can the transpose give a wrong normality test?

Test yourself: Linear Algebra final exam practice

Ready to check what you know? Work a full-length practice final, then open the complete worked solutions.

Want all 1,114 practice problems?

This tutorial covers one topic. Calc Pass covers the full Calculus AB & BC curriculum (the core of Calculus I and II): full review, every practice question, and step-by-step solutions to each one.

QR code to download the Calc Pass app

Get the Calc Pass App

Scan to download and practice on the go, with instant scoring.

Calc Pass Volume 1: 1,114 Practice Questions book cover

Vol. 1: 1,114 Questions

The complete practice question bank, organized by AB & BC topic.

View on Amazon
Calc Pass Volume 2: Solutions Manual book cover

Vol. 2: Solutions Manual

Fully worked solutions to every question in Volume 1.

View on Amazon