\usepackageamsmath\usepackageamsfonts\usepackageamssymb\usepackagefancybox\usepackagegraphics\usepackagemathrsfs\usepackage[all]xy\usepackagepgfplots\pgfplotssetcompat=newest\usetikzlibraryintersections,calc,arrows.meta

対角化の基本定理|正方行列の固有値の個数と対角化可能性

線形代数学の基本
線形代数学の基本

正方行列対角化はできればとても便利ですが,残念ながらどんな正方行列でも対角化できるわけではありません.

そのため,正方行列が対角化可能であるかが判定できる方法があると嬉しいですね.

実はn次正方行列が異なるn個の固有値を持てば,その時点で対角化可能であることが従います.また,どのように対角化されるかが計算しなくても分かります.

この記事では

  • 対角化可能性の基本定理と具体例
  • 対角化可能性の基本定理の証明

を順に説明します.

なお,この記事では特に断らない限り複素行列・複素ベクトルを扱います.

線形代数学の参考文献

以下は線形代数学に関するオススメの教科書です.

大学教養 線形代数(加藤文元 著)

数学科など理論系の学生向けの線形代数の入門書です.平易な例から丁寧に説明されています.

手を動かしてまなぶ 線形代数(藤岡敦 著)

理論と演習のバランスをとりながら勉強したい人にオススメの入門書です.

対角化可能性の基本定理

最初に冒頭で紹介した正方行列対角化可能であることが簡単に分かる場合の定理を紹介し,具体例を考えましょう.

定理

次の定理がこの記事の主定理です.

n次正方行列A異なるn個の固有値λ1,λ2,,λnをもつなら,Aは対角化可能である.

また,Aの固有値λ1,λ2,,λnそれぞれに属する固有ベクトルを1つずつ取りp1,,pnとすると,

P=[p1,p2,,pn]

正則行列となり,このPによりAが対角化できる:

P1AP=[λ1000λ2000λn]

つまり,n次正方行列Aが異なる固有値をn個もてば,その時点でAは対角化可能であること分かるわけですね.

さらに,それぞれの固有値に属する固有ベクトルを並べてできる正方行列PによりAは対角化でき,P1APの対角行列の成分は全て固有値になるということも併せて理解しておきましょう.

n次正方行列が異なる固有値をn個もたなくても対角化可能であることもありますが,そのような場合については次の記事で解説しています.

具体例1(2次正方行列)

正方行列A=[1221]を対角化せよ.また,Aを対角化する正則行列Pを1つ求めよ.

この正方行列A固有値1,3

  • 固有値1に属する固有ベクトルc[11]c0
  • 固有値3に属する固有ベクトルはd[11]d0

であることは前回の記事で示した.

よって,正方行列Aは異なる固有値を2個もつ2次正方行列だから,上で示した対角化可能性の基本定理より対角化可能である.

さらに,例えばc=d=1とすると

  • 固有値1に属する固有ベクトル[11]
  • 固有値3に属する固有ベクトル[11]

が得られ,これらの固有ベクトルを並べてP=[1111]とすると,上の定理より

P1AP=[1003]

と対角化可能である.

わざわざP1APを計算しなくても,固有値が対角成分に並ぶことが分かるのは強力ですね.

Pに並べる列ベクトルは固有ベクトルであればいいので,例えばP=[2323]P=[2π2π]などとしても

P1AP=[1003]

と対角化可能されます.

また,Pに並べる固有ベクトルの順番とP1APの対角成分の順番は対応するので,例えばP=[1111]と固有ベクトルを逆に並べた場合は対角化も

P1AP=[3001]

と対角成分の固有値は逆になります.

具体例2(3次正方行列)

3次正方行列に対して,固有値固有ベクトルを求めるところから対角化まで通してみましょう.

固有値・固有ベクトルの求め方は固有方程式固有方程式)を用いる方法が基本的だったことを思い出しておきましょう.

正方行列A=[220010301]を対角化せよ.また,Aを対角化する正則行列Pを1つ求めよ.

正方行列Aの固有多項式は

|xIA|=|x2200x+1030x1|=(x2)(x+1)(x1)

