行列指数関数

線型代数学における行列の指数関数は、正方行列に対して定義される行列値関数で、通常の指数関数に対応するものである。より抽象的には、行列リー群とその行列リー代数の間の対応関係を行列の指数函数が記述する。

n}}次 実または複素正方行列 X}} の指数関数 esup または exp}} は、冪級数

で定義される n}}次正方行列である。この級数は任意の X}} に対して収束するから、行列 X}} の指数関数は well-defined である。

X}} が 1}}次正方行列のとき、X}}乗 esup は 1}}次正方行列であり、その唯一の成分は X}} の唯一の成分に対する通常の指数関数に一致する。これらはしばしば同一視される。この意味において行列の指数函数は、通常の指数函数の一般化である。

1. 性質
-----

//X//, //Y//}} を n}}次複素正方行列、//a//, //b//}} を複素数とし、n}}次単位行列を I}}、n}}次正方零行列を O}} でそれぞれ表すことにする。また、X}} の転置を //X//sup、共役転置を //X//sup と表すことにする。行列の指数関数は以下の性質を満たす：

//XY// =}} //YX//}} ならば
Y}} が正則ならば
exp(//X//sup. このことから X}} が対称行列ならばその行列乗 esup もまた対称であり、X}} が歪対称であるなら esup は直交行列になる。
exp(//X//sup. このことから X}} がエルミートならば esup もまたエルミートであり、X}} が歪エルミートならば esup はユニタリ行列になる。

1.1. 線型微分方程式
------------

行列の指数関数が重要であることの一つの理由として、常微分方程式系の解を求める際に使うことができることが挙げられる。以下の方程式

の解は、A}} を定行列として、次のように与えられる。

行列の指数関数はまた以下の様な非等質微分方程式に対しても有効である。

A'}} が定行列でないとき、

の形の微分方程式は解を閉じた形の式として陽に表すことはできないが、マグヌス級数が無限和の形で解を与える。

1.2. 和に対する指数函数
--------------

実数//x//, //y//}} について、通常の指数関数が //e//sup を満たすことはよく知られている。同じことは可換な行列に対しても成り立つ。即ち、行列 //X//, //Y//}} が交換可能ならば

が成り立つ。しかし可換でない行列については上記の関係は成り立たない。この場合、ベイカー＝キャンベル＝ハウスドルフの公式が //e//sup の計算に利用できる。

逆は一般には成り立たない。即ち、等式 //e//sup は X}} と Y}} が可換であることを意味しない。

エルミート行列について、行列指数関数の跡に関係する2つの注目すべき定理を挙げる。ゴールデン–トンプソン不等式 は以下の定理である。

定理 //A//, //H//}} がエルミートであるとき、次の不等式が成り立つ。ここで可換性は要求されないことに注意する。ゴールデン–トンプソン不等式を 3つの行列に対するものに拡張できないことを示す反例が知られている。そもそもエルミート行列 A, B, C}} に対して trexpexp)}} が実になること自体が保証されないのだが、次に示すリーブの定理はある意味でそのような保証を与える：

定理 固定されたエルミート行列 H}} について、関数は正定値行列錐上の凹関数であるLiebEpstein。

1.3. 指数写像
---------

複素行列の指数関数が常に正則行列であるということに注意する。これは複素変数の指数関数が常に零でないことに対応する事実である。ゆえに、行列の指数関数 n}}次正方行列の全体の成す空間から n}}次元の一般線型群への写像

を定めている。実は、この写像は全射、すなわちどんな正則行列も何らかの行列乗として書くことができる。

2つの行列 //X//, //Y//}} について

が成り立つ。ここで !!}} · !!}}}} は任意の行列ノルムである。ここから、指数写像はコンパクト部分集合 //Msub 上で連続かつリプシッツ連続であることが従う。//

写像

は //t// =}} 0}} で単位元を通る、一般線型群内の滑らかな曲線を定義する。実は

が成り立つから、これらは一般線型群の1パラメータ部分群を与えている。

この曲線の t}} 上の微分係数は で与えられる。//t// =}} 0}} での微分係数はまさに行列 X}} であり、これはつまり X}} がこの一径数部分群を生成することを示している。

より一般にWilcox、t}} に依存する生成的指数 //X//}} に対して

となる。右辺の //e//sup を積分記号の外へ出して、残った被積分関数をアダマールの補題を使って展開すれば、以下の有用な行列乗の微分係数の表示

が得られる。この式における係数はもとの指数函数の成分に現れているものとは異なることに注意せよ。また閉じた形の式は指数写像の微分を参照。

1.4. 行列の指数関数の行列式
----------------

