Asymptotic expansion of Poisson averages
TheoremLogarithmic asymptotics of the square-root Poisson series#✓fully proved
For real \(x\) let \(S(x) = \sum_{k=1}^{\infty} x^k \sqrt{k} / k!\). Then, as \(x \to \infty\),
theorem ramanujan_part1_ch3_entry10_example1_sqrt_poisson_asymptotic :
(fun x : ℝ =>
Real.log (∑' (j : ℕ), x ^ (j + 1) * Real.sqrt (↑(j + 1) : ℝ) /
(↑(Nat.factorial (j + 1)) : ℝ)) - (x + (1 / 2 : ℝ) * Real.log x - 1 / (8 * x) - 1 /
(16 * x ^ 2))) =O[Filter.atTop] (fun x : ℝ =>
1 / x ^ 3)
Sources
- Berndt, Ramanujan's Notebooks Part I, Chapter 3, Entry 10, Example 1, printed p. 64 / PDF p. 74.
TheoremEventual summability of the square-root Poisson series#✓fully proved
For all sufficiently large real \(x\), the series \(\sum_{k=1}^{\infty} x^k \sqrt{k} / k!\) is summable.
theorem ramanujan_part1_ch3_entry10_example1_summable_eventually :
∀ᶠ (x : ℝ) in Filter.atTop,
Summable (fun j : ℕ => x ^ (j + 1) * Real.sqrt (↑(j + 1) : ℝ) /
(↑(Nat.factorial (j + 1)) : ℝ))
Sources
- Berndt, Ramanujan's Notebooks Part I, Chapter 3, Entry 10, Example 1, printed p. 64 / PDF p. 74.
Used by 1
TheoremEventual positivity of the square-root Poisson series#✓fully proved
For all sufficiently large real \(x\), the sum \(\sum_{k=1}^{\infty} x^k \sqrt{k} / k!\) is strictly positive.
theorem ramanujan_part1_ch3_entry10_example1_sum_positive_eventually :
∀ᶠ (x : ℝ) in Filter.atTop, 0 < ∑' (j : ℕ), x ^ (j + 1) * Real.sqrt (↑(j + 1)) /
(↑(Nat.factorial (j + 1)) : ℝ)
Sources
- Berndt, Ramanujan's Notebooks Part I, Chapter 3, Entry 10, Example 1, printed p. 64 / PDF p. 74.
Uses 1
TheoremAsymptotics of the logarithmic Poisson series#✓fully proved
There are real constants \(C > 0\) and \(x_0\) such that for every real \(x \ge x_0\) the series \(\sum_{k=1}^{\infty} x^k \log(k + 1) / k!\) has a sum \(S\) satisfying
theorem ramanujan_part1_ch3_entry10_example2_log_poisson_asymptotic :
∃ (C : ℝ) (x₀ : ℝ), 0 < C ∧ ∀ (x : ℝ), x₀ ≤ x → ∃ (S : ℝ),
HasSum (fun j : ℕ => x ^ (j + 1) * Real.log ((j + 2 : ℝ)) / (Nat.factorial (j + 1) : ℝ)) S ∧
|Real.exp (-x) * S - Real.log x - 1 / (2 * x) - 1 / (12 * x ^ 2)| ≤ C / x ^ 3
Sources
- Berndt, Ramanujan's Notebooks Part I, Chapter 3, Entry 10, Example 2, printed p. 65 / PDF p. 75.
TheoremSummability of the logarithmic Poisson series#✓fully proved
For every real \(x\), the series \(\sum_{k=1}^{\infty} e^{-x} x^k \log(k + 1) / k!\) is summable.
theorem ramanujan_part1_ch3_entry10_example2_summable :
∀ (x : ℝ),
Summable (fun j : ℕ => Real.exp (-x) * x ^ (j + 1) * Real.log ((j : ℝ) + 2) /
(Nat.factorial (j + 1) : ℝ))
Sources
- Berndt, Ramanujan's Notebooks Part I, Chapter 3, Entry 10, Example 2, printed p. 65 / PDF p. 75.
TheoremSummability of the harmonic Poisson series#✓fully proved
Let \(H_k = \sum_{j=1}^{k} 1/j\) be the \(k\)-th harmonic number (so \(H_0 = 0\)). For every real \(x\), the series \(\sum_{k=0}^{\infty} e^{-x} x^k H_k / k!\) is summable.
theorem ramanujan_part1_ch3_entry10_example4_summable (x : ℝ) :
Summable (fun k : ℕ => Real.exp (-x) * x ^ k *
(∑ j ∈ Finset.Icc (1 : ℕ) k, (1 : ℝ) / (j : ℝ)) / (Nat.factorial k : ℝ))
Sources
- Berndt, Ramanujan's Notebooks Part I, Chapter 3, Entry 10, Example 4, printed pp. 65-66 / PDF pp. 75-76.
Used by 1
TheoremAsymptotics of the harmonic Poisson series#✓fully proved
Let \(H_k = \sum_{j=1}^{k} 1/j\) be the \(k\)-th harmonic number (so \(H_0 = 0\)) and let \(\gamma\) be the Euler–Mascheroni constant. Then, as \(x \to \infty\),
theorem ramanujan_part1_ch3_entry10_example4_harmonic_poisson_asymptotic :
Asymptotics.IsBigO Filter.atTop (fun x : ℝ => Real.exp (-x) *
(∑' (k : ℕ),
(x ^ k * (∑ j ∈ Finset.Icc 1 k, (1 : ℝ) / (j : ℝ)) /
(k.factorial : ℝ))) - (Real.log x + Real.eulerMascheroniConstant)) (fun x : ℝ =>
(1 : ℝ) / x)
Sources
- Berndt, Ramanujan's Notebooks Part I, Chapter 3, Entry 10, Example 4, printed pp. 65-66 / PDF pp. 75-76.
Uses 1
Definition2-associated Stirling numbers in Ramanujan's Entry 10#✓fully proved
For natural numbers \(n\) and \(k\), \(S_2(n, k)\) is the integer \(\sum_{j=0}^{k} (-1)^j \binom{n}{j} S(n - j, k - j)\), where \(S\) denotes the Stirling numbers of the second kind.
\(S_2(n, k)\) is the number of partitions of an \(n\)-element set into \(k\) blocks, each containing at least two elements.
def twoAssocStirlingSecond (n k : ℕ) : ℤ :=
∑ j ∈ Finset.range (k + 1), ((-1 : ℤ) ^ j) * (Nat.choose n j : ℤ) *
(Nat.stirlingSecond (n - j) (k - j) : ℤ)
TheoremThe Entry 10 Stirling numbers are the 2-associated Stirling numbers#✓fully proved
For all natural numbers \(n\) and \(k\), the integer \(\sum_{j=0}^{k} (-1)^j \binom{n}{j} S(n - j, k - j)\) used in Ramanujan's Entry 10, where \(S\) denotes the Stirling numbers of the second kind, equals the 2-associated Stirling number of the second kind \(S_2(n, k)\).
theorem twoAssocStirlingSecond_eq_canonical (n k : ℕ) :
twoAssocStirlingSecond n k = MetaMathlibExt.twoAssocStirlingSecond n k
def twoAssocStirlingSecond (n k : ℕ) : ℤ :=
∑ j ∈ Finset.range (k + 1), ((-1 : ℤ) ^ j) * (Nat.choose n j : ℤ) *
(Nat.stirlingSecond (n - j) (k - j) : ℤ)
MathlibExt/NumberTheory/Ramanujan/Part1Ch3Entry10PoissonAsymptotic.lean:53
def twoAssocStirlingSecond (n k : ℕ) : ℤ :=
∑ j ∈ Finset.range (k + 1), ((-1 : ℤ) ^ j) * (Nat.choose n j : ℤ) *
(Nat.stirlingSecond (n - j) (k - j) : ℤ)
MathlibExt/Combinatorics/AssociatedStirlingSecond.lean:30
Uses 2 · Used by 1
TheoremRamanujan's asymptotic expansion of Poisson averages#✓fully proved
Let \(M \ge 1\) be an integer, \(\varphi : \mathbb{R} \to \mathbb{C}\) and \(G : \mathbb{R} \to \mathbb{R}\) functions, \(A \ge 1\) a real number and \(F\) a finite set of natural numbers. Write \(\varphi^{(n)}\) for the \(n\)-th iterated derivative of \(\varphi\) and \(S_2(n, k)\) for the 2-associated Stirling numbers of the second kind. Assume that there are \(C \ge 0\) and \(p \in \mathbb{N}\) such that \(|\varphi(x)| \le C x^p\) and \(|G(x)| \le C x^p\) for all sufficiently large \(x\); that \(G(x) \ge 1\) for all sufficiently large \(x\); and that for every integer \(m \ge 1\) and all sufficiently large \(x\), the derivatives \(\varphi^{(0)}, \dots, \varphi^{(m-1)}\) are differentiable at \(x\) and \(|\varphi^{(m)}(x)| / m! \le G(x) (A / x)^m\). Then, as \(x \to \infty\),
Proof: Taylor expansion on the central Poisson range is combined with the associated-Stirling central-moment identity from Comtet; Berndt's route and even-moment tails control the remainder.
MathlibExt.NumberTheory.Ramanujan.Part1Ch3.Entry10PoissonAsymptotic.ramanujan_part1_ch3_entry10_poisson_asymptotic · 19 lines
theorem ramanujan_part1_ch3_entry10_poisson_asymptotic
(M : ℕ) (hM : 0 < M)
(φ : ℝ → ℂ) (G : ℝ → ℝ) (A : ℝ) (hA : 1 ≤ A)
(F : Finset ℕ)
(h_poly : ∃ C : ℝ, 0 ≤ C ∧ ∃ p : ℕ, ∀ᶠ x in Filter.atTop, ‖φ x‖ ≤ C * x ^ p ∧ ‖G x‖ ≤ C * x ^ p)
(hG_ge : ∀ᶠ x in Filter.atTop, (1 : ℝ) ≤ G x)
(h_deriv : ∀ m : ℕ, 0 < m → ∀ᶠ x in Filter.atTop,
(∀ k ∈ Finset.Icc 1 m, DifferentiableAt ℝ (iteratedDeriv (k - 1) φ) x) ∧
‖iteratedDeriv m φ x / (Nat.factorial m : ℂ)‖ ≤ G x * (A / x) ^ m) :
Asymptotics.IsBigO Filter.atTop
(fun x : ℝ =>
((Real.exp (-x) : ℂ)) *
(∑' j : ℕ, if j ∈ F then (0 : ℂ) else
(((x ^ j / (Nat.factorial j : ℝ) : ℝ) : ℂ) * φ (j : ℝ)))
- φ x
- ∑ k ∈ Finset.Icc 2 M, ∑ n ∈ Finset.Icc k (2 * k - 2),
(((twoAssocStirlingSecond n (n + 1 - k) : ℤ) : ℂ) * ((x : ℂ) ^ (n - k + 1)) *
(iteratedDeriv n φ x / (Nat.factorial n : ℂ))))
(fun x : ℝ => G x / (x : ℝ) ^ M)
def twoAssocStirlingSecond (n k : ℕ) : ℤ :=
∑ j ∈ Finset.range (k + 1), ((-1 : ℤ) ^ j) * (Nat.choose n j : ℤ) *
(Nat.stirlingSecond (n - j) (k - j) : ℤ)
MathlibExt/NumberTheory/Ramanujan/Part1Ch3Entry10PoissonAsymptotic.lean:53
Sources
- Berndt, Ramanujan's Notebooks Part I, Chapter 3, Entry 10, printed pp. 57-58 / PDF pp. 67-68.
TheoremEventual summability of the series in Entry 10#✓fully proved
Let \(\varphi : \mathbb{R} \to \mathbb{C}\) be polynomially bounded at infinity, that is, there are a real number \(C \ge 0\) and a natural number \(p\) such that \(|\varphi(x)| \le C x^p\) for all sufficiently large real \(x\). Let \(F\) be a finite set of natural numbers. Then for all sufficiently large real \(x\), the family indexed by the natural numbers \(j\) whose \(j\)-th term is \(0\) for \(j \in F\) and \(x^j \varphi(j) / j!\) for \(j \notin F\) is summable.
theorem ramanujan_part1_ch3_entry10_poisson_summable
(φ : ℝ → ℂ) (F : Finset ℕ)
(hφ : ∃ (C : ℝ) (p : ℕ), 0 ≤ C ∧ ∀ᶠ (x : ℝ) in Filter.atTop, ‖φ x‖ ≤ C * x ^ p) :
∀ᶠ (x : ℝ) in Filter.atTop,
Summable (fun j : ℕ => if j ∈ F then 0 else ((x : ℂ) ^ j * φ (j : ℝ) / (Nat.factorial j : ℂ)))
Sources
- Berndt, Ramanujan's Notebooks Part I, Chapter 3, Entry 10, printed pp. 57-58 / PDF pp. 67-68.
Loading dependency explorer…
The dependency explorer needs JavaScript and its graph data.