行列の分解

ウィキペディア

行列の分解

出典: フリー百科事典『ウィキペディア（Wikipedia）』 (2022/05/17 08:25 UTC 版)

線型代数学という数学の分野において、行列の分解（ぎょうれつのぶんかい、英: matrix decomposition, matrix factorization）とは、行列の行列の積への因数分解である．多くの異なった行列の分解があり、それぞれがある問題のために利用される。リー群の分解はこれらのより本質的な視点を与える。

例

数値解析において、異なる分解が効率的な行列アルゴリズムを実装するために用いられる。

例えば、線型方程式系（連立一次方程式）Ax = b を解くとき、行列 A はLU分解により分解できる。LU分解は行列を下三角行列 L と上三角行列 U の積に分解する。系 L(Ux) = b と Ux = L⁻¹b は、もとの系 Ax = b と比べて解くのに必要な加法や乗法が少ないが、浮動小数点のような不正確な計算ではかなりの桁数を必要とし得る。

同様に、QR分解は A を直交行列 Q と上三角行列 R の積 QR として表す。系 Q(Rx) = b は Rx = ^tQb = c によって解かれ、系 Rx = c は '後退代入' によって解かれる。必要な加法と乗法の回数はLU分解のときの約2倍だが、QR分解は数値的に安定なため不正確な計算においてより多くの桁数が必要とならない。

線型方程式系を解くことと関係する分解

LU分解

• 適用：正方行列 A
• 分解：A = LU,、ただし L は下三角行列で U は上三角行列
• 関連：LDU分解は A = LDU である、ただし L は下三角行列で対角線に 1 が並び、U は上三角行列で対角線に 1 が並び、D は対角行列である．
• 関連：LUP分解は A = LUP である、ただし L は下三角行列で、U は上三角行列で、P は置換行列である．
• 存在： LUP 分解は任意の正方行列 A に対して存在する。P が単位行列のとき、LUP分解はLU分解となる。LU分解が存在すればLDU分解も存在する。
• コメント：LUP 分解と LU 分解は n × n の線型方程式系 Ax = b を解く際に有用である。これらの分解はガウスの消去法の過程を行列の形にまとめたものである。行列 P はガウスの消去法の過程で行われる任意の行の交換を表す。ガウスの消去法で行の交換なしに行階段形になれば P = I であり、したがって LU 分解は存在する。

階数因数分解

• 適用：階数 r の m × n 行列 A
• 分解：A = CF、ただし C は m × r の full column rank matrix であり、F は r × n の full row rank matrix である。
• コメント：階数因数分解は A のムーア・ペンローズ擬逆行列を計算するのに使え、適用して線型系 Ax = b の全ての解を得ることができる。

コレスキー分解

• 適用：正方，対称，正定値行列 A
• 分解：A = ^tUU,、ただし U は上三角行列で対角成分は正
• コメント：実対称正定値行列のコレスキー分解は一意にできる（上三角行列Uの対角要素を正にとる）。
• コメント：実対称と複素対称にも一応適用ができるが、行列が正則でも分解が存在しない（破綻する）可能性がある。
• コメント：コレスキー分解は複素エルミート行列にも適用できる（その場合にはUの転置はUのエルミート転置に読み替える）。必ずしも分解が存在しないが、行列が正定値なら必ず分解できる。
• コメント：代替はLDL分解であり，平方根を引き出すことを避けられる．

QR分解

• 適用：m × n 行列 A
• 分解：A = QR,、ただし Q は m 次直交行列であり，R は m × n の上三角行列である。
• コメント：QR分解は方程式系 Ax = b を A の逆行列を求めずに解く別の方法を提供する。Q が直交行列であることは ^tQQ = I を意味するので、Ax = b は Rx = ^tQb と同値であり、後者は R が三角行列だから解きやすい。

固有値や関連概念に基づく分解

固有値分解

• スペクトル分解とも呼ばれる。
• 適用：相異なる固有ベクトルを持つ正方行列 A（固有値は同じものがあってもよい）。
• 分解：A = VDV⁻¹、 ただし D は A の固有値からなる対角行列で，V の行は対応する A の固有ベクトル。
• 存在：n × n 行列 A はつねに（重複を込めて） n 個の固有値を持ち、それらを並べて n × n の対角行列 D と対応する零でない行の行列 V を作ることができ、固有値方程式 AV = VD を満たす。n 個の固有ベクトルが相異なるとき、V は可逆であり、分解 A = VDV⁻¹ ができる。
• コメント：固有ベクトルの長さが 1 であるように正規化することがいつでもできる。A が実対称行列であれば、V はいつでも可逆であり正規化された列を持つようにできる。すると等式 V^tV = I が成り立つ、なぜならば各固有ベクトルは互いに直交するからである。したがって、分解は A = VD^tV となる。
• コメント：n 個の相異なる固有値をもつという条件は十分ではあるが必要ではない。必要十分条件は各固有値の幾何学的重複度がその代数的重複度に等しいことである。
• コメント：固有分解は線型常微分方程式系あるいは線型差分方程式系の解の理解に有用である。例えば，初期条件 x₀ = c から始まる差分方程式 x_{k + 1} = Ax_k は x_k = ^kAc によって解かれ、これは x_k = VD^kV⁻¹c に同値であり、ここで V と D は A の固有ベクトルと固有値から作られる行列である。D は対角行列だから、冪 D^k は単に各対角成分を k 乗するだけである。A は普通対角でないから A を k 乗するよりもはるかに容易である。

