行列の平方根

数学のおもに線型代数学および函数解析学における行列の平方根（ぎょうれつのへいほうこん、）は、数に対する通常の平方根の概念を行列に対して拡張するものである。すなわち、行列 が行列 の平方根であるとは、行列の積に関して が に等しいときに言う。
「実数の平方根は必ずしも実数にならないが、複素数は必ず複素数の範囲で平方根を持つ」ことに対応する事実として、実行列の平方根は（存在しても）必ずしも実行列にならないが、複素行列が平方根を持てばそれは必ず複素行列の範囲で取れる。
平方根を持たない行列も存在する。
また一般に、ひとつの行列が複数の平方根を持ち得る。実際、 単位行列は次のように無数の平方根を持つ。formula_1 
このように行列の平方根は無数に存在しうるが、半正定値行列の範疇で行列の主平方根 ("principal square root") の概念が定義できて「半正定値行列の主平方根はただ一つ」である（これは「非負実数が非負の平方根（主平方根）をただ一つだけ持つ」という事実に対応する）。
 行列が、相異なる二つの非零固有値を持つならば、それは四つの平方根を持つ（より一般に、相異なる 個の非零固有値を持つ 行列は 個の平方根を持つ）。実際に、そのような仮定を満たす行列 は の固有ベクトルを列ベクトルに持つ行列 とそれに対応する固有値を対角成分に持つ対角行列 を用いて と固有値分解できるから、 の平方根は で与えられることがわかる。ただし、 は の任意の平方根で、それは の対角成分の任意の平方根を同じ位置の対角成分として持つ対角行列であり、その選び方は 通りある。同じ理由で、上で述べた「半正定値行列の主平方根がただ一つに定まる」ことも言える—半正定値行列 の全ての非負固有値の主平方根を対角成分に持つ対角行列を とする行列 はただ一つしかない。
適当な冪零行列 を用いて の形に書ける行列の平方根 は、二項級数に対する汎函数計算で求められる。同様に、行列の指数函数 , 対数函数 が既知ならば、 を の（主）平方根とすることができる（収束性に注意せよ）。
定義.
「非負実数が非負の平方根（主平方根）をただ一つだけ持つ」という事実に対応して
命題
主平方根をとる操作は行列全体の成す集合上で連続である。このとき、考えている行列が実行列ならば、その主平方根もまた実行列になる。主平方根に関する性質は、行列に対するの帰結として得られる。あるいは主平方根の存在と一意性はジョルダン標準形を用いて直截に示せる（後述）。
計算法.
明示公式.
 行列の場合は、すべての成分を明示的に計算することによって平方根を求めることはそう難しくない。固有値が退化していない場合の平方根はとして記述できる。
すなわち、formula_2とし、その行列式をformula_3、
特性方程式formula_4の判別式を
formula_5としたとき、
formula_6ならば、
formula_7の平方根は、
formula_8、formula_9、formula_10、formula_11と明示的に表記できる。
平方根となることは、実際に2乗を計算すればformula_12から容易にわかる。
あるいは、2次のケイリー・ハミルトンの定理formula_13からformula_14、formula_15としても良い。
これら以外に平方根が存在しないことについては、formula_16とした場合、formula_6よりformula_7は2つの相異なる固有値formula_19、formula_20と、独立な固有ベクトルformula_21、formula_22を持つが、任意の2次列ベクトルは、formula_23、formula_24の1次結合で表せるので、formula_25、formula_26とすると、
formula_27、
formula_28すなわち、
formula_29であるが、
formula_30のため、解はformula_31、formula_32、formula_33に定まる。これにより任意の2次列ベクトルformula_34がformula_35によりどう変換されるかが定まるが、これはformula_35が定まることを意味する。formula_7が固有値ゼロを持たない場合は解が4組、固有値ゼロを持つ場合（formula_38の場合）は解が2組であるが、これは上記の明示公式で尽くされているので、これら以外には、平方根は存在しない。
formula_39の場合は、複雑になる。
 が 対角行列ならば、 の対角成分の任意の平方根を対応する位置の対角成分に持つ対角行列 を作れば平方根が得られる。 の対角成分が非負の実数ならば、先の対角行列 で各成分の符号を全て正としたものは の主平方根である。