ヤコビの公式から、任意の複素正方行列について次のトレース恒等式が成り立つ：

計算に役立つだけでなく、上記の等式の右辺は常に非零であるから、左辺の行列式は非零 det(//esup であり、したがって行列指数関数 esup は常に正則であることが分かる。//

実行列の場合、上記の公式から写像

が全射ではないことも分かる。なぜならば、実行列について公式の右辺は常に正であるが、行列式が負の正則行列は存在するからである。このことは先に触れた複素行列の場合とは対照的である。

2. 指数函数の計算
----------

一般の行列乗の計算を確度と精度を以って行うことは非常に難しく、現在においても数学、特に数値解析において重要な研究トピックの一つである。MATLABやGNU Octaveはパデ近似を使っている。

いくつかの行列のクラスに関しては、比較的容易に計算ができる。

2.1. 対角行列の場合
------------

対角行列

に対して、行列 A}}乗は単に主対角成分のそれぞれを肩に載せた

で与えられる。これは対角行列同士の行列の積は単に成分ごとの積に等しいということからの帰結である。特に通常の指数函数は「一次元」の場合の対角行列の指数函数とみなせる。

これを利用すれば対角化可能行列乗も計算できる。つまり //A// =}} //UDU//sup かつ D}} が対角行列ならば

である。シルベスターの公式を応用しても同じ結果が得られる。

2.2. 正射影行列の場合
-------------

考える行列が射影行列ならば、これは冪等だから、行列乗は

となることが指数函数の定義より容易に分かる。実際、冪等性により //Psup だから、//

である。

2.3. 冪零行列の場合
------------

冪零行列 N}} は適当な正整数 q}} に対して //Nsup を満たす。N}}乗 esup は指数函数の定義級数から直接に//

と計算できる。

2.4. より一般の場合
------------

行列 X}} に対してその最小多項式が一次式の積に分解されるとき、行列 X}} は

A}}：対角化可能
N}}：冪零
A}} と N}} は可換 なる形に書くことができる。このとき X}}乗の計算は

により、先の対角化可能行列および冪零行列の計算に帰着される。後の等号で A}} と N}} との可換性が必要であることに注意せよ。

同様の方法は、代数閉体上の行列に対してジョルダン標準形を取ることで与えられる。即ち J}} が X}} のジョルダン標準形で //X// =}} //PJP//sup と書くとき、

である。ジョルダン細胞の直和として

と書けば、

となるから、後はジョルダン細胞乗が計算できればよい。各ジョルダン細胞は特別な形をした冪零行列 N}} を用いて

なる形に書けるのだから、

が得られる。

2.5. ローラン級数による評価
----------------

ケイリー・ハミルトンの定理を考えれば、n}}次正方行列乗はその行列の高々次数 //n// − 1}} の多項式として表示できるはずである。

非零な一変数多項式 P}} および Qsub は //P// =}} 0}} なるものとする。有理型函数

が整函数ならば

が成り立つ。これを示すには上記等式において //P//}} を掛けて z}} を A}} で置き換えればよい。

さてこのような多項式 //Qsub は以下のように見つけることができる。a}} は P}} の根として、 //Q//sub は P}} に f}} の a}} におけるローラン級数の主要部を掛けることで得られる。これは関連するフロベニウス共変行列に比例する。a}} が P}} の根を亙るときの //Q//sub 全ての和 Ssub が所期の Qsub として取れる。他全ての Qsub は //Ssub に P}} の定数倍を加えることで得られる。特に、ラグランジュ–シルヴェスター多項式 //Ssub は P}} より次数が低くなる唯一の Qsub である。//

3. 行列の行列乗
---------

行列の指数函数と行列の対数函数が既知であるならば、正規かつ正則な n}}次正方行列 X}} と n}}次複素正方行列 Y}} に対して、行列の行列乗 を

と定義することができる。ここに、行列の乗法は非可換であるから、行列の行列乗も左冪 sup と右冪 Xsup の別が生じることに注意せよ。さらに言えば、

X}} が正規かつ正則ならば、Xsup と sup は固有値集合が一致する。
X}} が正規かつ正則で、Y}} が正規であり、かつ //XY// =}} //YX//}} が成り立つならば、//Xsup が成り立つ。//
X}} が正規かつ正則で、//X//, //Y//, //Z//}} がどの2つも可換ならば、//X//sup が成り立つ。

4. 応用
-----

連立常微分方程式の数値解法であるexponential integratorの研究においては、行列指数関数は重要視されている。

5. 脚注
-----

6. 参考文献
-------

7. 関連項目
-------

行列値関数
行列の対数函数
指数関数
リー環の指数写像
リー・トロッター積公式

8. 外部リンク
--------

[1]