MATLABを使って行列の指数関数を計算

MATLABを使って前回の記事で考察した行列の指数関数を計算する。簡単に復習をすると行列の指数関数は

\begin{align}
e^A=\sum_{k=0}^{\infty}\frac{(P^{-1} A P)^k}{k!}=P^{-1} e^{A} P
\end{align}

で得られる。MATLABのコードはこれ

syms t
A = [-3 1;2 -2];
expm((-A*t))

結果

[    (2*exp(4*t))/3 + exp(t)/3,     exp(t)/3 - exp(4*t)/3]
[(2*exp(t))/3 - (2*exp(4*t))/3, exp(4*t)/3 + (2*exp(t))/3]

見ずらいが計算できている。