■因数分解の算法（その２）
（その１）では，
ａ^2＋ｂ^2＋２ａｂ＝（ａ＋ｂ）^2
ａ^3＋ｂ^3＋ｃ^3−３ａｂｃ＝（ａ＋ｂ＋ｃ）（ａ^2＋ｂ^2＋ｃ^2−ａｂ−ｂｃ−ｃａ）
は因数分解可能，しかるに，
ａ^4＋ｂ^4＋ｃ^4＋ｄ^4＋４ａｂｃｄ
は不可能という問題を取り上げました．
ところで，整数係数・有理数係数という制限を取り払うと，因数分解可能になる式があります．たとえば，
ｘ^2−３ｘｙ＋ｙ^2＋ｘ＋ｙ−１
は，無理数である黄金比：τ＝（√５＋１）／２，τ^(-1)＝（√５−１）／２を用いると
（τｘ−τ^(-1)ｙ−１）（τ^(-1)ｘ−τ^(-1)ｙ＋１）
と因数分解できます．
さらに実数という制限も外してみることにしましょう．
ｘ^2＋ｙ^2
を実数の世界で因数分解することはできませんが，複素数を使うと
ｘ^2＋ｙ^2＝（ｘ＋ｙｉ）（ｘ−ｙｉ）
のように因数分解できます．また，
ａ^2＋ｂ^2＋ｃ^2−ａｂ−ｂｃ−ｃａ
は，ω（１の虚立方根）＝（ｉ√３−１）／２を用いると
（ａ＋ｂω＋ｃω^2）（ａ＋ｂω^2＋ｃω）
と書くことができます．
ｘ^2＋ｙ^2＋ｚ^2
は複素数の助けを借りても因数分解できませんが，４元数を用いると，
ｘ^2＋ｙ^2＋ｚ^2＋ｗ^2＝（ｘ＋ｙｉ＋ｚｊ＋ｗｋ）（ｘ−ｙｉ−ｚｊ−ｗｋ）
ですから，
ｘ^2＋ｙ^2＋ｚ^2＝−（ｘｉ＋ｙｊ＋ｚｋ）^2
のように因数分解できます．
今回のコラムのネタは複素数を使っても因数分解できなかった
ｘ^2＋ｙ^2＋ｚ^2
ｘ^2＋ｙ^2＋ｚ^2＋ｗ^2
をトリックを用いて分解する方法を紹介したいと思います．
［参］佐武一郎「線型代数学」裳華房
［参］村上雅人「なるほど線形代数」海鳴社
からの受け売りですが，なるほど・・・のタイトルに違わぬ面白い本でした．
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
【１】複素数と行列
前回のコラム「４次元の特殊性」で取り扱った平面の回転行列は
［ｃｏｓθ，−ｓｉｎθ］
［ｓｉｎθ， ｃｏｓθ］
の形に書けました．
ある複素数に虚数単位ｉをかけると
ｚｉ＝（ｘ＋ｙｉ）ｉ＝−ｙ＋ｘｉ
となり，この操作は９０°回転に対応することがわかります．そこで，回転行列にθ＝π／２を代入すると
Ｊ＝［０，−１］
［１， ０］
となります．
複素数平面でｉが果たす役割と行列Ｊが果たす役割は等しいのですが，実際にこの行列を２乗すると
Ｊ^2＝［１，０］＝−Ｅ
［０，１］
となって，虚数のもっている性質を備えていることがわかります．
このことを踏まえると，複素数に対応した行列を導入することができます．
Ｚ＝ｘＥ＋ｙＪ＝［ｘ，−ｙ］
［ｙ， ｘ］
ここで，
Ｅ＝［１，０］ Ｊ＝［０，−１］
［０，１］ ［１， ０］
Ｚ’＝ｘＥ−ｙＪ＝［ ｘ，ｙ］
［−ｙ，ｘ］
Ｚ・Ｚ’＝［ｘ^2＋ｙ^2，０］＝（ｘ^2＋ｙ^2）Ｅ
［０，ｘ^2＋ｙ^2］
ですから，（ｘ^2＋ｙ^2）Ｅという行列が（虚数単位ｉを陽に用いることなしに）行列Ｚと行列Ｚ’の積に分解できたことになります．
また，オイラーの公式
ｅｘｐ（ｉθ）＝ｃｏｓθ＋ｉｓｉｎθ
を行列で表現すると
ｅｘｐ（Ｊθ）＝（ｃｏｓθ）Ｅ＋（ｓｉｎθ）Ｊ
＝［ｃｏｓθ，−ｓｉｎθ］
［ｓｉｎθ， ｃｏｓθ］
となって，確かに回転行列になっていることがわかります．
［問］２次の正方行列でＸ^2＝−Ｅを満たすものは，複素数の範囲では
［ｉ，０］
［０，ｉ］
などがある．しかし，実数の範囲でも解は無数にある．その解をすべて求めよ．
［補］Ｅ，Ｊはそれぞれ対称行列（Ｂ’＝Ｂ），交代行列（Ｂ’＝−Ｂ）になっています．任意のｎ次正方行列Ａに対して，
Ｂ＝（Ａ＋Ａ’）／２，Ｃ＝（Ａ−Ａ’）／２
とおけば，Ｂは対称行列，Ｃは交代行列であって，
Ａ＝Ｂ＋Ｃ
のような分解は一意で与えられます．
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
【２】パウリ行列
前節で述べたことにより，行列を使うと因数分解ができるようになる可能性が開けてきます．円に相当するｘ^2＋ｙ^2ができたわけですから，球に相当するｘ^2＋ｙ^2＋ｚ^2も分解してみたい・・・．そして実際に，
ｘ^2＋ｙ^2＋ｚ^2
の因数分解を可能にするのが「パウリ行列」です．
パウリ行列は
σx＝［０，１］ σy＝［０，−ｉ］ σz＝［１， ０］
［１，０］ ［ｉ， ０］ ［０，−１］
の３組の２×２行列で与えられるのですが，いずれも２乗すると単位行列になります．
σx^2＝Ｅ，σy^2＝Ｅ，σz^2＝Ｅ
また，行列のかけ算は非可換なのですが，パウリ行列では，
σxσy＝ｉσz，σyσx＝−ｉσz
のように符号が逆となり，
σxσy＋σyσx＝Ｏ（ゼロ行列）
σxσy−σyσx＝２ｉσz
のような関係が成立します．
ここで，行列
ｘσx＋ｙσy＋ｚσz＝［ ｚ，ｘ−ｙｉ］
［ｘ＋ｙｉ， −ｚ］
を考え，この行列を２乗してみます．すると，
（ｘσx＋ｙσy＋ｚσz）^2＝［ｘ^2＋ｙ^2＋ｚ^2，０］
［０，ｘ^2＋ｙ^2＋ｚ^2］
＝（ｘ^2＋ｙ^2＋ｚ^2）Ｅ
結局，（ｘ^2＋ｙ^2＋ｚ^2）Ｅという行列は，（ｘσx＋ｙσy＋ｚσz）^2に分解できたことになります．４元数を使わないとできなかった因数分解が，行列を利用すると分解できるトリックは，行列の成分として虚数単位を含んでいるうえに，行列自体にも虚数の働きがあり，普通の数にはない機能を２重に使っているからと考えられます．
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
【３】ディラック行列
パウリの行列において
ｘ^2＋ｙ^2＋ｚ^2＝ｒ^2
は３次元空間での球に相当するわけですが，歴史的にはパウリが基本粒子のスピンを数学的に表現するために考案したものです．この考えがヒントになって，ディラックが４次元時空における時間の項を加えた
ｘ^2＋ｙ^2＋ｚ^2＋ｔ^2
を因数分解するために「ディラック行列」を使いました．
αx＝［０，０，０，１］ αy＝［０， ０，０，−ｉ］
［０，０，１，０］ ［０， ０，ｉ， ０］
［０，１，０，０］ ［０，−ｉ，０， ０］
［１，０，０，０］ ［ｉ， ０，０， ０］
αz＝［０， ０，１， ０］ β＝［１，０， ０， ０］
［０， ０，０，−１］ ［０，１， ０， ０］
［１， ０，０， ０］ ［０，０，−１， ０］
［０，−１，０， ０］ ［０，０， ０，−１］
これらの４組の２×２行列をディラック行列と呼ぶのですが，前項と同様に，
ｘαx＋ｙαy＋ｚαz＋ｔβ
＝［ ｔ， ０， ｚ， ｘ−ｙｉ］
［ ０， ｔ， ｘ＋ｙｉ，−ｚ ］
［ ｚ， ｘ−ｙｉ，−ｔ， ０ ］
［ｘ＋ｙｉ，−ｚ ，０， −ｔ ］
（ｘαx＋ｙαy＋ｚαz＋ｔβ）^2
＝［ｘ^2＋ｙ^2＋ｚ^2＋ｔ^2，０，０，０］
［０，ｘ^2＋ｙ^2＋ｚ^2＋ｔ^2，０，０］
［０，０，ｘ^2＋ｙ^2＋ｚ^2＋ｔ^2，０］
［０，０，０，ｘ^2＋ｙ^2＋ｚ^2＋ｔ^2］
＝（ｘ^2＋ｙ^2＋ｚ^2＋ｔ^2）Ｅ
という関係が確かめられます．これで，（ｘ^2＋ｙ^2＋ｚ^2＋ｔ^2）Ｅという行列は，（ｘαx＋ｙαy＋ｚαz＋ｔβ）^2に分解できることがわかりました．
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
ちなみに，ディラック行列をパウリ行列で表現すると，
αx＝［０，σx］ αy＝［０，σy］ αz＝［０，σz］
［σx，０］ ［σy，０］ ［σz，０］
β＝［Ｅ， ０］
［０，−Ｅ］
ｘαx＋ｙαy＋ｚαz＋ｔβ＝［ｔＥ， ｘσx＋ｙσy＋ｚσz］
［ｘσx＋ｙσy＋ｚσz，−ｔＥ］
と表すことができます．
このことから
αx^2＝Ｅ，αy^2＝Ｅ，αz^2＝Ｅ，β^2＝Ｅ
αxαy＋αyαx＝Ｏ（ゼロ行列）
αxαy−αyαx＝２ｉ［σz，０］
［０，σz］
αxβ＋βαx＝Ｏ（ゼロ行列）
αxβ−βαx＝２［０，−σx］
［σx， ０］
と計算されます．ディラック行列がパウリ行列に一工夫加えた様子を窺い知ることができるでしょう．
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
【４】４元数と行列
積の交換法則が成り立たない代数として「行列」があります．
Ｅ＝［１，０］ Ｊ＝［０，−１］ Ｊ^2＝−Ｅ
［０，１］ ［１， ０］
とおけば，
Ａ＝［ａ1，−ａ2］
［ａ2， ａ1］
は
Ａ＝ａ1Ｅ＋ａ2Ｊ
と表されます．
Ａ＝ａ1Ｅ＋ａ2Ｊ，Ｂ＝ｂ1Ｅ＋ｂ2Ｊ
の形の行列全体は加法および乗法に関して閉じています．
Ａ＋Ｂ＝（ａ1＋ｂ1）Ｅ＋（ａ2＋ｂ2）Ｊ
ＡＢ＝（ａ1ｂ1−ａ2ｂ2）Ｅ＋（ａ1ｂ2＋ａ2ｂ1）Ｊ
乗法の可換性は成立しません．すなわち，【１】節では行列による表現を利用して複素数を導入したわけですが，類似の方法で４元数を行列の中に実現させる方法もあります．
Ｅ＝［１，０，０，０］
［０，１，０，０］
［０，０，１，０］
［０，０，０，１］
ｉ＝［０，−１，０， ０］ ｊ＝［０， ０，−１，０］
［１， ０，０， ０］ ［０， ０， ０，１］
［０， ０，０，−１］ ［１， ０， ０，０］
［０， ０，１， ０］ ［０，−１， ０，０］
ｋ＝［０，０， ０，−１］ Ａ＝［ａ1，−ａ2，−ａ3，−ａ4］
［０，０，−１， ０］ ［ａ2， ａ1，−ａ4， ａ3］
［０，１， ０， ０］ ［ａ3， ａ4， ａ1，−ａ2］
［１，０， ０， ０］ ［ａ4，−ａ3， ａ2， ａ1］
とおけば
Ａ＝ａ1Ｅ＋ａ2ｉ＋ａ3ｊ＋ａ4ｋ
と書くことができます．
この体系では，４元数同様，
ｉ^2＝−Ｅ，ｊ^2＝−Ｅ，ｋ^2＝−Ｅ，
ｉｊ＝ｋ，ｊｋ＝ｉ，ｋｉ＝ｊ，
ｊｉ＝−ｋ，ｋｊ＝−ｉ，ｉｋ＝−ｊ
なる性質をもっていて，加法および乗法に関して閉じています．また，乗法の可換性は成立しません．
この体系を用いると，
（ｘ^2＋ｙ^2＋ｚ^2）Ｅ＝−（ｘｉ＋ｙｊ＋ｚｋ）^2
（ｘ^2＋ｙ^2＋ｚ^2＋ｗ^2）Ｅ＝（ｘ＋ｙｉ＋ｚｊ＋ｗｋ）（ｘ−ｙｉ−ｚｊ−ｗｋ）
のように，虚数単位ｉを陽に用いることなしに２つの行列の積に分解できますが，それでは４元数そのままであって，虚数単位ｉを使ったパウリ行列やディラック行列よりも面白味に欠けるかもしれません．
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
行列ｊを変更して
ｉ＝［０，−１，０， ０］ ｊ＝［０，０，−１， ０］
［１， ０，０， ０］ ［０，０， ０，−１］
［０， ０，０，−１］ ［１，０， ０， ０］
［０， ０，１， ０］ ［０，１， ０， ０］
ｋ＝［０，０， ０，−１］ Ａ＝［ａ1，−ａ2，−ａ3，−ａ4］
［０，０，−１， ０］ ［ａ2， ａ1，−ａ4，−ａ3］
［０，１， ０， ０］ ［ａ3， ａ4， ａ1，−ａ2］
［１，０， ０， ０］ ［ａ4， ａ3， ａ2， ａ1］
とおくと，Ａの上三角部分，下三角部分が歪対称（ｇij＝−ｇji）になるので，形がスッキリすると思われます．
ところが，実際にやってみると
ｊ^2＝−１
は成り立つものの
ｉｊ＝ｋ，ｊｋ＝ｉ，ｋｉ＝ｊ，
ｊｉ＝−ｋ，ｋｊ＝−ｉ，ｉｋ＝−ｊ
が成立しません．なかなかうまくはいかないものです．
なお，８元数：
ｉ^2＝ｊ^2＝ｋ^2＝ｌ^2＝ｍ^2＝ｎ^2＝ｏ^2＝−１，
ｉ＝ｊｋ＝ｌｍ＝ｏｎ＝−ｋｊ＝−ｍｌ＝−ｎｏ，
ｊ＝ｋｉ＝ｌｎ＝ｍｏ＝−ｉｋ＝−ｎｌ＝−ｏｍ，
ｋ＝ｉｊ＝ｌｏ＝ｎｍ＝−ｊｉ＝−ｏｌ＝−ｍｎ，
ｌ＝ｍｉ＝ｎｊ＝ｏｋ＝−ｉｍ＝−ｊｎ＝−ｋｏ，
ｍ＝ｉｌ＝ｏｊ＝ｋｎ＝−ｌｉ＝−ｊｏ＝−ｎｋ，
ｎ＝ｊｌ＝ｉｏ＝ｍｋ＝−ｌｊ＝−ｏｉ＝−ｋｍ，
ｏ＝ｎｉ＝ｊｍ＝ｋｌ＝−ｉｎ＝−ｍｊ＝−ｌｋ
では，乗法の結合法則も破れていて（ａ（ｂｃ）≠（ａｂ）ｃ），積の交換法則も結合法則も成り立ちませんが，それでも分配法則は成り立っています．行列は結合法則を満たすので，８元数は行列の一部とはみなせないのです．なお，結合法則が成り立たない数の体系（非結合的な体）としては，８元数，リー代数，ジョルダン代数の３つが代表的です．
＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝＝
