Orthogonal Matrices
When a matrix has orthonormal columns, its transpose is its inverse — and it preserves every length and angle in space.
QTQ = I
A square matrix Q is orthogonal if its columns are orthonormal. The condition QTQ = I encodes two facts: unit lengths on the diagonal, zero dot products off the diagonal.
Q−1 = QT
Inverting a general matrix costs O(n³). For orthogonal Q, the inverse is just the transpose — an O(n²) read. Solving Qx = b is simply x = QTb.
Rigid Motions
Orthogonal matrices preserve lengths and angles. Multiplying by Q rotates or reflects space — but never stretches or squishes it.
- ‖Qx‖ = ‖x‖ — length preserved
- (Qx)·(Qy) = x·y — dot product preserved
- Angle between Qx and Qy = angle between x and y
- det(Q) = ±1 always
Rotation by θ
Householder Reflection
Reflect across the hyperplane perpendicular to unit vector u: H = I − 2uuᵀ.
H is symmetric, H² = I, det(H) = −1. Used in numerically stable QR factorization.
det = +1 or −1
det(Q)² = 1 because det(QᵀQ) = det(I) = 1. Two families emerge:
Orthogonal Transforms
- DFT matrix — unitary (complex orthogonal), F⁻¹ = F* free
- DCT (JPEG) — orthogonal, energy-preserving image transform
- MIMO precoding — Q rotates signal into channel eigenvectors
- Beamforming — QᵀQy = y exact matched filtering
- Householder QR — numerically stable, used in LAPACK
M5-L4: Least Squares
You've mastered orthogonal matrices — the rigid transformations that preserve shape. Next: what to do when a system has no exact solution at all — project onto the column space and solve the normal equations.