Alphabeta Math
TheoremStatement: Literature-sourcedProof: AI-adaptedPipeline-generatedprecheck pass
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 A∈Fm×n, b∈Fm, and τ≥0, where F=R or C. Let A=UΣV∗ be an exact SVD. Define Στ by retaining σi when σi>τ and replacing each σi≤τ by zero, and set Aτ=UΣτV∗. Put q=#{i:σi>τ}, the spectral numerical rank of A at threshold τ. The vector xτ=Aτ+b is the unique minimum-Euclidean-norm least-squares solution for the truncated model Aτx≈b. 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 xτ; a failure or exhausted budget returns inconclusive. No convergence within an arbitrary budget is asserted. Moreover, ∥A−Aτ∥2=σq+1≤τ,∥Axτ−b∥2≤∥Aτxτ−b∥2+τ∥xτ∥2, with σp+1=0 when q=p=min⁡(m,n). If every discarded singular value is zero, then Aτ=A and xτ 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 b, and the associated truncated matrix and pseudoinverse.

[F1]

Every linear map between finite-dimensional real or complex inner product spaces admits a singular value decomposition: A=UΣV∗ with orthonormal singular vectors and nonnegative diagonal singular values.

[F2]

The truncated singular-value pseudoinverse obtained by discarding singular values below a declared numerical-rank threshold: Aτ+=VΣτ+U∗, where diagonal entries are 1/σi for σi>τ and zero otherwise.

[F3]

For every right-hand side b, A+b is the unique least-squares solution of minimum Euclidean norm: for any matrix M and right-hand side b, M+b is the unique least-squares solution of minimum Euclidean norm.

[F4]

The best rank-at-most-k approximation in operator norm is the rank-k truncation of a singular value decomposition: for 0≤q<n, the SVD truncation to rank at most q has spectral residual σq+1, including the stated zero-tail convention.

[F5]

Spectral numerical rank counts singular values above threshold over R or C: for m,n≥1, the spectral numerical rank at threshold τ counts precisely the singular values strictly greater than τ.

[F6]

The Moore--Penrose pseudoinverse A+ 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.

[F7]

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 ∥Mx∥2 over unit vectors otherwise. Thus ∥Mx∥2≤∥M∥2∥x∥2 for every vector x by scaling when x≠0, and directly when x=0.

[F8]

The inner-product norm is definite, homogeneous, and satisfies the triangle inequality: the Euclidean vector norm satisfies ∥u+v∥2≤∥u∥2+∥v∥2.

[F9]

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.

[F10]

Rank-thresholded least squares by the singular value decomposition: on exact-arithmetic convergence with the returned SVD record the procedure computes xτ=Aτ+b; a failure or exhausted sweep budget returns inconclusive.

Proof

1.1F1F2F6algebra

Let I={i≤p:σi>τ}. By [F1]–[F2], Aτ=∑i∈Iσiuivi∗ and Aτ+=∑i∈Iσi−1viui∗. Orthonormality gives AτAτ+=∑i∈Iuiui∗ and Aτ+Aτ=∑i∈Ivivi∗, which are self-adjoint. Using ui∗uj=vi∗vj=δij in those sums verifies the two Penrose product equations. Together with self-adjointness, these are the four equations; by [F6], Aτ+ is the unique Moore--Penrose pseudoinverse of Aτ. The calculation includes I=∅, when both matrices and sums are zero.

1.2F1F4F5F7F9algebra

If p=0, the matrix has no singular-value terms, so Aτ=A, q=0, and its spectral numerical rank is zero because A has rank zero and distance zero from itself by [F9]. If p>0, [F5] identifies q with the spectral numerical rank retained by the strict threshold. When q<p, one has q<p≤n, so [F4] applies with k=q and gives ∥A−Aτ∥2=σq+1≤τ. When q=p, [F1] shows directly that all positive singular terms are retained, hence Aτ=A and the residual norm is zero, equal to the stated convention σp+1=0. This last branch also covers p=0 without applying [F4] or [F5] outside their dimension hypotheses.

2.1step 1.1F3F10

Apply [F3] to M=Aτ. Its unique minimum-norm least-squares solution is Aτ+b=xτ, 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.

2.2step 1.2F7F8algebra

By [F7], ∥(A−Aτ)xτ∥2≤∥A−Aτ∥2∥xτ∥2; [F8] gives ∥Axτ−b∥2≤∥Aτxτ−b∥2+∥(A−Aτ)xτ∥2. Combining with step 1.2 proves the stated residual bound.

3.1step 2.1step 1.2algebra

If every discarded singular value is zero, the SVD formula gives A=Aτ, so step 2.1 is the original-matrix solution. If any positive singular value is discarded, the proof has established optimality only for Aτ; it supplies no normal-equation or minimum-residual claim for A.

4.1step 1.1step 1.2step 2.1step 2.2step 3.1∎

Steps 1.1, 1.2, 2.1, 2.2, and 3.1 prove all stated claims and the limitation.

Depends on

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