For AI agents: use /llms.txt for the Nakafa content index.
Home Orthogonal and Unitary Matrices Copy Content
Open
Orthogonal and unitary matrices are very special types of matrices. Imagine them as "clean" transformations that don't change distances and angles in space, only rotating or reflecting objects.
The difference is simple. Orthogonal matrices work with real numbers, while unitary matrices work with complex numbers. Both have the same properties, just different versions.
A square real matrix A ∈ R n × n A \in \mathbb{R}^{n \times n} A ∈ R n × n is called orthogonal if:
This means, to get the inverse of this matrix, we just transpose it. Very practical, right?
A square complex matrix A ∈ C n × n A \in \mathbb{C}^{n \times n} A ∈ C n × n is called unitary if:
Here A H A^H A H is the conjugate transpose of A A A . The concept is similar, just for complex numbers.
This is also equivalent to:
Real orthogonal matrices are actually a special case of unitary matrices, since R n × n ⊂ C n × n \mathbb{R}^{n \times n} \subset \mathbb{C}^{n \times n} R n × n ⊂ C n × n .
What's interesting about orthogonal and unitary matrices is that their determinants always have absolute value 1 1 1 . Why is this?
For a unitary matrix A A A , we have A H A = I A^H A = I A H A = I . If we calculate its determinant:
So ∣ det A ∣ = 1 |\det A| = 1 ∣ det A ∣ = 1 . For orthogonal matrices, the proof is the same, just using A T A = I A^T A = I A T A = I .
Eigenvalues of orthogonal and unitary matrices also have special properties. Every eigenvalue λ \lambda λ always satisfies:
Why is this? Suppose A ⋅ v = λ ⋅ v A \cdot v = \lambda \cdot v A ⋅ v = λ ⋅ v for an eigenvector v ≠ 0 v \neq 0 v = 0 . For the complex case, we can calculate:
Since v H v ≠ 0 v^H v \neq 0 v H v = 0 , then ∣ λ ∣ 2 = 1 |\lambda|^2 = 1 ∣ λ ∣ 2 = 1 , so ∣ λ ∣ = 1 |\lambda| = 1 ∣ λ ∣ = 1 .
For real orthogonal matrices, the eigenvalues can be 1 1 1 or − 1 -1 − 1 if real. But if complex, they can be written as:
This means complex eigenvalues lie on the unit circle in the complex plane.
Let's look at a familiar example, the rotation matrix:
We can check that this is an orthogonal matrix:
The characteristic polynomial is:
The result is λ 1 , 2 = e ± i α \lambda_{1,2} = e^{\pm i\alpha} λ 1 , 2 = e ± i α .
The transformation R 2 → R 2 : x ↦ A ⋅ x \mathbb{R}^2 \to \mathbb{R}^2 : x \mapsto A \cdot x R 2 → R 2 : x ↦ A ⋅ x represents a rotation by angle α \alpha α . For α ≠ 0 \alpha \neq 0 α = 0 and α ≠ π \alpha \neq \pi α = π , this matrix has no real eigenvalues, but has two complex eigenvalues.