Case of Edelman's bound: is the integral (equality)
ProvedGaussianMatrix.wishart_lambda_min_cdf_rank_onechi-square-distributiongaussian-vectorpolar-coordinatesprobabilitysmallest-singular-value
Let and . Let be a standard Gaussian vector, viewed as a Gaussian matrix , so that . Then
The two integrands agree by Legendre's duplication formula.
This is the law of . It is exactly Edelman's density bound (wishart_lambda_min_cdf_density_bound) in the case , where the bound holds with equality. It serves as a proved sanity check that the constant in that statement is the right one; it is not needed by the reduction of wishart_lambda_min_tail.
Formalization Note. gaussianMatrix 1 k is the product law on Fin 1 → Fin k → ℝ. The right-hand side is a lintegral over . The constant is written exactly as the instance of the general constant , after simplifying and .
Preamble
import Definitions.Def_GaussianMatrix_basic open MeasureTheory ProbabilityTheory open scoped Matrix
Formal statement
namespace GaussianMatrix
theorem wishart_lambda_min_cdf_rank_one {k : ℕ} (hk : 1 ≤ k) (t : ℝ) (ht : 0 ≤ t) :
(gaussianMatrix 1 k) {G | sMin (Matrix.of G)ᵀ ^ 2 ≤ t}
= ∫⁻ x in Set.Ioc 0 t, ENNReal.ofReal
(2 ^ (((k : ℝ) - 2) / 2) * Real.Gamma (((k : ℝ) + 1) / 2)
/ (Real.Gamma (1 / 2) * Real.Gamma (k : ℝ))
* x ^ (((k : ℝ) - 2) / 2) * Real.exp (-x / 2)) := by sorry
end GaussianMatrixSource
standard fact: if g ~ N(0, I_k) then ‖g‖² ~ χ²_k, with density λ^{k/2−1}e^{−λ/2}/(2^{k/2}Γ(k/2)) on (0, ∞) (e.g. N. L. Johnson, S. Kotz, N. Balakrishnan, Continuous Univariate Distributions, Vol. 1, 2nd ed., Ch. 18). This is the case r = 1 of Edelman 1988, Prop. 5.1 / Tropp–Webber (B.5), where the bound is an equality.