行列指数関数・行列微分の概要
線形代数と解析の接点として重要なテーマの一つに 行列指数関数(matrix exponential) と 行列微分(matrix differentiation) があります。これらは、線形微分方程式の解法、制御理論、量子力学、マルコフ過程などの応用に広く現れます。
1. 行列指数関数
定義
スカラーの指数関数と同様に、行列指数関数は冪級数で定義されます。正方行列 A に対して、
eA=k=0∑∞k!1Ak
と定義されます。
主な性質
-
線形微分方程式の解
系 dtdx(t)=Ax(t) の解は
x(t)=eAtx(0)
で与えられる。
-
可換性の重要性
一般に eAeB=eA+B は AB=BA が成り立つ場合のみ成立。
-
対角化との関係
A=PDP−1 が成り立つ場合(対角化可能な場合)、
eA=PeDP−1,eD=diag(eλ1,…,eλn)
-
ジョルダン標準形との関係
一般の行列もジョルダン分解を通じて指数計算が可能。
2. 行列微分
基本的な考え方
行列関数 F(A) を行列 A に関して微分する場合、要素ごとの偏微分ではなく、行列全体に作用する線形写像としての導関数(フレシェ微分)を用いる。
-
スカラー値関数の場合
f:Rm×n→R のとき、
df=tr((∂A∂f)TdA)
で定義される。ここで ∂A∂f は 行列の勾配 を表す。
-
行列値関数の場合
F:Rm×n→Rp×q のとき、微分は dF=L(dA) として線形作用素で表す。
例
-
二次形式の微分
f(A)=xTAx⇒∂A∂f=xxT
-
トレースの微分
f(A)=tr(ATB)⇒∂A∂f=B
-
行列指数関数の微分
dtdeAt=AeAt=eAtA
3. 応用例
-
線形微分方程式系の解法:
dtdx(t)=Ax(t) の一般解を eAt で表現。
-
確率過程(マルコフ連鎖の連続時間版):
遷移行列 P(t)=eQt (ただし Q は生成行列)。
-
量子力学:
シュレーディンガー方程式の解にユニタリ行列 e−iHt が現れる。
-
最適化や機械学習:
損失関数の行列に関する勾配計算で行列微分を利用。
まとめ
-
行列指数関数は、スカラー指数関数を行列に拡張したものであり、特に線形微分方程式の解に現れる。
-
行列微分は、行列を変数とする関数に対して微分を拡張したものであり、フレシェ微分やトレースを通じて表現される。
-
両者は、制御理論、量子力学、確率論、機械学習など幅広い分野で重要な役割を果たす。