Mathwords logoMathwords

Matrix Exponential — Definition, Formula & Examples

The matrix exponential eAe^A extends the ordinary exponential function to square matrices. It is defined as an infinite series and is the key tool for solving systems of linear differential equations with constant coefficients.

For an n×nn \times n matrix AA, the matrix exponential is defined by the convergent power series eA=k=0Akk!e^A = \sum_{k=0}^{\infty} \frac{A^k}{k!}, where A0=IA^0 = I (the identity matrix). This series converges for every square matrix AA over R\mathbb{R} or C\mathbb{C}.

Key Formula

eA=k=0Akk!=I+A+A22!+A33!+e^{A} = \sum_{k=0}^{\infty} \frac{A^k}{k!} = I + A + \frac{A^2}{2!} + \frac{A^3}{3!} + \cdots
Where:
  • AA = An n × n square matrix
  • II = The n × n identity matrix (= A⁰)
  • kk = Summation index running from 0 to infinity

How It Works

To compute eAe^A in practice, you rarely sum the series term by term. For diagonalizable matrices, write A=PDP1A = PDP^{-1} where DD is diagonal; then eA=PeDP1e^A = P\, e^D\, P^{-1}, and eDe^D is simply the diagonal matrix with entries edie^{d_i}. For non-diagonalizable matrices, you use the Jordan normal form instead. The matrix exponential satisfies ddtetA=AetA\frac{d}{dt} e^{tA} = A\, e^{tA}, which is why the general solution to x=Ax\mathbf{x}' = A\mathbf{x} with initial condition x(0)=x0\mathbf{x}(0) = \mathbf{x}_0 is x(t)=etAx0\mathbf{x}(t) = e^{tA}\mathbf{x}_0.

Worked Example

Problem: Compute the matrix exponential eAe^A where A=(0000)A = \begin{pmatrix} 0 & 0 \\ 0 & 0 \end{pmatrix} replaced by a more interesting case: A=(1002)A = \begin{pmatrix} 1 & 0 \\ 0 & 2 \end{pmatrix}.
Recognize the structure: A is already a diagonal matrix with diagonal entries 1 and 2.
A=(1002)A = \begin{pmatrix} 1 & 0 \\ 0 & 2 \end{pmatrix}
Apply the exponential to each diagonal entry: For a diagonal matrix, the matrix exponential is the diagonal matrix of the exponentials of the entries.
eA=(e100e2)e^A = \begin{pmatrix} e^1 & 0 \\ 0 & e^2 \end{pmatrix}
Write the result: Evaluate the scalar exponentials.
eA=(e00e2)(2.718007.389)e^A = \begin{pmatrix} e & 0 \\ 0 & e^2 \end{pmatrix} \approx \begin{pmatrix} 2.718 & 0 \\ 0 & 7.389 \end{pmatrix}
Answer: eA=(e00e2)e^A = \begin{pmatrix} e & 0 \\ 0 & e^2 \end{pmatrix}

Why It Matters

The matrix exponential is central to solving systems of linear ODEs such as x=Ax\mathbf{x}' = A\mathbf{x}, which arise in physics, engineering, and control theory. It also appears in Lie group theory, quantum mechanics (the time-evolution operator eiHt/e^{-iHt/\hbar}), and network analysis for computing diffusion on graphs.

Common Mistakes

Mistake: Assuming eA+B=eAeBe^{A+B} = e^A e^B for all matrices.
Correction: This identity holds only when AA and BB commute (AB=BAAB = BA). In general, eA+BeAeBe^{A+B} \neq e^A e^B. Use the Baker–Campbell–Hausdorff formula or other methods for non-commuting matrices.