Theorem 6.1, proof — Hoeffding tail bound for at every
ProvedTraceEstimation.Rayleigh.hoeffding_tailLet be a nonzero symmetric positive semi-definite matrix, , and let be a normalized Rayleigh-quotient trace estimator of : are independent random vectors on a probability space with almost surely and . Let be the ratio between the largest and smallest nonzero eigenvalue of . Then for every ,
This is Hoeffding's inequality for the independent summands , each confined to an interval of length .
Formalization Note The model is a general probability space with random vectors z : Fin M → Ω → Fin n → ℝ satisfying IsNormalizedRayleighSample P A z (measurable, mutually independent, not necessarily identically distributed, normalized almost surely, unbiased for this ). The hypotheses are satisfiable, for instance by with independent uniform indices . The denominator is positive under the hypotheses (), so no division by zero occurs.
import Mathlib import Definitions.Def_TraceEstimation_Rayleigh_kappaF import Definitions.Def_TraceEstimation_Rayleigh_rayleighEstimator
namespace TraceEstimation.Rayleigh
open MeasureTheory ProbabilityTheory Matrix
/-- Avron–Toledo, proof of Theorem 6.1 (p. 8:10), third display (Hoeffding's inequality): let
`A ≠ 0` be symmetric positive semi-definite and `R_M` a normalized Rayleigh-quotient trace
estimator of `A` with `M ≥ 1` samples (Definition 3.2: independent `z_i` with `z_iᵀz_i = n` a.s.
and `E(z_iᵀAz_i) = trace(A)`). For every `t > 0`,
`Pr(|R_M − trace(A)| ≥ t) ≤ 2 exp(−2M² rank²(A) t² / (M n² trace²(A) κ_f²(A)))`. -/
theorem hoeffding_tail {Ω : Type*} [MeasurableSpace Ω] (P : Measure Ω) [IsProbabilityMeasure P]
{n M : ℕ} (A : Matrix (Fin n) (Fin n) ℝ) (hA : A.PosSemidef) (hA0 : A ≠ 0) (hM : 0 < M)
(z : Fin M → Ω → Fin n → ℝ) (hz : IsNormalizedRayleighSample P A z) (t : ℝ) (ht : 0 < t) :
P.real {ω | t ≤ |rayleighEstimator A z ω - A.trace|} ≤
2 * Real.exp (-(2 * (M : ℝ) ^ 2 * (A.rank : ℝ) ^ 2 * t ^ 2) /
((M : ℝ) * (n : ℝ) ^ 2 * A.trace ^ 2 * kappaF hA.isHermitian ^ 2)) := by sorry
end TraceEstimation.Rayleigh
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.