...

固有空間やJordan標準形の説明

by user

on
Category: Documents
11

views

Report

Comments

Transcript

固有空間やJordan標準形の説明
Contents
1. 多項式に関する基本的な扱い
2
2. 固有多項式
4
3. 正方行列の三角化
5
4. 固有空間 E(A, α), 広義固有空間 V (A, α)
7
5. 正方行列の対角化
9
6. 最小多項式
11
1
JORDAN 標準形に関する基本事項
学習院大学 澤野嘉宏
1. 多項式に関する基本的な扱い
ここでは,多項式に関して必要な基礎事項をまとめる.初めに,割り算とは何かから復習する.
多項式 A(X) を B(X) で割り算するとは,B(X) よりも次数が低い剰余と呼ばれる多項式 R(X)
と商と呼ばれる多項式 Q(X) を用いて A(X) = B(X)Q(X) + R(X) と表されることである.この
Q(X) と R(X) は整数の割り算と同じように筆算で求められるが,その方法には深入りしない.
次に,剰余の定理を復習しよう.
定理 1. α を定数とする.P (X) を (X − α) で割り算したときの剰余 (あまり) は P (α) で与えら
れる.
証明. あまりは 0 次でないといけないから,定数多項式 R を用いて P (X) = (X − α)Q(X) + R
と表される.この恒等式で X = α とすると,P (α) = R が得られる.したがって,P (α) が余り
である.
1 次式の割り算を特に考えたいのであるが,係数が R だと,X 2 + 1 のように因数分解が出来
ないものが出てくる.C にするとこのような問題が解消される.
定理 2 (代数学の基本定理). P (X) が複素数係数の多項式のとき,
P (X) = β(X − α1 )(X − α2 ) · · · (X − αn ), β ̸= 0
という因数分解が出来る.ここで,n は P (X) の次数.
もしくは,同一の因子をまとめてしまって次のようにして表わしてもよい.
定理 3 (代数学の基本定理). 複素係数多項式 P (X) は複素数 α と相異なる複素数 β1 , · · · , βk を
用いて P (X) = α(X − β1 )n1 · · · (X − βk )nk と因数分解する.
例 1.1 (複素数係数の多項式の因数分解の例).
(1) X 3 − 3iX 2 − 3X + i = (X − i)3 である.
(2) X 2 + 1 は実数係数の範囲ではこれ以上因数分解できないが,複素数係数では,X 2 + 1 =
(X − i)(X + i) と因数分解される.
これは有名な定理なので認めてしまう.20通りの証明があるらしい.
定理 4. 定数ではない多項式 P (X) と Q(X) が複素数の範囲で共通根を持たないとする.このと
き,ある多項式 R(X) と S(X) が存在して
(1.1)
P (X)R(X) + Q(X)S(X) = 1
が成り立つ.
2
JORDAN 標準形に関する基本事項
3
証明. (1.1) のようにある多項式 R(X) と S(X) が存在して T (X) = P (X)R(X) + Q(X)S(X) と
表される多項式全体の集まりを M とおく.M の中で 0 でない多項式のうち次数が最低のものを
V (X) とおく.M の中の多項式だから,当然ある多項式 R(X) と S(X) が存在して
(1.2)
P (X)R(X) + Q(X)S(X) = V (X)
が成り立つ.
主張 1. P (X) は V (X) で割り切れる.つまり,P (X) を V (X) で割ると余りが出ない.
主張 1 の証明. (背理法で) もしそうではないと仮定しよう.P (X) を V (X) で割り算して P (X) =
V (X)W (X) + U (X) と割り算すると U (X) の次数は V (X) より低くなってしまう.ここで,U (X)
は余りで,W (X) は商である.すると,
(1.3)
U (X) = −V (X)W (X) + P (X) = P (X)(1 − R(X)W (X)) + (−S(X)W (X))Q(X)
となる.R∗ (X) = 1 − R(X)W (X), S ∗ (X) = −S(X)W (X) とおく.先ほどの計算結果 (1.3) か
ら U (X) = R∗ (X)P (X) + S ∗ (X)Q(X) ∈ M であるから V (X) が M の中で次数が最低であるこ
とに反する.
主張 2. Q(X) も V (X) で割り切れる.
主張 2 の証明. P (X) のときと同じである.
したがって,主張 1 と 2 より,V (X) の根があると P (X), Q(X) の共通根となる.仮定によっ
て,これはありえないから,V (X) = c ̸= 0 となる.P (X)[c−1 R(X)] + Q(X)[c−1 S(X)] = 1 で
あるから,R(X), S(X) をおきなおせば定理にあるような多項式が取れる.
【注意】主張とは証明中などに現れる証明の方針を示したものである.主張を証明しないと定
理の証明が完成しない.
定義 1.2. 複素数係数の多項式 P (X) と Q(X) の最大公約式とは P (X) と Q(X) を割ることがで
きる複素数係数の多項式の中で最高次数のものをさす.
先ほどの定理は汎用性があるので,共通因数をもつ場合に書き改めておこう.
定理 5. P (X) と Q(X) が 0 でない複素係数多項式とする.S(X) を最大公約式とすると,
S(X) = P (X)T (X) + Q(X)V (X)
と多項式 T (X), V (X) を用いて表される.
証明. P (X), Q(X) を S(X) で割った多項式に先ほどの定理を適用すればよい.
帰納法で次が成り立つ.
定理 6. P1 (X), P2 (X), · · · , Pn (X) を複素係数多項式とする.これらの最大公約式を P (X) とす
るとき,
Q1 (X)P1 (X) + Q2 (X)P2 (X) + · · · + Qn (Q)Pn (X) = P (X)
が成り立つ.
学習院大学 澤野嘉宏
4
2. 固有多項式
定義 2.1. A を複素 n × n-行列とする.A の (複素) 固有値とは A x = α x を満たしている 0 でな
い (複素) ベクトル x が存在するような複素数 α のことである.
形式的に,実数の固有値も定義できる.当然次のように定義する.
定義 2.2. A を実 n × n-行列とする.A の (実) 固有値とは A x = α x を満たしている 0 でない
(実) ベクトル x が存在するような実数 α のことである.
ここで,
A x = α x ⇐⇒ A x = α E x ⇐⇒ (A − α E) x = 0 ⇐⇒ (α E − A) x = 0
に注意すると,連立方程式の理論から次のことが導き出される.
定理 7. A を複素 n × n-行列とする.α ∈ C が A の固有値であるための必要十分条件は
det(α E − A) = 0
で与えられる.
定義 2.3. A を複素 n × n-行列とする.λ に関する方程式
det(λ E − A) = 0
を固有方程式という.
定義を与えたので,先ほどの定理を言い換えると次のようになる.
定理 8. A を複素 n × n-行列とする.α ∈ C が A の固有値であるための必要十分条件は α が A
の固有方程式の解であることである.
行列式の定義に戻って計算するとわかるように,
det(λ E − A) = λn + · · ·
の形をした n 次多項式であるから,代数学の基本定理より次のことがいえる.
定理 9. A を複素 n × n-行列とする.A には 1 個から n 個の (複素) 固有値が存在する.
実数に関しても事情は同じである.ただし,代数学の基本定理が使えないので,固有値がない
こともある.
定理 10. A を複素 n × n-行列とする.A に 0 個から n 個の (複素) 固有値が存在する.
代数学の基本定理が使えないと固有値が得られないかもしれないので,複素数を使うことが多
い.以後,
(代数学の基本定理を活用したいがために)複素数のときのみに説明をしていく.
JORDAN 標準形に関する基本事項
5
3. 正方行列の三角化
線形変換 f が与えられたとする.ある基底を取ったときに,A という行列で f を記述できた
とする.基底を取り替えれば,それに応じた可逆行列 P で P −1 AP で f は記述されることにな
る.したがって,P −1 AP が“ きれいな ”行列で表されれば,それだけ f を捕らえるのが楽にな
る.何を以って“ きれい ”というかわからないが,複素正方行列に対しては次のことが言える.
定義 3.1. 上三角行列とは,A = {aij }i,j=1,··· ,n と表したときに,aij = 0, i > j が成立すること
である.
定理 11. A を複素正方行列とするとき,適当な対角行列 P を以って P −1 AP を上三角行列にで
きる.
証明. 帰納法によって証明する.n = 1 のときは明らかである.n = k 次正方行列 B はすべて適
当な対角行列 Q を以って Q−1 AQ を上三角行列にできると仮定する.n = k + 1 次正方行列 A が
与えられたとする.A の固有値 α をひとつ取る.A の α に固有ベクトル v0 をとる.v1 , v2 , · · · , vk
を v0 , v1 , v2 , · · · , vk が Cn の基底になるように取る.
Av0 = α v0 , Avj =
k
∑
αjl vl
l=0
とする.行列 P0 , X を

