Matrix decomposition

<h2 id="example">Example</h2>
In <a href="/facts/Numerical_analysis/QT9cYa9z">numerical analysis</a>, different decompositions are used to implement efficient matrix <a href="/facts/Algorithm/fnl5NmRt">algorithms</a>.
For example, when solving a <a href="/facts/System_of_linear_equations/fk9CFdgp">system of linear equations</a> 
 
 
 
 A
 
 x
 
 =
 
 b
 
 
 
 {\displaystyle A\mathbf {x} =\mathbf {b} }
 
, the matrix A can be decomposed via the <a href="/facts/LU_decomposition/M13Ew8RD">LU decomposition</a>. The LU decomposition factorizes a matrix into a <a href="/facts/Lower_triangular_matrix/qjXSFMM3">lower triangular matrix</a> L and an <a href="/facts/Upper_triangular_matrix/qjXSFMM3">upper triangular matrix</a> U. The systems 
 
 
 
 L
 (
 U
 
 x
 
 )
 =
 
 b
 
 
 
 {\displaystyle L(U\mathbf {x} )=\mathbf {b} }
 
 and 
 
 
 
 U
 
 x
 
 =
 
 L
 
 −
 1
 
 
 
 b
 
 
 
 {\displaystyle U\mathbf {x} =L^{-1}\mathbf {b} }
 
 require fewer additions and multiplications to solve, compared with the original system 
 
 
 
 A
 
 x
 
 =
 
 b
 
 
 
 {\displaystyle A\mathbf {x} =\mathbf {b} }
 
, though one might require significantly more digits in inexact arithmetic such as <a href="/facts/Floating_point/eIckahxe">floating point</a>.
Similarly, the <a href="/facts/QR_decomposition/NswfKLGY">QR decomposition</a> expresses A as QR with Q an <a href="/facts/Orthogonal_matrix/WDXBXYB5">orthogonal matrix</a> and R an upper triangular matrix. The system Q(Rx) = b is solved by Rx = QTb = c, and the system Rx = c is solved by '<a href="/facts/Triangular_matrix/qjXSFMM3">back substitution</a>'. The number of additions and multiplications required is about twice that of using the LU solver, but no more digits are required in inexact arithmetic because the QR decomposition is <a href="/facts/Numerically_stable/jRDdo8zV">numerically stable</a>.

<h2 id="decompositions-related-to-solving-systems-of-linear-equations">Decompositions related to solving systems of linear equations</h2>
<h3>LU decomposition</h3>
Main article: <a href="/facts/LU_decomposition/M13Ew8RD">LU decomposition</a>
<ul><li>Traditionally applicable to: <a href="/facts/Square_matrix/5TP9mNNn">square matrix</a> A, although rectangular matrices can be applicable.<a class="footnote-ref" id="fnref:1" href="#fn:1">1</a><a class="footnote-ref" id="fnref:2" href="#fn:2">2</a></li>
<li>Decomposition: 
 
 
 
 A
 =
 L
 U
 
 
 {\displaystyle A=LU}
 
, where L is <a href="/facts/Triangular_matrix/qjXSFMM3">lower triangular</a> and U is <a href="/facts/Triangular_matrix/qjXSFMM3">upper triangular</a>.</li>
<li>Related: the <a href="/facts/LDU_decomposition/M13Ew8RD">LDU decomposition</a> is 
 
 
 
 A
 =
 L
 D
 U
 
 
 {\displaystyle A=LDU}
 
, where L is <a href="/facts/Triangular_matrix/qjXSFMM3">lower triangular</a> with ones on the diagonal, U is <a href="/facts/Triangular_matrix/qjXSFMM3">upper triangular</a> with ones on the diagonal, and D is a <a href="/facts/Diagonal_matrix/tjIAHrE6">diagonal matrix</a>.</li>
<li>Related: the <a href="/facts/LUP_decomposition/M13Ew8RD">LUP decomposition</a> is 
 
 
 
 P
 A
 =
 L
 U
 
 
 {\displaystyle PA=LU}
 
, where L is <a href="/facts/Triangular_matrix/qjXSFMM3">lower triangular</a>, U is <a href="/facts/Triangular_matrix/qjXSFMM3">upper triangular</a>, and P is a <a href="/facts/Permutation_matrix/xC1ran5s">permutation matrix</a>.</li>
<li>Existence: An LUP decomposition exists for any square matrix A. When P is an <a href="/facts/Identity_matrix/Glbcf6ol">identity matrix</a>, the LUP decomposition reduces to the LU decomposition.</li>
<li>Comments: The LUP and LU decompositions are useful in solving an n-by-n system of linear equations 
 
 
 
 A
 
 x
 
 =
 
 b
 
 
 
 {\displaystyle A\mathbf {x} =\mathbf {b} }
 
. These decompositions summarize the process of <a href="/facts/Gaussian_elimination/23Rvf4ht">Gaussian elimination</a> in matrix form. Matrix P represents any row interchanges carried out in the process of Gaussian elimination. If Gaussian elimination produces the <a href="/facts/Row_echelon_form/U9LQHBYq">row echelon form</a> without requiring any row interchanges, then P = I, so an LU decomposition exists.</li></ul>
<h3>LU reduction</h3>
Main article: <a href="/facts/LU_reduction/BvBlNfYR">LU reduction</a>
<h3>Block LU decomposition</h3>
Main article: <a href="/facts/Block_LU_decomposition/XnbqvSM1">Block LU decomposition</a>
<h3>Rank factorization</h3>
Main article: <a href="/facts/Rank_factorization/7pbe4yfK">Rank factorization</a>
<ul><li>Applicable to: m-by-n matrix A of rank r</li>
<li>Decomposition: 
 
 
 
 A
 =
 C
 F
 
 
 {\displaystyle A=CF}
 
 where C is an m-by-r full column rank matrix and F is an r-by-n full row rank matrix</li>
<li>Comment: The rank factorization can be used to <a href="/facts/Moore%25E2%2580%2593Penrose_pseudoinverse/CU73tntq">compute the Moore–Penrose pseudoinverse</a> of A,<a class="footnote-ref" id="fnref:3" href="#fn:3">3</a> which one can apply to <a href="/facts/Moore%25E2%2580%2593Penrose_pseudoinverse/CU73tntq">obtain all solutions of the linear system</a> 
 
 
 
 A
 
 x
 
 =
 
 b
 
 
 
 {\displaystyle A\mathbf {x} =\mathbf {b} }
 