冪等行列の平方根は、自身を平方根に持つ。
対角化の利用.
対角化可能行列 に対し、適当な行列 と対角行列 が存在して と書ける。これは が を張る 個の固有値を持つことと同値である。このとき はその列ベクトルが 個の固有ベクトルであるように選べる。そうして の平方根は の任意の平方根を用いて formula_40 と書ける。実際、formula_41 である。
 が対称行列ならば対角化に用いる行列 は固有ベクトルを適当に選んで直交行列となるようにとれる。この場合、 の逆行列はたんに転置をとるだけであるから、formula_42 と書ける。
ジョルダン分解の利用.
正方行列 formula_43のジョルダン標準形 をformula_44とすると、次が言える。
このため、ジョルダン標準形formula_44の全ての平方根formula_45を知ることができれば、formula_48により、formula_43の全ての平方根formula_51を知ることができる。
formula_63とし、formula_64とすれば、formula_65は、formula_46の平方根のうちの一つである。
逆に、formula_67、ただしformula_68はジョルダン標準形で、formula_69とformula_70は共通の固有値を持たないとすると、formula_46の平方根は、formula_72ただし、formula_73に限られる。
これは、formula_74とすると、
formula_76となるが、formula_77、formula_70の対角成分（固有値）をformula_79と置き、第１列に注目すれば、formula_80だが、formula_69とformula_70は共通の固有値を持たないため、formula_83が言え、順次、第２列、第３列に注目すればformula_84が言え、formula_85が言える。
formula_86からも同様に、formula_87と置き、第k行に注目すればformula_88だが、formula_69とformula_70は共通の固有値を持たないため、formula_91が言え、順次、第k-1行、第k-2行に注目すればformula_92が言え、formula_93が言える。このため、上記が言える。
ジョルダン標準形の平方根には、ジョルダン細胞の平方根であるものと、
のようにジョルダン細胞の平方根ではないもの（同じ固有値のジョルダン細胞が複数あるときに発生する）があるので、注意が必要である。
ジョルダン細胞の平方根.
ジョルダン細胞formula_95とは次正方行列で、formula_96、formula_97、formula_98、formula_99のときformula_100となるものを言う。
formula_101のとき、ジョルダン細胞formula_95の平方根は、下記の行列formula_45およびformula_104である。
formula_109のとき、ジョルダン細胞formula_110は、
例 formula_113は平方根を持たない。
formula_101のとき、ジョルダン細胞formula_95の平方根が２つしかないことは、次から言える。
formula_116となる行列が存在したとし、formula_117の成分を考える。
formula_120を比較すると、formula_121このためformula_122
formula_123を比較すると、formula_124このためformula_125
このため、formula_45は上三角行列で、斜めに同じ値が並ばなければならない。formula_116のformula_128成分を比較することにより、formula_129が言え、以下formula_130成分formula_131を比較することにより、formula_45の全ての成分が順番に１次方程式で定まるため、平方根が２つしかないことが言える。
英語版からの直訳.
対角化可能でない行列の場合にはジョルダン標準形が利用できる。 
すべての固有値が正の実数であるような任意の複素行列が、同じ条件の平方根を持つことを見るには、ジョルダンブロックの場合に証明すれば十分である。そのようなブロックは実数 および冪零行列 を用いて の形に書ける。平方根の二項級数展開 （収束域は ）に対し、形式冪級数としての平方は に等しい。 を に置き換えれば、冪零性により有限個を除く全ての項は零となり、 が固有値 に属するジョルダンブロックの平方根を与える。
一意性を見るには の場合に確認すれば十分である。上で構成した平方根を の形に書けば、 は定数項を持たない の多項式である。固有値が正の実数となる他の任意の平方根 は の形で が冪零かつ と（したがって と）可換となるようにとれる。しかしこのとき であり、また と の可換性により は冪零ゆえ は可逆（逆行列はノイマン級数で与えられる）となるから、したがって .
すべての固有値が正の実数であるような行列 の最小多項式を とするとき、 の一般固有空間へのジョルダン分解は の部分分数分解から導かれる。すなわち、対応する一般固有空間の上への射影は の実係数多項式として与えられ、各固有空間上で は上記の通り の形をしている。固有空間上での平方根の冪級数展開は、 の主平方根が実係数多項式 に対する の形をしていることを示すものである。
現実的な計算法.
「対角化」の方法でも「ジョルダン分解」の方法でも、すべての固有値を算出することが必要となるが、それは行列の特性方程式（あるいは最小方程式）のすべての解を求めることと同じであり、行列の次数が大きくなれば非現実的となる。このため、現実的な平方根の求め方が必要となる。
行列対数関数、行列指数関数による求め方.
実数formula_133の平方根formula_134がformula_135で求まることと同様に、
n次実数値正方行列formula_43の全ての特性根の実数部分が正である場合、
行列対数関数をformula_137と定義し（cは任意の正数、級数が収束すればformula_138に係らず同じ値に収束する）
行列指数関数をformula_139と定義すれば（formula_140はn次実数値正方行列）、
2乗するとformula_43となり、かつ全ての特性根の実数部分が正となる行列formula_142は、
formula_143により計算でき、かつこの行列に一意に定まる。
この方法は、固有値を全て求める必要がないこと（「全ての特性根の実数部分が正」という条件は、特性根を全て求めなくても、十分条件がいくつか知られている）、収束計算が速いこと、対称行列に限らず一般の行列に利用可能であることなど、現実的かつ速い計算方法になっている。
また、行列の平方根に限らず、ｎ乗根も同様に計算することができる。
ニュートン法.
実数の方程式formula_144をニュートン法で解く方法を、行列にそのまま適用して求める方法である。
ｎ次正方行列formula_7に対し、ｎ次正方行列の列formula_146を次の漸化式で定める
formula_147
この列が適当な初期値formula_148について収束すれば、収束値formula_149について、formula_150となる。
このことは、収束すればformula_151が成り立つことから言える。
対称行列（エルミート行列）に限定した議論.
以下では、対称行列（あるいはエルミート行列）に限定した行列の平方根についての性質を示す。
「正定値行列」とは、対称行列（あるいはエルミート行列）で、その全ての固有値が正の実数であるものをいう。「半正定値行列」とは、対称行列（あるいはエルミート行列）で、その全ての固有値がゼロまたは正の実数であるものをいう。
定義.
転置あるいはエルミート共軛を用いれば、より一般に非対称あるいは非エルミートな矩形行列の範疇で「平方根」をとることができる。
 がエルミート（実係数の場合は対称）ならば、 は上で述べたの平方根と一致する。任意の正定値エルミート行列 に対し、それ自身正定値エルミートとなる平方根は一意であり、これを主平方根 (unique square root, principal square root)と呼ぶ。
非対称平方根のユニタリ自由度.
正実数の平方根は、主平方根に を掛けたものですべて与えられた。これに対応するように、正定値エルミート行列の任意の非エルミート平方根は、ユニタリ変換によって関連付けられる:
実際、主平方根を と書けば、 が正定値のとき は可逆で、がユニタリであることは formula_152 からわかる。 が正定値でない半正定値行列のときは逆行列の代わりにムーア・ペンローズ擬逆行列 が取れて、作用素 は部分等長だから、 の核の上で自明となるように拡張して が得られる。
応用.
平方根およびそのユニタリ自由度は線型代数学および函数解析学の全般に応用を持つ。
極分解.
可逆行列 に対して、ユニタリ行列 および正定値行列 が一意に存在して と書ける。これを の極分解と呼ぶ。この正定値行列 は正定値行列 の主平方根であり、 は で求まる。
 が可逆でないときでも、適当な方法で が定まれば（それは一意であり）極分解が定義される。極分解におけるユニタリ作用素 は一意ではないが、以下のようにして「自然な」ユニタリ行列は求められる: は の値域からそれ自身への作用素であり、これは の核上自明に延長してユニタリ作用素 にできるから、この を極分解に用いればよい。