markdown
特異値分解の入口md 8af0916
lecture/math/linear-algebra/introduction-to-singular-value-decomposition.lecture.n.md
Download PDF

特異値分解とくいちぶんかいsingular value decompositionSVDの入口いりぐち

date2026-07-14document_iddoc_0c18ace358b22600fb65290b07a5b343description特異値分解を、長方行列を直交変換と軸方向の伸縮に分ける方法として導入し、階数・最小二乗法・擬似逆行列との関係を整理する講義である。prerequisites対称行列と直交対角化 / 階数の基本 / 最小二乗法の基本type講義content_typelecturestatusactiverelateddata/lecture/math/linear-algebra/linear-algebra-portal.lecture.n.md / data/lecture/math/linear-algebra/symmetric-matrices-and-orthogonal-diagonalization.lecture.n.md / data/lecture/math/linear-algebra/complex-inner-products-and-unitary-matrices.lecture.n.md / data/lecture/math/linear-algebra/rank-basics.lecture.n.md / data/lecture/math/linear-algebra/least-squares-basics.lecture.n.md / data/lecture/math/linear-algebra/pseudoinverse-basics.lecture.n.md / data/exercise/math/linear-algebra/eigenvalues-diagonalization-and-extensions.exercise.n.md
mathlinear-algebraundergraduatelecture

Introduction to singular value decomposition特異値分解とくいちぶんかいSVD

1導入どうにゅう

この講義こうぎで重要じゅうようなのは、特異値分解とくいちぶんかいsingular value decompositionSVDは、正方行列せいほうぎょうれつsquare matrixに限定げんていされない行列ぎょうれつmatrixを、直交変換ちょっこうへんかんと軸方向じくほうこうの伸縮しんしゅくへ分解ぶんかいする方法ほうほうであるということである。

固有値こゆうちeigenvalueと対角化たいかくかdiagonalizationは正方行列せいほうぎょうれつsquare matrixを主対象しゅたいしょうにする。一方いっぽう、実務じつむや応用おうようでは、縦横たてよこの大おおきさが異ことなる行列ぎょうれつmatrixが頻出ひんしゅつする。特異値分解とくいちぶんかいsingular value decompositionSVDは、そのような長方行列ちょうほうぎょうれつにも適用てきようできる標準的ひょうじゅんてきな分解ぶんかいである。

SVD は、一般いっぱんの行列ぎょうれつmatrixに対たいする「直交基底ちょっこうきていつきの最もっとも安定あんていした分解ぶんかい」と位置付いちづけられる。入力側にゅうりょくがわの正規直交基底せいきちょっこうきていorthonormal basis、出力側しゅつりょくがわの正規直交基底せいきちょっこうきていorthonormal basis、軸方向じくほうこうの非負ひふの伸縮率しんしゅくりつを分離ぶんりするため、階数かいすうrank、最小二乗法さいしょうにじょうほうleast squares method、擬似逆行列ぎじぎゃくぎょうれつpseudoinverse matrix、低階数近似ていかいすうきんじを同一どういつの枠組わくぐみで記述きじゅつできる。

1Introduction

The key point of this lecture is that singular value decomposition特異値分解とくいちぶんかいSVD decomposes a matrix, not necessarily square, into orthogonal changes of coordinates and stretching along coordinate axes.

eigenvalues固有値こゆうち and diagonalization mainly concern square matrices. In applications, rectangular matrices are common. SVD is the standard decomposition that works for every real or complex matrix.

SVD separates the input-side orthonormal basis基底きてい, the output-side orthonormal basis基底きてい, and the nonnegative scale factors. This lets rank階数かいすう, least squares, pseudoinverses, and low-rank階数かいすう approximation be described in one framework.

2用語ようごと定義ていぎ

特異値分解とくいちぶんかいSingular value decomposition とは、実じつrealの m×n 行列ぎょうれつmatrix A を

A=UΣVT

と表示ひょうじすることである。ここで U は m×m 直交行列ちょっこうぎょうれつorthogonal matrix、V は n×n 直交行列ちょっこうぎょうれつorthogonal matrix、Σ は対角方向たいかくほうこうに非負数ひふすうを並ならべた m×n 行列ぎょうれつmatrixである。

Σ の対角成分たいかくせいぶん