α
0

0

P0 = (v0 v1 · · · vn ), X =  .
 ..

0
0
α01
α11
α21
..
.
α02
α12
α22
..
.
···
···
···
..
.
α0n
α1n
α2n
..
.
αn−1 1
αn1
αn−1 2
αn2
···
···
αn−1
αnn









n
となるようにすると,P0−1 AP0 = X となる.ここで,帰納法の仮定を使うと k 次正方行列 Q0 で




α11
α12
···
α1n
β11 β12 · · ·
β1n−1
β1n
 α21
 0 β22 · · ·
α22
···
α2n 
β2n−1
β2n 







.
.
.
.
.
.
−1 
.
.
..
..
..  Q0 =  ..
..
..
..
Q0  ..





αn−1 1 αn−1 2 · · · αn−1 n 
 0
0 · · · βn−1 n−1 βn−1 n 
αn1
αn2
···
αnn
0
0 ···
0
βnn
(
)
1 0
となる.ここで,P1 =
とおくと,
0 Q0


α γ01 γ02 · · ·
γ0n−1
γ0n
 0 β11 β12 · · ·
β1n−1
β1n 


0
0 β22 · · ·
β2n−1
β2n 


−1
−1
P1 P0 AP0 P1 =  .

..
..
..
..
 ..

