NETMaths
The bookUnit 1 · Inner Product Spaces and Forms27 / 83

Gram–Schmidt, orthogonal/unitary/normal matrices, spectral theorem

Why this is asked: Spectral theorem: real symmetric ⇒ orthogonally diagonalisable with real eigenvalues; complex normal ⇒ unitarily diagonalisable. Know which matrix classes are normal.

In one sentence

Normal is exactly the class the spectral theorem applies to, Hermitian and unitary and skew-Hermitian are the three special cases whose eigenvalues sit on the real axis, the unit circle and the imaginary axis, and every question here is asking which of those you are in.

Why the exam asks it

Because "diagonalisable" and "unitarily diagonalisable" are different conditions, and the second is the useful one — it gives an orthonormal eigenbasis, which is what makes the spectral theorem a tool rather than a classification. And because the word "symmetric" carries real weight over R\mathbb{R} and almost none over C\mathbb{C}.

The idea

An inner product gives v=v,v\|v\| = \sqrt{\langle v,v\rangle}, and

u,vuv(Cauchy–Schwarz),|\langle u,v\rangle| \le \|u\|\,\|v\| \qquad \text{(Cauchy–Schwarz)},

with equality exactly when uu and vv are dependent.

Gram–Schmidt turns any independent list into an orthonormal one spanning the same subspaces step by step. Two things worth noting: it needs the inputs to be independent (a dependent vector produces the zero vector at that step), and it gives the QRQR factorisation directly.

Orthogonal projection onto a subspace WW minimises distance: vPWvvw\|v - P_W v\| \le \|v - w\| for every wWw \in W, with equality only at w=PWvw = P_W v. That is the best-approximation property, and it is what least squares is.

V=WWV = W \oplus W^{\perp} for finite-dimensional WW, and (W)=W(W^{\perp})^{\perp} = W.

The adjoint and the three families

AA^* is the conjugate transpose. Then:

  • Hermitian (self-adjoint): A=AA^* = A.
  • Unitary: AA=IA^*A = I; over R\mathbb{R}, orthogonal, ATA=IA^TA = I.
  • Skew-Hermitian: A=AA^* = -A.
  • Normal: AA=AAAA^* = A^*A.

All three of the first families are normal, and normal is strictly larger than their union.

normal: AA* = A*A exactly the unitarily diagonalisable ones Hermitian A* = A unitary A*A = I skew-Hermitian A* = −A these three do not exhaust the normal ones where the eigenvalues sit real axis Hermitian imaginary axis skew-Hermitian unit circle unitary A general normal matrix has eigenvalues anywhere in ℂ.
Hermitian, unitary and skew-Hermitian all sit inside the normal matrices, and each family is pinned to a different part of the complex plane.

The spectral theorem

Over C\mathbb{C}: AA is normal     \iff AA is unitarily diagonalisable, that is A=UDUA = U D U^* with UU unitary.

Over R\mathbb{R}: AA is symmetric     \iff AA is orthogonally diagonalisable, A=QDQTA = Q D Q^T with QQ orthogonal and DD real.

Consequences that get examined directly:

  • Hermitian \Rightarrow all eigenvalues real.
  • Unitary \Rightarrow all eigenvalues have modulus 1.
  • Skew-Hermitian \Rightarrow all eigenvalues purely imaginary.
  • For a normal matrix, eigenvectors belonging to distinct eigenvalues are orthogonal — not merely independent, which is all you get in general.

Each converse holds too, given normality: a normal matrix with all-real eigenvalues is Hermitian, and so on. Without normality the converses fail — (1101)\begin{pmatrix}1&1\\0&1\end{pmatrix} has a real eigenvalue and is not Hermitian.

Real orthogonal matrices

A rotation (0110)\begin{pmatrix}0&-1\\1&0\end{pmatrix} is orthogonal, hence normal, hence unitarily diagonalisable over C\mathbb{C} with eigenvalues ±i\pm i — both of modulus 11, as required. Over R\mathbb{R} it is not diagonalisable at all, because x2+1x^2+1 has no real root.

So "normal" does not mean "diagonalisable over the reals". The real spectral theorem needs symmetric, and orthogonal is a different condition entirely. Note also det=±1\det = \pm 1 for orthogonal matrices, not +1+1: reflections have determinant 1-1.