σ1[PARSE ERROR: Undefined("Command(\"ge\")")]σ2[PARSE ERROR: Undefined("Command(\"ge\")")]⋯[PARSE ERROR: Undefined("Command(\"ge\")")]0

を A の特異値とくいちSingular valuesという。複素行列ふくそぎょうれつでは

A=UΣV*

と表示ひょうじする。

2Terms and definitions

For a real m×n matrix A, a singular value decomposition特異値分解とくいちぶんかい is

A=UΣVT,

where U is an m×m orthogonal matrix, V is an n×n orthogonal matrix, and Σ is an m×n matrix with nonnegative diagonal entries.

The diagonal entries

σ1[PARSE ERROR: Undefined("Command(\"ge\")")]σ2[PARSE ERROR: Undefined("Command(\"ge\")")]⋯[PARSE ERROR: Undefined("Command(\"ge\")")]0

are the singular values特異値とくいち of A.

For complex matrices,

A=UΣV*,

with unitary U and V.

3方針ほうしん

方針ほうしんは、実行列じつぎょうれつでは ATA、複素行列ふくそぎょうれつでは A*A を解析かいせきすることである。ATA は実対称じつたいしょうかつ非負定値ひふていちであり、A*A はエルミート行列ぎょうれつmatrixかつ非負定値ひふていちであるため、直交ちょっこうorthogonalまたはユニタリに対角化たいかくかdiagonalizationできる。

なぜ A*A を見みるかというと、任意にんいの x に対たいして

⟨A*Ax,x⟩=⟨Ax,Ax⟩=∥Ax∥2[PARSE ERROR: Undefined("Command(\"ge\")")]0

となるからである。したがって A*A はエルミート行列ぎょうれつHermitian matrixかつ非負定値ひふていちであり、固有値こゆうちeigenvalueは 0 以上いじょうの実数じっすうreal numberになる。この固有値こゆうちeigenvalueの平方根へいほうこんが特異値とくいちsingular valueである。

ATA=VΛVT

複素ふくそcomplexの場合ばあいは

A*A=VΛV*

である。この固有値こゆうちeigenvalue λi は非負ひふであり、σi=λi が特異値とくいちsingular valueになる。

data/lecture/math/linear-algebra/symmetric-matrices-and-orthogonal-diagonalization.lecture.n.md data/lecture/math/linear-algebra/complex-inner-products-and-unitary-matrices.lecture.n.md

3Plan

Study ATA in the real case and A*A in the complex case. These matrices are symmetric or Hermitian and positive semidefinite. Therefore they can be orthogonally or unitarily diagonalized.

The reason for looking at A*A is that, for every x,

⟨A*Ax,x⟩=⟨Ax,Ax⟩=∥Ax∥2[PARSE ERROR: Undefined("Command(\"ge\")")]0.

Thus A*A is Hermitian and positive semidefinite, so its eigenvalues固有値こゆうち are nonnegative real numbers. In the real case we write

ATA=VΛVT,

and in the complex case

A*A=VΛV*.

The singular values特異値とくいち of A are the square roots of these nonnegative eigenvalues:

σi=λi(A*A).
data/lecture/math/linear-algebra/symmetric-matrices-and-orthogonal-diagonalization.lecture.n.md data/lecture/math/linear-algebra/complex-inner-products-and-unitary-matrices.lecture.n.md

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

特異値分解とくいちぶんかいsingular value decompositionSVDは、行列ぎょうれつmatrix A の作用さようを

直交変換⟶軸方向の伸縮⟶直交変換

に分解ぶんかいする。直交変換ちょっこうへんかんは長ながさと角度かくどを保存ほぞんし、Σ だけが長ながさを変更へんこうする。したがって特異値とくいちsingular valueは、行列ぎょうれつmatrixが各主方向かくしゅほうこうをどれだけ伸縮しんしゅくするかを表あらわす。

4Intuitive explanation

SVD decomposes the action of A as

orthogonalchangeofinputcoordinates⟶axiswisestretching[PARSE ERROR: Undefined("RBrace")]⟶orthogonalplacementintheoutputspace.

Orthogonal or unitary factors preserve lengths and angles. The matrix Σ is the only part that changes lengths. Thus singular values特異値とくいち measure how much A stretches each principal input direction.

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

5.11. ATA または A*A から右特異みぎとくいベクトルを得える

