Prove2Me
Navigate
DiscoverFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

rademacher_sampled_matrix_moment_from_row_column_energy

Proved

by Shuze Chen · Jun 14, 2026 · Mathlib 0df444a (Lean v4.33.1)

candes-rechtconvex-optimizationlean4matrix-completionmoment-boundsprobabilityrademacher

Role. It is part of the symmetrization and matrix-moment machinery behind the spectral norm concentration estimates.

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 M∈Rn1×n2M\in\mathbb R^{n_1\times n_2}M∈Rn1​×n2​ has rank rrr, mmm entries are observed, and n=max⁡(n1,n2)n=\max(n_1,n_2)n=max(n1​,n2​). Recovery means nuclear-norm minimization: minimize ∥X∥∗\|X\|_*∥X∥∗​ among matrices XXX agreeing with MMM on the observed entries. Probability notation. successProb⁡(m,M)\operatorname{successProb}(m,M)successProb(m,M) is the fixed-cardinality success probability: Ω\OmegaΩ is chosen uniformly among all subsets of n1n2n_1n_2n1​n2​ entries with ∣Ω∣=m|\Omega|=m∣Ω∣=m, and the event is that the convex program uniquely returns MMM. In Bernoulli nodes, Pp(E)\mathbb P_p(E)Pp​(E) or bernoulliEventProb⁡(p,E)\operatorname{bernoulliEventProb}(p,E)bernoulliEventProb(p,E) means each entry is sampled independently with probability ppp, usually p=m/(n1n2)p=m/(n_1n_2)p=m/(n1​n2​). Coherence notation. The object SSS records SVD/singular-vector data for MMM. The hypotheses A0(S,μ0)A0(S,\mu_0)A0(S,μ0​) and A1(S,μ1)A1(S,\mu_1)A1(S,μ1​) are the Candes-Recht incoherence assumptions: μ0\mu_0μ0​ measures how spread out the singular vector spaces are, and μ1\mu_1μ1​ measures the largest entry of the sign matrix UV⊤UV^\topUV⊤. The parameter β>2\beta>2β>2 controls polynomial failure probabilities such as n−βn^{-\beta}n−β. For sampled row/column nodes, Ni(Ω)N_i(\Omega)Ni​(Ω) counts observed entries in row iii, Nj(Ω)N^j(\Omega)Nj(Ω) counts observed entries in column jjj, and the corresponding energies sum Xij2X_{ij}^2Xij2​ over sampled entries. These estimates feed the noncommutative Khintchine and spectral-norm concentration bounds.

Claim. Noncommutative Khintchine step from Section 6.1: conditional on the sampled row/column energies, the Rademacher-signed sampled matrix has the expected spectral moment at scale (qn/p) ∥X∥∞\sqrt(q n/p)\,\lVert X\rVert_{\infty}(​qn/p)∥X∥∞​.

Lecture-note formulation:

EΩ,ε∥p−1∑(i,j)∈ΩεijXijeiej⊤∥q≤(Cqnp∥X∥∞)q.\mathbb E_{\Omega,\varepsilon} \left\|p^{-1}\sum_{(i,j)\in\Omega}\varepsilon_{ij}X_{ij}e_ie_j^\top\right\|^q \le \left(C\sqrt{\frac{qn}{p}}\|X\|_\infty\right)^q .EΩ,ε​​p−1(i,j)∈Ω∑​εij​Xij​ei​ej⊤​​q≤(Cpqn​​∥X∥∞​)q.

The constants in this node are universal existential constants; the theorem asserts that some positive constants with these roles exist.

Decomposition status. A corresponding proof sketch reduces this node to smaller mathematical subclaims. The checked reduction uses 3 subclaims: Rademacher sampled matrix conditional Khintchine bound; Bernoulli Rademacher moment bound by conditional Khintchine scales; conditional Khintchine scale moment from row column energy moment.

Preamble
import Definitions.Def_matrix_completion_rademacher
open MatrixCompletion
Formal statement
theorem rademacher_sampled_matrix_moment_from_row_column_energy
    (Cenergy : ℝ) :
    0 < Cenergy →
    ∃ Crad : ℝ, 0 < Crad ∧
      ∀ (β : ℝ), 2 < β →
      ∀ (n₁ n₂ m q : ℕ) (X : Matrix (Fin n₁) (Fin n₂) ℝ),
        0 < n₁ → 0 < n₂ → m ≤ n₁ * n₂ →
        (m : ℝ) ≥ β * (↑(max n₁ n₂)) *
          Real.log (↑(max n₁ n₂)) →
        1 ≤ q →
        (q : ℝ) ≥ β * Real.log (↑(max n₁ n₂)) →
        (q : ℝ) ≤ 2 * (β * Real.log (↑(max n₁ n₂))) →
        (q : ℝ) ≤
          ((m : ℝ) / ((n₁ : ℝ) * (n₂ : ℝ))) *
            (↑(max n₁ n₂)) →
        bernoulliExpectation ((m : ℝ) / ((n₁ : ℝ) * (n₂ : ℝ)))
            (fun Omega =>
              (max (sampledRowEnergyMax Omega X)
                (sampledColumnEnergyMax Omega X)) ^ q) ≤
          (Cenergy * ((m : ℝ) / ((n₁ : ℝ) * (n₂ : ℝ))) *
            (↑(max n₁ n₂)) * entrySupNorm X ^ 2) ^ q →
        bernoulliExpectation ((m : ℝ) / ((n₁ : ℝ) * (n₂ : ℝ)))
            (fun Omega =>
              rademacherExpectation
                (fun eps =>
                  spectralNorm
                    (rademacherSampledMatrix Omega eps
                      ((m : ℝ) / ((n₁ : ℝ) * (n₂ : ℝ))) X) ^ q)) ≤
          (Crad * Real.sqrt
            (((q : ℝ) * (↑(max n₁ n₂))) /
              ((m : ℝ) / ((n₁ : ℝ) * (n₂ : ℝ)))) *
            entrySupNorm X) ^ q := by
  sorry
Source
Candes, Emmanuel, and Benjamin Recht. "Exact matrix completion via convex optimization." Communications of the ACM 55.6 (2012): 111-119.

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me