Skip to content

Chapter 56. Hermitian and Normal Matrices

Hermitian and normal matrices are central objects in complex linear algebra. Hermitian matrices are the complex analogue of real symmetric matrices. Normal matrices form a broader class that includes Hermitian, unitary, and skew-Hermitian matrices.

Their importance comes from the spectral theorem. Hermitian matrices can be diagonalized by unitary matrices with real eigenvalues. More generally, every normal matrix can be unitarily diagonalized. This means that after a change to orthonormal coordinates, the matrix becomes diagonal. The geometry of the space is preserved because the change of basis is unitary. (en.wikipedia.org)

56.1 Hermitian Matrices

A complex square matrix ACn×nA\in\mathbb{C}^{n\times n} is Hermitian if

A=A, A^*=A,

where

A=AT A^*=\overline{A}^{\,T}

is the conjugate transpose.

Entrywise, this means

aij=aji a_{ij}=\overline{a_{ji}}

for all i,ji,j.

For example,

A=[21+i1i3] A= \begin{bmatrix} 2 & 1+i\\ 1-i & 3 \end{bmatrix}

is Hermitian because

1+i=1i. \overline{1+i}=1-i.

The diagonal entries of a Hermitian matrix are always real. Indeed,

aii=aii, a_{ii}=\overline{a_{ii}},

which implies

aiiR. a_{ii}\in\mathbb{R}.

Hermitian matrices are the complex analogue of symmetric matrices:

Real caseComplex case
AT=AA^T=AA=AA^*=A
SymmetricHermitian

56.2 Quadratic Forms

For a Hermitian matrix AA, the quantity

xAx x^*Ax

is always real.

Indeed,

xAx=(xAx)=xAx. \overline{x^*Ax} = (x^*Ax)^* = x^*A^*x.

Since A=AA=A^*,

xAx=xAx. x^*A^*x=x^*Ax.

Thus

xAx=xAx, \overline{x^*Ax}=x^*Ax,

so xAxRx^*Ax\in\mathbb{R}.

This property is fundamental. Hermitian matrices represent real-valued quadratic energy expressions in complex coordinates.

56.3 Positive Definite Hermitian Matrices

A Hermitian matrix AA is positive definite if

xAx>0for all nonzero xCn. x^*Ax>0 \quad \text{for all nonzero }x\in\mathbb{C}^n.

It is positive semidefinite if

xAx0 x^*Ax\ge 0

for all xx.

For example,

A=[2003] A= \begin{bmatrix} 2&0\\ 0&3 \end{bmatrix}

is positive definite because

xAx=2x12+3x22 x^*Ax = 2|x_1|^2+3|x_2|^2

is strictly positive for every nonzero vector.

Positive definite Hermitian matrices define Hermitian inner products:

x,yA=xAy. \langle x,y\rangle_A=x^*Ay.

They appear throughout optimization, statistics, differential equations, and numerical linear algebra.

56.4 Eigenvalues of Hermitian Matrices

Hermitian matrices have real eigenvalues.

Suppose

Av=λv Av=\lambda v

for some nonzero vector vv. Then

vAv=v(λv). v^*Av=v^*(\lambda v).

Using scalar multiplication,

vAv=λvv. v^*Av=\lambda v^*v.

Since AA is Hermitian, the quantity vAvv^*Av is real. Also,

vv>0. v^*v>0.

Therefore

λ=vAvvv \lambda = \frac{v^*Av}{v^*v}

is real.

Thus every eigenvalue of a Hermitian matrix lies in R\mathbb{R}. This is one of the most important structural properties of Hermitian matrices. (mathworld.wolfram.com)

56.5 Orthogonality of Eigenvectors

Eigenvectors corresponding to distinct eigenvalues of a Hermitian matrix are orthogonal.

Suppose

Av=λv,Aw=μw, Av=\lambda v, \qquad Aw=\mu w,

with