実行列じつぎょうれつでは ATA が実対称行列じつたいしょうぎょうれつである。したがって正規直交せいきちょっこう固有ベクトルeigenvector v1,[PARSE ERROR: Undefined("Command(\"dots\")")],vn を選択せんたくできる。

ATAvi=σi2vi

である。この vi が右特異みぎとくいベクトルである。複素行列ふくそぎょうれつでは A*A を用もちい、

A*Avi=σi2vi

と記述きじゅつする。この変更へんこうにより、定義ていぎの A=UΣV* と導出どうしゅつの記法きほうが整合せいごうする。

5.22. 左特異ひだりとくいベクトルを構成こうせいする

σi>0 のとき

ui=Aviσi

と定義ていぎする。この ui は互たがいに直交ちょっこうorthogonalし、長ながさ 1 である。これを補完ほかんして Rm の正規直交基底せいきちょっこうきていorthonormal basisにすると、U が得えられる。

左特異ひだりとくいベクトルは AAT からも特徴付とくちょうづけられる。実行列じつぎょうれつでは

AATui=σi2ui

であり、複素ふくそcomplex行列ぎょうれつmatrixでは

AA*ui=σi2ui

である。つまり、ATA または A*A は入力側にゅうりょくがわの主方向しゅほうこうを与あたえ、AAT または AA* は出力側しゅつりょくがわの主方向しゅほうこうを与あたえる。

5.33. 階数かいすうrankとの関係かんけい

非零特異値ひぜろとくいちの個数こすうは A の階数かいすうrankに等ひとしい。

[PARSE ERROR: Undefined("Command(\"boxed\")")]rank(A)=#{i∣σi>0}

小ちいさい特異値とくいちsingular valueは、行列ぎょうれつmatrixが情報じょうほうを強つよく圧縮あっしゅくする方向ほうこうを示しめす。

5.44. コンパクト SVD と基本部分空間きほんぶぶんくうかん

A の階数かいすうrankを r とし、正せいの特異値とくいちsingular valueだけを残のこすと

A=UrΣrVrT

と表示ひょうじできる。これをコンパクト SVDCompact SVDという。複素ふくそcomplexの場合ばあいは Vr* を用もちいる。

この表示ひょうじでは、Ur の列れつcolumnが Col(A) の正規直交基底せいきちょっこうきていorthonormal basisを与あたえ、Vr の列れつcolumnが Row(A) の正規直交基底せいきちょっこうきていorthonormal basisを与あたえる。さらに、零特異値ぜろとくいちに対応たいおうする右特異みぎとくいベクトルは ker(A) を、零特異値ぜろとくいちに対応たいおうする左特異ひだりとくいベクトルは ker(AT) または ker(A*) を記述きじゅつする。

5.55. 低階数近似ていかいすうきんじへの接続せつぞく

特異値とくいちsingular valueを大おおきい順じゅんに並ならべ、上位じょうい k 個こだけを残のこすと

Ak=∑i=1kσiuiviT

を得える。複素ふくそcomplexの場合ばあいは vi* を用もちいる。Eckart-Young の定理ていりにより、Ak は階数かいすうrank k 以下いかの行列ぎょうれつmatrixの中なかで A に最もっとも近ちかい近似きんじを与あたえる。

∥A-Ak∥F=σk+12+⋯+σr2,∥A-Ak∥2=σk+1

である。この事実じじつにより、SVD は分解ぶんかいの存在そんざいだけでなく、圧縮あっしゅくや雑音除去ざつおんじょきょの基礎きそにもなる。

5Precise explanation

5.11. Right singular vectors from ATA or A*A

In the real case, ATA is symmetric. Choose an orthonormal eigenbasis v1,…,vn:

ATAvi=σi2vi.

These vi are the right singular vectors. In the complex case,

A*Avi=σi2vi.

5.22. Constructing left singular vectors

For σi>0, define

ui=Aviσi.

Then the ui are orthonormal:

⟨ui,uj⟩=1σiσj⟨Avi,Avj⟩=1σiσj⟨A*Avi,vj⟩.

Using A*Avi=σi2vi and the orthonormality of the vi, this equals δij. Extend these ui to an orthonormal basis基底きてい of the output space to obtain U.

