MATLABを使って前回の記事で考察した行列の指数関数を計算する。簡単に復習をすると行列の指数関数は
\begin{align} e^A=\sum_{k=0}^{\infty}\frac{(P^{-1} A P)^k}{k!}=P^{-1} e^{A} P \end{align}
で得られる。MATLABのコードはこれ
1 2 3 | syms t A = [-3 1;2 -2]; expm((-A*t)) |
結果
1 2 | [ (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] |
見ずらいが計算できている。
コメント