スポンサードリンク

こんにちは、ももやまです。
今回は、対角化の中でも、直交行列を用いた対角化についてのまとめです。

前回の記事(第16羽 行列の対角化)はこちら↓

www.momoyama-usagi.com

スポンサードリンク

1.直交行列で対角化できる条件

直交行列の対角化を行うためには通常の対角化に比べ、対角化を行う行列 A に対し、より厳しい条件が課せられます。

直交行列の対角化

行列 A実対称行列(A=tAを満たすとき、直交行列*1 P を用いて、P1AP と対角化をすることができる。

もちろん直交行列で対角化できるような行列(つまり実対称行列)は普通に対角化を行うこともできます

実対称行列についてはこちらの記事を、直交行列についてはこちらの記事をご覧ください。

スポンサードリンク

2.直交行列の対角化(重解なしの場合)

では、まずは固有値の重解がない場合の例題を解きながら直交行列の対角化の流れを理解しましょう。

例題1

行列A=(5222)が実対称行列であることを確かめ、直交行列 P を用いて対角化しなさい。

解説1

行列 A と転置行列 tA は等しくなりますね。
なので、直交行列を用いた対角化を行うことができます。

固有値を t とすると、固有方程式は、|AtE|=|5t222t|=(t5)(t2)4=t27t+6=(t1)(t6)=0より固有値は1, 6となる。

(固有値を出すまでは通常の対角化と同じ流れです。)

つぎに、固有ベクトルを求めていきます。

前回は正則行列 P での対角化なので固有ベクトルの大きさは考えなくてもOKでしたが、今回は直交行列 P を用いた対角化のため、それぞれの固有ベクトルを正規直交基底にする必要があります。そのため、それぞれの固有ベクトルの大きさは1でなければなりません。

まず、固有値が1のときの固有ベクトルを求めましょう。

固有値が1のときの固有ベクトルは、(A1E)= (4221) (2100)となる。2x+y=0を解くと、任意定数 k を用いて(xy)=k(12)と表せます。

固有ベクトル p1 は大きさを1に正規化したベクトルなので、p1=15(12)となる*2

つぎに固有値が6のときの固有ベクトルを求めましょう。

固有値が6のときの固有ベクトルは、(A6E)= (1224) (1200)となる。x2y=0を解くと、任意定数 k を用いて(xy)=k(21)と表せます。

大きさを1に正規化した固有ベクトル p2 は、p2=15(21)となる。

ここで p1, p2 が正規直交基底になることを確認しましょう。

それぞれのベクトルは正規化済なので大きさが1である確認は不要です。

あとは p1p2 が直交することを確認するだけです。p1p2=(12)(21)=0となり、p1, p2 は直交しますね。

よって、p1, p2 は正規直交基底となるので、直交行列P=(p1,p2)=15(1221)を用いて、P1AP=(1006)と対角化することができます。*3

さて、先ほど固有値の異なる固有ベクトル p1, p2 は互いに直交することを確認しましたね。しかしこれは偶然なのでしょうか?

実は、実対称行列の場合、異なる固有値に対する固有ベクトルは必ず直交します。簡単に証明してみましょう。

[証明] (余裕ない人は飛ばしてOKです)

実対称行列 A の固有値 α に対する固有ベクトルを p1、固有値 β に対する固有ベクトルを p2とします(ただし αβ とする)。すると、Ap1=αp1   Ap2=βp2が成立しますね。

ここで内積には、(Ap1)p2=p1(tAp2)=k(p1p2)=kp1p2が成立します。

さらに実対象行列なので、A=tA ですね。つまり、(Ap1)p2=p1(Ap2)となります。

さらに、(Ap1p2)=(αp1)p2=α(p1p2) p1(Ap2)=p1(βp2)=β(p1p2)

なので、α(p1p2)=(Ap1p2)=p1(Ap2)=β(p1p2)となりますね。つまり、(αβ)(p1p2)=0となります。

ここで、固有値 α, β は互いにことなるので、両辺を αβ で割り、p1p2=0となるので異なる固有値に属する固有ベクトル p1, p2 が直交することが示せましたね*4

[証明終わり]

異なる固有値に対する固有ベクトルは必ず直交する実対称行列と固有ベクトルの直交性

スポンサードリンク

3.直交行列の対角化(重解ありの場合)

つぎに固有値が重解となる場合の直交行列を用いた対角化について説明していきたいとおもいます。

固有値が重解の場合、同じ固有値に対する固有ベクトルは直交していないため、グラムシュミットの直交化法を用いて同じ固有値同士の固有ベクトルを直交化させる必要があります。

グラムシュミットの直交化を忘れてしまった(orよくわかんない)という人はこちらの記事で確認をお願いします↓

www.momoyama-usagi.com

例題2

行列A=(322232223)が実対称行列であることを確かめ、直交行列 P を用いて対角化しなさい。

解説2

まずは、いつもどおり固有値を求めていきます。

固有値を t とすると、固有方程式は、|AtE|=|3t+2+22+3t+22+2+3t23t2223t|=|7t7t7t23t2223t|=(7t)|111223t22222223t2|=(7t)|11101t0001t|=(1t)2(7t)=0となる*5。ので、固有値は1(2重解)と7になる。

つぎに固有ベクトルを求めていく。

(1) 固有値が1(2重解)のとき(2重解)

(A1E)= (222222222) (111000000)となる。x+y+z=0を解くと任意定数 s, t を用いて(xyz)=s(110)+t(101)となる。

ここで、ベクトルa1=(110)   a2=(101)は固有ベクトルだが、互いに直交していない。なので、グラムシュミットを用いることで a1, a2 を正規直交化し、正規直交するような固有ベクトルを求める。

a1, a2 を正規直交化したベクトルを p1, p2 とする。p1=1|a1|a1=12(110)となる。

a2u1=12=12なので、b2=a2(a2u1)u1=(101)1212(110)=12(202)12(110)=12(112)となる。

よって、p2=1|b2|p2=16(112)となる。

よって、正規直交化された固有ベクトルは、p1=12(110),    p2=16(112)となる。

(2) 固有値が7のとき

重解ではないのでただ正規化するだけでOK。

(A7E)= (4+4282+422422+4224) (066066224) (06606+6662224+2) (011000202) (101011000)となる。{xz=0yz=0を解くと任意定数 k を用いて(xyz)=k(111)となる。

大きさを1に正規化した固有ベクトル p3 は、p3=13(111)となる。

よって、p1, p2, p3 は正規直交基底となるので、直交行列P=(p1,p2,p3)=16(312312022)を用いて、P1AP=(100010007)と対角化することができます。

[検算] AP=PD=16(317231720272)

4.練習問題

では、2問ほど練習をしましょう。

1問は2次正方行列、もう1問は3次正方行列です。

練習1

行列A=(3113)を直交行列 P を用いて対角化しなさい。

練習2

行列A=(321202123)を直交行列 P を用いて対角化しなさい。

5.練習問題の答え

対角化したあとは必ず AP=PD を計算して正しく対角化できているかを検算しましょう!*6P がルートになるので若干検算がめんどくさくなりますが……)

解説1

固有値を t とすると、固有方程式は、|AtE|=|3t113t|=(t3)21=t26t+8=(t2)(t4)=0より固有値は2, 4となる。

つぎに、p1, p2 が正規直交基底となるような固有ベクトルをもとめる。

(1) 固有値が2のとき(A2E)= (1111) (1100)となる。x+y=0を解くと、任意定数 k を用いて(xy)=k(11)と表せる。

固有ベクトル p1 は大きさを1に正規化したベクトルなので、p1=12(11)となる。

(2) 固有値が4のとき(A4E)= (1111) (1100)となる。xy=0を解くと、任意定数 k を用いて(xy)=k(11)と表せる。

固有ベクトル p1 は大きさを1に正規化したベクトルなので、p1=12(11)となる。

よって、p1, p2 は正規直交基底となるので、直交行列P=(p1,p2)=12(1111)を用いて、P1AP=(2004)と対角化することができます。

[検算] AP=PD=12(2424)となるので対角化が正しいことがわかる。

解説2

固有値を t とすると、固有方程式は、|AtE|=|3t212t2123t|=|3t+1221t32t2123t|=|3t+1221t32t2123t|=|4t04+t2t2123t|=(4t)|10122t2+21+123t1|=(4t)|1010t4022t|=(4t)|t422t|となる。

ここで、|t422t|=t(t2)8=t22t8=(t+2)(t4)=0となるので、|AtE|=(4t)(t4)(t+2)=(t4)2(t+2)=0を満たす t が固有値となり、固有値は4(2重解)と-2となる。

つぎに、p1, p2, p3 が正規直交基底となるような固有ベクトルをもとめる。

(1) 固有値が4(2重解)のとき

(A1E)= (121242121) (121000000)となる。x2y+z=0を解くと任意定数 s, t を用いて(xyz)=s(101)+t(210)となる。

ここで、ベクトルa1=(101)   a2=(210)をグラムシュミットを用いることで a1, a2 を正規直交化する。

a1, a2 を正規直交化した固有ベクトルを p1, p2 とする。p1=1|a1|a1=12(101)となる。

a2u1=12=22なので、b2=a2(a2u1)u1=(210)2212(101)=(210)(101)=(111)となる。

よって、p2=1|b2|,   p2=13(111)となる。

よって、正規直交化された固有ベクトルは、p1=12(101),   p2=13(111)となる。

(2) 固有値が-2のとき

重解ではないのでただ正規化するだけでOK。

(A7E)= (521222125) (552+101+25222+42+10125) (01224061212254) (012000101) (101012000)となる。{xz=0y+2z=0を解くと任意定数 k を用いて(xyz)=k(121)となる。

大きさを1に正規化した固有ベクトル p3 は、p3=16(121)となる。

よって、p1, p2, p3 は正規直交基底となるので、直交行列P=(p1,p2,p3)=16(321022321)を用いて、P1AP=(400040002)と対角化することができます。

[検算] AP=PD=16(43422042443422)

6.さいごに

今回は、直交行列を用いた対角化について説明しました。

直交行列の対角化は、2次形式を標準形にする際に必ず使うので覚えておきましょう。

また、直交行列 P で対角化を行うと、tP=P1 が成立するため、行列の n 乗などを計算する際に逆行列を簡単に求められるなどのメリットがあります。

次回は行列の対角化の応用として行列の n 乗を求める方法についてまとめたいと思います。

*1:直交行列は、tPP=E   tP=P1を満たすような行列のこと。

*2:正規化をする際には任意定数 k を考えずにb1=(12)を正規化することを考えればよい。この場合、b1|b1|(12)=15(12)と正規化できる。(任意定数は定数倍を表しているだけなので正規化するときに消えるため。他の問題においても任意定数部分を無視して正規化するとよい。)

*3: [検算] AP=PD=15(11226)となるので対角化が正しいことがわかる。

*4:逆にいうと同じ固有値同士の場合、αβ=0 となり、p1p2=0 が成立するとは限らないことがわかる。

*5:上三角行列、もしくは下三角行列の行列式は対角成分の積となるのを利用した。

*6:もちろん Ap=tp で検算してもOK(行列と固有ベクトルの積は固有値倍されたベクトルに等しい)。

関連広告・スポンサードリンク

おすすめの記事