The left singular vectors are also eigenvectors固有こゆうベクトル of AAT in the real case:

AATui=σi2ui,

and of AA* in the complex case.

5.33. rank階数かいすう

The rank階数かいすう of A is the number of nonzero singular values特異値とくいち:

[PARSE ERROR: Undefined("Command(\"boxed\")")]rank(A)=#{i∣σi>0}.

Small singular values特異値とくいち indicate directions that are strongly compressed and often contribute to numerical instability.

5.44. Compact SVD and fundamental subspaces

If rank(A)=r, keeping only positive singular values特異値とくいち gives the compact SVD

A=UrΣrVrT

in the real case, and A=UrΣrVr* in the complex case.

The columns of Ur form an orthonormal basis基底きてい of Col(A). The columns of Vr form an orthonormal basis基底きてい of the row space行空間ぎょうくうかん. Right singular vectors for σi=0 span ker(A), and left singular vectors for zero singular values特異値とくいち span ker(AT) or ker(A*).

SpaceSVD description
Col(A)spanned by left singular vectors for σi>0
Row(A)spanned by right singular vectors for σi>0
ker(A)spanned by right singular vectors for σi=0
ker(AT) or ker(A*)spanned by left singular vectors for σi=0

5.55. Low-rank階数かいすう approximation

Keeping only the largest k singular values特異値とくいち gives

Ak=∑i=1kσiuiviT

in the real case, and vi* in the complex case. By the Eckart-Young theorem, Ak is the closest matrix to A among matrices of rank階数かいすう at most k:

∥A-Ak∥F=σk+12+⋯+σr2,∥A-Ak∥2=σk+1.

This supports compression and denoising, but discarding small singular values特異値とくいち also discards information. If a small-singular-value direction contains meaningful signal, truncation must be used carefully.

6具体例ぐたいていれい

A=(3001)

では、すでに軸方向じくほうこうの伸縮しんしゅくが分離ぶんりされている。したがって

U=I,Σ=(3001),V=I

とでき、特異値とくいちsingular valueは 3,1 である。

一般いっぱんの行列ぎょうれつmatrixでは、まず VT が入力空間にゅうりょくくうかんinput spaceを主方向しゅほうこうへ回転かいてんし、Σ が伸縮しんしゅくし、U が出力空間しゅつりょくくうかんへ配置はいちする。

非対角ひたいかくな例れいとして、

A=(1101)

を考かんがえる。このとき

ATA=(1112)

であり、その固有値こゆうちeigenvalueは

λ±=3±52

である。したがって特異値とくいちsingular valueは

σ±=3±52

となる。A 自身じしんの固有値こゆうちeigenvalueは 1 だけであるが、特異値とくいちsingular valueは A が長ながさをどれだけ伸のばすかを記録きろくする。SVD は固有値こゆうちeigenvalueの代替だいたいではなく、別べつの、より幾何的きかてきで安定あんていした量りょうquantityを測はかる。

6Concrete examples

For

A=(3001),

we may take

U=I,Σ=(3001),V=I.

The singular values特異値とくいち are 3 and 1.

For a general matrix, VT first rotates the input space into principal directions, Σ stretches those directions, and U places the result in the output space.

For a non-diagonal example,

A=(1101).

Then

ATA=(1112),

whose eigenvalues固有値こゆうち are

λ±=3±52.

Therefore the singular values特異値とくいち are

σ±=3±52.

Although A has only the eigenvalue固有値こゆうち 1, its singular values特異値とくいち record how lengths are stretched by A. SVD is not a replacement for eigenvalues固有値こゆうち; it measures a different and often more stable geometric quantity.

7別べつの観点かんてん

特異値分解とくいちぶんかいsingular value decompositionSVDは、最小二乗法さいしょうにじょうほうleast squares methodと擬似逆行列ぎじぎゃくぎょうれつpseudoinverse matrixの基礎きそである。A が正則せいそくでなくても、非零特異値ひぜろとくいちの方向ほうこうだけを反転はんてんすることで、最小さいしょうノルム解かいを定義ていぎできる。

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

  • 特異値とくいちsingular valueは ATA の固有値こゆうちeigenvalueの平方根へいほうこんである。
  • 非零特異値ひぜろとくいちの個数こすうは階数かいすうrankである。
  • 特異値分解とくいちぶんかいsingular value decompositionSVDは長方行列ちょうほうぎょうれつにも適用てきようできる。
  • 小ちいさい特異値とくいちsingular valueは情報じょうほうの圧縮あっしゅくや数値的不安定性すうちてきふあんていせいと関係かんけいする。
  • コンパクト SVD は、列空間れつくうかんcolumn space・行空間ぎょうくうかんrow space・核かくkernelを同時どうじに記述きじゅつする。

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