.
.
.
.


0
0
0 · · · βn−1 n−1 βn−1 n 
0
0
0 ···
0
βnn
となり,k + 1 次正方行列 A も P = P0 P1 とおけば,P −1 AP = P1 −1 P0 −1 AP0 P1 が正方行列に
なることがわかった.
この定理の大事な応用としてケーリーハミルトンの定理がある.
系 3.2 (ケーリーハミルトンの定理). ΦA (λ) = det(λ E − A) とする.すると,ΦA (A) = 0 が成
り立つ.
6
学習院大学 澤野嘉宏
証明. 定理 11 によって A が上三角行列であるとしてよい.もし,


a00 a01 a02 · · ·
a0n−1
a0n
 0 a11 a12 · · ·
a1n−1
a1n 


 0
0
a
·
·
·
a
a2n 
22
2n−1


A= .

..
..
..
..
 ..

.
.
.
.


 0
0
0 · · · an−1 n−1 an−1 n 
0
0
0 ···
0
ann
のときは,
ΦA (λ) = (λ − a00 )(λ − a11 ) · · · (λ − ann )
となるから,具体的に計算できる.
JORDAN 標準形に関する基本事項
7
4. 固有空間 E(A, α), 広義固有空間 V (A, α)
(A の固有値)α と n ∈ N に対して,線形空間 Vn (A, α) を
Vn (A, α) = {x ∈ Cn : (A − αI)n x = 0}
と定める.Vn (A, α) をもとにして
E(A, α) = V1 (A, α), V (A, α) =
∞
∪
Vn (A, α)
n=1
と定める.
【注意】
(1) V (A, α), Vn (A, α), E(A, α) はみな一般的な記号ではなくここだけで通用する記号と考え
たほうがよい.
(2) E(A, α) ⊂ V (A, α) が成り立つ.
V (A, α) と E(A, α) の違いを一番簡単で顕著な例を見てみよう.
(
)
0 1
例 4.1. A =
の固有値は 0 のみで,
0 0
{( )
(
)( )
}
x1
0 1
x1
E(A, 0) = V1 (A, 0) =
∈ C2 :
=0
x2
0 0
x2
{( )
( )
}
x1
x
=
∈ C2 : A 1 = 0
x2
x2
{( )
}
x1
=
∈ C 2 : x1 ∈ C
0
である一方,n ≥ 2 のとき,An = 0 であるから,
{( )
(
)n ( )
}
x1
0 1
x1
Vn (A, 0) =
∈ C2 :
=0
x2
0 0
x2
{( )
(
)( )
}
x1
0 0
x1
=
∈ C2 :
=0
x2
0 0
x2
{( )
}
x1
=
∈ C2 : x1 , x2 ∈ C = C2
x2
となる.したがって,V (A, 0) =
∞
∪
Vn (A, 0) = C2 が得られる.以上より,固有空間 E(A, 0) と
n=1
広義固有空間 V (A, 0) は確かに異なって,E(A, 0) ( V (A, 0) となる.
次の命題は,広義固有空間の定義を覚えやすくしてくれる.
命題 4.2. α ∈ C を n × n-行列 A の固有値とするとき,次の等式が成り立つ.
E(A, α) = {x ∈ Cn : (A − α E)x = 0}
V (A, α) = {x ∈ Cn : (A − α E)n x = 0}
証明. E(A, α) のほうは定義そのものだから,V (A, α) の等式を証明する.一方の包含関係 V (A, α) ⊂
∞
∪
{x ∈ Cn : (A − α E)n x = 0} は明らかである.そこで,x ∈ V (A, α) とする.すると,
の定
k=1
学習院大学 澤野嘉宏
8
義によって x ∈ Vk (A, α) となる k が存在する.(λ − α)k と ΦA (λ) の最大公約式は (λ − α)l , l ≤ k
である.したがって,
P (λ)(λ − α)k + Q(λ)ΦA (λ) = (λ − α)l
したがって,
(A − α E)n x = (P (A)(A − α E)k + Q(A)ΦA (A))(A − α E)n−l x
= (A − α E)n−l P (A){(A − α E)k x} + Q(A)(A − α E)n−l {ΦA (A)x} = 0
となる.
次の定理は重宝する.
定理 12. A の固有値を α1 , · · · , αr とする.このとき,
Cn = V (A, α1 ) ⊕ V (A, α2 ) ⊕ · · · ⊕ V (A, αr )
が成立する.
証明. ΦA (λ) = (λ − α1 )m1 · · · (λ − αr )mr と表す.j = 1, 2, · · · , r に対して φj (λ) =
とおく.多項式 ψ1 , ψ2 , · · · , ψr を
ΦA (λ)
(λ − αj )mj
ψ1 (λ)φ1 (λ) + ψ2 (λ)φ2 (λ) + · · · + ψr (λ)φr (λ) = 1
となるように取る.
j = 1, 2, · · · , r と x ∈ Cn に対して,(A − αj E)mj ψj (A)φj (A)x = ψj (A)ΦA (A)x = 0 となる
から,ψj (A)φj (A)x ∈ V (A, αj ) が成り立つ.つまり,
Cn = V (A, α1 ) + V (A, α2 ) + · · · + V (A, αr )
が証明された.
最後に,この等式の + は ⊕ で置き換えられることを示そう.仮に
x1 + x2 + · · · + xr = 0, x1 ∈ V (A, α1 ), x2 ∈ V (A, α2 ), · · · , xr ∈ V (A, αr )
を満たしていたとする.xj ∈ V (A, αj ) よりある k が存在して,(A − αj E)k xj = 0 が成り立つ.
このような k に対して,Pj (λ)(λ − αj )k + Qj (λ)φj (λ) = 1 となるような多項式を取る.すると,
∑
xj = Pj (A){(A − αj E)k xj } + Qj (A)φj (A)xj = Qj (A)φj (A)xj = −
Qj (A)φj (A)xl = 0.
l̸=j
したがって,+ は ⊕ で置き換えられる.
次の事実は実際に計算するときに使う大事なので,きちんとまとめておくことにする.
系 4.3. A の固有値を α1 , · · · , αr とする.j = 1, 2, · · · , r に対して Xj = {xj,l }l∈Lj ⊂ V (A, αj )
r
∪
が与えられたとする.このとき,各 Xj = {xj,l }l∈Lj が一次独立ならば,
Xj も一次独立である.
j=1
証明. 仮に,係数の組 λj,l ∈ C を与えて
r ∑
∑
λj,l xj,l = 0
j=1 l∈Lj
が成り立ったとすると,Cn = V (A, α1 ) ⊕ V (A, α2 ) ⊕ · · · ⊕ V (A, αr ) なので,
∑
λj,l xj,l = 0 と
l∈Lj
なる.Xj = {xj,l }l∈Lj が一次独立なので,λj,l = 0 となる.
JORDAN 標準形に関する基本事項
9
5. 正方行列の対角化
定義 5.1. 上三角行列とは,A = {aij }i,j=1,··· ,n と表したときに,aij = 0, i ̸= j が成立すること
である.
命題 5.2. A を n 次正方行列とする.A の固有値を α1 , α2 , · · · , αk と書くとき,A が対角化可能
であることと,
k
∑
dimC E(A, αj ) = n
j=1
は同値である.
証明. A が対角化可能であろうがなかろうが,
k
∑
dimC E(A, αj ) ≤ n
j=1
に注意する.
A が対角化可能であるとする.すると,一次独立なベクトル v1 , v2 , · · · , vn が存在して,P =
(v1 v2 · · · vn ) とおくとき,


β1 0
0 ··· 0
 0 β2 0 · · · 0 



0 β3 · · · 0 
P −1 AP =  0

 ..
..
..
.. 
..
.
.
.
.
. 
0
0
0
···
βn
が成り立つ.ここで,P を取り替えて
α1 = β1 , β2 , · · · , βm1
α2 = βm1 +1 , βm1 +2 , · · · , βm1
..
.
αr = βm1 +···+mr−1 +1 , βm1 +···+mr−1 +2 , · · · , βm1 +···+mr−1 +mr
としてよい.すると,P を構成する列ベクトルは一次独立であるから,
k
∑
dimC E(A, αj ) ≥ n
j=1
が成り立つ.したがって,
k
∑
dimC E(A, αj ) = n
j=1
がいえた.
逆に
k
∑
j=1
dimC E(A, αj ) = n
10
学習院大学 澤野嘉宏
とする.mj = dimCE(A, αj ) と書く.また,β1 , β2 , · · · , βn を
β1 , β2 , · · · , βm1 = α1
βm1 +1 , βm1 +2 , · · · , βm1 +m2 = α2
..
.
βm1 +···+mr−1 +1 , βm1 +···+mr−1 +2 , · · · , βm1 +···+mr−1 +mr = αr
と定める.仮定によって一次独立なベクトル v1
, v2 , · · · , vn を A vj = αj
vj となるようにとる.そ
β1 0
0 ··· 0
 0 β2 0 · · · 0 



−1
0 β3 · · · 0 
こで,P = (v1 v2 · · · vn ) とおけば,P AP =  0
 が成り立つ.
 ..
..
..
.. 
..

.
.
.
.
.
0
0
0 · · · βn
JORDAN 標準形に関する基本事項
11
6. 最小多項式
P (A) = 0 となる 0 ではない次数最小の多項式 P (t) を A の最小多項式という.φA (t) と書く.
固有多項式はここでは,ΦA (t) と表すことにする.
命題 6.1.
(1) 固有多項式 ΦA (t) は最小多項式 φA (t) で割り切れる.
(2) A の固有値が α のとき,最小多項式は t − α で割り切れる.
命題 6.2. 固有値に対応する固有ベクトルは 1 つは存在する.
定理 13. 行列 A が対角化可能であるための必要十分条件は φA (λ) = 0 が重解を持たないことで
ある.
証明. φA (λ) = 0 が重解を持たないとすると,φA (λ) = (λ − α1 ) · · · (λ − αr ) と因数分解できる.
φA (λ)
φj (λ) =
とすると,多項式 P1 (λ), P2 (λ), · · · , Pr (λ) で
λ − αj
P1 (λ)φ1 (λ) + P2 (λ)φ2 (λ) + · · · + Pr (λ)φr (λ) = 1
が得られる.j = 1, 2, · · · , r と x ∈ Cn に対して,(A − αj E)Pj (A)φj (A)x = Pj (A)ΦA (A)x = 0
となるから,Pj (A)φj (A)x ∈ V (A, αj ) が成り立つ.つまり,
Cn = E(A, α1 ) + E(A, α2 ) + · · · + E(A, αr )
が証明された.ここで,
Cn = V (A, α1 ) ⊕ V (A, α2 ) ⊕ · · · ⊕ V (A, αr )
n
∑
であったから,E(A, αj ) = V (A, αj ) となり,
dimE(A, αj ) = n つまり A が対角化可能であ
j=1
ることが得られた.

α1
0


A が対角化可能であるとすると,P −1 AP =  0
 ..
 .
0
α2
0
..
.
0
0
α3
..
.
···
···
···
..
.
0
0
0
..
.




 となる対角行列 P が存


0
0
0 · · · αn
在する.α1 , α2 , · · · , αn のうち,重複しているものを除いて αk1 , αk2 , · · · , αkr と表すと,P (λ) =
(λ − αk1 ) · · · (λ − αkr ) に対して,P (A) = 0 となることは対角行列の性質より明らかである.最
小多項式は P (λ) の約数でないといけないから,φA (λ) = 0 は重解を持たない.
Fly UP