How statement and proof provenance work
The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.
- Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
- AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
- AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.
These labels describe origin, not correctness: citations and verification chips remain separate evidence.
The SVD threshold solver gives the unique minimum-norm fit to its truncated model
Statement
Let , , and , where or . Let be an exact SVD. Define by retaining when and replacing each by zero, and set . Put , the spectral numerical rank of at threshold . The vector is the unique minimum-Euclidean-norm least-squares solution for the truncated model . An exact-arithmetic run of Rank-thresholded least squares by the singular value decomposition that converges within its supplied finite sweep budget and supplies this SVD record returns ; a failure or exhausted budget returns inconclusive. No convergence within an arbitrary budget is asserted. Moreover,
with when . If every discarded singular value is zero, then and is the unique minimum-norm least-squares solution for the original matrix. If a positive singular value is discarded, no original-matrix least-squares optimality is asserted.
Facts & Assumptions
Given: An exact SVD, threshold , right-hand side , and the associated truncated matrix and pseudoinverse.
Every linear map between finite-dimensional real or complex inner product spaces admits a singular value decomposition: with orthonormal singular vectors and nonnegative diagonal singular values.
The truncated singular-value pseudoinverse obtained by discarding singular values below a declared numerical-rank threshold: , where diagonal entries are for and zero otherwise.
For every right-hand side , is the unique least-squares solution of minimum Euclidean norm: for any matrix and right-hand side , is the unique least-squares solution of minimum Euclidean norm.
The best rank-at-most-k approximation in operator norm is the rank-k truncation of a singular value decomposition: for , the SVD truncation to rank at most has spectral residual , including the stated zero-tail convention.
Spectral numerical rank counts singular values above threshold over R or C: for , the spectral numerical rank at threshold counts precisely the singular values strictly greater than .
The Moore--Penrose pseudoinverse as the solution of the four Penrose equations and Every finite real or complex matrix has a unique Moore--Penrose pseudoinverse: the four Penrose equations define the Moore--Penrose pseudoinverse, and every finite real or complex matrix has a unique matrix satisfying them.
The operator norm is zero on the zero domain and otherwise is max_{||v||=1} ||Tv||: the matrix operator norm is zero on a zero domain and is the maximum of over unit vectors otherwise. Thus for every vector by scaling when , and directly when .
The inner-product norm is definite, homogeneous, and satisfies the triangle inequality: the Euclidean vector norm satisfies .
Numerical rank relative to a declared norm, scale, and tolerance defines numerical rank as the least rank of a matrix within the declared norm tolerance, for all finite matrix dimensions, including zero.
Rank-thresholded least squares by the singular value decomposition: on exact-arithmetic convergence with the returned SVD record the procedure computes ; a failure or exhausted sweep budget returns inconclusive.
Proof
Let . By [F1]–[F2], and . Orthonormality gives and , which are self-adjoint. Using in those sums verifies the two Penrose product equations. Together with self-adjointness, these are the four equations; by [F6], is the unique Moore--Penrose pseudoinverse of . The calculation includes , when both matrices and sums are zero.
If , the matrix has no singular-value terms, so , , and its spectral numerical rank is zero because has rank zero and distance zero from itself by [F9]. If , [F5] identifies with the spectral numerical rank retained by the strict threshold. When , one has , so [F4] applies with and gives . When , [F1] shows directly that all positive singular terms are retained, hence and the residual norm is zero, equal to the stated convention . This last branch also covers without applying [F4] or [F5] outside their dimension hypotheses.
Apply [F3] to . Its unique minimum-norm least-squares solution is , proving the truncated-model claim independently of any iteration budget. For an exact-arithmetic run converging within its supplied budget and returning the stated SVD record, [F10] identifies this vector as the returned solution; the failure branches return inconclusive.
By [F7], ; [F8] gives . Combining with step 1.2 proves the stated residual bound.
If every discarded singular value is zero, the SVD formula gives , so step 2.1 is the original-matrix solution. If any positive singular value is discarded, the proof has established optimality only for ; it supplies no normal-equation or minimum-residual claim for .
Steps 1.1, 1.2, 2.1, 2.2, and 3.1 prove all stated claims and the limitation.
Depends on
- Numerical rank relative to a declared norm, scale, and tolerance
- Rank-thresholded least squares by the singular value decomposition
- The truncated singular-value pseudoinverse obtained by discarding singular values below a declared numerical-rank threshold
- The Moore--Penrose pseudoinverse $A^+$ as the solution of the four Penrose equations
- Every finite real or complex matrix has a unique Moore--Penrose pseudoinverse
- Every linear map between finite-dimensional real or complex inner product spaces admits a singular value decomposition
- For every right-hand side $b$, $A^+b$ is the unique least-squares solution of minimum Euclidean norm
- The best rank-at-most-k approximation in operator norm is the rank-k truncation of a singular value decomposition
- Spectral numerical rank counts singular values above threshold over R or C
- The operator norm is zero on the zero domain and otherwise is max_{||v||=1} ||Tv||
- The inner-product norm is definite, homogeneous, and satisfies the triangle inequality
Used by
Dependency tree · two levels
35 results within two dependency steps of this one, each drawn at its shortest distance from it. An arrow runs from a result to what uses it, so the chart reads left to right and ends at this result, which carries a heavier outline. Every node is a link to that result. Click elsewhere on the chart to enlarge it.
Sources
- Wen-Wei Lin, Lecture Notes of Matrix Computations, §3.2.2 (standard reference, not scraped)