| 9-2 正規行列の直交分解 | ||
| f-denshi.com 最終更新日: | ||
[1] T をベクトル空間 V 上の対角化可能な線形演算子(正方行列)で,その固有値を λ1,λ2,・・・λm それに対応する固有空間を E (λ1),E (λ2),・・・,E (λm) とすれば,ベクトル空間V は,
V = E(λ1) E(λ2) ・・・ E(λm)
と,固有空間の直和で表せることを4章で[#] 述べました。ここでは Tとその随伴演算子 T* とが,
TT*=T*T
を満たすならば,T を用いた V の直和分解は,V の直交分解[#]にもなっていることを示します。つまり,異なる固有空間に属するベクトルどおしはお互いに直交している,すなわち,
V = E (λ1)⊥E (λ2)⊥・・・・・⊥E (λm) [ 直交分解 ]
であることを示します。
[2] きちんと述べておくと,
|
定理1 [正規行列による直交分解] V 上の線形演算子 T の固有空間: E (λj) によって V が直交分解: V = E (λ1)⊥E (λ2)⊥・・・・・⊥E (λm) [ 直交分解 ] するための必要十分条件は,T が正規演算子(正規行列)であることである。 |
[必要性] T の固有空間 E (λj) によって V が直交分解しているとします。すると,任意のベクトル x ∈V は,
x = x1+x2+・・・+xm ; xj∈E (λj)
と一意的に表され,互いに直交している固有空間,E(λj)への射影演算子 Pm [#]を用いれば,
Tx =Tx1+Tx2+・・・+Txm
=λ1x1+λ2x2+・・・+λmxm
↓ Pjx = xj
=λ1P1x+λ2P2x+・・・+λmPmx
=(λ1P1+λ2P2+・・・+λmPm)x
すなわち,
(1) T = λ1P1 + λ2P2 + ・・・ +λmPm
と射影演算子を用いて展開できることがわかります。この表現を演算子Tのスペクトル分解といいます。
[3] さらに射影演算子の基本的性質 [#] より次の結果,(2),(3)も得られます。
正規演算子T の射影演算子を用いた展開 (1) T = λ1P1 + λ2P2 + ・・・ +λmPm(2) T* =λ1*P1 +λ2*P2 + ・・・+λm*Pm (Pk =Pk*より )つまり,T,T*,T-1 について,
(3) T-1= 1 P1+ 1 P2+・・・+ 1 Pm (PjPk=0 (i≠j), Pj2 =Pj
P1+・・・+Pm+・・・+Pk = E より )λ1 λ2 λm となります。
− T T* T-1 固有空間 E(λm)で共通 固有値 λm λm* 1/λm
すると,(1),(2)を用いて,PjPk=0 (i≠j),Pj2 =Pj に注意して計算すれば,
TT* = T*T (=Σ|λ|2Pk)
であることがわかります[#]。 (必要性の証明終)
[4] [十分性] 今度は逆にTT* = T*T であったとして,V がE (λj) によって直交分解されることを示します。 そのために次の補題が必要です。
補題
線形演算子AとBが可換,すなわち,AB=BAであるとき,AとBに共通な固有ベクトル(=同時固有ケット)が存在する。
証明:
Aの固有値の一つをμ,その固有空間をE(μ),固有ベクトルの一つをxとします。AB=BAであるとき,
A(Bx)=B(Ax)=B(μx)=μ(Bx)
これは,Bx がAの固有ベクトルであって固有値μをもつこと,すなわち,Bx∈E(μ)であることを示しています。したがって,Bの固有値の一つλとそれに対応する固有ベクトルをy,すなわち,
By=λy
を考えると,λy,すなわち,y∈E(μ)であることが分かります。そして,固有空間 E (μ)に属する任意のベクトルの固有値はμである[#]ことから,
Ay=μy
でなければいけません。つまり,このy はAとBに共通な固有ベクトルです。(補題の証明終わり)
[5] この補題からV上の線形演算子TとT*が正規行列(TT*=T*T)であるとき,共通の固有ベクトルz1∈V が存在して,
Tz1=μz1
T*z1=λz1
この固有ベクトルの張る1次元ベクトル空間,
V1={cz1|c∈C}
に対する直交補空間W1を,
W1={w1|(w1,z1)=0} ⇔ V=V1⊥W1
と定義します。すると,
(Tw1,z1)=(w1,T*z1)=(w1,λz1)=λ*(w1,z1)=0
(T*w1,z1)=(w1,Tz1)=(w1,μz1)=μ*(w1,z1)=0
と計算できる[#]ので,「w1∈W1ならば,Tw1,T*w1∈W1である」ことが分かります。つまり,互いに可換なTとT*はV上の線形演算子であるだけでなく,W1上の線形演算子[#]でもあることが分かりました。
[6] すると,W1上の線形演算子TとT*に共通の固有ベクトルz2∈W1が存在して,
V2={cz2|c∈C}
に対する直交補空間W2を考えることができ,W1=V2⊥W2,すなわち,
V=V1⊥V2⊥W2
とすることができます。そして,TとT*はW2上の線形演算子であることを先ほどと同様に示すことができます。
これを繰り返していけば,結局,ベクトル空間Vは,
V=V1⊥V2⊥・・・⊥Vn
と互いに直交する1次元ベクトル空間の直和に分解されます。
[7] ここで,固有値が等しい空間Vj どおしを集めて分類し,
E(λk)=Vk1⊥・・・⊥Vkr ; k=1,2,・・・,m
とすれば,
V=E(λ1)⊥・・・⊥E(λm)
となり,Vがその固有空間によって直交分解されることが分かります。 (証明終わり)
[1] 正規行列に関するもっとも重要な定理をもう一度(再掲)。
| 定理2 [ 正規行列の対角化 ] n次元正規行列Tは適当なユニタリ行列 U によって,対角行列 D に変換できる。すなわち, U-1TU = D なる U が存在する。逆に複素n次元正方行列 T が対角化可能であるならば,T は正規行列である。 |
[証明]
[必要性] T が正規行列 ⇒ T が対角化可能
定理1よりベクトル空間 V 上の線形演算子 T が正規行列であるなら,ベクトル空間 V の任意のベクトルx は一意的に,
x =x1+x2+・・・+xk
ただし,
x1=x1e1 +・・・+xses ∈E(λ1),基底:Σ1={e1, ・・・ ,es} ; dimE(λ1)=s
x2=xs+1es+1+・・・+xs+tes+t ∈E(λ2),基底:Σ2={es+1,・・・,es+t} ; dimE(λ2)=t
・・・
xk=xv+1ev+1+・・・+xv+uev+u ∈E(λm),基底:Σm={ev+1,・・・,ev+u} ; dimE(λm)=u
と書くことができ,x1,x2,・・・,xk はそれぞれ互いに直交しており,n=dimV =dimE(λ1)+dimE(λ2)+・・・+dimE(λm) を満足します。
特に各固有空間 E(λk) の基底Σk を正規直交系になるように選んでおけば[#],これらの基底すべてを並べた,
Σ={e1, ・・・,es,es+1,・・・,es+t,ev+1,・・・,ev+u}
はV の正規直交基底となります。この基底におけるTの表現行列[#]は対角線上に固有値が並んだ対角行列Dです。すなわち,もとの正規直交基底からこの基底Σへの基底変換[#]を表すユニタリ行列をUとすれば,
U-1TU = D =
λ1 O : λ1 λ2 O : λ2 : λm
が成り立ちます。
[2] [十分性] T が対角化可能 ⇒ T が正規行列 を示します。
もし,T があるユニタリ行列で,
U-1TU = D (つまり,T = UDU-1=UDU*, T* =UD*U*= UD*U-1) [#]
と対角化されるならば,
TT* = (UDU-1)(UD*U-1) = UDD*U-1
T*T = (UD*U-1)(UDU-1) = UD*DU-1
となります。ここで対角行列は,DD*=D*D なので結局,TT* = T*T がいえます。 (終)
以下,補遺
特に,Tの固有ベクトルが,e1,e2,・・・en であるときは,
T|ek>=λk|ek>
なので,
T=ΣΣ|ej><ej|T|ek><ek|
=ΣΣ|ej><ej|λk|ek><ek|
=ΣΣλk|ej><ej|ek><ek|
=Σλk|ek><ek|
=ΣλkΛk
T=ΣλkΛk
CopyRight フジエダ電子出版