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:

    “The vulgar call good fortune that which really is produced by the calculations of genius.”
    —Ralph Waldo Emerson (1803–1882)

    “Now, since our condition accommodates things to itself, and transforms them according to itself, we no longer know things in their reality; for nothing comes to us that is not altered and falsified by our Senses. When the compass, the square, and the rule are untrue, all the calculations drawn from them, all the buildings erected by their measure, are of necessity also defective and out of plumb. The uncertainty of our senses renders uncertain everything that they produce.”
    —Michel de Montaigne (1533–1592)

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