13

2次の正方行列のn乗根ってn個なの?

1283
0

この記事 is 何?

一般に複素数範囲の数であれば n 乗根は重解含めて n 通り存在するわけですが、行列の場合にもあてはまるのでしょうか? そこで、最も単純な 2 次の単位行列の n 乗根について考えてみました(*´∀`*)

2次の単位行列の自明なn乗根

2次の単位行列といえば E=(1001) ですよね。
E=(1001) は、この空間が (1,0)(0,1) を基底の元として構成されてるよってことを表した数で、この空間からみた (1,0)(0,1) を新たな基底の元とする異空間を想定すると、全く同じクローン空間ができます。

(1001)(1001)=(1001)

このことは、 E を整数乗しても E と同じだよってことを意味しているわけですが、では整数乗根はどうなるでしょうか?

En乗根として自明的に思いつくのは、基底の姿勢を回転させてその角度で対応させるもの。

例えば基底の各元を空間内で原点を中心に 2πn rad 回転させた (cos2πn,sin2πn)(sin2πn,cos2πn) を新たな基底として想定すれば、同じ操作を n 回繰り返すことで1周して最初と同じ姿勢に戻りますから、これは n 乗根といえますね。

(例)(cos2πnsin2πnsin2πncos2πn)n=(1001)

少し反則的ではありますが、回転方向は同一空間上に限る必要はありません。各元に対して垂直な次元方向へ回転させても、(想定している2次元の空間からは飛び出してしまいますが)n 乗したら 1周して元の空間の単位行列へ戻ってきます。

(例)(e2πni00e2πni)n=(1001)

このような異次元への飛び出しを認めるなら回転面は無限に存在しますので n 乗根は無限に存在するといえるでしょう。ただ先程も書きましたようにこれは少し反則的です。といいますのも、この例では実質的に2次元×2次元=4次元を対象空間としていますので、通常の数(という表現もヘンですが)でいうところの「複素数範囲に n 個」ではなく「四元数範囲に無限個」と対比させないと割に合いません(←?)

そんなわけで、今回は実数行列の範囲に絞ったときの n 乗根について考察していきたいと思います(*´ω`*)/

2次の単位行列の非自明なn乗根 その1

行列の累乗ときたら、まずは対角化でしょう。早速、En 乗根の行列 R=(x0x1y0y1)x0,y0,x1,y1R)が対角化可能であるようなケースから考えてみます。行列に関する詳細解説は本題ではないため、導出の流れをざっくり追いかける感じで証明などは省略させてくださいませ。

R の固有値を λ とすると λ の総和は trR=x0+y1、総積は detR=x0y1x1y0 に等しいため2次方程式の解と係数の関係より λ2trAλ+detA=0 を解いて

λ=trA±(trA)24detA2=x0+y12±(x0+y12)2(x0y1+x1y0)=x0+y12±(x0y12)2+x1y0

固有ベクトルを (vx,vy) と表すと (x0λx1y0y1λ)(vxvy)=(00) であることから

vx:vy=x1:(x0λ)=(y1λ):y0 より導出される固有ベクトルの定数倍からなる行列を P、対角行列を D として
R=PDP1Rn=(PDP1)n=PDnP1=EP1RnP=P1(PDnP1)P=Dn=E

これを見ると、Dn の対角成分はいずれも 1、すなわち D の対角成分である R2つの固有値 λ0λ1 はいずれも 1n乗根であることが分かります。であれば、固有値の各項との関係性より mZnZ+ を用いて

λ=x0+y12±(x0y12)2+x1y0=cos2mπn±isin2mπn

と表わせますので、各項で連立方程式を立てて解いてましょう。