λμ. \lambda\ne\mu.

Compute

vAw. v^*Aw.

Since Aw=μwAw=\mu w,

vAw=μvw. v^*Aw=\mu v^*w.

On the other hand,

vAw=(Av)w. v^*Aw=(A^*v)^*w.

Because A=AA=A^*,

Av=Av=λv. A^*v=Av=\lambda v.

Thus

vAw=(λv)w=λvw. v^*Aw = (\lambda v)^*w = \overline{\lambda}v^*w.

Since λ\lambda is real,

λ=λ. \overline{\lambda}=\lambda.

Therefore

λvw=μvw. \lambda v^*w=\mu v^*w.

Hence

(λμ)vw=0. (\lambda-\mu)v^*w=0.

Since λμ\lambda\ne\mu,

vw=0. v^*w=0.

Thus vwv\perp w.

56.6 Spectral Theorem for Hermitian Matrices

Every Hermitian matrix can be unitarily diagonalized.

More precisely, if ACn×nA\in\mathbb{C}^{n\times n} is Hermitian, then there exists a unitary matrix UU and a real diagonal matrix Λ\Lambda such that

A=UΛU. A=U\Lambda U^*.

Equivalently,

UAU=Λ. U^*AU=\Lambda.

The columns of UU are orthonormal eigenvectors of AA, and the diagonal entries of Λ\Lambda are the eigenvalues of AA.

This is the complex spectral theorem for Hermitian matrices. It is one of the central theorems of linear algebra. Hermitian matrices admit an orthonormal eigenbasis and are therefore unitarily diagonalizable. (mathworld.wolfram.com)

56.7 Example of Hermitian Diagonalization

Consider

A=[2112]. A= \begin{bmatrix} 2&1\\ 1&2 \end{bmatrix}.

This matrix is Hermitian because it is real symmetric.

The characteristic polynomial is

det(AλI)=2λ112λ=(2λ)21. \det(A-\lambda I) = \begin{vmatrix} 2-\lambda & 1\\ 1 & 2-\lambda \end{vmatrix} = (2-\lambda)^2-1.

Thus

λ24λ+3=0. \lambda^2-4\lambda+3=0.

The eigenvalues are

λ1=3,λ2=1. \lambda_1=3, \qquad \lambda_2=1.

For λ1=3\lambda_1=3,

(A3I)v=0 (A-3I)v=0

gives eigenvectors proportional to

[11]. \begin{bmatrix} 1\\ 1 \end{bmatrix}.

For λ2=1\lambda_2=1,

(AI)v=0 (A-I)v=0

gives eigenvectors proportional to

[11]. \begin{bmatrix} 1\\ -1 \end{bmatrix}.

Normalize them:

u1=12[11],u2=12[11]. u_1= \frac{1}{\sqrt2} \begin{bmatrix} 1\\ 1 \end{bmatrix}, \qquad u_2= \frac{1}{\sqrt2} \begin{bmatrix} 1\\ -1 \end{bmatrix}.

Then

U=12[1111] U= \frac{1}{\sqrt2} \begin{bmatrix} 1&1\\ 1&-1 \end{bmatrix}

is unitary, and

A=U[3001]U. A = U \begin{bmatrix} 3&0\\ 0&1 \end{bmatrix} U^*.

56.8 Normal Matrices

A matrix ACn×nA\in\mathbb{C}^{n\times n} is normal if

AA=AA. A^*A=AA^*.

This condition means that AA commutes with its adjoint.

Every Hermitian matrix is normal because

A=A A^*=A

implies

AA=A2=AA. A^*A=A^2=AA^*.

Every unitary matrix is normal because

AA=I=AA. A^*A=I=AA^*.

Every skew-Hermitian matrix is also normal.

The class of normal matrices includes several important families:

Matrix typeDefining condition
HermitianA=AA^*=A
Skew-HermitianA=AA^*=-A
UnitaryAA=IA^*A=I
NormalAA=AAA^*A=AA^*

