logo

行列の指数関数と対数関数 📂行列代数

行列の指数関数と対数関数

定義

指数関数 exp\exp対数関数 log\log行列 に対して一般化しようとしている。

行列の指数

指数関数を行列に対して一般化した exp:Cn×nCn×n\exp : \mathbb{C}^{n \times n} \to \mathbb{C}^{n \times n} は行列 ACn×nA \in \mathbb{C}^{n \times n} に対して次のように定義される。 expA:=k=0Akk! \exp A := \sum_{k=0}^{\infty} \frac{A^{k}}{k!} expA\exp A は簡単に eAe^{A} とも表され、AA行列指数matrix exponential という。

行列の対数

対数関数を行列に対して一般化した log:Cn×nCn×n\log : \mathbb{C}^{n \times n} \to \mathbb{C}^{n \times n}expA=B\exp A = B を満たす行列 AABB が存在する時に AAlogB\log B と表し、 AABB行列対数matrix logarithm という。

定理 1

もし二つの行列 A,BCn×nA, B \in \mathbb{C}^{n \times n}AB=BAAB = BA を満たせば、次が成立する。 logAB=logA+logB2πiU(logA+logB) \log AB = \log A + \log B - 2 \pi i \mathcal{U} \left( \log A + \log B \right) 特に AA固有値 μk\mu_{k}BB の固有値 νk\nu_{k} に対して次の系を得られる。 logAB=logA+logB    k=1,,n:argμk+argνk(π,π] \log AB = \log A + \log B \iff \forall k = 1 , \cdots , n : \arg \mu_{k} + \arg \nu_{k} \in ( - \pi , \pi ] ここで arg\arg複素数の偏角U\mathcal{U} は次のように定義された 行列アンワインディング関数matrix unwinding function である。 U(A):=AlogeA2πi,ACn×n \mathcal{U} (A) := {\frac{ A - \log e^{A} }{ 2 \pi i }} \qquad , A \in \mathbb{C}^{n \times n}

解説

行列の指数と対数は 複素数 で定義される指数と対数を形式的に拡張したものである。見ての通り、行列の指数は行列の冪級数を用いて定義される。

エルミートおよび正定値行列の空間において

エルミートおよび正定値行列の空間:

  • エルミート行列の空間: サイズが n×nn \times nエルミート行列集合 は次のように表される。スカラー体 R\mathbb{R} に対して Hn\mathbb{H}_{n}ベクトル空間 である。 Hn:={ACn×n:A=A} \mathbb{H}_{n} := \left\{ A \in \mathbb{C}^{n \times n} : A = A^{\ast} \right\}
  • 正定値行列の集合: サイズが n×nn \times n正定値行列 の集合を Pn\mathbb{P}_{n} と表す。 PnHn\mathbb{P}_{n} \subset \mathbb{H}_{n}Hn\mathbb{H}_{n}凸錐 である。

すべての行列の空間ではなく、エルミート行列の空間 Hn\mathbb{H}_{n} が定義域で正定値行列の凸錐 Pn\mathbb{P}_{n} が値域である exp:HnPn\exp : \mathbb{H}_{n} \to \mathbb{P}_{n}全単射 であり、その逆関数は log:PnHn\log : \mathbb{P}_{n} \to \mathbb{H}_{n} である。特にこれらは 微分同型写像 である。

HP.png

関連項目


  1. Aprahamian, M., & Higham, N. J. (2014). The matrix unwinding function, with an application to computing the matrix exponential. SIAM Journal on Matrix Analysis and Applications, 35(1), 88-109. https://doi.org/10.1137/130920137 Lemma 3.12 ↩︎