{x0+y12=cos2mπny1=2cos2mπnx0±(x0y12)2+x1y0=isin2mπn(x0y12)2+x1y0=sin22mπn(x0cos2mπn)2+x1y0=sin22mπn=cos22mπn1x022(cos2mπn)x0+x1y0+1=0x0=cos2mπn±(cos2mπn)2(x1y0+1)
よって、
R=(cos2mπn±(cos2mπn)2(x1y0+1)x1y0cos2mπn(cos2mπn)2(x1y0+1))

と表せることが判明しました。y0x1 が遊んでますが、これは任意の数でおっけーってことで!

さて、2次の正方行列における対角化可能な必要十分条件は2つの固有値が異なる値を取ることです。
今回の場合、±sin22mπn=0 のとき重解 cos22mπn を取ることが自明ですから、n=1n=2 のときのみですね。

n0n1n2mZnZ+x1,y0R において
(cos2mπn±(cos2mπn)2(x1y0+1)x1y0cos2mπn(cos2mπn)2(x1y0+1))n=(1001)

また、対角化不可能というのは正則な P を作れなくて P1 を求められないってことなのですが、R が最初から単位行列の定数倍の場合は P は任意となるため、結果として P を一意に求められないというケースを考慮する必要があります。上の例でいえば x1=y0=0、かつ、cos2mπn=±1 (n=1,2) の場合。代入してみると (1001)1(1001)2 となり、いずれも単位行列になることが分かりました。こちらも追記しておきましょう。

mZn=1n=2 において
(cos2mπn00cos2mπn)n=(1001)\
(この形なら、nZ+ でいけますね。)

ところで余談ですが、y0=1x1=1 とすると、Rn=(0112cos2mπn)n=(1001) を得られます。実はコレ、+1e2mπni を基底の元とする斜交座標系の基本的な回転を表しており、この数式から非常に興味深い考察を得ることができます。よろしければ ガラパゴ累乗定理 ガラパゴ三角関数 を参照してみてくださいませ(*´ω`*)

2次の単位行列の非自明なn乗根 その2

続いて、対角化不可能かつ R が対角行列ではないようなケースを考察してみます。

R=(x0x1y0y1) の固有値が重解、つまり (x0y12)2+x1y0=0 で、かつ、R は単位行列のn乗根かつ実数行列の範囲ということから detR=x0y1x1y0=1 を踏まえると

(x0y12)2+x0y11=(x0+y12)21=0x0+y1=±2trRy1=x0±2x0y1=x0(x0±2)=x02±2x0x1y0=(x022x0+1)=(x01)2

0 ではない任意定数 k を用いれば
y0=k(x01)x1=x01k
と表現することができますね。
よって、この条件を満たす行列は R=(x0x01kk(x01)x0±2) となります。
もちろん、対角行列ではないという条件付きなので x0±1 なことはお忘れなく。

固有値を算出してみると、λ2(±2)λ+1=(λ1)2=0 より λ=±1 と期待通りに重解となっていますね。

固有ベクトル (vx,vy)
{(x01)vx+(x01k)vy=0k(x01)vx+(x0±1)vy=0 より
λ±1 のいずれであっても (1,k) となります。

さて、固有値が重解のため対角化はできませんが、対角化がダメならジョルダン標準形を使えばいいじゃない? というわけでサクっと試してみましょう。

(λ10λ)=(±110±1)

(x0x01kk(x01)x0±2)(1pxkpy)=(1pxkpy)(±110±1)(±1pxx0(x01)pyk±kk(x01)px(x02)py)=(±11±px±kk±py){pxx0(x01)pyk=1±pXk(x01)px(x02)py=k±py

x0±1 という条件下なので {px=1x01py=0

(x0x01kk(x01)x0±2)=(11x01k0)(±110±1)(11x01k0)1=(11x01k0)(±110±1)(01kx01x01k)(x0x01kk(x01)x0±2)n=(11x01k0)(±110±1)n(01kx01x01k)

x0+y1=1 のとき
(x0x01kk(x01)x0+2)n=(11x01k0)(1n01)(01kx01x01k)=(1n(x01)n(x01)kkn(x01)n(x01)+1)

x0+y1=1 のとき
(x0x0+1kk(x0+1)x02)n=(11x01k0)((1)n(1)n+1n0(1)n)(01kx01x01k)=±(1n(x1)n(x0+1)kkn(x0+1)n(x0+1)+1)

x0±1 という縛りがあるため、これらが単位行列になるには少なくとも n=0 でないといけませんね。
要するに、自身が対角行列ではなく、かつ、対角行列を作れないときは単位行列の n 乗根にはなり得ないのです。

結論

2 次の正方行列の n乗根は、実数行列の範囲において以下の通り無数に存在します。(*´ω`*)

n0n1n2mZnZ+x1,y0R において
(cos2mπn±(cos2mπn)2(x1y0+1)x1y0cos2mπn(cos2mπn)2(x1y0+1))n=(1001)

mZnZ+ において
(cos2mπn00cos2mπn)n=(1001)

最後に、この記事の検算・検証にご協力下さった @nayuta_ito さんに感謝します。
どこかおかしな点がありましたら、やさしく教えてね!

投稿日:20201116
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。
バッチを贈って投稿者を応援しよう

バッチを贈ると投稿者に現金やAmazonのギフトカードが還元されます。

投稿者

https://mathlog.info/articles/323         数学を愛する会 副会長 CCO / ガラパゴ数学 開拓者 / 猫舌・甘党・薄味派

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中
  1. 2次の単位行列の自明なn乗根
  2. 2次の単位行列の非自明なn乗根 その1
  3. 2次の単位行列の非自明なn乗根 その2
  4. 結論