markdown
Symmetric matrices and orthogonal diagonalizationmd 9e6fa92
lecture/math/linear-algebra/symmetric-matrices-and-orthogonal-diagonalization.lecture.n.md
Download PDF

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

mathlinear-algebraundergraduatelecture

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.

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.

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

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.

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 u0 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.

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.

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.

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.

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.

10Final forms

[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μ

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.

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.

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 xq1, then

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

Thus Axq1, 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

14Non-example: diagonalizable does not imply orthogonally diagonalizable

The matrix

A=(1102)

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

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