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

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

mathlinear-algebraundergraduatelecture

Symmetric matrices対称行列たいしょうぎょうれつ and orthogonal diagonalization

1導入どうにゅう

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

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

1Introduction

The important point of this lecture is that a symmetric matrix対称行列たいしょうぎょうれつ is special not only because of its eigenvalues固有値こゆうち, but also because it guarantees orthogonality of eigenvectors in the right sense.

A general matrix may fail to be diagonalizable. Even if it is diagonalizable, its eigenvector basis need not be orthogonal. By contrast, every real symmetric matrix is diagonalizable by an orthonormal basis正規直交基底せいきちょっこうきてい. This is the finite-dimensional 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*

と表示ひょうじする。

2Terms and definitions

A real square matrix A is symmetric対称たいしょう if

AT=A.

A complex square matrix A is Hermitianエルミート if

A*=A.

A unitary matrixユニタリ行列ぎょうれつ satisfies

U*U=UU*=I.

It preserves the complex inner product, lengths, and orthogonality.

A normal matrix正規行列せいきぎょうれつ satisfies

A*A=AA*.

Hermitian and unitary matrices are basic examples of normal matrices.

Orthogonal diagonalization直交対角化ちょっこうたいかくか of a real matrix means writing

A=QDQT

with Q orthogonal and D diagonal. In the complex case one writes

A=UDU*

with U unitary.

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

3Plan

First, symmetry is interpreted through the inner-product identity

⟨Au,v⟩=⟨u,Av⟩.

Then this identity is used to prove that eigenvectors belonging to distinct eigenvalues are 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である。

4Intuitive explanation

A symmetric matrix is a transformation whose input and output contributions are balanced in the inner product. Such a transformation can be decomposed into stretches along mutually orthogonal principal directions. The principal directions are eigenvectors, and the stretch factors are eigenvalues.

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

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

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

⟨Au,v⟩=⟨u,Av⟩

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

λ⟨u,v⟩=μ⟨u,v⟩

である。したがって

(λ-μ)⟨u,v⟩=0

であり、λ≠μ であるため ⟨u,v⟩=0 である。よって u と v は直交ちょっこう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⟩

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

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

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

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

ここで Q は直交行列ちょっこうぎょうれつorthogonal matrix、D は実固有値じつこゆうちを並ならべた対角行列たいかくぎょうれつdiagonal matrixである。Q の列れつcolumnは A の正規直交せいきちょっこう固有ベクトル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できることは、内積構造ないせきこうぞうと両立りょうりつすることを意味いみする。

5Precise explanation

5.11. Eigenvectors for distinct eigenvalues are orthogonal

Let A be a real symmetric matrix, and suppose

Au=λu,Av=μv,λ≠μ.

Because A is symmetric,

⟨Au,v⟩=⟨u,Av⟩.

Substituting the eigenvalue equations gives

λ⟨u,v⟩=μ⟨u,v⟩.

Hence

(λ-μ)⟨u,v⟩=0.

Since λ≠μ, ⟨u,v⟩=0, so u and v are orthogonal.

The same argument works for complex Hermitian matrices. If Au=λu, then

λ⟨u,u⟩=⟨Au,u⟩=⟨u,Au⟩=λ_⟨u,u⟩.

Since u≠0 implies ⟨u,u⟩>0, we get λ=λ_, so Hermitian eigenvalues are real.

5.22. Spectral theorem for real symmetric matrices

For a real symmetric matrix A in finite dimensions,

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

Here Q is an orthogonal matrix直交行列ちょっこうぎょうれつ, D is diagonal with real eigenvalues, and the columns of Q are orthonormal eigenvectors of A.

In the complex Hermitian case,

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

where U is unitary and the diagonal entries of D are real.

5.33. Relation with normal matrices

Over finite-dimensional complex inner product spaces,

A*A=AA*⟺Aisunitarilydiagonalizable.

Hermitian matrices are the normal matrices whose eigenvalues are real. Real symmetric matrices are the real version of this self-adjoint situation.

6具体例ぐたいれい

A=(2112)

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

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

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

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

により

A=QDQT

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

6Concrete example

