Matrix Exponential - Calculations

Calculations

Suppose that we want to compute the exponential of

B=\begin{bmatrix} 21 & 17 & 6 \\ -5 & -1 & -6 \\ 4 & 4 & 16 \end{bmatrix}.

Its Jordan form is

J = P^{-1}BP = \begin{bmatrix} 4 & 0 & 0 \\ 0 & 16 & 1 \\ 0 & 0 & 16 \end{bmatrix},

where the matrix P is given by

P=\begin{bmatrix} -\frac14 & 2 & \frac54 \\ \frac14 & -2 & -\frac14 \\ 0 & 4 & 0 \end{bmatrix}.

Let us first calculate exp(J). We have

The exponential of a 1×1 matrix is just the exponential of the one entry of the matrix, so exp(J1(4)) = . The exponential of J2(16) can be calculated by the formula exp(λ + N) = eλ exp(N) mentioned above; this yields

\begin{align} \exp \left( \begin{bmatrix} 16 & 1 \\ 0 & 16 \end{bmatrix} \right) & = e^{16} \exp \left( \begin{bmatrix} 0 & 1 \\ 0 & 0 \end{bmatrix} \right) \\ & = e^{16} \left(\begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix} + \begin{bmatrix} 0 & 1 \\ 0 & 0 \end{bmatrix} + {1 \over 2!}\begin{bmatrix} 0 & 0 \\ 0 & 0 \end{bmatrix} + \cdots \right) = \begin{bmatrix} e^{16} & e^{16} \\ 0 & e^{16} \end{bmatrix}. \end{align}

Therefore, the exponential of the original matrix B is

\begin{align} \exp(B) & = P \exp(J) P^{-1} = P \begin{bmatrix} e^4 & 0 & 0 \\ 0 & e^{16} & e^{16} \\ 0 & 0 & e^{16} \end{bmatrix} P^{-1} \\ & = {1\over 4} \begin{bmatrix} 13e^{16} - e^4 & 13e^{16} - 5e^4 & 2e^{16} - 2e^4 \\ -9e^{16} + e^4 & -9e^{16} + 5e^4 & -2e^{16} + 2e^4 \\ 16e^{16} & 16e^{16} & 4e^{16} \end{bmatrix}. \end{align}

Read more about this topic:  Matrix Exponential

Famous quotes containing the word calculations:

    “Nowhere are our calculations more frequently upset than in war.”
    —Titus Livius (Livy)

    “What, really, is wanted from a neighborhood? Convenience, certainly, an absence of major aggravation, to be sure. But perhaps most of all, ideally, what is wanted is a comfortable background, a breathing space of intermission between the intensities of private life and the calculations of public life.”
    —Joseph Epstein (b. 1937)

    “Heaven’s calculations don’t follow man’s calculations.”
    —Chinese proverb.