行列の平方根

数学のおもに線型代数学および函数解析学における行列の平方根は、数に対する通常の平方根の概念を行列に対して拡張するものである。すなわち、行列 B}} が行列 A}} の平方根であるとは、行列の積に関して //B//exp が A}} に等しいときに言う。

「実数の平方根は必ずしも実数にならないが、複素数は必ず複素数の範囲で平方根を持つ」ことに対応する事実として、実行列の平方根は必ずしも実行列にならないが、複素行列が平方根を持てばそれは必ず複素行列の範囲で取れる。

平方根を持たない行列も存在する例えば }}。

また一般に、ひとつの行列が複数の平方根を持ち得るたとえば、行列 は行列 およびこれらの[[加法逆元。実際、2 × 2}} 単位行列は次のように無数の平方根を持つ。

このように行列の平方根は無数に存在しうるが、半正定値行列これはふつう、[[対称行列の範疇で行列の主平方根  の概念が定義できて「半正定値行列の主平方根はただ一つ」である。

2 × 2}} 行列が、相異なる二つの非零固有値を持つならば、それは四つの平方根を持つ。実際に、そのような仮定を満たす行列 A}} は A}} の固有ベクトルを列ベクトルに持つ行列 V}} とそれに対応する固有値を対角成分に持つ対角行列 D}} を用いて 1=//A// = //VDV//−1}} と固有値分解できるから、A}} の平方根は //VD½ V//−1}} で与えられることがわかる。ただし、//D//exp は D}} の任意の平方根で、それは D}} の対角成分の任意の平方根を同じ位置の対角成分として持つ対角行列であり、その選び方は 2exp 通りある。同じ理由で、上で述べた「半正定値行列の主平方根がただ一つに定まる」ことも言える—半正定値行列正定値行列となるための必要十分条件はそのすべての固有値が正となることであった}} A}} の全ての非負固有値の主平方根を対角成分に持つ対角行列を //D//exp とする行列 //VDsup はただ一つしかない。//

適当な冪零行列 N}} を用いて //I// + //N//}} の形に書ける行列の平方根 sup は、二項級数に対する汎函数計算で求められる。同様に、行列の指数函数 exp}}, 対数函数 log}} が既知ならば、 exp)}} を A}} の平方根とすることができる。

1. 定義
-----

定義 : 行列 B}} が行列 A}} の平方根であるとは、//B//exp を満たすときに言う。 定義 「非負実数が非負の平方根をただ一つだけ持つ」という事実に対応して

命題

半正定値行列は、それ自身が半正定値となるような平方根をただ一つ持つ。
一般に、すべての固有値が正の実数となる複素行列はすべての固有値が正の実数となる平方根をただ一つ持つ。 が成り立つ。そのように定まるただ一つの  平方根は主平方根  と呼ばれる。主平方根をとる操作は行列全体の成す集合上で連続である。このとき、考えている行列が実行列ならば、その主平方根もまた実行列になる。主平方根に関する性質は、行列に対する正則汎函数計算の帰結として得られる。あるいは主平方根の存在と一意性はジョルダン標準形を用いて直截に示せる。

注意: 記号 sqrt や •sup は、主平方根を表すために用いる場合や、平方根の任意の一つを表すために用いる場合などがあるので、文脈に注意すべきである。

2. 計算法
------

2.1. 明示公式
---------

2 × 2}} 行列の場合は、すべての成分を明示的に計算することによって平方根を求めることはそう難しくない。固有値が退化していない場合の平方根は二次正方行列の平方根として記述できる。

すなわち、とし、その行列式を、 特性方程式の判別式をとしたとき、

ならば、の平方根は、

、、、と明示的に表記できる。

平方根となることは、実際に2乗を計算すればから容易にわかる。

あるいは、2次のケイリー・ハミルトンの定理から、としても良い。

これら以外に平方根が存在しないことについては、とした場合、よりは2つの相異なる固有値、と、独立な固有ベクトル、を持つが、任意の2次列ベクトルは、、の1次結合で表せるので、、とすると、、すなわち、であるが、のため、解は、、に定まる。これにより任意の2次列ベクトルがによりどう変換されるかが定まるが、これはが定まることを意味する。が固有値ゼロを持たない場合は解が4組、固有値ゼロを持つ場合は解が2組であるが、これは上記の明示公式で尽くされているので、これら以外には、平方根は存在しない。

の場合は、複雑になる。

D}} が //n// × //n//}} 対角行列ならば、D}} の対角成分の任意の平方根を対応する位置の対角成分に持つ対角行列 R}} を作れば平方根が得られる。D}} の対角成分が非負の実数ならば、先の対角行列 R}} で各成分の符号を全て正としたものは D}} の主平方根である。

冪等行列の平方根は、自身を平方根に持つ。

2.2. 対角化の利用
-----------