Let

A=(2112).

This matrix is symmetric. Take

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

The corresponding eigenvalues are 3 and 1. Therefore

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

and

A=QDQT.

7別べつの観点かんてん

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

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

7Another viewpoint

Ordinary diagonalization chooses a basis that simplifies a matrix. Orthogonal diagonalization makes the stronger claim that this basis can be orthonormal. Because orthogonal and unitary matrices preserve lengths and angles, the change of coordinates does not destroy the inner-product structure.

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

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

8Criteria

  • A real symmetric matrix is orthogonally diagonalizable.
  • A Hermitian matrix is unitarily diagonalizable.
  • Over complex finite-dimensional spaces, a matrix is unitarily diagonalizable exactly when it is normal.
  • Eigenvectors for distinct eigenvalues of a symmetric or Hermitian matrix are orthogonal.
  • Eigenvalues describe stretch factors along principal directions.

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

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

9Scope and limitations

Orthogonal diagonalization does not hold for every diagonalizable matrix. It requires symmetry compatible with the inner product. In complex linear algebra, the correct analogue of a real symmetric matrix is a Hermitian matrix, not merely a matrix satisfying AT=A.

10最終形さいしゅうけい

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

10Final forms

[PARSE ERROR: Undefined("Command(\"boxed\")")]AT=A⟹A=QDQT
[PARSE ERROR: Undefined("Command(\"boxed\")")]A*=A⟹A=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 matrix、D は対角行列たいかくぎょうれつdiagonal matrixである。

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

11In one sentence

A symmetric matrix対称行列たいしょうぎょうれつ decomposes space into mutually orthogonal eigendirections, each with its own stretch factor.

A real symmetric matrix対称行列たいしょうぎょうれつ has enough mutually orthogonal eigenvectors to form an orthonormal basis. Thus it can be written as A=QDQT, or equivalently QTAQ=D, where Q is orthogonal and D is diagonal.

The repeated-eigenvalue case is the main place where mistakes occur. Eigenvectors for different eigenvalues are automatically orthogonal, but inside one eigenspace one may still need Gram-Schmidt to choose an orthonormal basis. Orthogonal diagonalization preserves lengths and angles, which is why it is the right form for quadratic forms and principal-axis calculations.

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 の対角成分たいかくせいぶんが実数じっすうになる。

12Theorem: real and complex spectral theorems

In a finite-dimensional real inner product space, every real symmetric matrix A=AT has an orthonormal basis q1,…,qn of eigenvectors:

Aqi=λiqi.

Equivalently,

A=QDQT.

In a finite-dimensional complex inner product space, every normal matrix正規行列せいきぎょうれつ satisfying A*A=AA* is unitarily diagonalizable:

A=UDU*.

Hermitian matrices are the special case in which the diagonal entries of D are real.

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

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

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

なので、Ax∈q1⊥ である。つまり q1⊥ は A で不変ふへんである。この不変部分空間ふへんぶぶんくうかんに帰納法きのうほうを適用てきようすると、残のこりの正規直交せいきちょっこう固有ベクトルeigenvectorを得える。

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

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

13Proof idea: from one eigenvector to its orthogonal complement

For a real symmetric matrix, one first chooses a unit eigenvector q1. Then consider the orthogonal complement q1⊥. If x∈q1⊥, then

⟨Ax,q1⟩=⟨x,Aq1⟩=λ1⟨x,q1⟩=0.

Thus Ax∈q1⊥, so q1⊥ is invariant under A. Applying induction to this invariant subspace gives the remaining orthonormal eigenvectors.

The crucial point is that symmetry preserves orthogonal complements. This is exactly what fails for a general matrix.

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

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

A=(1102)

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

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

14Non-example: diagonalizable does not imply orthogonally diagonalizable

The matrix

A=(1102)

has distinct eigenvalues 1 and 2, so it is diagonalizable. But AT≠A, and it cannot be orthogonally diagonalized. If QTAQ=D with Q orthogonal and D diagonal, then A=QDQT, which is symmetric. That contradicts AT≠A.

This example separates diagonalizability from orthogonal diagonalizability.

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
タブを全て閉じる