特異値分解とくいちぶんかいsingular value decompositionSVDは任意にんいの実行列じつぎょうれつ・複素行列ふくそぎょうれつに成立せいりつする。固有値分解こゆうちぶんかいとは異ことなり、正方行列せいほうぎょうれつsquare matrixや対角化可能性たいかくかかのうせいを前提ぜんていにしない。

10存在定理そんざいていりと証明しょうめいの骨格こっかく

任意にんいの m×n 実行列じつぎょうれつ A について、直交行列ちょっこうぎょうれつorthogonal matrix U,V と非負ひふの対角成分たいかくせいぶんをもつ m×n 行列ぎょうれつ Σ が存在そんざいして

A=UΣVT

と書かける。複素行列ふくそぎょうれつcomplex matrixでは

A=UΣV*

である。

証明しょうめいの骨格こっかくは ATA または A*A のスペクトル定理ていりspectral theoremである。A*A はエルミート行列ぎょうれつHermitian matrixであり、

⟨A*Ax,x⟩=∥Ax∥2[PARSE ERROR: Undefined("Command(\"ge\")")]0

なので非負定値ひふていちである。したがって正規直交せいきちょっこう固有ベクトルeigenvector vi と非負ひふの固有値こゆうちeigenvalue λi を持もつ。σi=λi と置おき、σi>0 について

ui=Aviσi

と定義ていぎする。すると

⟨ui,uj⟩=1σiσj⟨Avi,Avj⟩=1σiσj⟨A*Avi,vj⟩=σi2σiσj⟨vi,vj⟩

より ui は互たがいに直交ちょっこうorthogonalし、長ながさ 1 になる。残のこりの ui を補完ほかんすれば U が得えられる。この構成こうせいが SVD の存在そんざいを説明せつめいする。

data/lecture/math/linear-algebra/symmetric-matrices-and-orthogonal-diagonalization.lecture.n.md

7Another viewpoint

Singular value decompositionSVD特異値分解とくいちぶんかい is the basis of the least squares method最小二乗法さいしょうにじょうほう and the pseudoinverse matrix擬似逆行列ぎじぎゃくぎょうれつ. Even when A is not invertible, one can invert only the directions with nonzero singular values and thereby define the minimum-norm solution.

8Criteria

  • Singular values特異値とくいち are the square roots of the eigenvalues固有値こゆうち of ATA or A*A.
  • The number of nonzero singular values is the rank階数かいすう.
  • SVD applies to rectangular matrices as well as square matrices.
  • Small singular values特異値とくいち are related to compression and numerical instability.
  • Compact SVD simultaneously describes the column space列空間れつくうかん, row space行空間ぎょうくうかん, and kernel核かく.

9Scope

SVD holds for every real or complex matrix. Unlike eigendecomposition, it does not require the matrix to be square or diagonalizable.

10Existence theorem and proof outline

For every real m×n matrix A, there exist orthogonal matrices直交行列ちょっこうぎょうれつ U,V and an m×n matrix Σ with nonnegative diagonal entries such that

A=UΣVT.

For a complex matrix,

A=UΣV*.

The proof is based on the spectral theorem for ATA or A*A. The matrix A*A is Hermitianエルミート and

⟨A*Ax,x⟩=∥Ax∥2[PARSE ERROR: Undefined("Command(\"ge\")")]0,

so it is positive semidefinite. Hence it has an orthonormal eigenbasis vi with nonnegative eigenvalues固有値こゆうち λi. Set σi=λi, and for σi>0 define

ui=Aviσi.

Then

⟨ui,uj⟩=1σiσj⟨Avi,Avj⟩=1σiσj⟨A*Avi,vj⟩=σi2σiσj⟨vi,vj⟩,

so the ui are mutually orthogonal直交ちょっこう and have length 1. Completing the remaining ui gives U. This construction explains the existence of the SVD.

