固有値問題・固有空間・スペクトル分解 のバックアップ(No.3)

更新


前の単元 <<<?                線形代数Ⅱ?                >>> 次の単元?

$n$ 次正方行列 $A$ の固有値問題

1年生でしっかりやったはずなので、ここでは簡単におさらい。

固有値問題

A\bm x=\lambda\bm x   ただし  \bm x\ne\bm 0    ( \lambda :固有値、 \bm x :固有ベクトル)

(A-\lambda E)\bm x=\bm 0

もし (A-\lambda E) が正則なら \bm x=(A-\lambda E)^{-1}\bm 0=\bm 0 となってしまってダメなので、

|A-\lambda E|=0   

これは \lambda に関する n 次方程式( A の固有方程式と呼ぶ) となるので、これを複素数の範囲で解けば、重複度を含めて n 個の \lambda が求まる(代数学の基本定理)

求まったそれぞれの \lambda について、

(A-\lambda E)\bm x=\bm 0 \bm x について解くことで、

1コ以上、重複度コ以下の一次独立な固有ベクトルの線形結合の形で \bm x が求まる

練習

&math( A=\begin{pmatrix} \cos\theta&-\sin\theta\\ \sin\theta&\cos\theta\\ \end{pmatrix} ); について固有値問題を解け。ただし \theta!=2m\pi

A は回転を表わす行列。一方で、固有値問題は A\bm x=\lambda\bm x だから、 A を掛けても向きの変わらないベクトルを探す問題。すなわち・・・解無しになるのでは? と思うのが正しい感覚。

実際、実数の範囲では固有値は存在しない。一方、複素数の範囲ではどんな行列にも必ず重複度を含めて n コの固有値が存在し、それぞれの固有値に対して最低1コの固有ベクトルが見付かる。

回転を表わす上記の A に対しても、複素数の範囲ならばちゃんと固有値・固有ベクトルが見付かることを確認せよ。

対角化

すべての固有値に対して、ちょうど重複度コの一次独立な固有ベクトルが求まる場合には (必ずしも見付からない)、それらをすべて集めると n コの一次独立な固有ベクトルが見つかる。

A\bm x_1=\lambda_1\bm x_1

A\bm x_2=\lambda_2\bm x_2

\phantom{A\bm x_2}\vdots

A\bm x_n=\lambda_n\bm x_n

これらをまとめて書くと、

&math( A\underbrace{\Big(\bm x_1\ \bm x_2\ \dots\ \bm x_n\Big)}_P &=\Big(\lambda_1\bm x_1\ \lambda_2\bm x_2\ \dots\ \lambda_n\bm x_n\Big)\\ &=\underbrace{\Big(\bm x_1\ \bm x_2\ \dots\ \bm x_n\Big)}_P \underbrace{\begin{pmatrix} \lambda_1\\ &\lambda_2\\ &&\ddots\\ &&&\lambda_n \end{pmatrix}}_\Lambda );

AP=P\Lambda

P の列ベクトルが一次独立であることから P は正則なので、 左から P^{-1} を掛ければ、

P^{-1}AP=\Lambda

として、正則行列 P により A を対角化可能。

対称行列の場合

  • 固有値 \lambda はすべて実数 \lambda_k\in\mathbb R
  • 直交行列 {}^t\!R=R^{-1} により対角化可能 R^{-1}AR={}^T\!RAR=\Lambda

1年生の時に証明した。

線形変換の固有値問題

線形変換 T:V\to V に対して、

   T\bm x=\lambda\bm x   ただし  \bm x\ne\bm 0

を満たすスカラー \lambda を固有値、ベクトル \bm x を固有ベクトルと呼ぶ。

この方程式の、ある基底 B に対する表現は

   T_B\bm x_B=\lambda\bm x_B   ただし  \bm x_B\ne\bm 0