なので,固有方程式|xIA|=0の解は2,1,1である.

よって,正方行列Aは異なる固有値を3個もつ3次正方行列だから,上で示した対角化可能性の基本定理より対角化可能である.

さらに,固有値2,1,1それぞれに属する固有ベクトルp1,p2,p3を任意にとり,P=[p1,p2,p3]とおくと

P1AP=[200010001]()

と対角化される.よって,あとはp1,p2,p3を求めればよい.

[1] 固有値2に属する固有ベクトルp1Ap1=2p1(A2I)p1=0を満たし,行基本変形より

A2I=[020030301][301010000]

となるから,掃き出し法よりp1=c1[103] (c10)である.

[2] 固有値1に属する固有ベクトルp2Ap2=1p2(AI)p2=0を満たし,行基本変形より

AI=[120020300][100010000]

となるから,掃き出し法よりp2=c2[001] (c20)である.

[3] 固有値1に属する固有ベクトルp3Ap3=(1)p3(A+I)p3=0を満たし,

A+I=[320000302]

となるから,p3=c3[233] (c30)である.

以上より,[1]-[3]で例えばc1=c2=c3=1として得られる固有ベクトルを並べたP=[102003313]によって,A()と対角化可能である.





対角化可能性の基本定理の証明

まずは定理の証明の鍵となる命題を示し,上の定理を証明しましょう.

異なる固有値に属する固有ベクトル

正方行列Aの異なる固有値λ1,,λrに対して,それぞれの固有値に属する固有ベクトルa1,,arとすると,a1,,ar線形独立である.

数学的帰納法により示す.一般にただひとつのベクトルは線形独立なのでa1は線形独立である.

a1,,amm<r)は線形独立であると仮定する.線形関係

c1a1++cm+1am+1=0()

を考える.線形関係()の両辺に左からTを施して

λ1c1a1++λm+1cm+1am+1=0

であり,線形関係()の両辺にλm+1をかけて

λm+1c1a1++λm+1cm+1am+1=0

である.これら2式の辺々引いて

(λm+1λ1)c1a1++(λm+1λm)cmam=0

となる.帰納法の仮定よりa1,,amは線形独立なので(λm+1λk)ck=0k=1,,m)が成り立つ.

固有値λ1,,λm+1は異なると仮定していたから,λr+1λk0なのでck=0を得る.

これをもとの線形関係()に代入するとcm+1am+1=0となるから,am+10よりcm+1=0を得る.

よって,a1,,ar,am+1も線形独立である.

この命題から次が成り立ちますね.

n次正方行列Aが異なるn個の固有値をもつとき,それぞれの固有値に属する固有ベクトルを1つずつ取りp1,,pnとすると,正方行列[p1,,pn]正則行列である.

一般にa1,,anCnが線形独立なら正方行列[a1,,an]は正則である.

また,上で示した命題からp1,,pnは線形独立だから,n次正方行列P=[p1,,pn]は正則行列である.

定理の証明

それでは本題の定理を証明しましょう.

n正方行列A異なるn個の固有値λ1,λ2,,λnをもつとする.

このとき,Aの固有値λ1,λ2,,λnそれぞれに属する固有ベクトルを1つずつ取りp1,,pnとすると,

P=[p1,p2,,pn]

正則行列となり,このPによってA

P1AP=[λ1000λ2000λn]

対角化可能である.

対角行列B

B=[λ1000λ2000λn]

で定める.固有値・固有ベクトルの定義よりλkpk=Apk (k=1,2,,n)が成り立つから,

PB=[p1,p2,,pn][λ1000λ2000λn]=[λ1p1,λ2p2,,λnpn]=[Ap1,Ap2,,Apn]=A[p1,p2,,pn]=AP

となる.いまλ1,,λnは全て異なるから,上で示した系よりP=[p1,,pn]は正則行列である.

よって,PB=APの両辺に左からP1をかけてB=P1APと対角化可能である.

証明を読んでみるとp1,p2,,pnが線形独立でない場合でもPB=APまでは成り立ちますね.

そこからB=P1APの形に変形するにはPが正則行列でなければならず,それはAが異なる固有値をn個もつことから従うわけですね.

コメント