56.9 Skew-Hermitian Matrices

A matrix AA is skew-Hermitian if

A=A. A^*=-A.

Its diagonal entries are purely imaginary or zero, because

aii=aii. a_{ii}=-\overline{a_{ii}}.

If

Av=λv, Av=\lambda v,

then

vAv=λvv. v^*Av=\lambda v^*v.

But

(vAv)=vAv=vAv. (v^*Av)^* = v^*A^*v = -v^*Av.

Thus

vAv v^*Av

is purely imaginary, so λ\lambda is purely imaginary.

Skew-Hermitian matrices are the complex analogue of skew-symmetric real matrices.

56.10 Spectral Theorem for Normal Matrices

A matrix AA is normal if and only if it is unitarily diagonalizable.

That is,

AA=AA A^*A=AA^*

if and only if there exists a unitary matrix UU and a diagonal matrix Λ\Lambda such that

A=UΛU. A=U\Lambda U^*.

Unlike the Hermitian case, the diagonal entries of Λ\Lambda need not be real.

Thus a normal matrix may have complex eigenvalues, but it still has an orthonormal eigenbasis.

This theorem characterizes normal matrices completely in terms of orthonormal diagonalization. A complex matrix is normal exactly when it can be diagonalized by a unitary similarity transformation. (en.wikipedia.org)

56.11 Example of a Normal Matrix

Consider

A=[100i]. A= \begin{bmatrix} 1&0\\ 0&i \end{bmatrix}.

Then

A=[100i]. A^*= \begin{bmatrix} 1&0\\ 0&-i \end{bmatrix}.

Compute

AA=[1001]=AA. A^*A = \begin{bmatrix} 1&0\\ 0&1 \end{bmatrix} = AA^*.

Thus AA is normal.

It is not Hermitian because

AA. A^*\ne A.

Its eigenvalues are

1andi. 1 \quad \text{and} \quad i.

The standard basis vectors form an orthonormal eigenbasis, so the matrix is already unitarily diagonalized.

56.12 Unitary Diagonalization

Suppose

A=UΛU A=U\Lambda U^*

with UU unitary and Λ\Lambda diagonal.

Then

A=(UΛU)=UΛU. A^* = (U\Lambda U^*)^* = U\Lambda^*U^*.

Thus

AA=UΛΛU, A^*A = U\Lambda^*\Lambda U^*,

and

AA=UΛΛU. AA^* = U\Lambda\Lambda^*U^*.

Since diagonal matrices commute,

ΛΛ=ΛΛ. \Lambda^*\Lambda = \Lambda\Lambda^*.

Therefore

AA=AA. A^*A=AA^*.

Thus every unitarily diagonalizable matrix is normal.

This proves one direction of the spectral theorem for normal matrices.

56.13 Real Symmetric Matrices

A real symmetric matrix is a Hermitian matrix with real entries.

Thus the complex spectral theorem contains the real symmetric spectral theorem as a special case.

If

AT=A, A^T=A,

then there exists an orthogonal matrix QQ such that

A=QΛQT. A=Q\Lambda Q^T.

The eigenvalues are real, and the eigenvectors may be chosen orthonormal.

This is the real form of the spectral theorem.

56.14 Unitary Matrices and Eigenvalues

If UU is unitary and

Uv=λv Uv=\lambda v

for some nonzero vector vv, then

Uv2=v2. \|Uv\|_2=\|v\|_2.

But

Uv=λv, Uv=\lambda v,

so

λv2=λv2. \|\lambda v\|_2=|\lambda|\|v\|_2.

Therefore

λ=1. |\lambda|=1.

Thus every eigenvalue of a unitary matrix lies on the complex unit circle.

This geometric condition reflects norm preservation.

56.15 Normality and Orthogonality

For arbitrary matrices, eigenvectors corresponding to distinct eigenvalues need not be orthogonal.