data/lecture/math/linear-algebra/symmetric-matrices-and-orthogonal-diagonalization.lecture.n.md

11基本部分空間きほんぶぶんくうかんの対応たいおう

SVD は、4 つの基本部分空間きほんぶぶんくうかんを同時どうじに整理せいりする。

空間くうかんSVD での記述きじゅつ
Col(A)正せいの特異値とくいちsingular valueに対応たいおうする左特異ひだりとくいベクトル ui が張はる
Row(A)正せいの特異値とくいちsingular valueに対応たいおうする右特異みぎとくいベクトル vi が張はる
ker(A)σi=0 に対応たいおうする vi が張はる
ker(AT) または ker(A*)σi=0 に対応たいおうする ui が張はる

このため SVD は、階数かいすうrankだけでなく、解かいが失うしなわれる方向ほうこう、列空間れつくうかんへの射影しゃえいprojection、最小二乗解さいしょうにじょうかいを同時どうじに記述きじゅつできる。

11Correspondence of fundamental subspaces

SVD organizes the four fundamental subspaces基本部分空間きほんぶぶんくうかん at the same time.

Space空間くうかんDescription by SVD
Col(A)spanned by the left singular vectors特異とくいベクトル ui corresponding to positive singular values特異値とくいち
Row(A)spanned by the right singular vectors特異とくいベクトル vi corresponding to positive singular values特異値とくいち
ker(A)spanned by the vi corresponding to σi=0
ker(AT) or ker(A*)spanned by the ui corresponding to σi=0

Therefore SVD describes not only the rank階数かいすう, but also the directions in which solutions are lost, the projection射影しゃえい onto the column space列空間れつくうかん, and least-squares solutions最小二乗解さいしょうにじょうかい.

12別べつの見方みかた

SVD は、最小二乗法さいしょうにじょうほうleast squares methodと擬似逆行列ぎじぎゃくぎょうれつpseudoinverse matrixの基礎きそでもある。A が可逆かぎゃくinvertibleでなくても、正せいの特異値とくいちsingular valueをもつ特異方向とくいほうこうsingular directionは逆向ぎゃくむきに戻もどせる一方いっぽう、σi=0 の方向ほうこうは失うしなわれた情報じょうほうとして扱あつかう。

12Another viewpoint

SVD is the foundation of the least squares method最小二乗法さいしょうにじょうほう and the pseudoinverse matrix擬似逆行列ぎじぎゃくぎょうれつ. Even when A is not invertible可逆かぎゃく, nonzero singular directions can be inverted while zero singular directions are treated as lost information.

13計算例けいさんれい:非対角ひたいかくな 2×2 行列ぎょうれつ

A=(1101)

を考かんがえる。このとき

ATA=(1112)

であり、固有値こゆうちeigenvalueは

λ±=3±52

である。したがって特異値とくいちsingular valueは

σ±=3±52

である。ここで A 自身じしんの固有値こゆうちeigenvalueは 1 の重複ちょうふくであるが、特異値とくいちsingular valueは ATA の平方根へいほうこんとして別べつの量りょうを測はかっている。つまり SVD は「固有値こゆうちeigenvalueの代替だいたい」ではなく、入力方向にゅうりょくほうこうの長ながさがどれだけ伸のびるかを測はかる分解ぶんかいである。

13Calculation example: a non-diagonal 2×2 matrix

Consider

A=(1101).

Then

ATA=(1112),

and the eigenvalues固有値こゆうち are

λ±=3±52.

Therefore the singular values特異値とくいち are

σ±=3±52.

Here the eigenvalue固有値こゆうち of A itself is 1 with multiplicity two, but the singular values特異値とくいち measure a different quantity as square roots of the eigenvalues of ATA. In other words, SVD is not a substitute for eigenvalues固有値こゆうち; it is a decomposition分解ぶんかい that measures how much lengths in input directions are stretched.

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

  • 特異値とくいちsingular valueは ATA または A*A の固有値こゆうちeigenvalueの平方根へいほうこんである。
  • 正せいの特異値とくいちsingular valueの個数こすうが階数かいすうrankである。
  • SVD は長方行列ちょうほうぎょうれつにも適用てきようできる。
  • 小ちいさい特異値とくいちsingular valueは、圧縮あっしゅくや数値的すうちてきな不安定性ふあんていせいを示しめすことがある。
  • コンパクト SVD は、列空間れつくうかんcolumn space、行空間ぎょうくうかんrow space、核かくを同時どうじに記述きじゅつする。

