ベルヌーイ数
ベルヌーイ数 (ベルヌーイすう、Bernoulli number) は数論における基本的な係数を与える数列であり、もともと、連続する整数のべき乗和を定式化する際の展開係数として1713年にヤコブ・ベルヌーイが著書 Ars Conjectandi (推測術) にて導入した[1]ことからこの名称がついた。 ベルヌーイ数は、べき乗和の展開係数にとどまらず、級数展開の係数や剰余項、リーマンゼータ関数においても登場する。また、ベルヌーイ数はすべてが有理数である。
定義
ベルヌーイ数 <math>B_n</math> は下に示すマクローリン展開 (テイラー展開) の展開係数として定義される。
- <math>
f(x) = \frac{x}{e^x-1} = \sum_{n=0}^{\infty} \frac{B_n}{n!} x^n
</math> この定義から、関数 <math>x/(e^x - 1)</math> を繰り返し微分していけばベルヌーイ数を得ることができるが、そのような手段でベルヌーイ数を得るのは容易ではない。 ベルヌーイ数を計算するには、マクローリン展開ではなく、次の漸化式を用いる。この漸化式から、ベルヌーイ数がすべて有理数であることがわかる。
- <math>
B_0 = 1,\quad{}B_n = -{1 \over n+1}\sum^{n-1}_{k=0}{n+1 \choose k}B_k
</math> ここで、<math>\left({\textstyle{n+1\atop k}}\right)</math>は二項係数 (二項定理 参照) である。 上の漸化式を用いて、ベルヌーイ数を第 29 項までを算出すると下表のようになる。 この表には、ベルヌーイ数が有理数であるので、分子と分母を記載している。
n | 分子 | 分母 | n | 分子 | 分母 | n | 分子 | 分母 | ||
---|---|---|---|---|---|---|---|---|---|---|
0 | 1 | 1 | 10 | 5 | 66 | 20 | −174611 | 330 | ||
1 | −1 | 2 | 11 | 0 | — | 21 | 0 | — | ||
2 | 1 | 6 | 12 | −691 | 2 730 | 22 | 85 4513 | 138 | ||
3 | 0 | — | 13 | 0 | — | 23 | 0 | — | ||
4 | −1 | 30 | 14 | 7 | 6 | 24 | −236 364 091 | 2 730 | ||
5 | 0 | — | 15 | 0 | — | 25 | 0 | — | ||
6 | 1 | 42 | 16 | −3 617 | 510 | 26 | 8 553 103 | 6 | ||
7 | 0 | — | 17 | 0 | — | 27 | 0 | — | ||
8 | −1 | 30 | 18 | 43 867 | 798 | 28 | −23 749 461 029 | 870 | ||
9 | 0 | — | 19 | 0 | — | 29 | 0 | — |
ベルヌーイ数の漸化式は、上記の関数 <math>f(x)=x/(e^x-1)</math> の逆数をテイラー展開し、その 2 つの積が 1 になることから導出できる。その漸化式は厳密な計算には有用であるが、<math>n</math> が大きくなると途中の式の値が非常に大きくなるため、浮動小数点数を使って計算する場合、精度が著しく悪くなる計算として知られている。
奇数番目のベルヌーイ数は <math>B_1</math> 以外はすべて 0 であり、偶数番目は <math>B_0</math> を除いて正の数と負の数が交互に並ぶ。 ベルヌーイ数の第 3 項以降の奇数項が 0 となることは、正接関数 (tangent) のマクローリン展開から証明[2]できる。
ベルヌーイ数の一般項
第 2 種スターリング数との関係 (スターリング数 参照) から、次のようなベルヌーイ数の一般項を算出する公式が存在する。
- <math>
B_n = \sum_{j=0}^n (-1)^j\,j^n\sum_{m=j}^n \frac{1}{m+1}{m\choose j}
</math> この公式は、総和記号が二重になっているため、上に示した漸化式ほど手軽にベルヌーイ数を計算する公式ではない。
漸近的性質
ベルヌーイ数とリーマンゼータ関数の関係から、ベルヌーイ数に次の漸近的性質があることが導かれる。
- <math>
|B_{2n}| \simeq \frac{2\,(2n)!}{(2\pi)^{2n}}\quad\quad(n\rightarrow\infty)
</math>
ベルヌーイ数を用いた級数展開
ベルヌーイ数は、いくつかの双曲線関数と三角関数の級数展開における展開係数となる。 ベルヌーイ数を展開係数とする関数とそのローラン級数による表現を挙げる。 まず、余接関数 (cotangent) のローラン級数展開は次のようになる。
- <math>
\begin{align}
& \coth z = \frac1z + \sum_{k=1}^{\infty} \frac{2^{2k}B_{2k}}{(2k)!}z^{2k-1},\\ & \cot z = \frac1z + \sum_{k=1}^{\infty} (-1)^k \frac{2^{2k}B_{2k}}{(2k)!}z^{2k-1},\\
\end{align} </math> 第 1 の関係式は、ベルヌーイ数が <math>x/(e^x-1)</math> の展開係数であることを利用して数式変形すれば得られる。 第 2 の関係式は <math>\cot z = i\,\coth z</math> であることを利用すれば、第 1 の関係式から導き出される。これらの級数の収束半径は <math>|z| < \pi</math> である。 次に正接関数 (tangent) のローラン級数展開は次のようになる。
- <math>
\tan z = \sum_{k=1}^{\infty} \frac{(-1)^k\,(2^{2k}-4^{2k})\,B_{2k}}{(2k)!}z^{2k-1}
</math> この関係式は、<math>\tan z = \cot z - 2\cot 2z</math> を利用して余接関数のローラン級数展開を変形すれば導出できる。 なお、この級数の収束半径は <math>|z| < \pi/2</math> である。 この正接関数のローラン級数展開の展開係数による数列はタンジェント数と呼ばれる。 一方、余割関数 (cosecant) は次のようにローラン級数展開される。
- <math>
\csc z = \frac{1}{\sin z} = \frac{1}{z} +\sum_{k=1}^\infin \frac{(-1)^k\,(2-2^{2k})\,B_{2k}}{(2k)!}z^{2k-1}
</math> この関係式は、<math>\csc 2z = (\tan z + \cot z)/2</math> を利用すれば導出できる。 なお、この級数の収束半径は<math>|z| < \pi</math> である。
べき乗和による導入
ベルヌーイ数は、もともと、連続する整数のべき乗和を定式化する際に、展開係数として導入された。 現代の表記法によって書くならば、定式化するべき乗和とは、
- <math>
S_k(n) \equiv \sum_{j=0}^{n-1} j^k = 0^k + 1^k + 2^k + \cdots + (n-1)^k
</math> なる級数である。この級数は、ベルヌーイ数を用いて、
- <math>
S_k(n) = \frac{1}{k+1} \sum_{j=0}^k {k+1\choose j}B_j\,n^{k-j+1}
</math> のように書くことができる。 ベルヌーイ数の漸化式は、べき乗和を定式化した際の考察から得られる。 さらに、ベルヌーイ数の指数型母関数 (母関数 参照) が <math>x/(e^x - 1)</math> となることから、その母関数を現在ではベルヌーイ数の定義とする。
ヤコブ・ベルヌーイは彼の著書『推測術』でベルヌーイ数を導入した際、べき乗和を上に書いたような <math>0</math> から <math>n-1</math> にわたる和でなく、<math>1</math> から <math>n</math> にわたる和:
- <math>
\hat{S}_k(n) \equiv \sum_{j=1}^{n} j^k = 1^k + 2^k + 3^k + \cdots + n^k
</math> として扱っていた。 ベルヌーイは、その著書で整数のべき乗 <math>n^c</math> の和を計算する公式として、次の数式を記している[3]。
- <math>
\begin{align}
\int n^c & = {} \frac{1}{c+1}n^{c+1} + \frac{1}2n^c + \frac{c}2 An^{c-1} + \frac{c.c-1.c-2}{2.3.4}Bn^{c-3}\\ & \quad\quad\quad{} + \frac{c.c-1.c-2.c-3.c-4}{2.3.4.5.6}Cn^{c-5} \\ & \quad\quad\quad{} + \frac{c.c-1.c-2.c-3.c-4.c-5.c-6}{2.3.4.5.6.7.8}Dn^{c-7}+ \ldots\ldots
\end{align} </math> この数式に記載されている展開係数 <math>A, B, C, D, \ldots</math> がベルヌーイ数 (<math>B_2</math> 以降) である。ベルヌーイが記した数式は、
- <math>
\hat{S}_k(n) = \frac{1}{k+1} \sum_{j=0}^k {k+1\choose j}\hat{B}_j\,n^{k-j+1}
</math> に相当する。この数式に用いた展開係数 <math>\hat{B}_j</math> は、
- <math>
\hat{B}_1 = 1/2\,(= -B_1),\quad \hat{B}_j = B_j\quad(j\neq1)
</math> のように、<math>j\neq1</math> においてベルヌーイ数と一致する。一部の文献[4][5]では <math>B_j</math> の代わりに <math>\hat{B}_j</math> をベルヌーイ数と呼んでいる。
一方、日本ではベルヌーイとほぼ同時期に関孝和がべき乗和を定式化し、ベルヌーイ数を発見していた[6]。 そのため、ベルヌーイ数を関・ベルヌーイ数と書いている文献[7]もある。
一般化されたベルヌーイ数
一般化されたベルヌーイ数は代数的数で、ベルヌーイ数がリーマンゼータ函数の特殊値に関連する方法と同じ方法で、ディリクレのL-函数の特殊値に関連して定義される。
χ を modulo f のディリクレ指標とすると、一般化されたベルヌーイ数は、
- <math>\sum_{a=1}^f\chi(a)\frac{te^{at}}{e^{ft}-1}=\sum_{k=0}^\infty B_{k,\chi}\frac{t^k}{k!}</math>
により定義される。B1,1=1/2 を除き、 χ(-1) ≠ (-1)k であれば、任意のディリクレ指標 χ に対し、Bk,χ = 0 である。
正でない整数でのリーマンゼータ函数の値とベルヌーイ数の間の関係を一般化し、全ての k ≥ 1 に対し、
- <math>L(1-k,\chi)=-\frac{B_{k,\chi}}{k}</math>
が成り立つ。ここに L(s, χ) は χ のディリクレのL-函数である。[8]
関連項目
脚注
- ↑ Julian Havil, "オイラーの定数 ガンマ," 新妻弘 訳, 共立出版, 初版, pp.97-99, 2009.
- ↑ 正接関数のマクローリン展開の結果において、実数変数を仮定した場合、ベルヌーイ数の第 3 項以降の奇数項は虚数項に対応する。 実数変数における正接関数が実数関数でなければならないので、そのマクローリン展開に虚数項に対応する項が存在してはならない。 よって、ベルヌーイ数の第3項以降の奇数項はゼロでなければならない。
- ↑ E. Hairer, G. Wanner, "解析教程 上," 蟹江幸博 訳, シュプリンガー・ジャパン, 新装版, p.18, 2006.
- ↑ 例えば、 荒木恒男, 伊吹山知義, 金子昌信, "ベルヌーイ数とゼータ関数," 牧野書店, 2001.
- ↑ Wikipedia ファウルハーバーの公式 もベルヌーイの記述に基づき、第 1 項を1/2とする記述で説明している。
- ↑ 小川束, "関孝和によるベルヌーイ数の発見," 数理解析研究所講究録, 第1583巻, 2008.
- ↑ 例えば、 桜井進, 中村義作, "天才たちが愛した美しい数式," PHP研究所, 第1版, p.205, 2008.
- ↑ テンプレート:Harvnb