For normal matrices, they are orthogonal.

Suppose

Av=λv,Aw=μw, Av=\lambda v, \qquad Aw=\mu w,

with

λμ. \lambda\ne\mu.

Then

Av=λv A^*v=\overline{\lambda}v

for eigenvectors of a normal matrix associated with a unitary diagonalization. Using this property and the commutation relation

AA=AA, A^*A=AA^*,

one proves

vw=0. v^*w=0.

Thus normality is precisely the condition that allows diagonalization by an orthonormal basis.

56.16 Functional Calculus

If

A=UΛU A=U\Lambda U^*

is unitarily diagonalized, then functions of AA are computed by applying the function to the diagonal entries:

f(A)=Uf(Λ)U, f(A)=Uf(\Lambda)U^*,

where

f(Λ)=diag(f(λ1),,f(λn)). f(\Lambda) = \operatorname{diag}(f(\lambda_1),\ldots,f(\lambda_n)).

For example,

Ak=UΛkU, A^k = U\Lambda^kU^*,

and

eA=UeΛU. e^A = Ue^\Lambda U^*.

This is important in differential equations, quantum mechanics, and matrix analysis.

56.17 Rayleigh Quotients

For a Hermitian matrix AA, the quantity

RA(x)=xAxxx R_A(x)=\frac{x^*Ax}{x^*x}

is called the Rayleigh quotient.

Since AA is Hermitian,

RA(x)R. R_A(x)\in\mathbb{R}.

If xx is an eigenvector with eigenvalue λ\lambda, then

RA(x)=λ. R_A(x)=\lambda.

The minimum and maximum eigenvalues satisfy

λminRA(x)λmax. \lambda_{\min} \le R_A(x) \le \lambda_{\max}.

Rayleigh quotients connect eigenvalues with optimization.

56.18 Positive Semidefinite Decomposition

A Hermitian matrix AA is positive semidefinite if and only if all its eigenvalues are nonnegative.

Indeed, if

A=UΛU, A=U\Lambda U^*,

then

xAx=xUΛUx. x^*Ax = x^*U\Lambda U^*x.

Let

y=Ux. y=U^*x.

Since UU is unitary,

xAx=yΛy=i=1nλiyi2. x^*Ax=y^*\Lambda y = \sum_{i=1}^n \lambda_i|y_i|^2.

This quantity is nonnegative for all xx exactly when

λi0 \lambda_i\ge 0

for all ii.

Thus positivity is determined entirely by the eigenvalues.

56.19 Numerical Importance

Hermitian and normal matrices are computationally important because their spectral structure is stable and geometrically well behaved.

Hermitian matrices have:

PropertyConsequence
Real eigenvaluesStable physical interpretation
Orthogonal eigenspacesStable decompositions
Unitary diagonalizationNumerically robust spectral methods
Positive definitenessStable optimization and factorization

Algorithms for Hermitian matrices are often faster and more stable than general matrix algorithms.

Examples include:

ProblemSpecialized method
Hermitian eigenvalue problemSymmetric QR, divide-and-conquer
Positive definite systemsCholesky factorization
Singular valuesHermitian reductions
OptimizationHessian analysis

Many scientific computing libraries exploit Hermitian structure explicitly.

56.20 Summary

A Hermitian matrix satisfies

A=A. A^*=A.

Its eigenvalues are real, and eigenvectors for distinct eigenvalues are orthogonal.

Every Hermitian matrix can be diagonalized by a unitary matrix:

A=UΛU. A=U\Lambda U^*.

A matrix is normal if

AA=AA. A^*A=AA^*.

Normal matrices are exactly the matrices that admit unitary diagonalization.

Hermitian, skew-Hermitian, and unitary matrices are all normal. Their structure is governed by orthonormal eigenbases and the geometry of Hermitian inner products.

The spectral theorem for Hermitian and normal matrices is one of the central organizing principles of linear algebra.