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 とは、じつrealm×n 行列ぎょうれつmatrix A

A=UΣVT

表示ひょうじすることである。ここで Um×m 直交行列ちょっこうぎょうれつorthogonal matrixVn×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=Ax2[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=Ax2[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

orthogonalchangeofinputcoordinatesaxiswisestretching[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階数かいすうrankr とし、せい特異値とくいちsingular valueだけをのこすと

A=UrΣrVrT

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

この表示ひょうじでは、UrれつcolumnCol(A)正規直交基底せいきちょっこうきていorthonormal basisあたえ、VrれつcolumnRow(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-AkF=σk+12++σr2,A-Ak2=σ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σjAvi,Avj=1σiσjA*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-AkF=σk+12++σr2,A-Ak2=σ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 value3,1 である。

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

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

A=(1101)

かんがえる。このとき

ATA=(1112)

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

λ±=3±52

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

σ±=3±52

となる。A 自身じしん固有値こゆうちeigenvalue1 だけであるが、特異値とくいちsingular valueAながさをどれだけばすかを記録きろくする。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 valueATA固有値こゆうち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=Ax2[PARSE ERROR: Undefined("Command(\"ge\")")]0

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

ui=Aviσi

定義ていぎする。すると

ui,uj=1σiσjAvi,Avj=1σiσjA*Avi,vj=σi2σiσjvi,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=Ax2[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σjAvi,Avj=1σiσjA*Avi,vj=σi2σiσjvi,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 自身じしん固有値こゆうちeigenvalue1重複ちょうふくであるが、特異値とくいちsingular valueATA平方根へいほうこんとしてべつりょうはかっている。つまり 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 valueATA または 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
タブを全て閉じる