【線形代数】行列指数関数の求め方(2パターン)

はじめに

本問、意外と出題されたりします。現代制御での出題が最も多いですが、線形代数でも何かしらの初期値を与えて、そこからの時間変化をグラフに描く問題が出題されることがあります。

知っていれば得点できますので、是非この記事で復習いただければと思います。

本記事で覚えたいこと

固有値を用いた行列指数関数の求め方

行列$A$を対角化する作業とあまり変わりません。行列Dの部分が、exp項に変わったところに注意すれば、最終的に答えに行き着くのではないでしょうか。

ラプラス変換を用いた状態遷移行列の求め方

ひょっとすると、こちらのやり方のほうがシンプルかもしれません。

${\displaystyle \frac{dx(t)}{dt}}=Ax(t)$をラプラス変換すると、$sX(s)-x(0)=AX(s)$となるため

$$\begin{array}{}\text{(1)}& (sE-A)X(s)=x(0)\iff X(s)=(sE-A{)}^{-1}x(0)\iff x(t)=e^{At}x(0)\end{array}$$

になることから、上記の手順(CHECK)が発生しています。

解答例

固有値を用いた場合

特性方程式$|A-\lambda E|$について

$$\begin{array}{}\text{(2)}& \begin{array}{rl}\left|\begin{array}{cc}-\lambda & 1\\ -3& -4-\lambda \end{array}\right|& ={\lambda }^2+4\lambda +3=0\\ \lambda & =-3,-1\end{array}\end{array}$$

(i)$\lambda =-3$のとき、固有ベクトルは、$v_1=\left(\begin{array}{c}-\frac{1}{3}\\ 1\end{array}\right)$

(ii)$\lambda =-3$のとき、固有ベクトルは、$v_2=\left(\begin{array}{c}-1\\ 1\end{array}\right)$

これより、$P=\left(\begin{array}{cc}-\frac{1}{3}& -1\\ 1& 1\end{array}\right),P^{-1}={\displaystyle \frac{1}{2}}\left(\begin{array}{cc}3& 3\\ -3& -1\end{array}\right)$を得る。

よって、求める状態遷移行列$e^{At}$は

$$\begin{array}{}\text{(3)}& e^{At}& ={\displaystyle \frac{1}{2}}\left(\begin{array}{cc}-\frac{1}{3}& -1\\ 1& 1\end{array}\right)\left(\begin{array}{cc}{e}^{-3t}& 0\\ 0& e^{-t}\end{array}\right)\left(\begin{array}{cc}3& 3\\ -3& -1\end{array}\right)\text{(4)}& & ={\displaystyle \frac{1}{2}}\left(\begin{array}{cc}-\frac{1}{3}{e}^{-3t}& -e^{-t}\\ e^{-3t}& e^{-t}\end{array}\right)\left(\begin{array}{cc}3& 3\\ -3& -1\end{array}\right)\text{(5)}& & ={\displaystyle \frac{1}{2}}\left(\begin{array}{cc}-e^{-3t}+3e^{-t}& -e^{-3t}+e^{-t}\\ 3e^{-3t}-3e^{-t}& 3e^{-3t}-e^{-t}\end{array}\right)\end{array}$$

ラプラス変換を用いた場合

$(sE-A)=\left(\begin{array}{cc}s& -1\\ 3& s+4\end{array}\right)$だから、これの逆行列は

$$\begin{array}{}\text{(6)}& (sE-A{)}^{-1}& ={\displaystyle \frac{1}{s^2+4s+3}}\left(\begin{array}{cc}s+4& 1\\ -3& s\end{array}\right)\text{(7)}& & ={\displaystyle \frac{1}{2}}\left(\begin{array}{cc}-\frac{1}{s+3}+\frac{3}{s+1}& -\frac{1}{s+3}+\frac{1}{s+1}\\ \frac{3}{s+3}-\frac{3}{s+1}& \frac{3}{s+3}-\frac{1}{s+1}\end{array}\right)\end{array}$$

よって、逆ラプラス変換することで、行列指数関数

$$\begin{array}{}\text{(8)}& e^{At}={\displaystyle \frac{1}{2}}\left(\begin{array}{cc}-e^{-3t}+3e^{-t}& -e^{-3t}+e^{-t}\\ 3e^{-3t}-3e^{-t}& 3e^{-3t}-e^{-t}\end{array}\right)\end{array}$$

を得る。

最後に

参考文献(P22)にあるように、微分方程式を分解し、行列指数関数を求める第3の方法もあります。ただ、手計算で行うには大変でしたので、今回は割愛します。興味がある方は是非チェックしてください。

また、ラプラス変換を使用した方法の途中で、部分分数分解を行った計算があります。効率的に部分分数分解する方法を下記のページで紹介しています。宜しければチェック下さると幸いです。

ヘビサイドの展開公式を利用した部分分数分解

制御の問題を解いていると、手計算では逆ラプラス変換できない複雑な伝達関数が良く出てきます。 部分分数分解によって変換できる形に式変形する。が常套手段ですが、計算ミスが多く発生します。 本問では、ヘビサイドの展開定理を用い、複雑な伝達関数の部分分数分解を迅速に行う方法を紹介します。

denjoforest.com

2023.12.27

参考文献

http://shiwasu.ee.ous.ac.jp/in/03_web.pdf