.</li></ul>
<h3>Cholesky decomposition</h3>
Main article: <a href="/facts/Cholesky_decomposition/Q3IjkuTs">Cholesky decomposition</a>
<ul><li>Applicable to: <a href="/facts/Square_matrix/5TP9mNNn">square</a>, <a href="/facts/Symmetric_matrix/08CGwjNl">hermitian</a>, <a href="/facts/Positive-definite_matrix/Hbr6AuPS">positive definite</a> matrix 
 
 
 
 A
 
 
 {\displaystyle A}
 
</li>
<li>Decomposition: 
 
 
 
 A
 =
 
 U
 
 ∗
 
 
 U
 
 
 {\displaystyle A=U^{*}U}
 
, where 
 
 
 
 U
 
 
 {\displaystyle U}
 
 is upper triangular with real positive diagonal entries</li>
<li>Comment: if the matrix 
 
 
 
 A
 
 
 {\displaystyle A}
 
 is Hermitian and positive semi-definite, then it has a decomposition of the form 
 
 
 
 A
 =
 
 U
 
 ∗
 
 
 U
 
 
 {\displaystyle A=U^{*}U}
 
 if the diagonal entries of 
 
 
 
 U
 
 
 {\displaystyle U}
 
 are allowed to be zero</li>
<li>Uniqueness: for positive definite matrices Cholesky decomposition is unique. However, it is not unique in the positive semi-definite case.</li>
<li>Comment: if 
 
 
 
 A
 
 
 {\displaystyle A}
 
 is real and symmetric, 
 
 
 
 U
 
 
 {\displaystyle U}
 
 has all real elements</li>
<li>Comment: An alternative is the <a href="/facts/LDL_decomposition/Q3IjkuTs">LDL decomposition</a>, which can avoid extracting square roots.</li></ul>
<h3>QR decomposition</h3>
Main article: <a href="/facts/QR_decomposition/NswfKLGY">QR decomposition</a>
<ul><li>Applicable to: m-by-n matrix A with linearly independent columns</li>
<li>Decomposition: 
 
 
 
 A
 =
 Q
 R
 
 
 {\displaystyle A=QR}
 
 where 
 
 
 
 Q
 
 
 {\displaystyle Q}
 
 is a <a href="/facts/Unitary_matrix/8fsxv3KD">unitary matrix</a> of size m-by-m, and 
 
 
 
 R
 
 
 {\displaystyle R}
 
 is an <a href="/facts/Triangular_matrix/qjXSFMM3">upper triangular</a> matrix of size m-by-n</li>
<li>Uniqueness: In general it is not unique, but if 
 
 
 
 A
 
 
 {\displaystyle A}
 
 is of full <a href="/facts/Matrix_rank/Yp9b5LK3">rank</a>, then there exists a single 
 
 
 
 R
 
 
 {\displaystyle R}
 
 that has all positive diagonal elements. If 
 
 
 
 A
 
 
 {\displaystyle A}
 
 is square, also 
 
 
 
 Q
 
 
 {\displaystyle Q}
 
 is unique.</li>
<li>Comment: The QR decomposition provides an effective way to solve the system of equations 
 
 
 
 A
 
 x
 
 =
 
 b
 
 
 
 {\displaystyle A\mathbf {x} =\mathbf {b} }
 
. The fact that 
 
 
 
 Q
 
 
 {\displaystyle Q}
 
 is <a href="/facts/Orthogonal_matrix/WDXBXYB5">orthogonal</a> means that 
 
 
 
 
 Q
 
 
 T
 
 
 
 Q
 =
 I
 
 
 {\displaystyle Q^{\mathrm {T} }Q=I}
 
, so that 
 
 
 
 A
 
 x
 
 =
 
 b
 
 
 
 {\displaystyle A\mathbf {x} =\mathbf {b} }
 
 is equivalent to 
 
 
 
 R
 
 x
 
 =
 
 Q
 
 
 T
 
 
 
 
 b
 
 
 
 {\displaystyle R\mathbf {x} =Q^{\mathsf {T}}\mathbf {b} }
 
, which is very easy to solve since 
 
 
 
 R
 
 
 {\displaystyle R}
 
 is <a href="/facts/Triangular_matrix/qjXSFMM3">triangular</a>.</li></ul>
<h3>RRQR factorization</h3>
Main article: <a href="/facts/RRQR_factorization/mvHv2QKE">RRQR factorization</a>
<h3>Interpolative decomposition</h3>
Main article: <a href="/facts/Interpolative_decomposition/hKlwbCf9">Interpolative decomposition</a>
<h2 id="decompositions-based-on-eigenvalues-and-related-concepts">Decompositions based on eigenvalues and related concepts</h2>
<h3>Eigendecomposition</h3>
Main article: <a href="/facts/Eigendecomposition_(matrix)/a2nfF7hJ">Eigendecomposition (matrix)</a>
<ul><li>Also called <a href="/facts/Spectral_decomposition_(Matrix)/a2nfF7hJ">spectral decomposition</a>.</li>
<li>Applicable to: <a href="/facts/Square_matrix/5TP9mNNn">square matrix</a> A with linearly independent eigenvectors (not necessarily distinct eigenvalues).</li>
<li>Decomposition: 
 
 
 
 A
 =
 V
 D
 
 V
 
 −
 1
 
 
 
 
 {\displaystyle A=VDV^{-1}}
 
, where D is a <a href="/facts/Diagonal_matrix/tjIAHrE6">diagonal matrix</a> formed from the <a href="/facts/Eigenvalue/8TjEoT8u">eigenvalues</a> of A, and the columns of V are the corresponding <a href="/facts/Eigenvector/8TjEoT8u">eigenvectors</a> of A.</li>
<li>Existence: An n-by-n matrix A always has n (complex) eigenvalues, which can be ordered (in more than one way) to form an n-by-n diagonal matrix D and a corresponding matrix of nonzero columns V that satisfies the <a href="/facts/Eigenvalue_equation/qAwPXmpR">eigenvalue equation</a> 
 
 
 
 A
 V
 =
 V
 D
 
 
 {\displaystyle AV=VD}
 
. 
 
 
 
 V
 
 
 {\displaystyle V}
 
 is invertible if and only if the n eigenvectors are <a href="/facts/Linear_independence/8JrHfIIa">linearly independent</a> (that is, each eigenvalue has <a href="/facts/Geometric_multiplicity/8TjEoT8u">geometric multiplicity</a> equal to its <a href="/facts/Algebraic_multiplicity/8TjEoT8u">algebraic multiplicity</a>). A sufficient (but not necessary) condition for this to happen is that all the eigenvalues are different (in this case geometric and algebraic multiplicity are equal to 1)</li>
