User:AugPi/sandbox

In mathematics, the matrix exponential is a matrix function on square matrices analogous to the ordinary exponential function. It is used to solve systems of linear differential equations. In the theory of Lie groups, the matrix exponential gives the connection between a matrix Lie algebra and the corresponding Lie group.

Let X be an n×n real or complex matrix. The exponential of X, denoted by e^X or exp(X), is the n×n matrix given by the power series

$e^{X}=\sum _{k=0}^{\infty }{1 \over k!}X^{k}$

where $X^{0}$ is defined to be the identity matrix $I$ with the same dimensions as $X$.

The above series always converges, so the exponential of X is well-defined. If X is a 1×1 matrix the matrix exponential of X is a 1×1 matrix whose single element is the ordinary exponential of the single element of X.

Properties

Elementary properties

Let X and Y be n×n complex matrices and let a and b be arbitrary complex numbers. We denote the n×n identity matrix by I and the zero matrix by 0. The matrix exponential satisfies the following properties.

We begin with the properties that are immediate consequences of the definition as a power series:

• e⁰ = I
• exp(X^T) = (exp X)^T, where X^T denotes the transpose of X.
• exp(X^∗) = (exp X)^∗, where X^∗ denotes the conjugate transpose of X.
• If Y is invertible then e^YXY−1 = Ye^XY⁻¹.

The next key result is this one:

• If $XY=YX$ then $e^{X}e^{Y}=e^{X+Y}$.

The proof of this identity is the same as the standard power-series argument for the corresponding identity for the exponential of real numbers. That is to say, as long as $X$ and $Y$ commute, it makes no difference to the argument whether $X$ and $Y$ are numbers or matrices. It is important to note that this identity typically does not hold if $X$ and $Y$ do not commute.

Consequences of the preceding identity are the following:

• e^aXe^bX = e^{(a + b)X}
• e^Xe^−X = I

Using the above results, we can easily verify the following claims. If X is symmetric then e^X is also symmetric, and if X is skew-symmetric then e^X is orthogonal. If X is Hermitian then e^X is also Hermitian, and if X is skew-Hermitian then e^X is unitary.

Finally, we have the following:

• A Laplace transform of matrix exponentials amounts to the resolvent,

$\int _{0}^{\infty }e^{-ts}e^{tX}\,dt=(sI-X)^{-1}$

for all sufficiently large positive values of s.

Linear differential equation systems

Main article: matrix differential equation

One of the reasons for the importance of the matrix exponential is that it can be used to solve systems of linear ordinary differential equations. The solution of

${\frac {d}{dt}}y(t)=Ay(t),\quad y(0)=y_{0},$

where A is a constant matrix, is given by

$y(t)=e^{At}y_{0}.\,$

The matrix exponential can also be used to solve the inhomogeneous equation

${\frac {d}{dt}}y(t)=Ay(t)+z(t),\quad y(0)=y_{0}.$

See the section on applications below for examples.

There is no closed-form solution for differential equations of the form

${\frac {d}{dt}}y(t)=A(t)\,y(t),\quad y(0)=y_{0},$

where A is not constant, but the Magnus series gives the solution as an infinite sum.

The determinant of the matrix exponential

By Jacobi's formula, for any complex square matrix the following trace identity holds:

$\det(e^{A})=e^{\operatorname {tr} (A)}~.$

These are the generators of the SU(3) group in the triplet representation, and they are normalized as

$\operatorname {tr} (\lambda _{j}\lambda _{k})=2\delta _{jk}.$

The Lie algebra structure constants of the group are given by the commutators of $\lambda _{k}$

$[\lambda _{j},\lambda _{k}]=2if_{jkl}\lambda _{l}~,$

where $f_{jkl}$ are the structure constants completely antisymmetric and are analogous to the Levi-Civita symbol $\epsilon _{jkl}$ of SU(2).

In general, they vanish, unless they contain an odd number of indices from the set {2,5,7}, corresponding to the antisymmetric λs. Note $f_{ljk}={\frac {-i}{2}}\mathrm {tr} ([\lambda _{l},\lambda _{j}]\lambda _{k})$.

Moreover,

$\{\lambda _{j},\lambda _{k}\}={\frac {4}{3}}\delta _{jk}+2d_{jkl}\lambda _{l}$

where $d_{jkl}$ are the completely symmetric coefficient constants. They vanish if the number of indices from the set {2,5,7} is odd.