markdown
対称行列と直交対角化md 9e6fa92
lecture/math/linear-algebra/symmetric-matrices-and-orthogonal-diagonalization.lecture.n.md
Download PDF

対称行列たいしょうぎょうれつsymmetric matrix直交対角化ちょっこうたいかくか

mathlinear-algebraundergraduatelecture

1導入どうにゅう

この講義こうぎ重要じゅうようなのは、対称行列たいしょうぎょうれつsymmetric matrixは、固有値こゆうちeigenvalueだけでなく固有ベクトルeigenvector直交性ちょっこうせいまで保証ほしょうする特別とくべつ行列ぎょうれつmatrixであるということである。

一般いっぱん行列ぎょうれつmatrixは、対角化たいかくかdiagonalizationできるとはかぎらない。対角化たいかくかdiagonalizationできる場合ばあいであっても、基底きていbasis構成こうせいする固有ベクトルeigenvectorたがいに直交ちょっこうorthogonalするとはかぎらない。これにたいして、実対称行列じつたいしょうぎょうれつ正規直交基底せいきちょっこうきていorthonormal basisによって対角化たいかくかdiagonalizationできる。この事実じじつ有限次元ゆうげんじげんスペクトル定理ていりSpectral theoremである。

2用語ようご定義ていぎ

対称行列たいしょうぎょうれつSymmetric matrix とは、実正方行列じつせいほうぎょうれつ A

AT=A

たす行列ぎょうれつmatrixである。

エルミート行列ぎょうれつHermitian matrix とは、複素正方行列ふくそせいほうぎょうれつ A

A*=A

たす行列ぎょうれつmatrixである。ただし A*共役転置きょうやくてんちである。

ユニタリ行列ぎょうれつUnitary matrix とは、複素正方行列ふくそせいほうぎょうれつ U

U*U=UU*=I

たす行列ぎょうれつmatrixである。これは複素ふくそcomplex内積ないせきinner productながさ、直交ちょっこうorthogonal保存ほぞんする行列ぎょうれつmatrixである。

正規行列せいきぎょうれつNormal matrix とは、複素正方行列ふくそせいほうぎょうれつ A

A*A=AA*

たす行列ぎょうれつmatrixである。エルミート行列ぎょうれつmatrixとユニタリ行列ぎょうれつmatrix正規行列せいきぎょうれつ代表例だいひょうれいである。

直交対角化ちょっこうたいかくかOrthogonal diagonalization とは、実行列じつぎょうれつ A直交行列ちょっこうぎょうれつorthogonal matrix Q対角行列たいかくぎょうれつdiagonal matrix D によって

A=QDQT

表示ひょうじすることである。複素ふくそcomplex場合ばあいは、ユニタリ行列ぎょうれつmatrix U により

A=UDU*

表示ひょうじする。

3方針ほうしん

方針ほうしん二段階にだんかいである。まず、対称性たいしょうせい内積ないせきinner productなか

Au,v=u,Av

意味いみすることを確認かくにんする。つぎに、この関係かんけいから、ことなる固有値こゆうちeigenvalue対応たいおうする固有ベクトルeigenvector直交ちょっこうorthogonalすることを導出どうしゅつする。

data/lecture/math/linear-algebra/inner-product-space-basics.lecture.n.md data/lecture/math/linear-algebra/complex-inner-products-and-unitary-matrices.lecture.n.md

4直感的ちょっかんてき説明せつめい

対称行列たいしょうぎょうれつsymmetric matrixは、入力方向にゅうりょくほうこう出力方向しゅつりょくほうこう寄与きよ内積ないせきinner product測定そくていしたとき、左右さゆう交換こうかんしてもあたい変化へんかしない行列ぎょうれつmatrixである。

この対称性たいしょうせいは、空間くうかんゆがみを生成せいせいする変換へんかんではなく、たがいに直交ちょっこうorthogonalする主方向しゅほうこうごとの伸縮しんしゅくとして分解ぶんかいできることを示唆しさする。主方向しゅほうこう固有ベクトルeigenvector伸縮率しんしゅくりつ固有値こゆうちeigenvalueである。

5厳密げんみつ説明せつめい

5.11. ことなる固有値こゆうちeigenvalue固有ベクトルeigenvector直交ちょっこうorthogonalする

A実対称行列じつたいしょうぎょうれつとし、Au=λuAv=μvλμ とする。内積ないせきinner product対称性たいしょうせいから

Au,v=u,Av

である。左辺さへんleft-hand side右辺うへんright-hand side固有値こゆうちeigenvalue表現ひょうげんすると

λu,v=μu,v

である。したがって

(λ-μ)u,v=0

であり、λμ であるため u,v=0 である。よって uv直交ちょっこうorthogonalする。