<li>Comment: One can always normalize the eigenvectors to have length one (see the definition of the eigenvalue equation)</li>
<li>Comment: Every <a href="/facts/Normal_matrix/WVKxrEK9">normal matrix</a> A (that is, matrix for which 
 
 
 
 A
 
 A
 
 ∗
 
 
 =
 
 A
 
 ∗
 
 
 A
 
 
 {\displaystyle AA^{*}=A^{*}A}
 
, where 
 
 
 
 
 A
 
 ∗
 
 
 
 
 {\displaystyle A^{*}}
 
 is a <a href="/facts/Conjugate_transpose/7eBWtqVG">conjugate transpose</a>) can be eigendecomposed. For a <a href="/facts/Normal_matrix/WVKxrEK9">normal matrix</a> A (and only for a normal matrix), the eigenvectors can also be made orthonormal (
 
 
 
 V
 
 V
 
 ∗
 
 
 =
 I
 
 
 {\displaystyle VV^{*}=I}
 
) and the eigendecomposition reads as 
 
 
 
 A
 =
 V
 D
 
 V
 
 ∗
 
 
 
 
 {\displaystyle A=VDV^{*}}
 
. In particular all <a href="/facts/Unitary_matrix/8fsxv3KD">unitary</a>, <a href="/facts/Hermitian_matrix/qArPblYo">Hermitian</a>, or <a href="/facts/Skew-Hermitian_matrix/wxjjnvCC">skew-Hermitian</a> (in the real-valued case, all <a href="/facts/Orthogonal_matrix/WDXBXYB5">orthogonal</a>, <a href="/facts/Symmetric_matrix/08CGwjNl">symmetric</a>, or <a href="/facts/Skew-symmetric_matrix/nw5CQeLN">skew-symmetric</a>, respectively) matrices are normal and therefore possess this property.</li>
<li>Comment: For any real <a href="/facts/Symmetric_matrix/08CGwjNl">symmetric matrix</a> A, the eigendecomposition always exists and can be written as 
 
 
 
 A
 =
 V
 D
 
 V
 
 
 T
 
 
 
 
 
 {\displaystyle A=VDV^{\mathsf {T}}}
 
, where both D and V are real-valued.</li>
<li>Comment: The eigendecomposition is useful for understanding the solution of a system of linear ordinary differential equations or linear difference equations. For example, the difference equation 
 
 
 
 
 x
 
 t
 +
 1
 
 
 =
 A
 
 x
 
 t
 
 
 
 
 {\displaystyle x_{t+1}=Ax_{t}}
 
 starting from the initial condition 
 
 
 
 
 x
 
 0
 
 
 =
 c
 
 
 {\displaystyle x_{0}=c}
 
 is solved by 
 
 
 
 
 x
 
 t
 
 
 =
 
 A
 
 t
 
 
 c
 
 
 {\displaystyle x_{t}=A^{t}c}
 
, which is equivalent to 
 
 
 
 
 x
 
 t
 
 
 =
 V
 
 D
 
 t
 
 
 
 V
 
 −
 1
 
 
 c
 
 
 {\displaystyle x_{t}=VD^{t}V^{-1}c}
 
, where V and D are the matrices formed from the eigenvectors and eigenvalues of A. Since D is diagonal, raising it to power 
 
 
 
 
 D
 
 t
 
 
 
 
 {\displaystyle D^{t}}
 
, just involves raising each element on the diagonal to the power t. This is much easier to do and understand than raising A to power t, since A is usually not diagonal.</li></ul>
<h3>Jordan decomposition</h3>
The <a href="/facts/Jordan_normal_form/M7ezYbnl">Jordan normal form</a> and the <a href="/facts/Jordan%25E2%2580%2593Chevalley_decomposition/w4PN88Xo">Jordan–Chevalley decomposition</a>

<ul><li>Applicable to: <a href="/facts/Square_matrix/5TP9mNNn">square matrix</a> A</li>
<li>Comment: the Jordan normal form generalizes the eigendecomposition to cases where there are repeated eigenvalues and cannot be diagonalized, the Jordan–Chevalley decomposition does this without choosing a basis.</li></ul>
<h3>Schur decomposition</h3>
Main article: <a href="/facts/Schur_decomposition/EYNrqHmh">Schur decomposition</a>
<ul><li>Applicable to: <a href="/facts/Square_matrix/5TP9mNNn">square matrix</a> A</li>
<li>Decomposition (complex version): 
 
 
 
 A
 =
 U
 T
 
 U
 
 ∗
 
 
 
 
 {\displaystyle A=UTU^{*}}
 
, where U is a <a href="/facts/Unitary_matrix/8fsxv3KD">unitary matrix</a>, 
 
 
 
 
 U
 
 ∗
 
 
 
 
 {\displaystyle U^{*}}
 
 is the <a href="/facts/Conjugate_transpose/7eBWtqVG">conjugate transpose</a> of U, and T is an <a href="/facts/Upper_triangular/qjXSFMM3">upper triangular</a> matrix called the complex <a href="/facts/Schur_form/EYNrqHmh">Schur form</a> which has the <a href="/facts/Eigenvalue/8TjEoT8u">eigenvalues</a> of A along its diagonal.</li>
<li>Comment: if A is a <a href="/facts/Normal_matrix/WVKxrEK9">normal matrix</a>, then T is diagonal and the Schur decomposition coincides with the spectral decomposition.</li></ul>
<h3>Real Schur decomposition</h3>
<ul><li>Applicable to: <a href="/facts/Square_matrix/5TP9mNNn">square matrix</a> A</li>
<li>Decomposition: This is a version of Schur decomposition where 
 
 
 
 V
 
 
 {\displaystyle V}
 
 and 
 
 
 
 S
 
 
 {\displaystyle S}
 
 only contain real numbers. One can always write 
 
 
 
 A
 =
 V
 S
 
 V
 
 
 T
 
 
 
 
 
 {\displaystyle A=VSV^{\mathsf {T}}}
 
 where V is a real <a href="/facts/Orthogonal_matrix/WDXBXYB5">orthogonal matrix</a>, 
 
 
 
 
 V
 
 
 T
 
 
 
 
 
 {\displaystyle V^{\mathsf {T}}}
 
 is the <a href="/facts/Matrix_transpose/8wmsagGS">transpose</a> of V, and S is a <a href="/facts/Block_matrix/fPI8HX9f">block upper triangular</a> matrix called the real <a href="/facts/Schur_form/EYNrqHmh">Schur form</a>. The blocks on the diagonal of S are of size 1×1 (in which case they represent real eigenvalues) or 2×2 (in which case they are derived from <a href="/facts/Complex_conjugate/7TEaoQDe">complex conjugate</a> eigenvalue pairs).</li></ul>
