Independent equal-block Gaussian limit in the bounded mixing case
ProvedMarkovChainCLT.exists_independent_blocks_of_bounded_of_summable_alphaLet be a measurable, centered, strictly stationary real sequence on a probability space with measure . Suppose almost surely for every and . Write
Assume . There exist integer block lengths and gap lengths such that, setting ,
and, for every real ,
The final limit concerns a triangular array of independent copies of the big-block sums. Together with the Bernstein comparison estimate, these conclusions give the bounded strongly mixing CLT. The statement isolates the remaining block-schedule and independent-array limit argument; it does not assume independence of the original sequence.
Formalization Note Independent copies are represented through the power of their common characteristic function, so no auxiliary product probability space is needed. The displayed limits concern ; at the inverse convention and zero exponent are harmless.
import Definitions.Def_MixingCoefficients import Mathlib open MeasureTheory ProbabilityTheory Filter MarkovChainCLT open scoped ENNReal NNReal Topology ProbabilityTheory
theorem MarkovChainCLT.exists_independent_blocks_of_bounded_of_summable_alpha
{Ω : Type*} [MeasurableSpace Ω]
(P : Measure Ω) [IsProbabilityMeasure P] (Y : ℕ → Ω → ℝ)
(hY : ∀ n, Measurable (Y n)) (hstat : IsStrictlyStationary P Y)
(hcent : ∫ ω, Y 0 ω ∂P = 0)
(B : ℝ) (hB : ∀ n, ∀ᵐ ω ∂P, |Y n ω| < B)
(hα : Summable (fun n => alphaMixingCoef P Y n))
(hvarlim : Tendsto (fun n : ℕ => Var[∑ i ∈ Finset.range n, Y i; P] / (n : ℝ))
atTop (𝓝 (seqAsymptoticVariance P Y)))
(hvar : 0 < seqAsymptoticVariance P Y) :
∃ p q : ℕ → ℕ,
(∀ n, 0 < p n) ∧
Tendsto p atTop atTop ∧ Tendsto q atTop atTop ∧
Tendsto (fun n : ℕ => (p n : ℝ) / n) atTop (𝓝 0) ∧
Tendsto (fun n : ℕ => ((n - (n / (p n + q n)) * p n : ℕ) : ℝ) / n)
atTop (𝓝 0) ∧
Tendsto (fun n : ℕ => ((n / (p n + q n) - 1 : ℕ) : ℝ) *
alphaMixingCoef P Y (q n + 1)) atTop (𝓝 0) ∧
∀ t : ℝ, Tendsto (fun n : ℕ =>
(charFun (P.map (fun ω => (Real.sqrt n)⁻¹ * ∑ i ∈ Finset.range (p n), Y i ω)) t) ^
(n / (p n + q n)))
atTop (𝓝 (charFun (gaussianReal 0 (seqAsymptoticVariance P Y).toNNReal) t)) := by sorry