Analysis MOC

Matrix exponential

The matrix exponential exp :ℂ𝑛×𝑛 →ℂ𝑛×𝑛 uses the power series definition of the exponential function on matrices. Let 𝐀 be a 𝑛 ×𝑛 real/complex matrix. Then 𝑒𝐀 is given by #m/def/anal/vec

𝑒𝐀=βˆžβˆ‘π‘˜=0π€π‘˜π‘˜!

This is convergent for all 𝐀 βˆˆβ„‚π‘›Γ—π‘› under any norm.

Proof of convergence

Let β€– β‹…β€– denote the Operator norm. Then

0β‰€β€–π€π‘˜π‘˜!β€–=β€–π€π‘˜β€–π‘˜!β‰€β€–π€β€–π‘˜π‘˜!

Since 𝑒‖𝐴‖ =βˆ‘βˆžπ‘˜=0β€–π΄β€–π‘˜π‘˜! converges, the sequence converges absolutely and uniformly by the Weierstraß M-test. ℂ𝑛×𝑛 is finite-dimensional: All finite dimensional normed vector spaces are Banach and All norms on a finite dimensional space are equivalent. Thus the series converges, and does so regardless of norm.

Properties

For any 𝐀 βˆˆβ„‚π‘›Γ—π‘›, the following properties hold: #m/thm/anal/vec

  1. For any invertible 𝐓 ∈GL(𝑛), π‘’π“π€π“βˆ’1 =π“π‘’π€π“βˆ’1.
  2. 𝑒𝑑𝐀 uniquely solves ˙𝐗(𝑑) =𝐀𝐗(𝑑) with initial condition 𝐗(0) =𝐈.
  3. 𝑒𝑑𝐀𝑒𝑠𝐁 =𝑒𝑑𝐀+𝑠𝐁 if 𝐀𝐁 =𝐁𝐀 for all 𝑑,𝑠 βˆˆβ„‚.
  4. (𝑒𝐀)† =𝑒(𝐀†)
  5. det𝑒𝐀 =𝑒tr⁑𝐀
  6. π‘’βˆ’π‘–πœ™βƒ—πœŽβ‹…βƒ—π§/2 =cos⁑(πœ™/2)𝐈 βˆ’π‘–sin⁑(πœ™/2)βƒ—πœŽ ⋅⃗𝐧 for ⃗𝐧 βˆˆπ•Š3 (see Pauli matrices)
Proof of properties 1–5

Let 𝐓 ∈GL(n). Then

π‘’π“π€π“βˆ’1=βˆžβˆ‘π‘˜=0(π“π€π“βˆ’1)π‘˜π‘˜!=βˆžβˆ‘π‘˜=0π“π€π‘˜π“βˆ’πŸπ‘˜!=π“π‘’π€π“βˆ’1

proving ^P1

Let

𝐗(𝑑)=𝑒𝑑𝐀=βˆžβˆ‘π‘˜=0(𝑑𝐀)π‘˜π‘˜!

Then

˙𝐗(𝑑)=βˆžβˆ‘π‘˜=1π‘‘π‘˜βˆ’1π€π‘˜(π‘˜βˆ’1)!=π€βˆžβˆ‘π‘˜=0(𝑑𝐀)π‘˜π‘˜!=𝐀𝐗(𝑑)

and by the Existence and uniqueness theorem for IVPs this is unique, proving ^P2.

By Binomial expansion

𝑒𝑑𝐀𝑒𝑠𝐁=(βˆžβˆ‘π‘˜=0(𝑑𝐀)π‘˜π‘˜!)(βˆžβˆ‘β„“=0(𝑠𝐁)β„“β„“!)=βˆžβˆ‘π‘˜=0π‘˜βˆ‘β„“=0π‘‘β„“π‘ π‘˜βˆ’β„“π€β„“ππ‘˜βˆ’β„“β„“!(π‘˜βˆ’β„“)!=βˆžβˆ‘π‘˜=0π‘˜βˆ‘β„“=0π‘˜!β„“!(π‘˜βˆ’β„“)!π‘‘β„“π‘ π‘˜βˆ’β„“π€β„“πβ„“βˆ’π‘˜π‘˜!=βˆžβˆ‘π‘˜=0π‘˜βˆ‘β„“=0(π‘˜β„“)π‘‘β„“π‘ π‘˜βˆ’β„“π€β„“πβ„“βˆ’π‘˜π‘˜!=βˆžβˆ‘π‘˜=0(𝑑𝐀+𝑠𝐁)π‘˜π‘˜!=𝑒(𝑑𝐀+𝑠𝐁)

proving ^P3

By basic properties of the Conjugate transpose

(𝑒𝐀)†=(βˆžβˆ‘π‘˜=0π€π‘˜π‘˜!)†=βˆžβˆ‘π‘˜=0(π€π‘˜)β€ π‘˜!=βˆžβˆ‘π‘˜=0(𝐀†)π‘˜π‘˜!=𝑒(𝐀†)

proving ^P4.

Let 𝐃 be the Jordan canonical form so 𝐀 =π’πƒπ’βˆ’1. Then

det𝑒𝐀=det𝑒𝐃=π‘›βˆπ‘˜=1π‘’π·π‘˜π‘˜=𝑒tr⁑𝐃=𝑒tr⁑𝐀

proving ^P5

Generalisations


#state/develop | #lang/en | #SemBr