対角化可能行列 A}} に対し、適当な行列 V}} と対角行列 D}} が存在して //A// =}} //VDV//−1}} と書ける。これは A}} が Csup を張る n}} 個の固有値を持つことと同値である。このとき V}} はその列ベクトルが n}} 個の固有ベクトルであるように選べる。そうして A}} の平方根は D}} の任意の平方根を用いて と書ける。実際、 である。 A}} が対称行列ならば対角化に用いる行列 V}} は固有ベクトルを適当に選んで直交行列となるようにとれる。この場合、V}} の逆行列はたんに転置をとるだけであるから、 と書ける。

2.3. ジョルダン分解の利用
---------------

正方行列 のジョルダン標準形 をとすると、次が言える。

をの平方根とすると、は、より、の平方根となる。逆にをの平方根とすると、は、より、の平方根であり、である。このため、ジョルダン標準形の全ての平方根を知ることができれば、により、の全ての平方根を知ることができる。

とし、とすれば、は、の平方根のうちの一つである。

逆に、、ただしはジョルダン標準形で、とは共通の固有値を持たないとすると、の平方根は、ただし、に限られる。

これは、とすると、

よりとなるが、、の対角成分をと置き、第１列に注目すれば、だが、とは共通の固有値を持たないため、が言え、順次、第２列、第３列に注目すればが言え、が言える。

からも同様に、と置き、第k行に注目すればだが、とは共通の固有値を持たないため、が言え、順次、第k-1行、第k-2行に注目すればが言え、が言える。このため、上記が言える。

ジョルダン標準形の平方根には、ジョルダン細胞の平方根であるものと、

のようにジョルダン細胞の平方根ではないものがあるので、注意が必要である。

2.3.1. ジョルダン細胞の平方根
------------------

ジョルダン細胞とはn}}次正方行列で、のとき、、、のときとなるものを言う。

のとき、ジョルダン細胞の平方根は、下記の行列およびである。

のとき、、のときのとき、ジョルダン細胞は、

の場合、平方根0を持つの場合、平方根を持たない例 は平方根を持たない。

のとき、ジョルダン細胞の平方根が２つしかないことは、次から言える。となる行列が存在したとし、の成分を考える。

を比較すると、このため

を比較すると、このため

このため、は上三角行列で、斜めに同じ値が並ばなければならない。の成分を比較することにより、が言え、以下成分を比較することにより、の全ての成分が順番に１次方程式で定まるため、平方根が２つしかないことが言える。

2.3.2. 英語版からの直訳
---------------

対角化可能でない行列の場合にはジョルダン標準形が利用できる。行列の対数函数#非対角化可能行列の対数の項と同様の級数展開を用いる方法}}

すべての固有値が正の実数であるような任意の複素行列が、同じ条件の平方根を持つことを見るには、ジョルダンブロックの場合に証明すれば十分である。そのようなブロックは実数 //λ// > 0}} および冪零行列 N}} を用いて //λ//}} の形に書ける。平方根の二項級数展開 1=1/2 = 1 + //a//1 //z// + //a//2 //z//2 + ⋯}}に対し、形式冪級数としての平方は 1 + //z//}} に等しい。z}} を N}} に置き換えれば、冪零性により有限個を除く全ての項は零となり、1=//S// = sqrt が固有値 1=sqrt に属するジョルダンブロックの平方根を与える。

一意性を見るには 1=//λ// = 1}} の場合に確認すれば十分である。上で構成した平方根を 1=//S// = //I// + //L//}} の形に書けば、L}} は定数項を持たない N}} の多項式である。固有値が正の実数となる他の任意の平方根 T}} は 1=//T// = //I// + //M//}} の形で M}} が冪零かつ N}} と可換となるようにとれる。しかしこのとき 1=0 = //S//2 − //T//2 = 2/2)}} であり、また L}} と M}} の可換性により 1=//L// + //M//}} は冪零ゆえ 1=//I// + /2}} は可逆となるから、したがって 1=//L// = //M//}}.

すべての固有値が正の実数であるような行列 A}} の最小多項式を //p//}} とするとき、A}} の一般固有空間へのジョルダン分解は //p//sup の部分分数分解から導かれる。すなわち、対応する一般固有空間の上への射影は A}} の実係数多項式として与えられ、各固有空間上で A}} は上記の通り //λ//}} の形をしている。固有空間上での平方根の冪級数展開は、A}} の主平方根が実係数多項式 //q//}} に対する //q//}} の形をしていることを示すものである。

2.4. 現実的な計算法
------------

「対角化」の方法でも「ジョルダン分解」の方法でも、すべての固有値を算出することが必要となるが、それは行列の特性方程式のすべての解を求めることと同じであり、行列の次数が大きくなれば非現実的となる。このため、現実的な平方根の求め方が必要となる。