<h3>QZ decomposition</h3>
Main article: <a href="/facts/QZ_decomposition/EYNrqHmh">QZ decomposition</a>
<ul><li>Also called: generalized Schur decomposition</li>
<li>Applicable to: <a href="/facts/Square_matrix/5TP9mNNn">square matrices</a> A and B</li>
<li>Comment: there are two versions of this decomposition: complex and real.</li>
<li>Decomposition (complex version): 
 
 
 
 A
 =
 Q
 S
 
 Z
 
 ∗
 
 
 
 
 {\displaystyle A=QSZ^{*}}
 
 and 
 
 
 
 B
 =
 Q
 T
 
 Z
 
 ∗
 
 
 
 
 {\displaystyle B=QTZ^{*}}
 
 where Q and Z are <a href="/facts/Unitary_matrix/8fsxv3KD">unitary matrices</a>, the * superscript represents <a href="/facts/Conjugate_transpose/7eBWtqVG">conjugate transpose</a>, and S and T are <a href="/facts/Upper_triangular/qjXSFMM3">upper triangular</a> matrices.</li>
<li>Comment: in the complex QZ decomposition, the ratios of the diagonal elements of S to the corresponding diagonal elements of T, 
 
 
 
 
 λ
 
 i
 
 
 =
 
 S
 
 i
 i
 
 
 
 /
 
 
 T
 
 i
 i
 
 
 
 
 {\displaystyle \lambda _{i}=S_{ii}/T_{ii}}
 
, are the generalized <a href="/facts/Eigenvalue/8TjEoT8u">eigenvalues</a> that solve the <a href="/facts/Eigendecomposition_of_a_matrix/a2nfF7hJ">generalized eigenvalue problem</a> 
 
 
 
 A
 
 v
 
 =
 λ
 B
 
 v
 
 
 
 {\displaystyle A\mathbf {v} =\lambda B\mathbf {v} }
 
 (where 
 
 
 
 λ
 
 
 {\displaystyle \lambda }
 
 is an unknown scalar and v is an unknown nonzero vector).</li>
<li>Decomposition (real version): 
 
 
 
 A
 =
 Q
 S
 
 Z
 
 
 T
 
 
 
 
 
 {\displaystyle A=QSZ^{\mathsf {T}}}
 
 and 
 
 
 
 B
 =
 Q
 T
 
 Z
 
 
 T
 
 
 
 
 
 {\displaystyle B=QTZ^{\mathsf {T}}}
 
 where A, B, Q, Z, S, and T are matrices containing real numbers only. In this case Q and Z are <a href="/facts/Orthogonal_matrix/WDXBXYB5">orthogonal matrices</a>, the T superscript represents <a href="/facts/Matrix_transpose/8wmsagGS">transposition</a>, and S and T are <a href="/facts/Block_matrix/fPI8HX9f">block upper triangular</a> matrices. The blocks on the diagonal of S and T are of size 1×1 or 2×2.</li></ul>
<h3>Takagi's factorization</h3>
<ul><li>Applicable to: square, complex, symmetric matrix A.</li>
<li>Decomposition: 
 
 
 
 A
 =
 V
 D
 
 V
 
 
 T
 
 
 
 
 
 {\displaystyle A=VDV^{\mathsf {T}}}
 
, where D is a real nonnegative <a href="/facts/Diagonal_matrix/tjIAHrE6">diagonal matrix</a>, and V is <a href="/facts/Unitary_matrix/8fsxv3KD">unitary</a>. 
 
 
 
 
 V
 
 
 T
 
 
 
 
 
 {\displaystyle V^{\mathsf {T}}}
 
 denotes the <a href="/facts/Matrix_transpose/8wmsagGS">matrix transpose</a> of V.</li>
<li>Comment: The diagonal elements of D are the nonnegative square roots of the eigenvalues of 
 
 
 
 A
 
 A
 
 ∗
 
 
 =
 V
 
 D
 
 2
 
 
 
 V
 
 −
 1
 
 
 
 
 {\displaystyle AA^{*}=VD^{2}V^{-1}}
 
.</li>
<li>Comment: V may be complex even if A is real.</li>
<li>Comment: This is not a special case of the eigendecomposition (see above), which uses 
 
 
 
 
 V
 
 −
 1
 
 
 
 
 {\displaystyle V^{-1}}
 
 instead of 
 
 
 
 
 V
 
 
 T
 
 
 
 
 
 {\displaystyle V^{\mathsf {T}}}
 
. Moreover, if A is not real, it is not Hermitian and the form using 
 
 
 
 
 V
 
 ∗
 
 
 
 
 {\displaystyle V^{*}}
 
 also does not apply.</li></ul>
<h3>Singular value decomposition</h3>
Main article: <a href="/facts/Singular_value_decomposition/8t3v6wk3">Singular value decomposition</a>
<ul><li>Applicable to: m-by-n matrix A.</li>
<li>Decomposition: 
 
 
 
 A
 =
 U
 D
 
 V
 
 ∗
 
 
 
 
 {\displaystyle A=UDV^{*}}
 
, where D is a nonnegative <a href="/facts/Diagonal_matrix/tjIAHrE6">diagonal matrix</a>, and U and V satisfy 
 
 
 
 
 U
 
 ∗
 
 
 U
 =
 I
 ,
 
 V
 
 ∗
 
 
 V
 =
 I
 
 
 {\displaystyle U^{*}U=I,V^{*}V=I}
 
. Here 
 
 
 
 
 V
 
 ∗
 
 
 
 
 {\displaystyle V^{*}}
 
 is the <a href="/facts/Conjugate_transpose/7eBWtqVG">conjugate transpose</a> of V (or simply the <a href="/facts/Matrix_transpose/8wmsagGS">transpose</a>, if V contains real numbers only), and I denotes the identity matrix (of some dimension).</li>
<li>Comment: The diagonal elements of D are called the <a href="/facts/Singular_value/0QWHnqGG">singular values</a> of A.</li>
<li>Comment: Like the eigendecomposition above, the singular value decomposition involves finding basis directions along which matrix multiplication is equivalent to scalar multiplication, but it has greater generality since the matrix under consideration need not be square.</li>
<li>Uniqueness: the singular values of 
 
 
 
 A
 
 
 {\displaystyle A}
 
 are always uniquely determined. 
 
 
 
 U
 
 
 {\displaystyle U}
 
 and 
 
 
 
 V
 
 
 {\displaystyle V}
 
 need not to be unique in general.</li></ul>
