• Nakafa

    Nakafa

    Learn free and with quality.
Subject
    • Bachelor
Exercises
Holy
  • Quran
Articles
  • Politics
  • Community
  • About

Command Palette

Search for a command to run...

Linear Methods of AI

Orthogonal and Unitary Matrices

Getting to Know Orthogonal and Unitary Matrices

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.

Mathematical Definitions

Orthogonal Matrices

A square real matrix A∈Rn×nA \in \mathbb{R}^{n \times n}A∈Rn×n is called orthogonal if:

A−1=ATA^{-1} = A^TA−1=AT

This means, to get the inverse of this matrix, we just transpose it. Very practical, right?

This is equivalent to:

ATA=AAT=IA^T A = A A^T = IATA=AAT=I

Unitary Matrices

A square complex matrix A∈Cn×nA \in \mathbb{C}^{n \times n}A∈Cn×n is called unitary if:

A−1=AHA^{-1} = A^HA−1=AH

Here AHA^HAH is the conjugate transpose of AAA. The concept is similar, just for complex numbers.

This is also equivalent to:

AHA=AAH=IA^H A = A A^H = IAHA=AAH=I

Real orthogonal matrices are actually a special case of unitary matrices, since Rn×n⊂Cn×n\mathbb{R}^{n \times n} \subset \mathbb{C}^{n \times n}Rn×n⊂Cn×n.

Interesting Determinant Properties

What's interesting about orthogonal and unitary matrices is that their determinants always have absolute value 1. Why is this?

For a unitary matrix AAA, we have AHA=IA^H A = IAHA=I. If we calculate its determinant:

1=det⁡I=det⁡(AHA)=det⁡AH⋅det⁡A=det⁡A‾⋅det⁡A=∣det⁡A∣21 = \det I = \det(A^H A) = \det A^H \cdot \det A = \overline{\det A} \cdot \det A = |\det A|^21=detI=det(AHA)=detAH⋅detA=detA⋅detA=∣detA∣2

So ∣det⁡A∣=1|\det A| = 1∣detA∣=1. For orthogonal matrices, the proof is the same, just using ATA=IA^T A = IATA=I.

Special Eigenvalues

Eigenvalues of orthogonal and unitary matrices also have special properties. Every eigenvalue λ\lambdaλ always satisfies:

∣λ∣=1|\lambda| = 1∣λ∣=1

Why is this? Suppose A⋅v=λ⋅vA \cdot v = \lambda \cdot vA⋅v=λ⋅v for an eigenvector v≠0v \neq 0v=0. For the complex case, we can calculate:

vHv=vHAHAv=(A⋅v)H(A⋅v)=(λ⋅v)H(λ⋅v)v^H v = v^H A^H A v = (A \cdot v)^H (A \cdot v) = (\lambda \cdot v)^H (\lambda \cdot v)vHv=vHAHAv=(A⋅v)H(A⋅v)=(λ⋅v)H(λ⋅v)
=λ‾⋅λ⋅vHv=∣λ∣2⋅vHv= \overline{\lambda} \cdot \lambda \cdot v^H v = |\lambda|^2 \cdot v^H v=λ⋅λ⋅vHv=∣λ∣2⋅vHv

Since vHv≠0v^H v \neq 0vHv=0, then ∣λ∣2=1|\lambda|^2 = 1∣λ∣2=1, so ∣λ∣=1|\lambda| = 1∣λ∣=1.

Forms of Eigenvalues

For real orthogonal matrices, the eigenvalues can be 1 or -1 if real. But if complex, they can be written as:

λ=exp⁡(iφ)=cos⁡φ+isin⁡φ\lambda = \exp(i\varphi) = \cos \varphi + i \sin \varphiλ=exp(iφ)=cosφ+isinφ

This means complex eigenvalues lie on the unit circle in the complex plane.

Concrete Example of Rotation Matrix

Let's look at a familiar example, the rotation matrix:

A=(cos⁡α−sin⁡αsin⁡αcos⁡α)A = \begin{pmatrix} \cos \alpha & -\sin \alpha \\ \sin \alpha & \cos \alpha \end{pmatrix}A=(cosαsinα​−sinαcosα​)

We can check that this is an orthogonal matrix:

ATA=(cos⁡αsin⁡α−sin⁡αcos⁡α)(cos⁡α−sin⁡αsin⁡αcos⁡α)A^T A = \begin{pmatrix} \cos \alpha & \sin \alpha \\ -\sin \alpha & \cos \alpha \end{pmatrix} \begin{pmatrix} \cos \alpha & -\sin \alpha \\ \sin \alpha & \cos \alpha \end{pmatrix}ATA=(cosα−sinα​sinαcosα​)(cosαsinα​−sinαcosα​)
=(cos⁡2α+sin⁡2αcos⁡α(−sin⁡α)+sin⁡αcos⁡α(−sin⁡α)cos⁡α+cos⁡αsin⁡α(−sin⁡α)(−sin⁡α)+cos⁡2α)=(1001)= \begin{pmatrix} \cos^2 \alpha + \sin^2 \alpha & \cos \alpha(-\sin \alpha) + \sin \alpha \cos \alpha \\ (-\sin \alpha) \cos \alpha + \cos \alpha \sin \alpha & (-\sin \alpha)(-\sin \alpha) + \cos^2 \alpha \end{pmatrix} = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}=(cos2α+sin2α(−sinα)cosα+cosαsinα​cosα(−sinα)+sinαcosα(−sinα)(−sinα)+cos2α​)=(10​01​)

Finding Eigenvalues

The characteristic polynomial is:

χA(t)=det⁡(cos⁡α−t−sin⁡αsin⁡αcos⁡α−t)\chi_A(t) = \det \begin{pmatrix} \cos \alpha - t & -\sin \alpha \\ \sin \alpha & \cos \alpha - t \end{pmatrix}χA​(t)=det(cosα−tsinα​−sinαcosα−t​)
=(cos⁡α−t)2+sin⁡2α=cos⁡2α+sin⁡2α−2tcos⁡α+t2= (\cos \alpha - t)^2 + \sin^2 \alpha = \cos^2 \alpha + \sin^2 \alpha - 2t \cos \alpha + t^2=(cosα−t)2+sin2α=cos2α+sin2α−2tcosα+t2
=1−2tcos⁡α+t2= 1 - 2t \cos \alpha + t^2=1−2tcosα+t2

The eigenvalues are:

λ1,2=cos⁡α±cos⁡2α−1=cos⁡α±−sin⁡2α=cos⁡α±isin⁡α\lambda_{1,2} = \cos \alpha \pm \sqrt{\cos^2 \alpha - 1} = \cos \alpha \pm \sqrt{-\sin^2 \alpha} = \cos \alpha \pm i \sin \alphaλ1,2​=cosα±cos2α−1​=cosα±−sin2α​=cosα±isinα

The result is λ1,2=e±iα\lambda_{1,2} = e^{\pm i\alpha}λ1,2​=e±iα.

The transformation R2→R2:x↦A⋅x\mathbb{R}^2 \to \mathbb{R}^2 : x \mapsto A \cdot xR2→R2: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.

Previous

Eigenvalues of Diagonal and Triangular Matrices

Next

Symmetric and Hermitian Matrices

  • Orthogonal and Unitary MatricesDiscover orthogonal and unitary matrices with determinant properties, eigenvalue analysis, rotation examples, and their applications in linear algebra.
On this page
  • Getting to Know Orthogonal and Unitary Matrices
  • Mathematical Definitions
    • Orthogonal Matrices
    • Unitary Matrices
  • Interesting Determinant Properties
  • Special Eigenvalues
  • Forms of Eigenvalues
  • Concrete Example of Rotation Matrix
    • Finding Eigenvalues
  • Comments
  • Report
  • Source code