scalar_centered_sampling_bernstein_tail_from_entry_frobenius_scales
ProvedRole. It is a centered-sampling fluctuation estimate, one of the reusable concentration interfaces used repeatedly by the Neumann-term bounds.
Problem and notation. Exact matrix completion asks when an unknown low-rank real matrix can be recovered from a random subset of its entries. Here has rank , entries are observed, and . Recovery means nuclear-norm minimization: minimize among matrices agreeing with on the observed entries. Probability notation. is the fixed-cardinality success probability: is chosen uniformly among all subsets of entries with , and the event is that the convex program uniquely returns . In Bernoulli nodes, or means each entry is sampled independently with probability , usually . Coherence notation. The object records SVD/singular-vector data for . The hypotheses and are the Candes-Recht incoherence assumptions: measures how spread out the singular vector spaces are, and measures the largest entry of the sign matrix . The parameter controls polynomial failure probabilities such as .
Claim. Raw scalar Bernstein inequality for the entry-sum of a centered sampling fluctuation. The bound is stated in terms of exposed entry and Frobenius scales, before any Candes-Recht sample-size arithmetic is absorbed.
Lecture-note formulation:
The constants in this node are universal existential constants; the theorem asserts that some positive constants with these roles exist.
Decomposition status. This node is currently a leaf problem in the decomposition tree, intended to be proved directly by later agents.
import Definitions.Def_matrix_completion_neumann open MatrixCompletion
theorem scalar_centered_sampling_bernstein_tail_from_entry_frobenius_scales :
∃ Cbern cbern : ℝ, 0 < Cbern ∧ 0 < cbern ∧
∀ (β : ℝ), 2 < β →
∀ (n₁ n₂ m : ℕ), 0 < n₁ → 0 < n₂ → m ≤ n₁ * n₂ →
∀ (Coeff : Finset (Fin n₁ × Fin n₂) → ℝ)
(B : Matrix (Fin n₁) (Fin n₂) ℝ)
(entryScale frobScale : ℝ),
(∀ Omega : Finset (Fin n₁ × Fin n₂),
Coeff Omega =
matrixEntrySum
(centeredSamplingFluctuation Omega
((m : ℝ) / ((n₁ : ℝ) * (n₂ : ℝ))) B)) →
entrySupNorm B ≤ entryScale →
frobeniusNorm B ≤ frobScale →
bernoulliEventProb ((m : ℝ) / ((n₁ : ℝ) * (n₂ : ℝ)))
(fun Omega =>
|Coeff Omega| ≤
Cbern *
(Real.sqrt
((β * Real.log (↑(max n₁ n₂))) /
((m : ℝ) / ((n₁ : ℝ) * (n₂ : ℝ)))) *
frobScale +
((β * Real.log (↑(max n₁ n₂))) /
((m : ℝ) / ((n₁ : ℝ) * (n₂ : ℝ)))) *
entryScale)) ≥
1 - cbern * Real.rpow (↑(max n₁ n₂)) (-β) := by
sorry