<h3>Scale-invariant decompositions</h3>
Refers to variants of existing matrix decompositions, such as the SVD, that are invariant with respect to diagonal scaling.

<ul><li>Applicable to: m-by-n matrix A.</li>
<li>Unit-Scale-Invariant Singular-Value Decomposition: 
 
 
 
 A
 =
 D
 U
 S
 
 V
 
 ∗
 
 
 E
 
 
 {\displaystyle A=DUSV^{*}E}
 
, where S is a unique nonnegative <a href="/facts/Diagonal_matrix/tjIAHrE6">diagonal matrix</a> of scale-invariant singular values, U and V are <a href="/facts/Unitary_matrix/8fsxv3KD">unitary matrices</a>, 
 
 
 
 
 V
 
 ∗
 
 
 
 
 {\displaystyle V^{*}}
 
 is the <a href="/facts/Conjugate_transpose/7eBWtqVG">conjugate transpose</a> of V, and positive diagonal matrices D and E.</li>
<li>Comment: Is analogous to the SVD except that the diagonal elements of S are invariant with respect to left and/or right multiplication of A by arbitrary nonsingular diagonal matrices, as opposed to the standard SVD for which the singular values are invariant with respect to left and/or right multiplication of A by arbitrary unitary matrices.</li>
<li>Comment: Is an alternative to the standard SVD when invariance is required with respect to diagonal rather than unitary transformations of A.</li>
<li>Uniqueness: The scale-invariant singular values of 
 
 
 
 A
 
 
 {\displaystyle A}
 
 (given by the diagonal elements of S) are always uniquely determined. Diagonal matrices D and E, and unitary U and V, are not necessarily unique in general.</li>
<li>Comment: U and V matrices are not the same as those from the SVD.</li></ul>
Analogous scale-invariant decompositions can be derived from other matrix decompositions; for example, to obtain scale-invariant eigenvalues.<a class="footnote-ref" id="fnref:4" href="#fn:4">4</a><a class="footnote-ref" id="fnref:5" href="#fn:5">5</a>

<h3>Hessenberg decomposition</h3>
<ul><li>Applicable to: <a href="/facts/Square_matrix/5TP9mNNn">square matrix</a> A.</li>
<li>Decomposition: 
 
 
 
 A
 =
 P
 H
 
 P
 
 ∗
 
 
 
 
 {\displaystyle A=PHP^{*}}
 
 where 
 
 
 
 H
 
 
 {\displaystyle H}
 
 is the <a href="/facts/Hessenberg_matrix/vy7p7s3Y">Hessenberg matrix</a> and 
 
 
 
 P
 
 
 {\displaystyle P}
 
 is a <a href="/facts/Unitary_matrix/8fsxv3KD">unitary matrix</a>.</li>
<li>Comment: often the first step in the Schur decomposition.</li></ul>
<h3>Complete orthogonal decomposition</h3>
Main article: Complete orthogonal decomposition
<ul><li>Also known as: UTV decomposition, ULV decomposition, URV decomposition.</li>
<li>Applicable to: m-by-n matrix A.</li>
<li>Decomposition: 
 
 
 
 A
 =
 U
 T
 
 V
 
 ∗
 
 
 
 
 {\displaystyle A=UTV^{*}}
 