複素ふくそcomplexエルミート行列ぎょうれつHermitian matrixでもおな議論ぎろん成立せいりつする。A*=A自己随伴性じこずいはんせいself-adjointnessであり、内積ないせきinner productなかA左右さゆう移動いどうしても

Au,v=u,Av

成立せいりつする。したがって固有値こゆうちeigenvalueことなれば固有ベクトルeigenvector直交ちょっこうorthogonalする。また Au=λu のとき

λu,u=Au,u=u,Au=λ_u,u

である。u0 なら u,u>0 なので、λ=λ_ となり、固有値こゆうちeigenvalue実数じっすうreal numberである。

5.22. 実対称行列じつたいしょうぎょうれつのスペクトル定理ていり

有限次元ゆうげんじげん実内積空間じつないせきくうかんで、実対称行列じつたいしょうぎょうれつ A についてつぎ成立せいりつする。

[PARSE ERROR: Undefined("Command(\"boxed\")")]A=QDQT

ここで Q直交行列ちょっこうぎょうれつorthogonal matrixD実固有値じつこゆうちならべた対角行列たいかくぎょうれつdiagonal matrixである。QれつcolumnA正規直交せいきちょっこう固有ベクトルeigenvectorである。

複素ふくそcomplex場合ばあいは、エルミート行列ぎょうれつmatrix A について

[PARSE ERROR: Undefined("Command(\"boxed\")")]A=UDU*

成立せいりつする。ここで U はユニタリ行列ぎょうれつmatrixであり、D対角成分たいかくせいぶん実数じっすうreal numberである。

5.33. 正規行列せいきぎょうれつとの関係かんけい

複素ふくそcomplex有限次元ゆうげんじげん内積空間ないせきくうかんinner product spaceでは、正規行列せいきぎょうれつはユニタリに対角化たいかくかdiagonalizationできる。すなわち

A*A=AA*A=UDU*

である。エルミート行列ぎょうれつmatrix正規行列せいきぎょうれつのうち、固有値こゆうちeigenvalue実数じっすうreal numberになる特別とくべつ場合ばあいである。実対称行列じつたいしょうぎょうれつは、その実数版じっすうばんとして位置いちづけられる。

この整理せいりにより、一般いっぱん対角化可能行列たいかくかかのうぎょうれつ正規行列せいきぎょうれつ、エルミート行列ぎょうれつmatrix実対称行列じつたいしょうぎょうれつ差異さい明確めいかくになる。直交ちょっこうorthogonalまたはユニタリな基底きていbasis対角化たいかくかdiagonalizationできることは、内積構造ないせきこうぞう両立りょうりつすることを意味いみする。

6具体例ぐたいれい

A=(2112)

対称行列たいしょうぎょうれつsymmetric matrixである。固有ベクトルeigenvectorとして

q1=12(11),q2=12(1-1)

採用さいようすると、対応たいおうする固有値こゆうちeigenvalue3,1 である。したがって

Q=12(111-1),D=(3001)

により

A=QDQT

直交対角化ちょっこうたいかくかされる。

7べつ観点かんてん

通常つうじょう対角化たいかくかdiagonalizationは、適切てきせつ基底きていbasis選択せんたくして行列ぎょうれつmatrix単純化たんじゅんかする操作そうさoperationである。直交対角化ちょっこうたいかくかは、さらにその基底きていbasis正規直交基底せいきちょっこうきていorthonormal basisにできるというつよ主張しゅちょうである。

この差異さい計算上けいさんじょう幾何学上きかがくじょう重要じゅうようである。直交行列ちょっこうぎょうれつorthogonal matrixながさと角度かくど保存ほぞんするため、座標変換ざひょうへんかんによって内積構造ないせきこうぞう破壊はかいされない。

8判定基準はんていきじゅん

  • 実対称行列じつたいしょうぎょうれつなら、直交対角化ちょっこうたいかくかできる。
  • エルミート行列ぎょうれつmatrixなら、ユニタリ対角化たいかくかdiagonalizationできる。
  • 複素ふくそcomplexでは、ユニタリ対角化たいかくかdiagonalizationできる行列ぎょうれつmatrix正規行列せいきぎょうれつである。
  • ことなる固有値こゆうちeigenvalue対応たいおうする固有ベクトルeigenvector直交ちょっこうorthogonalする。
  • 固有値こゆうちeigenvalueは、主方向しゅほうこうごとの伸縮率しんしゅくりつとして解釈かいしゃくできる。

9どこまで成立せいりつするか

直交対角化ちょっこうたいかくか任意にんい対角化可能行列たいかくかかのうぎょうれつ成立せいりつするわけではない。直交対角化ちょっこうたいかくかには、内積ないせきinner product両立りょうりつする対称性たいしょうせい必要ひつようである。複素ふくそcomplexでは対称行列たいしょうぎょうれつsymmetric matrixではなくエルミート行列ぎょうれつmatrix基本対象きほんたいしょうにする。

