1
大学数学基礎解説
文献あり

ライプニッツ級数の連分数による加速

184
0
$$\newcommand{acoloneqq}[0]{\ &\hspace-2pt\coloneqq} \newcommand{ar}[1]{\operatorname{ar{#1}}} \newcommand{arc}[1]{\operatorname{arc{#1}}} \newcommand{asupplement}[1]{&\hspace{#1}\textsf} \newcommand{beginend}[2]{{\begin{#1}#2\end{#1}}} \newcommand{bm}[0]{\boldsymbol} \newcommand{bscolor}[1]{\color{var(--bs-#1)}} \newcommand{bscolorbox}[1]{\colorbox{var(--bs-#1)}} \newcommand{bsrowcolor}[1]{\rowcolor{var(--bs-#1)}} \newcommand{C}[0]{\mathbb C} \newcommand{Defarrow}[0]{\xLeftrightarrow{\textrm{def}}} \newcommand{DefOp}[1]{\DeclareMathOperator{#1}{#1}} \newcommand{fqty}[0]{\!\qty} \newcommand{hen}[1]{{(\textrm{{#1}辺})}} \newcommand{hygeo}[6]{{{}_{#1}{#2}_{#3}{\qty[\beginend{matrix*}{#4\\ #5}\ ;{#6}]}}} \newcommand{In}[0]{\in\mathbb} \newcommand{kfrac}[0]{\mathop{\Large\raise-.8pt{\textrm K}}} \newcommand{Kfrac}[0]{\mathop{\huge\raise-2.2pt{\textrm K}}} \newcommand{kome}[0]{\textreferencemark} \newcommand{leftshiftarrow}[0]{{\substack{\curvearrowright\\ \leftharpoondown}}} \newcommand{lr}[3]{{\left#1{#2}\right#3}} \newcommand{lvvr}[3]{\lr{#1}{\negmedspace\lr|{#2}|\negmedspace}{#3}} \newcommand{N}[0]{\mathbb N} \newcommand{ot}[0]{\leftarrow} \newcommand{P}[0]{\mathbb P} \newcommand{Q}[0]{\mathbb Q} \newcommand{R}[0]{\mathbb R} \newcommand{RANGE}[0]{}\newcommand{rangeex}[6][,]{{#2{#3}_{#5}#4#1\cdots#1#2{#3}_{#6}#4}}\newcommand{range}[2][,]{\rangeex[#1]{}{#2}{}}{} \newcommand{REQUIRE}[0]{}\require{physics}{} \newcommand{rightshiftarrow}[0]{{\substack{\curvearrowleft\\ \rightharpoondown}}} \newcommand{rmIm}[0]{\operatorname{Im}} \newcommand{rmRe}[0]{\operatorname{Re}} \newcommand{rprod}[0]{\mathop{\prod\!\llap\coprod}} \newcommand{sahen}[0]{\hen左} \newcommand{STIRLING}[0]{}\newcommand{stirling}[3][]{{\qty[\beginend{matrix}{{#2}\\ {#3}}{#1}]}}\newcommand{Stirling}[3][]{{\qty{\beginend{matrix}{{#2}\\ {#3}}{#1}}}}{} \newcommand{uhen}[0]{\hen右} \newcommand{vbin}[1]{\mathbin{{#1}\!\llap|\ }} \newcommand{Z}[0]{\mathbb Z} $$
$ \beginend{array}{l N\ge1 \\ \beginend{aligned}{ \frac\pi4 &= \sum_{n=0}^{N-1} \frac{(-1)^n}{2n+1} + \frac{(-1)^N}2\Kfrac_{n=0}^\infty \frac{\delta_{n,0}+n^2}{2N} \\&= \sum_{n=0}^{N-1} \frac{(-1)^n}{2n+1} + \cfrac{\lr.{(-1)^N}/2} {2N+\cfrac{1^2} {2N+\cfrac{2^2} {2N+\cfrac{3^2} {2N+\ddots}}}} } \\ N\ge0 \\ \beginend{aligned}{ \ln2 &= \sum_{n=1}^N \frac{(-1)^{n-1}}n + (-1)^N\Kfrac_{n=0}^\infty \frac{\delta_{n,0}+n^2}{2N+1} \\&= \sum_{n=1}^N \frac{(-1)^{n-1}}n + \cfrac{(-1)^N} {2N+1+\cfrac{1^2} {2N+1+\cfrac{2^2} {2N+1+\cfrac{3^2} {2N+1+\ddots}}}} }}$

定義

$\delta_{a,b}$ クロネッカーのデルタ
$\displaystyle \sum_{n=0}^{-1} f(n) \coloneqq 0$ 空和
$\displaystyle \Kfrac_{n=0}^\infty \frac{a_n}{b_n}$ 連分数
$B_n(z),E_n(z)$ ベルヌーイ多項式、オイラー多項式

一般的な式

以前書いた『 ゼータ関数を加速させよ 』でこのような式を導出しました。
$\displaystyle \sum_{n=0}^\infty \frac{(-1)^n}{z+n} = \sum_{n=0}^{N-1} \frac{(-1)^n}{z+n} + (-1)^N\sum_{n=0}^\infty E_{2n}\cdot(2z+2N-1)^{-2n-1}$
右辺の2つ目の級数は発散級数で、$N$を大きくし分母を大きくした上で
適切な項数の部分和をとることにより一つ目の部分和の誤差が打ち消されるというものです。
この発散級数は 連分数展開 できます。
数値計算によればこの連分数は恐らく収束します。
$N$を大きくしていくことで連分数の収束は急速に速まるため、
数値計算に有用です。
$\displaystyle 4\qty(\frac11-\frac13+\cdots-\frac1{15} + \cfrac{\frac12} {16+\cfrac{1^2} {16+\cfrac{2^2} {\phantom{16+{}}\cfrac{\small\ddots} {16+\cfrac{6^2} {16}}}}}) = 3.141592653\bscolor{secondary}60\cdots$

$N\ge0$
$\beginend{align}{ \sum_{n=0}^\infty \frac{(-1)^n}{z+n} &= \sum_{n=0}^{N-1} \frac{(-1)^n}{z+n} + \frac{(-1)^N}2\Kfrac_{n=0}^\infty \frac{\delta_{n,0}+n^2}{2z+2N-1} \\&= \sum_{n=0}^{N-1} \frac{(-1)^n}{z+n} + \cfrac{(-1)^N} {2z+2N-1+\cfrac{1^2} {2z+2N-1+\cfrac{2^2} {2z+2N-1+\cfrac{3^2} {2z+2N-1+\scriptsize\ddots}}}} }$

オイラーの定数

先程の式の左辺$\displaystyle \sum_{n=0}^\infty \frac{(-1)^n}{z+n}$$ \dfrac12\qty(\psi\fqty(\dfrac{z+1}2)-\psi\fqty(\dfrac z2))$と等しく、
ディガンマ関数が出現しますが、打ち消し合うためここからオイラーの定数$\gamma$を導くことはできません。
『ゼータ関数を加速させよ』で導いた式により、ディガンマ関数は次のように表されます。
$\displaystyle \psi(z) = \ln(z-a+N) - \sum_{n=0}^{N-1} \frac1{z+n} + \sum_{n=1}^\infty \frac{(-1)^{n-1}B_n(a)}n(z-a+N)^{-n}$
よって、$\displaystyle \sum_{n=1}^\infty \frac{B_n(a)}nz^n$を連分数展開できればオイラーの定数を効率的に求めることができそうです。
また、別のアプローチとして 指数積分 ガウスの連分数 ${}_2F_0$への一般化を用いる方法があります。

$ \beginend{align}{ x&>0\ ? \\ \gamma &= -\sum_{n=1}^\infty \frac{(-x)^n}{n\cdot n!} - \ln x - \cfrac{e^{-x}} {x+\cfrac1 {1+\cfrac1 {x+\cfrac2 {1+\cfrac2 {x+\cfrac3 {1+\cfrac3 {x+\scriptsize\ddots}}}}}}} \\&= -\sum_{n=1}^\infty \frac{(-x)^n}{n\cdot n!} - \ln x - e^{-x}\Kfrac_{n=0}^\infty \frac{\delta_{n,0}-n^2}{2n+1+x} \\&= -\sum_{n=1}^\infty \frac{(-x)^n}{n\cdot n!} - \ln x - \cfrac{e^{-x}} {1+x-\cfrac{1^2} {3+x-\cfrac{2^2} {5+x-\cfrac{3^2} {7+x-\ddots}}}} }$

1つ目の連分数からの 式変形 により2つ目の連分数が得られます。
$x$を大きくすれば連分数は速くなりますがその分級数が遅くなります。

参考文献

投稿日:20231130
更新日:823
OptHub AI Competition

この記事を高評価した人

高評価したユーザはいません

この記事に送られたバッジ

バッジはありません。

投稿者

著者の記事における命題は大半が自分で発見したものであり、 何かしらの論文などに基づいたものではありません。

コメント

他の人のコメント

コメントはありません。
読み込み中...
読み込み中