Unitary matrix

A unitary/orthogonal matrix is a complex square matrix with orthonormal column vectors.
A matrix is unitary if and only if:

where is the conjugate transpose (adjoint) of , which equals the inverse matrix .

By definition, , any orthogonal matrix is also orthonormal.

Geometric meaning

Unitary transformations preserve geometry in vector spaces:
Inner products are preserved:
They act as pure rotations (or reflections if the determinant is : 1) of vectors, preserving lengths and angles between vectors.
So just undoes the rotation.

Properties

Let be a unitary matrix. Then:

for any vector : Unitary transformations preserve lengths.
for any vectors : Unitary transformations preserve angles – they just rotate.
→ The absolute value of the determinant , since the basis vectors are orthonormal (on the unit circle/cube/… hence unitary) and this attribute is preserved (only the sign might change). → That’s equivalent to saying the eigenvectors of also lie on the unit circle in the complex plane.
→ The product of two unitary matrices is unitary.

Rotation matrix in

is a unitary matrix that rotates vectors in by an angle counterclockwise.
Clockwise rotation is achieved by using , or swapping the signs of the off-diagonal elements.
Intuitively this works because the first column is where the basis vector ends up after rotation by , and the second column is where the basis vector ends up after rotation by .

EXAMPLE

is a matrix that rotates vectors in by counterclockwise.


Common unitary matrices

Applications

Unitary matrices are fundamental in:

signal processing (preserving signal energy)
PCA and SVD (rotations in high dimensions)
error-correcting codes (preserving hamming distance)

References

Visual Kernel
SVD - Data-Driven Science and engineering - Steve Brunton

Footnotes

  1. A composition between a pure rotation and a reflection is called improper rotation. So strictly speaking, orthogonal matrices always produce improper rotations. But for many cases this does not matter, e.g. in eigendecomposition, where contains the normalized eigenvectors of , then a is just as much a normalized eigenvector as , so can be used interchangeably, and you can even pick the signs to ensure that produces a pure rotation.