10最終形さいしゅうけい

[PARSE ERROR: Undefined("Command(\"boxed\")")]AT=AA=QDQT
[PARSE ERROR: Undefined("Command(\"boxed\")")]A*=AA=UDU*
[PARSE ERROR: Undefined("Command(\"boxed\")")]A*A=AA*Aisunitarilydiagonalizable
[PARSE ERROR: Undefined("Command(\"boxed\")")]λμEλEμ

11一言ひとことでいうと

  • 対称行列たいしょうぎょうれつsymmetric matrixは、空間くうかん直交ちょっこうorthogonalする固有方向こゆうほうこうごとの伸縮しんしゅく分解ぶんかいできる行列ぎょうれつmatrixである。

じつreal対称行列たいしょうぎょうれつsymmetric matrixは、たがいに直交ちょっこうorthogonalする固有ベクトルeigenvector十分じゅうぶんち、それらから正規直交基底せいきちょっこうきていorthonormal basisつくれる。したがって A=QDQT同値どうちQTAQ=Dける。ここで Q直交行列ちょっこうぎょうれつorthogonal matrixD対角行列たいかくぎょうれつdiagonal matrixである。

重複固有値ちょうふくこゆうち場合ばあいあやまりやすい。ことなる固有値こゆうちeigenvalueぞくする固有ベクトルeigenvector自動的じどうてき直交ちょっこうorthogonalするが、おな固有空間こゆうくうかんeigenspaceなかでは、さらにグラム・シュミットの直交化ちょっこうか正規直交基底せいきちょっこうきていorthonormal basisえらぶことがある。直交対角化ちょっこうたいかくかながさと角度かくどたもつため、二次形式にじけいしきquadratic form主軸しゅじく計算けいさんてきしている。

12定理ていり:スペクトル定理ていり実数版じっすうばん複素数版ふくそすうばん

実数じっすう有限次元ゆうげんじげん内積空間ないせきくうかんinner product spaceでは、実対称行列じつたいしょうぎょうれつ A=AT について、正規直交基底せいきちょっこうきていorthonormal basis q1,[PARSE ERROR: Undefined("Command(\"dots\")")],qn存在そんざいし、それぞれ

Aqi=λiqi

たす。したがって

A=QDQT

である。

複素ふくそ有限次元ゆうげんじげん内積空間ないせきくうかんinner product spaceでは、正規行列せいきぎょうれつnormal matrix A*A=AA* について、ユニタリ行列ぎょうれつ (unitary matrix) U対角行列たいかくぎょうれつdiagonal matrix D存在そんざいして

A=UDU*

である。エルミート行列ぎょうれつHermitian matrixはこの特別とくべつ場合ばあいで、D対角成分たいかくせいぶん実数じっすうになる。

13証明しょうめい方針ほうしん:1 ほん固有こゆうベクトルから直交補空間ちょっこうほくうかん

実対称行列じつたいしょうぎょうれつ場合ばあい特性多項式とくせいたこうしき実根じつこんつことをもちいて、まず単位たんい固有ベクトルeigenvector q1る。つぎに q1かんがえる。xq1 なら

Ax,q1=x,Aq1=λ1x,q1=0

なので、Axq1 である。つまり q1A不変ふへんである。この不変部分空間ふへんぶぶんくうかん帰納法きのうほう適用てきようすると、のこりの正規直交せいきちょっこう固有ベクトルeigenvectorる。

この証明しょうめい要点ようてんは、対称性たいしょうせい直交補空間ちょっこうほくうかんorthogonal complementたもつことである。ここが一般いっぱん行列ぎょうれつでは失敗しっぱいする。

data/lecture/math/linear-algebra/orthogonal-complements-and-projections.lecture.n.md

14反例はんれい対角化可能たいかくかかのうdiagonalizableでも直交対角化ちょっこうたいかくかできない場合ばあい

A=(1102)

固有値こゆうちeigenvalue 1,2つので対角化可能たいかくかかのうである。しかし ATA であり、直交対角化ちょっこうたいかくかはできない。実際じっさい直交行列ちょっこうぎょうれつorthogonal matrix QQTAQ=Dけるなら A=QDQT となり、右辺うへん対称行列たいしょうぎょうれつになるため A対称たいしょうでなければならない。これは矛盾むじゅんである。

この反例はんれいは、「対角化可能たいかくかかのう」と「直交対角化可能ちょっこうたいかくかかのう」がことなる条件じょうけんであることをしめしている。

raw .n.md をコピー
loc をコピー (filepath:line ~ line)
copy share link
copy encoded share link
path をコピー
copy share link
copy encoded share link
copy share link
copy encoded share link
タブを全て閉じる