ジョルダン分解

ジョルダン標準形とジョルダン・シュヴァレー分解

• 適用：正方行列 A
• コメント：ジョルダン標準形は固有分解を固有値に重複があり対角化できない場合に一般化し、ジョルダン・シュヴァレー分解はこれを基底を選ばずに行う。

シューア分解

• 適用：正方行列 A
• コメント：この分解には2つのバージョンがある．複素シューア分解と実シューア分解である．複素行列は必ず複素シューア分解を持つ．
• 分解（複素バージョン）：A = UTU^∗、 ただし U がユニタリ行列で、U^∗ は U の共役転置で、T は A の固有値を対角線に持つ複素シューア標準形と呼ばれる上三角行列である。
• 分解（実バージョン）：A = VS^tV、ただし A, V, S は実数のみからなる行列である。V は直交行列で，^tV は V の転置で、S は実シューア標準形と呼ばれるブロック上三角行列である。S の対角にあるブロックのサイズは 1 × 1（実固有値を表す）かまたは 2 × 2（複素共役な固有値の対から導かれる）である。

QZ分解

• 別名：一般シューア分解
• 適用：正方行列 A と B
• コメント：この分解には2つのバージョンがある：複素と実．
• 分解（複素バージョン）：$A=QSZ^{*}$
  この節の加筆が望まれています。
  主に: examples and additional citations （2014年12月）

  quasimatrix（準行列）および cmatrix あるいは continuous matrix（連続行列）に対して SVD, QR, LU, コレスキー分解の類似がある。‘quasimatrix’は，行列のように、長方形の体系で、元は添え字付けられているが、1つの離散的な添え字が連続的な添え字に置き換えられる。同様に，‘cmatrix’は両方の添え字が連続である。cmatrix の例として，積分作用素の核を考えることができる。

  これらの分解は Fredholm (1903), Hilbert (1904), Schmidt (1907) による初期の研究に基づいている。これら独創的な論文の説明と英訳は Stewart (2011) を参照

  関連項目

  • Matrix splitting
  • Non-negative matrix factorization
  • Proper orthogonal decomposition

  参考文献

  • Choudhury, Dipa; Horn, Roger A. (April 1987). “A Complex Orthogonal-Symmetric Analog of the Polar Decomposition”. SIAM Journal on Algebraic Discrete Methods 8 (2). doi:10.1137/0608019. 
  • Fredholm, I. (1903), “Sur une classe d’´equations fonctionnelles” (French), Acta Mathematica 27: 365–390, doi:10.1007/bf02421317 
  • Hilbert, D. (1904), “Grundzüge einer allgemeinen Theorie der linearen Integralgleichungen” (German), Nachr. Königl. Ges. Gött 1904: 49–91 
  • Horn, Roger A.; Merino, Dennis I. (January 1995). “Contragredient equivalence: A canonical form and some applications”. Linear Algebra and its Applications 214. doi:10.1016/0024-3795(93)00056-6. 
  • Meyer, C. D. (2000), Matrix Analysis and Applied Linear Algebra, SIAM, ISBN 978-0-89871-454-8, http://www.matrixanalysis.com/ 
  • Schmidt, E. (1907), “Zur Theorie der linearen und nichtlinearen Integralgleichungen. I Teil. Entwicklung willkürlichen Funktionen nach System vorgeschriebener” (German), Mathematische Annalen 63: 433–476, doi:10.1007/bf01449770 
  • Simon, C.; Blume, L. (1994). Mathematics for Economists. Norton. ISBN 0-393-95733-0 
  • Stewart, G. W. (2011), Fredholm, Hilbert, Schmidt: three fundamental papers on integral equations, http://www.cs.umd.edu/~stewart/FHS.pdf 2015年1月6日閲覧。 
  • Townsend, A.; Trefethen, L. N. (2015), “Continuous analogues of matrix factorizations”, Proc. R. Soc. A 471 (2173), doi:10.1098/rspa.2014.0585 

  外部リンク

  • Online Matrix Calculator
  • Wolfram Alpha Matrix Decomposition Computation » LU and QR Decomposition
  • Springer Encyclopaedia of Mathematics » Matrix factorization
  • GraphLab en:GraphLab collaborative filtering library, large scale parallel implementation of matrix decomposition methods (in C++) for multicore.

ウィキペディア小見出し辞書

行列の分解

出典: フリー百科事典『ウィキペディア（Wikipedia）』 (2022/06/01 07:43 UTC 版)

「行列」の記事における「行列の分解」の解説

詳細は「行列の分解」を参照 行列を2つあるいは3つの行列の積に因数分解するには以下の方法が知られている。 LU分解 - 正方行列Aを下三角行列と上三角行列の積に分解。 A = LU コレスキー分解 - 正値対称行列(またはエルミート行列)Aを下三角行列と上三角行列の積に分解。 A = U*U QR分解 - (m,n)行列を直交行列(またはユニタリ行列)Qと上三角行列Rに分解 A = QR 固有値分解 - 特異値分解 - (m,n)行列を直交行列(またはユニタリ行列)U,Vと対角行列Dに分解 A = UDV*