となるから、行列 T_B の固有値問題を解けば T の固有値問題が解ける。

   (Tの固有値)&=(T_Bの固有値)\\&= (T_{B'}の固有値)

ただし基底 B' は別の基底であり、固有値は基底の取り方に依存しない。

  基底 B から基底 B' への基底の変換行列を P_{B\to B'} とすると、

     T_{B'}=P_{B\to B'}^{-1}T_BP_{B\to B'}

  であり、異なる基底に対する T の表現 T_B,T_{B'}
  互いに相似の関係にある。固有値が等しいのは当然!

一方、固有ベクトル \bm x に対する「表現」 \bm x_B,\bm x_B' は、 基底の取り方に依存する

練習

x の2次以下の多項式からなる線形空間 P^2[x] で定義された線形変換

   T:f(x)\mapsto f(x+1)

の固有値、固有ベクトルを求めよ。

  &math( T(ax^2+bx+c) &=a(x+1)^2+b(x+1)+c\\ &=ax^2+(2a+b)x+(a+b+c) );

基底を B=\langle x^2,x,1\rangle と取れば、

  &math( T_B=\begin{pmatrix} 1&0&0\\ 2&1&0\\ 1&1&1\\ \end{pmatrix} );

であるから、

  &math( |T_B-\lambda E|= \begin{vmatrix} 1-\lambda&0&0\\ 2&1-\lambda&0\\ 1&1&1-\lambda\\ \end{vmatrix}=(1-\lambda)^3=0 );

\therefore \lambda=1

  &math( (T_B-E)\bm x=\begin{pmatrix} 0&0&0\\ 2&0&0\\ 1&1&0\\ \end{pmatrix}\begin{pmatrix} a\\b\\c \end{pmatrix}=\bm 0 );

   a=b=0

   c=s と置けば、

  &math( \begin{pmatrix} a\\b\\c \end{pmatrix}= s \begin{pmatrix} 0\\0\\1 \end{pmatrix} );

&math( \bm x_B=\begin{pmatrix} 0\\0\\1 \end{pmatrix} ); に対応するベクトルは \bm x=0x^2+0x+1=1

すなわち、固有値 \lambda=1 、固有ベクトル 1 が答えになる。

固有空間

ある固有値 \lambda に属する固有ベクトルは部分空間 V(\lambda)=\set{\bm x|A\bm x=\lambda\bm x} を作る。

\because A\bm x=\lambda\bm x,A\bm x=\lambda\bm x ならば A(a\bm x+b\bm y)=\lambda(a\bm x+b\bm y) より、 V(\lambda) はベクトルの和とスカラー倍に対して閉じている。

この部分空間 V(\lambda) A \lambda に属する固有空間と呼ぶ。

固有値 \lambda に属する一次独立な固有ベクトルが r コあって \bm x_1,\bm x_2,\dots,\bm x_r であるならば、 V(\lambda) は これらのベクトルにより張られる空間に他ならない。

V(\lambda)=[\bm x_1,\bm x_2,\dots,\bm x_r]=\set{\bm x|\bm x=\sum_{k=1}^rc_k\bm x_k}

対角化可能な場合

特殊な行列の固有値

エルミート変換・エルミート行列

エルミート変換: (T\bm x,\bm y)=(\bm x,T\bm y)
エルミート行列: T_B^\dagger=T_B   ただし B は正規直交基底

エルミート変換の正規直交基底に対する表現行列はエルミート行列になる

これらの固有値は実数になる。

∵&math( &(\bm x,T\bm x)= \lambda(\bm x,\bm x)= \lambda|\bm x|^2=\\ &(T\bm x,\bm x)=\overline\lambda(\bm x,\bm x)=\overline\lambda|\bm x|^2\\ );

(\lambda-\overline\lambda)|\bm x|^2=0

\lambda=\overline\lambda

\lambda\in\mathbb R

注:エルミート行列は実行列では対称行列に対応する

ユニタリ変換・ユニタリ行列

ユニタリ変換: (T\bm x,T\bm y)=(\bm x,\bm y)    内積を保存する変換のこと
ユニタリ行列: T_B^\dagger=T_B^{-1}   ただし B は正規直交基底

ユニタリ変換の正規直交基底に対する表現行列はユニタリ行列になる

これらの固有値の絶対値は1になる。

∵&math( &(T\bm x,T\bm x)=\lambda\overline\lambda(\bm x,\bm x)=|\lambda|^2|\bm x|^2=\\ &(\bm x,\bm x)=|\bm x|^2\\ );

(|\lambda|^2-1)|\bm x|^2=0

|\lambda|=1

注:ユニタリ行列は実行列では直交行列に対応する

&math( A=\begin{pmatrix} \cos\theta&-\sin\theta\\ \sin\theta&\cos\theta\\ \end{pmatrix} ); はユニタリ行列。

正規行列 $A^\dagger A=AA^\dagger$ はユニタリ行列で対角化可能

正規行列であること A^\dagger A=AA^\dagger

ユニタリ行列 U^\dagger=U^{-1} で対角化可能 U^{-1}AU=U^\dagger AU=\Lambda は同値である。

詳しい証明は 線形代数Ⅱ/正規行列の対角化可能性? を参照のこと。


前の単元 <<<?                線形代数Ⅱ?                >>> 次の単元?

質問・コメント





Counter: 143438 (from 2010/06/03), today: 8, yesterday: 0