, where T is a <a href="/facts/Triangular_matrix/qjXSFMM3">triangular matrix</a>, and U and V are <a href="/facts/Unitary_matrix/8fsxv3KD">unitary matrices</a>.</li>
<li>Comment: Similar to the singular value decomposition and to the Schur decomposition.</li></ul>
<h2 id="other-decompositions">Other decompositions</h2>
<h3>Polar decomposition</h3>
Main article: <a href="/facts/Polar_decomposition/XWQLdp9B">Polar decomposition</a>
<ul><li>Applicable to: any square complex matrix A.</li>
<li>Decomposition: 
 
 
 
 A
 =
 U
 P
 
 
 {\displaystyle A=UP}
 
 (right polar decomposition) or 
 
 
 
 A
 =
 
 P
 ′
 
 U
 
 
 {\displaystyle A=P'U}
 
 (left polar decomposition), where U is a <a href="/facts/Unitary_matrix/8fsxv3KD">unitary matrix</a> and P and P' are <a href="/facts/Positive_semidefinite_matrix/Hbr6AuPS">positive semidefinite</a> <a href="/facts/Hermitian_matrices/qArPblYo">Hermitian matrices</a>.</li>
<li>Uniqueness: 
 
 
 
 P
 
 
 {\displaystyle P}
 
 is always unique and equal to 
 
 
 
 
 
 
 A
 
 ∗
 
 
 A
 
 
 
 
 {\displaystyle {\sqrt {A^{*}A}}}
 
 (which is always hermitian and positive semidefinite). If 
 
 
 
 A
 
 
 {\displaystyle A}
 
 is invertible, then 
 
 
 
 U
 
 
 {\displaystyle U}
 
 is unique.</li>
<li>Comment: Since any Hermitian matrix admits a spectral decomposition with a unitary matrix, 
 
 
 
 P
 
 
 {\displaystyle P}
 
 can be written as 
 
 
 
 P
 =
 V
 D
 
 V
 
 ∗
 
 
 
 
 {\displaystyle P=VDV^{*}}
 
. Since 
 
 
 
 P
 
 
 {\displaystyle P}
 
 is positive semidefinite, all elements in 
 
 
 
 D
 
 
 {\displaystyle D}
 
 are non-negative. Since the product of two unitary matrices is unitary, taking 
 
 
 
 W
 =
 U
 V
 
 
 {\displaystyle W=UV}
 
one can write 
 
 
 
 A
 =
 U
 (
 V
 D
 
 V
 
 ∗
 
 
 )
 =
 W
 D
 
 V
 
 ∗
 
 
 
 
 {\displaystyle A=U(VDV^{*})=WDV^{*}}
 
 which is the singular value decomposition. Hence, the existence of the polar decomposition is equivalent to the existence of the singular value decomposition.</li></ul>
<h3>Algebraic polar decomposition</h3>
<ul><li>Applicable to: square, complex, non-singular matrix A.<a class="footnote-ref" id="fnref:6" href="#fn:6">6</a></li>
<li>Decomposition: 
 
 
 
 A
 =
 Q
 S
 
 
 {\displaystyle A=QS}
 
, where Q is a complex orthogonal matrix and S is complex symmetric matrix.</li>
<li>Uniqueness: If 
 
 
 
 
 A
 
 
 T
 
 
 
 A
 
 
 {\displaystyle A^{\mathsf {T}}A}
 
 has no negative real eigenvalues, then the decomposition is unique.<a class="footnote-ref" id="fnref:7" href="#fn:7">7</a></li>
<li>Comment: The existence of this decomposition is equivalent to 
 
 
 
 A
 
 A
 
 
 T
 
 
 
 
 
 {\displaystyle AA^{\mathsf {T}}}
 
 being similar to 
 
 
 
 
 A
 
 
 T
 
 
 
 A
 
 
 {\displaystyle A^{\mathsf {T}}A}
 
.<a class="footnote-ref" id="fnref:8" href="#fn:8">8</a></li>
<li>Comment: A variant of this decomposition is 
 
 
 
 A
 =
 R
 C
 
 
 {\displaystyle A=RC}
 
, where R is a real matrix and C is a circular matrix.<a class="footnote-ref" id="fnref:9" href="#fn:9">9</a></li></ul>
<h3>Mostow's decomposition</h3>
<ul><li>Applicable to: square, complex, non-singular matrix A.<a class="footnote-ref" id="fnref:10" href="#fn:10">10</a><a class="footnote-ref" id="fnref:11" href="#fn:11">11</a></li>
<li>Decomposition: 
 
 
 
 A
 =
 U
 
 e
 
 i
 M
 
 
 
 e
 
 S
 
 
 
 
 {\displaystyle A=Ue^{iM}e^{S}}
 
, where U is unitary, M is real anti-symmetric and S is real symmetric.</li>
<li>Comment: The matrix A can also be decomposed as 
 
 
 
 A
 =
 
 U
 
 2
 
 
 
 e
 
 
 S
 
 2
 
 
 
 
 
 e
 
 i
 
 M
 
 2
 
 
 
 
 
 
 {\displaystyle A=U_{2}e^{S_{2}}e^{iM_{2}}}
 
, where U2 is unitary, M2 is real anti-symmetric and S2 is real symmetric.<a class="footnote-ref" id="fnref:12" href="#fn:12">12</a></li></ul>
<h3>Sinkhorn normal form</h3>
Main article: <a href="/facts/Sinkhorn%2527s_theorem/kgSm7y0h">Sinkhorn's theorem</a>
<ul><li>Applicable to: square real matrix A with strictly positive elements.</li>
<li>Decomposition: 
 
 
 
 A
 =
 
 D
 
 1
 
 
 S
 
 D
 
 2
 
 
 
 
 {\displaystyle A=D_{1}SD_{2}}
 
, where S is <a href="/facts/Doubly_stochastic_matrix/IwqKyx05">doubly stochastic</a> and D1 and D2 are real diagonal matrices with strictly positive elements.</li></ul>
<h3>Sectoral decomposition</h3>
<ul><li>Applicable to: square, complex matrix A with <a href="/facts/Numerical_range/4gCh6ShW">numerical range</a> contained in the sector 
 
 
 
 
 S
 
 α
 
 
 =
 
 {
 
 r
 
 e
 
 i
 θ
 
 
 ∈
 
 C
 
 ∣
 r
 >
 0
 ,
 
 |
 
 θ
 
 |
 
 ≤
 α
 <
 
 
 π
 2
 
 
 
 }
 
 
 
 {\displaystyle S_{\alpha }=\left\{re^{i\theta }\in \mathbb {C} \mid r>0,|\theta |\leq \alpha <{\frac {\pi }{2}}\right\}}
 
.</li>
<li>Decomposition: 
 
 
 
 A
 =
 C
 Z
 
 C
 
 ∗
 
 
 
 
 {\displaystyle A=CZC^{*}}
 
, where C is an invertible complex matrix and 
 
 
 
 Z
 =
 diag
 ⁡
 
 (
 
 
 e
 
 i
 
 θ
 
 1
 
 
 
 
 ,
 …
 ,
 
 e
 
 i
 
 θ
 
 n
 
 
 
 
 
 )
 
 
 
 {\displaystyle Z=\operatorname {diag} \left(e^{i\theta _{1}},\ldots ,e^{i\theta _{n}}\right)}
 
 with all 
 
 
 
 
 |
 
 θ
 
 j
 
 
 |
 
 ≤
 α
 
 
 {\displaystyle \left|\theta _{j}\right|\leq \alpha }
 
.<a class="footnote-ref" id="fnref:13" href="#fn:13">13</a><a class="footnote-ref" id="fnref:14" href="#fn:14">14</a></li></ul>
<h3>Williamson's normal form</h3>
<ul><li>Applicable to: square, <a href="/facts/Positive-definite_matrix/Hbr6AuPS">positive-definite</a> real matrix A with order 2n×2n.</li>
<li>Decomposition: 
 
 
 
 A
 =
 
 S
 
 
 T
 
 
 
 diag
 ⁡
 (
 D
 ,
 D
 )
 S
 
 
 {\displaystyle A=S^{\mathsf {T}}\operatorname {diag} (D,D)S}
 
, where 
 
 
 
 S
 ∈
 
 Sp
 
 (
 2
 n
 )
 
 
 {\displaystyle S\in {\text{Sp}}(2n)}
 
 is a <a href="/facts/Symplectic_matrix/kkWvzo5M">symplectic matrix</a> and D is a nonnegative n-by-n diagonal matrix.<a class="footnote-ref" id="fnref:15" href="#fn:15">15</a></li></ul>
<h3>Matrix square root</h3>
Main article: <a href="/facts/Square_root_of_a_matrix/uWGkQQqN">Square root of a matrix</a>
<ul><li>Decomposition: 
 
 
 
 A
 =
 B
 B
 
 
 {\displaystyle A=BB}
 
, not unique in general.</li>
<li>In the case of positive semidefinite 
 
 
 
 A
 
 
 {\displaystyle A}
 
, there is a unique positive semidefinite 
 
 
 
 B
 
 
 {\displaystyle B}
 
 such that 
 
 
 
 A
 =
 
 B
 
 ∗
 
 
 B
 =
 B
 B
 
 
 {\displaystyle A=B^{*}B=BB}
 
.</li></ul>
<h2 id="generalizations">Generalizations</h2>

There exist analogues of the SVD, QR, LU and Cholesky factorizations for quasimatrices and cmatrices or continuous matrices.<a class="footnote-ref" id="fnref:16" href="#fn:16">16</a> A ‘quasimatrix’ is, like a matrix, a rectangular scheme whose elements are indexed, but one discrete index is replaced by a continuous index. Likewise, a ‘cmatrix’, is continuous in both indices. As an example of a cmatrix, one can think of the kernel of an <a href="/facts/Integral_operator/AL6u7Wrj">integral operator</a>.
These factorizations are based on early work by Fredholm (1903), Hilbert (1904) and Schmidt (1907). For an account, and a translation to English of the seminal papers, see Stewart (2011).

<h2 id="see-also">See also</h2>
<ul><li><a href="/facts/Matrix_splitting/lSExxJBj">Matrix splitting</a></li>
<li><a href="/facts/Non-negative_matrix_factorization/DIBLH8B0">Non-negative matrix factorization</a></li>
<li><a href="/facts/Principal_component_analysis/pCbVTqdR">Principal component analysis</a></li></ul>

<h3>Notes</h3>

<h3>Citations</h3>

<h3>Bibliography</h3>
<ul><li>Choudhury, Dipa; Horn, Roger A. (April 1987). "A Complex Orthogonal-Symmetric Analog of the Polar Decomposition". SIAM Journal on Algebraic and Discrete Methods. 8 (2): 219–225. <a href="/facts/Doi_(identifier)/muM9Etpq">doi</a>:<a href="https://doi.org/10.1137%2F0608019">10.1137/0608019</a>.</li>
<li><a href="/facts/Ivar_Fredholm/x4YvpYLv">Fredholm, I.</a> (1903), "Sur une classe d'´equations fonctionnelles", Acta Mathematica (in French), 27: 365–390, <a href="/facts/Doi_(identifier)/muM9Etpq">doi</a>:<a href="https://doi.org/10.1007%2Fbf02421317">10.1007/bf02421317</a></li>
<li><a href="/facts/David_Hilbert/1vhlHn4b">Hilbert, D.</a> (1904), "Grundzüge einer allgemeinen Theorie der linearen Integralgleichungen", Nachr. Königl. Ges. Gött (in German), 1904: 49–91</li>
<li>Horn, Roger A.; Merino, Dennis I. (January 1995). <a href="https://doi.org/10.1016%2F0024-3795%2893%2900056-6">"Contragredient equivalence: A canonical form and some applications"</a>. Linear Algebra and Its Applications. 214: 43–92. <a href="/facts/Doi_(identifier)/muM9Etpq">doi</a>:<a href="https://doi.org/10.1016%2F0024-3795%2893%2900056-6">10.1016/0024-3795(93)00056-6</a>.</li>
<li>Meyer, C. D. (2000), <a href="http://www.matrixanalysis.com/">Matrix Analysis and Applied Linear Algebra</a>, <a href="/facts/Society_for_Industrial_and_Applied_Mathematics/2HY4T5WY">SIAM</a>, <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-0-89871-454-8</li>
<li><a href="/facts/Erhard_Schmidt/5UGgwTah">Schmidt, E.</a> (1907), <a href="https://zenodo.org/record/1428258">"Zur Theorie der linearen und nichtlinearen Integralgleichungen. I Teil. Entwicklung willkürlichen Funktionen nach System vorgeschriebener"</a>, Mathematische Annalen (in German), 63 (4): 433–476, <a href="/facts/Doi_(identifier)/muM9Etpq">doi</a>:<a href="https://doi.org/10.1007%2Fbf01449770">10.1007/bf01449770</a></li>
<li>Simon, C.; Blume, L. (1994). Mathematics for Economists. Norton. <a href="/facts/ISBN_(identifier)/15AdSPa9">ISBN</a> 978-0-393-95733-4.</li>
<li>Stewart, G. W. (2011), <a href="http://www.cs.umd.edu/~stewart/FHS.pdf">Fredholm, Hilbert, Schmidt: three fundamental papers on integral equations</a> (PDF), retrieved 2015-01-06</li>
<li>Townsend, A.; Trefethen, L. N. (2015), "Continuous analogues of matrix factorizations", <a href="/facts/Proceedings_of_the_Royal_Society/qyn2m2ja">Proc. R. Soc. A</a>, 471 (2173): 20140585, <a href="/facts/Bibcode_(identifier)/9HtdQSGB">Bibcode</a>:<a href="https://ui.adsabs.harvard.edu/abs/2014RSPSA.47140585T">2014RSPSA.47140585T</a>, <a href="/facts/Doi_(identifier)/muM9Etpq">doi</a>:<a href="https://doi.org/10.1098%2Frspa.2014.0585">10.1098/rspa.2014.0585</a>, <a href="/facts/PMC_(identifier)/dX1zMt71">PMC</a> <a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC4277194">4277194</a>, <a href="/facts/PMID_(identifier)/JlHAvMHt">PMID</a> <a href="https://pubmed.ncbi.nlm.nih.gov/25568618">25568618</a></li>
<li>Jun, Lu (2021), Numerical matrix decomposition and its modern applications: A rigorous first course, <a href="/facts/ArXiv_(identifier)/H6EtgnBe">arXiv</a>:<a href="https://arxiv.org/abs/2107.02579">2107.02579</a></li></ul>

<h2 id="external-links">External links</h2>
<ul><li><a href="http://www.bluebit.gr/matrix-calculator/">Online Matrix Calculator</a></li>
<li><a href="http://www.wolframalpha.com/input/?i=matrix+decomposition&rawformassumption={%22C%22,+%22matrix+decomposition%22}+-%3E+{%22Calculator%22}&rawformassumption={%22MC%22,%22%22}-%3E{%22Formula%22}">Wolfram Alpha Matrix Decomposition Computation » LU and QR Decomposition</a></li>
<li><a href="https://encyclopediaofmath.org/wiki/Matrix_factorization">Springer Encyclopaedia of Mathematics » Matrix factorization</a></li>
<li><a href="https://web.archive.org/web/20110314171151/http://www.graphlab.ml.cmu.edu/pmf.html">GraphLab</a> <a href="/facts/GraphLab/NEWTJ1e6">GraphLab</a> <a href="/facts/Collaborative_filtering/NhSFypjN">collaborative filtering</a> library, large scale parallel implementation of matrix decomposition methods (in C++) for multicore.</li></ul>

<h2 id="references">References</h2>

<ol>
<li id="fn:1">Lay, David C. (2016). Linear algebra and its applications. Steven R. Lay, Judith McDonald (Fifth Global ed.). Harlow. p. 142. ISBN 978-1-292-09223-2. OCLC 920463015.{{cite book}}: CS1 maint: location missing publisher (link) <a href="978-1-292-09223-2" target="_blank">978-1-292-09223-2</a> <a href="#fnref:1" class="footnote-back-ref">↩</a></li>
<li id="fn:2">If a non-square matrix is used, however, then the matrix U will also have the same rectangular shape as the original matrix A. And so, calling the matrix U upper triangular would be incorrect as the correct term would be that U is the 'row echelon form' of A. Other than this, there are no differences in LU factorization for square and non-square matrices. <a href="#fnref:2" class="footnote-back-ref">↩</a></li>
<li id="fn:3">Piziak, R.; Odell, P. L. (1 June 1999). "Full Rank Factorization of Matrices". Mathematics Magazine. 72 (3): 193. doi:10.2307/2690882. JSTOR 2690882. <a href="/wiki/Doi_(identifier)" target="_blank">/wiki/Doi_(identifier)</a> <a href="#fnref:3" class="footnote-back-ref">↩</a></li>
<li id="fn:4">Uhlmann, J.K. (2018), "A Generalized Matrix Inverse that is Consistent with Respect to Diagonal Transformations", SIAM Journal on Matrix Analysis and Applications, 239 (2): 781–800, doi:10.1137/17M113890X <a href="/wiki/Doi_(identifier)" target="_blank">/wiki/Doi_(identifier)</a> <a href="#fnref:4" class="footnote-back-ref">↩</a></li>
<li id="fn:5">Uhlmann, J.K. (2018), "A Rank-Preserving Generalized Matrix Inverse for Consistency with Respect to Similarity", IEEE Control Systems Letters, 3: 91–95, arXiv:1804.07334, doi:10.1109/LCSYS.2018.2854240, ISSN 2475-1456, S2CID 5031440 <a href="/wiki/ArXiv_(identifier)" target="_blank">/wiki/ArXiv_(identifier)</a> <a href="#fnref:5" class="footnote-back-ref">↩</a></li>
<li id="fn:6">Choudhury & Horn 1987, pp. 219–225 - Choudhury, Dipa; Horn, Roger A. (April 1987). "A Complex Orthogonal-Symmetric Analog of the Polar Decomposition". SIAM Journal on Algebraic and Discrete Methods. 8 (2): 219–225. doi:10.1137/0608019. <a href="https://doi.org/10.1137%2F0608019" target="_blank">https://doi.org/10.1137%2F0608019</a> <a href="#fnref:6" class="footnote-back-ref">↩</a></li>
<li id="fn:7">Bhatia, Rajendra (2013-11-15). "The bipolar decomposition". Linear Algebra and Its Applications. 439 (10): 3031–3037. doi:10.1016/j.laa.2013.09.006. <a href="/wiki/Doi_(identifier)" target="_blank">/wiki/Doi_(identifier)</a> <a href="#fnref:7" class="footnote-back-ref">↩</a></li>
<li id="fn:8">Horn & Merino 1995, pp. 43–92 - Horn, Roger A.; Merino, Dennis I. (January 1995). "Contragredient equivalence: A canonical form and some applications". Linear Algebra and Its Applications. 214: 43–92. doi:10.1016/0024-3795(93)00056-6. <a href="https://doi.org/10.1016%2F0024-3795%2893%2900056-6" target="_blank">https://doi.org/10.1016%2F0024-3795%2893%2900056-6</a> <a href="#fnref:8" class="footnote-back-ref">↩</a></li>
<li id="fn:9">Bhatia, Rajendra (2013-11-15). "The bipolar decomposition". Linear Algebra and Its Applications. 439 (10): 3031–3037. doi:10.1016/j.laa.2013.09.006. <a href="/wiki/Doi_(identifier)" target="_blank">/wiki/Doi_(identifier)</a> <a href="#fnref:9" class="footnote-back-ref">↩</a></li>
<li id="fn:10">Mostow, G. D. (1955), Some new decomposition theorems for semi-simple groups, Mem. Amer. Math. Soc., vol. 14, American Mathematical Society, pp. 31–54 <a href="https://archive.org/details/liealgebrasandli029541mbp" target="_blank">https://archive.org/details/liealgebrasandli029541mbp</a> <a href="#fnref:10" class="footnote-back-ref">↩</a></li>
<li id="fn:11">Nielsen, Frank; Bhatia, Rajendra (2012). Matrix Information Geometry. Springer. p. 224. arXiv:1007.4402. doi:10.1007/978-3-642-30232-9. ISBN 9783642302329. S2CID 118466496. <a href="9783642302329" target="_blank">9783642302329</a> <a href="#fnref:11" class="footnote-back-ref">↩</a></li>
<li id="fn:12">Bhatia, Rajendra (2013-11-15). "The bipolar decomposition". Linear Algebra and Its Applications. 439 (10): 3031–3037. doi:10.1016/j.laa.2013.09.006. <a href="/wiki/Doi_(identifier)" target="_blank">/wiki/Doi_(identifier)</a> <a href="#fnref:12" class="footnote-back-ref">↩</a></li>
<li id="fn:13">Zhang, Fuzhen (30 June 2014). "A matrix decomposition and its applications". Linear and Multilinear Algebra. 63 (10): 2033–2042. doi:10.1080/03081087.2014.933219. S2CID 19437967. <a href="https://zenodo.org/record/851661" target="_blank">https://zenodo.org/record/851661</a> <a href="#fnref:13" class="footnote-back-ref">↩</a></li>
<li id="fn:14">Drury, S.W. (November 2013). "Fischer determinantal inequalities and Highamʼs Conjecture". Linear Algebra and Its Applications. 439 (10): 3129–3133. doi:10.1016/j.laa.2013.08.031. <a href="https://doi.org/10.1016%2Fj.laa.2013.08.031" target="_blank">https://doi.org/10.1016%2Fj.laa.2013.08.031</a> <a href="#fnref:14" class="footnote-back-ref">↩</a></li>
<li id="fn:15">Idel, Martin; Soto Gaona, Sebastián; Wolf, Michael M. (2017-07-15). "Perturbation bounds for Williamson's symplectic normal form". Linear Algebra and Its Applications. 525: 45–58. arXiv:1609.01338. doi:10.1016/j.laa.2017.03.013. S2CID 119578994. <a href="/wiki/ArXiv_(identifier)" target="_blank">/wiki/ArXiv_(identifier)</a> <a href="#fnref:15" class="footnote-back-ref">↩</a></li>
<li id="fn:16">Townsend & Trefethen 2015 - Townsend, A.; Trefethen, L. N. (2015), "Continuous analogues of matrix factorizations", Proc. R. Soc. A, 471 (2173): 20140585, Bibcode:2014RSPSA.47140585T, doi:10.1098/rspa.2014.0585, PMC 4277194, PMID 25568618 <a href="https://ui.adsabs.harvard.edu/abs/2014RSPSA.47140585T" target="_blank">https://ui.adsabs.harvard.edu/abs/2014RSPSA.47140585T</a> <a href="#fnref:16" class="footnote-back-ref">↩</a></li>
</ol>

Matrix decomposition open-in-new

Matrix decomposition