matrix exponential

The exponential of a real valued square matrix $A$, denoted by $e^A$, is defined as

$e^A$ $=$ ${\displaystyle \sum _{k=0}^{\mathrm{\infty }}}{\displaystyle \frac{1}{k!}}{A}^k$

$=$ $I+A+{\displaystyle \frac{1}{2}}{A}^2+\mathrm{\cdots }$

Let us check that $e^A$ is a real valued square matrix. Suppose $M$ is a real number such for all entries $A_{ij}$ of $A$. Then for all entries in $A^2$, where $n$ is the order of $A$. (Alternatively, one could argue using matrix norms: We have $||e^A||\le e^{||A||}$ for the 2-norm, and hence the entries of $e^A$ are bounded by $M=||e^A||$.) Thus, in general, we have . Since ${\sum }_{k=0}^{\mathrm{\infty }}\frac{n^k}{k!}{M}^{k+1}$ converges, we see that $e^A$ converges to real valued $n\times n$ matrix.

Example 1. Suppose $A$ is nilpotent, i.e., $A^r=0$ for some natural number $r$. Then

$e^A$ $=$ $I+A+{\displaystyle \frac{1}{2!}}{A}^2+\mathrm{\cdots }+{\displaystyle \frac{1}{(r-1)!}}{A}^{r-1}.$

Example 2. If $A$ is diagonalizable, i.e., of the form $A=LD{L}^{-1}$, where $D$ is a diagonal matrix, then

$e^A$ $=$ ${\displaystyle \sum _{k=0}^{\mathrm{\infty }}}{\displaystyle \frac{1}{k!}}{(LD{L}^{-1})}^k$

$=$ ${\displaystyle \sum _{k=0}^{\mathrm{\infty }}}{\displaystyle \frac{1}{k!}}L{D}^k{L}^{-1}$

$=$ $L{e}^D{L}^{-1}.$

Further, if $D=diag\{a_1,\mathrm{\cdots },a_n\}$, then $D^k=diag\{a_1^k,\mathrm{\cdots },a_n^k\}$ whence

$e^A$ $=$ $Ldiag\{e^{a_1},\mathrm{\cdots },e^{a_n}\}{L}^{-1}.$

For diagonalizable matrix $A$, it follows that $det{e}^A=e^{trA}$. However, this formula is, in fact, valid for all $A$.

Let $A$ be a square $n\times n$ real valued matrix. Then the matrix exponential satisfies the following properties

1. For the $n\times n$ zero matrix $O$, $e^O=I$, where $I$ is the $n\times n$ identity matrix.

2. If $A=Ldiag\{a_1,\mathrm{\cdots },a_n\}{L}^{-1}$ for an invertible $n\times n$ matrix $L$, then

   $$e^A=Ldiag\{e^{a_1},\mathrm{\cdots },e^{a_n}\}{L}^{-1}.$$

3. If $A$ and $B$ commute, then $e^{A+B}=e^A{e}^B$.

4. The trace of $A$ and the determinant of $e^A$ are related by the formula

   $$det{e}^A=e^{trA}.$$

   In effect, $e^A$ is always invertible. The inverse is given by

   $${(e^A)}^{-1}=e^{-A}.$$

5. If $e^A$ is a rotational matrix, then $trA=0$.

A relevant example on property 3.
We report an interesting example where the cited property is valid. In the field of complex numbers consider the complex matrix

$$C=A+iB,$$ (1)

being $C$ hermitian, i.e. $C^{⊺}=\overline{C}$ (here ”$⊺$” and overline ”$-$” stand for tranposition and conjugation, respectively) and orthogonal, i.e $C^{-1}=C^{⊺}$ . From (1),

$$C^{⊺}=A^{⊺}+i{B}^{⊺}.$$

Since $C$ is orthogonal, from the complex equation $C{C}^{⊺}=I$ ($I$ is the identity matrix), we have

$$C{C}^{⊺}=(A+iB)(A^{⊺}+i{B}^{⊺})=(A{A}^{⊺}-B{B}^{⊺})+i(B{A}^{⊺}+A{B}^{⊺})=I,$$

whence the imaginary part leads to the equation

$$B{A}^{⊺}+A{B}^{⊺}=0.$$ (2)

But $C$ is also hermitian, so that

$$C^{⊺}=A^{⊺}+i{B}^{⊺}=\overline{C}=A-iB,$$

therefore $A^{⊺}=A$ is symmetric, and $B^{⊺}=-B$ is skew-symmetric. From these and (2), $BA=AB$, and this implies that $\exp (A)\cdot \exp (B)=\exp (A+B)$. So that, the real and imaginary parts of an orthogonal and hermitian matrix verifies the property. Likewise, it is easy to show that if the complex matrix is symmetric and unitary, its real an imaginary components also verify this property.

Title	matrix exponential
Canonical name	MatrixExponential
Date of creation	2013-03-22 13:33:27
Last modified on	2013-03-22 13:33:27
Owner	mathcam (2727)
Last modified by	mathcam (2727)
Numerical id	13
Author	mathcam (2727)
Entry type	Definition
Classification	msc 15A15
Classification	msc 15-00
Related topic	ProofOfEquivalenceOfFormulasForExp