2.4.1. 行列対数関数、行列指数関数による求め方
--------------------------

実数の平方根がで求まることと同様に、

n次実数値正方行列の全ての特性根の実数部分が正である場合、

行列対数関数をと定義し

行列指数関数をと定義すれば、

2乗するととなり、かつ全ての特性根の実数部分が正となる行列は、

により計算でき、かつこの行列に一意に定まる。

この方法は、固有値を全て求める必要がないこと、収束計算が速いこと、対称行列に限らず一般の行列に利用可能であることなど、現実的かつ速い計算方法になっている。

また、行列の平方根に限らず、ｎ乗根も同様に計算することができる。

2.4.2. ニュートン法
-------------

実数の方程式をニュートン法で解く方法を、行列にそのまま適用して求める方法である。

ｎ次正方行列に対し、ｎ次正方行列の列を次の漸化式で定める

この列が適当な初期値について収束すれば、収束値について、となる。

このことは、収束すればが成り立つことから言える。

3. 対称行列に限定した議論
--------------

以下では、対称行列に限定した行列の平方根についての性質を示す。 「正定値行列」とは、対称行列で、その全ての固有値が正の実数であるものをいう。「半正定値行列」とは、対称行列で、その全ての固有値がゼロまたは正の実数であるものをいう。

3.1. 定義
-------

転置あるいはエルミート共軛を用いれば、より一般に非対称あるいは非エルミートな矩形行列の範疇で「平方根」をとることができる。

定義: 半正定値実正方行列 A}} に対して、1=//A// = //B   supを満たす任意の矩形行列 B}} を A}} の非対称平方根  と呼ぶ。// 定義: 半正定値複素正方行列 A}} に対して、1=//A// = //BB//*}}を満たす任意の矩形行列 B}} を A}} の非エルミート平方根  と呼ぶ。B}} がエルミートならば、B}} は上で述べたA}}の平方根と一致する。任意の正定値エルミート行列 A}} に対し、それ自身正定値エルミートとなる平方根は一意であり、これを主平方根 と呼ぶ。

注: コレスキー分解からも平方根の例が得られるが、コレスキー因子と平方根とを混同してはならない。

3.2. 非対称平方根のユニタリ自由度
-------------------

正実数の平方根は、主平方根に ±1}} を掛けたものですべて与えられた。これに対応するように、正定値エルミート行列の任意の非エルミート平方根は、ユニタリ変換によって関連付けられる:

主張: 半正定値行列 T}} に対し、1=//T// = //A*A// = //B*B//}} ならばユニタリ行列 U}} が存在して 1=//A// = //UB//}} と書ける。実際、主平方根を //B// coloneqq}} //T//sup と書けば、T}} が正定値のとき B}} は可逆で、1=//U// = //AB//supがユニタリであることは からわかる。T}} が正定値でない半正定値行列のときは逆行列の代わりにムーア・ペンローズ擬逆行列 //B//sup が取れて、作用素 1=//B//+//A//}} は部分等長だから、T}} の核の上で自明となるように拡張して U}} が得られる。

3.2.1. 応用
---------

平方根およびそのユニタリ自由度は線型代数学および函数解析学の全般に応用を持つ。

3.2.2. 極分解
----------

可逆行列 A}} に対して、ユニタリ行列 U}} および正定値行列 P}} が一意に存在して 1=//A// = //UP//}} と書ける。これを A}} の極分解と呼ぶ。この正定値行列 P}} は正定値行列 A*A}} の主平方根であり、U}} は 1=//U// = //AP//−1}} で求まる。

A}} が可逆でないときでも、適当な方法で P}} が定まれば極分解が定義される。極分解におけるユニタリ作用素 U}} は一意ではないが、以下のようにして「自然な」ユニタリ行列は求められる: //AP//+}} は A}} の値域からそれ自身への作用素であり、これは A*}} の核上自明に延長してユニタリ作用素 U}} にできるから、この U}} を極分解に用いればよい。

3.3. 一般化
--------

有限次元数空間上で行列を考える代わりに、任意のヒルベルト空間上の有界作用素に対して、その平方根を考えることができる。とくに有界半正定値作用素に対して、半正定値な平方根としての主平方根は一意に決まる。あるいは非エルミート平方根に関しても同様に考えることができる。無限次元の場合には、平方根がユニタリ作用素を施す違いを除いて決まるという事実は、作用素が閉値域ならば正しい。非有界作用素に対しては、閉かつ稠密に定義された二つの平方根 A, B}} に対し部分等方な U}} で 1=//A = UB//}} とできることなどは言える。

4. 関連項目
-------

行列値関数
行列の対数
実二次正方行列#行列変数の函数

5. 脚注
-----

5.1. 注釈
-------

5.2. 出典
-------

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

title=A Course in Functional Analysis, Chapter IV, Reisz functional calculus