Where intuition breaks

"Symmetric, so diagonalisable." Over R\mathbb{R} yes. Over C\mathbb{C}, (1ii1)\begin{pmatrix}1&i\\i&-1\end{pmatrix} is symmetric, nonzero and squares to zero. The complex analogue of real symmetry is Hermitian, and the option will say "symmetric" and let you supply "real" yourself.

"Normal, so diagonalisable over R\mathbb{R}." The rotation matrix. Normal gives unitary diagonalisability over C\mathbb{C}, which says nothing about R\mathbb{R}.

"Orthogonal, so det=1\det = 1." det=±1\det = \pm1. The +1+1 ones form the rotations; the 1-1 ones include every reflection.

"Hermitian, so positive definite." I-I is Hermitian with eigenvalues 1-1. Hermitian fixes the eigenvalues to be real, not positive.

"Eigenvectors for distinct eigenvalues are orthogonal." True for normal matrices; false in general, where they are only independent. Take (1102)\begin{pmatrix}1&1\\0&2\end{pmatrix}: eigenvectors (1,0)(1,0) and (1,1)(1,1) are not orthogonal.

"AA and AA^* have the same eigenvalues." They are conjugates of each other. For Hermitian matrices this is the same thing, which is why the error survives.

"Gram–Schmidt works on any spanning set." It needs independence. On a dependent list it produces a zero vector, which cannot be normalised — the fix is to discard it, but the naive algorithm divides by zero.

"Unitarily diagonalisable is the same as diagonalisable." (1102)\begin{pmatrix}1&1\\0&2\end{pmatrix} is diagonalisable — distinct eigenvalues — and not normal, so not unitarily diagonalisable. The eigenbasis exists; an orthonormal one does not.

The exam's angle

  1. Compute AAAA^* and AAA^*A first. Normality is a two-line check and it decides whether the spectral theorem is available at all.
  2. Match eigenvalue location to family. Real \to Hermitian, modulus one \to unitary, purely imaginary \to skew-Hermitian — but only once you know the matrix is normal.
  3. Read R\mathbb{R} or C\mathbb{C}. Symmetric matters over R\mathbb{R}; over C\mathbb{C} it is nearly worthless and Hermitian is the condition.
  4. For "orthonormal basis of eigenvectors", the question is normality, not diagonalisability.
  5. Remember orthogonal matrices are normal but rarely real-diagonalisable. Rotations are the standard counterexample and appear constantly.

The night before

  • Normal     \iff unitarily diagonalisable, over C\mathbb{C}. Real symmetric     \iff orthogonally diagonalisable, over R\mathbb{R}.
  • Hermitian: real eigenvalues. Unitary: modulus 11. Skew-Hermitian: purely imaginary. Converses hold given normality.
  • Normal eigenvectors for distinct eigenvalues are orthogonal; in general only independent.
  • Complex symmetric is worthless: (1ii1)\begin{pmatrix}1&i\\i&-1\end{pmatrix} is symmetric and nilpotent. Use Hermitian.
  • Orthogonal det=±1\Rightarrow \det = \pm1. Rotation: normal, not diagonalisable over R\mathbb{R}, eigenvalues ±i\pm i.
  • Cauchy–Schwarz with equality     \iff dependence.
  • Gram–Schmidt needs independence; gives QRQR.
  • Orthogonal projection is the best approximation; V=WWV = W\oplus W^\perp.
  • Diagonalisable \ne unitarily diagonalisable.

See it move

Gram–Schmidt, one projection at a timeinteractive

Steer v₂ and watch its shadow on e₁ get subtracted — and see the algorithm genuinely halt when the vectors go dependent.

Unlock interactive visuals

The trap here

“A real matrix with all real eigenvalues is orthogonally diagonalisable” — false

[[1,1],[0,1]]

Eigenvalue 1 twice but only one eigenvector; orthogonal diagonalisability requires symmetry.

More on this →

Check yourself — select all that apply

Let U be a real orthogonal n × n matrix (UᵀU = I). Which of the following are true?

Next: Quadratic forms, positive definiteness, Sylvester's law

Create a free account to keep your place and have this feed your study plan.

Open this in the full syllabus view · Unit 1