14Criteria

  • singular values特異値とくいち are square roots of eigenvalues固有値こゆうち of ATA or A*A.
  • The number of nonzero singular values特異値とくいち is the rank階数かいすう.
  • SVD applies to rectangular matrices.
  • Small singular values特異値とくいち indicate compression and possible numerical instability.
  • Compact SVD describes column space列空間れつくうかん, row space行空間ぎょうくうかん, and kernels simultaneously.

15低階数近似ていかいすうきんじの意味いみと注意ちゅうい

Eckart-Young の定理ていりは、上位じょうい k 個この特異値とくいちsingular valueだけを残のこした

Ak=∑i=1kσiuivi*

が、階数かいすう k 以下いかの行列ぎょうれつmatrixの中なかで A に最もっとも近ちかいことを主張しゅちょうする。ただし、小ちいさい特異値とくいちsingular valueを切きることは情報じょうほうを捨すてることでもある。ノイズ除去じょきょや圧縮あっしゅくでは有効ゆうこうだが、小ちいさい特異値とくいちsingular valueの方向ほうこうが意味いみのある信号しんごうを含ふくむ場合ばあいには注意ちゅういが必要ひつようである。

15Meaning and cautions of low-rank approximation

The Eckart-Young theorem定理ていり states that

Ak=∑i=1kσiuivi*

obtained by keeping only the top k singular values特異値とくいち is the closest matrix to A among all matrices of rank階数かいすう at most k. However, cutting off small singular values特異値とくいち also means throwing away information. This is useful for denoising and compression, but caution is necessary when the directions of small singular values特異値とくいち contain meaningful signal.

16適用範囲てきようはんいと限界げんかい

SVD は、実じつまたは複素ふくそのすべての行列ぎょうれつmatrixに存在そんざいする。固有値分解こゆうちぶんかいeigenvalue decompositionと違ちがい、正方せいほうであることや対角化可能たいかくかかのうdiagonalizableであることを仮定かていしない。ただし、特異値とくいちsingular valueだけでは符号ふごう、位相いそう、固有こゆうベクトルeigenvectorの振ふる舞まいは捉とらえられない。特異値とくいちsingular valueが記述きじゅつするのは長ながさの伸のびである。

16Scope and limitations

SVD exists for every real or complex matrix. Unlike eigenvalue decomposition固有値分解こゆうちぶんかい, it does not assume squareness or diagonalizability. However, the singular values特異値とくいち alone do not capture sign, phase, or eigenvector固有こゆうベクトル behavior; they describe length stretching.

17最終形さいしゅうけい

[PARSE ERROR: Undefined("Command(\"boxed\")")]A=UΣVT
[PARSE ERROR: Undefined("Command(\"boxed\")")]A=UΣV*(complexcase)
[PARSE ERROR: Undefined("Command(\"boxed\")")]σi=λi(ATA)
[PARSE ERROR: Undefined("Command(\"boxed\")")]rank(A)=#{i:σi>0}
[PARSE ERROR: Undefined("Command(\"boxed\")")]Ak=∑i=1kσiuiviT

17Final forms

[PARSE ERROR: Undefined("Command(\"boxed\")")]A=UΣVT
[PARSE ERROR: Undefined("Command(\"boxed\")")]A=UΣV*(complexcase)
[PARSE ERROR: Undefined("Command(\"boxed\")")]σi=λi(ATA)
[PARSE ERROR: Undefined("Command(\"boxed\")")]rank(A)=#{i:σi>0}
[PARSE ERROR: Undefined("Command(\"boxed\")")]Ak=∑i=1kσiuiviT

18一言ひとことでいうと

  • 特異値分解とくいちぶんかいsingular value decompositionSVDは、任意にんいの行列ぎょうれつmatrixを直交変換ちょっこうへんかんと軸方向じくほうこうの伸縮しんしゅくへ分解ぶんかいする標準形ひょうじゅんけいである。

18In one sentence

Singular value decomposition特異値分解とくいちぶんかいSVD is the standard form that decomposes any matrix into orthogonal coordinate changes and axis-wise stretching.

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