調和数 (発散列)

出典: フリー百科事典『ウィキペディア(Wikipedia)』
移動: 案内検索
オアの調和数調和数の項を参照。また調和数の列は調和数列とは異なる。
n = ⌊x⌋ に対する調和数 Hn,1 のグラフ(赤)。これは γ + ln(x)(青)に漸近収斂する。

数学において、n-番目の調和数(ちょうわすう、: harmonic number)は 1 から n までの自然数逆数和

H_n= 1+\frac{1}{2}+\frac{1}{3}+\cdots+\frac{1}{n}=\sum_{k=1}^n \frac{1}{k}

である。これはまた、1 から n までの自然数の調和平均の逆数の n-倍に等しい。

調和数は遥か昔から研究され、数論の各分野において重要である。調和数の極限は、調和級数と呼ばれ(しばしば調和数もひっくるめて一口に調和級数と呼ぶこともある)、リーマンゼータ函数と近しい関係にあり、また種々の特殊函数のさまざまな表示に現れる。

十分大きな数の標本について、その出現頻度がジップの法則に従って分布するとき、全体の中で n-番目の頻度で現れる標本の総頻度は n-番目の調和数である。このことは長い尻尾およびネットワーク値英語版の驚くべき帰結の一種を導く。

調和数の計算法[編集]

調和数の積分表示

 H_n = \int_0^1 \frac{1 - x^n}{1 - x}\,dx

オイラーによる。この等式は簡単な代数的等式

\frac{1-x^n}{1-x}=1+x+\ldots +x^{n-1}

を使えば明らかである。また、積分の変数を単純に x = 1 − u と変換すれば、Hn のきれいな組合せ論的展開

 \begin{align} H_n 
 & = \int_{0}^{1} \frac{1-(1-u)^n}{u}du = \int_{0}^{1} \left[\sum_{k=1}^n(-1)^{k-1}\binom{n}{k} u^{k-1}\right]du\\
 & =\sum_{k=1}^{n} (-1)^{k-1}\binom{n}{k} \int_{0}^{1} u^{k-1}du = \sum_{k=1}^{n} (-1)^{k-1}\frac{1}{k}\binom{n}{k}
\end{align}

が得られる。同じ表現は、第三レトケシュ恒等式英語版x1 = 1, ..., xn = 1 とおき、

\Pi_k(1,\ldots,n)=(-1)^{n-k}(k-1)!(n-k)!

なる事実を用いることでも得られる。すなわち

H_n=H_{n,1}=\sum_{k=1}^n\frac{1}{k}=(-1)^{n-1}n!\sum_{k=1}^n\frac{1}{k^2\Pi_k(1,\ldots,n)}=\sum_{k=1}^n(-1)^{k-1}\frac{1}{k}\binom nk

が成り立つ。また、レトケシュ恒等式を x1 = 12, ..., xn = n2 に対して用いれば、この場合

\Pi_k(1^2,2^2,\ldots,n^2)=(-1)^{n-k}\frac{(n-k)!(n+k)!}{2k^2}

となるので、ζ(2) の第 n-部分和についての類似の公式

H_{n,2}=\sum_{k=1}^n\frac{1}{k^2}=2\sum_{k=1}^n(-1)^{k-1}\frac{1}{k^2}\frac{\binom nk}{\binom {n+k} k}

を得る。Hn の増大度は n自然対数 ln(n) と同程度の速さである。このことは、Hn を積分

\int_1^n {dx \over x}\quad(= \ln(n))

で近似することによって確認できる。数列 (Hn - ln(n)) は単調に減少して、

 \lim_{n \to \infty} H_n - \ln(n) = \gamma

なる定数(この定数 γ はオイラー・マスケローニ定数と呼ばれ、その値は 0.5772156649... である)を極限にもち、これに対応する漸近展開

H_n = \ln{n} + \gamma + \frac{1}{2}n^{-1} - \frac{1}{12}n^{-2} + \frac{1}{120}n^{-4} + \mathcal{O}(n^{-6})

で与えられる。

分数パラメータに対する特殊値[編集]

調和数 Hn のパラメータ n を積分

H_\alpha = \int_0^1\frac{1-x^\alpha}{1-x}dx

によって拡張すれば、0 と 1 の間の分数値をもつパラメータ α に対する解析的な特殊値を定めることができる。あるいはさらに漸化式

 H_\alpha = H_{\alpha-1}+\frac{1}{\alpha}

によって拡張することもでき、結局は任意の x > 0 に対して(x が整数のときもそうでないときも)

 H_{x} =  x \sum_{k=1}^\infty \frac{1}{k(x+k)}

が成立する。いくつかの特殊値について計算すれば、以下のようになる。

\begin{align} 
 H_{3/4} & = \frac{4}{3}-3\ln{2}+\frac{\pi}{2},\\
 H_{2/3} & = \frac{3}{2}(1-\ln{3})+\frac{\pi}{6}\sqrt{3},\\
 H_{1/2} & = 2 -2\ln{2},\\
 H_{1/3} & = 3-\frac{\pi}{2\sqrt{3}} -\frac{3}{2}\ln{3},\\ 
 H_{1/4} & = 4-\frac{\pi}{2} - 3\ln{2},\\
 H_{1/6} & = 6-\frac{\pi}{2}\sqrt{3} -2\ln{2} -\frac{3}{2}\ln(3),\\
 H_{1/8} & = 8-\frac{\pi}{2} - 4\ln{2} - \frac{1}{\sqrt{2}}\{\pi + \ln(2 + \sqrt{2}) - \ln(2 - \sqrt{2})\},\\
 H_{1/12} & = 12-3\!\left(\ln{2}+\frac{\ln{3}}{2}\right)-\pi\!\left(1+\frac{\sqrt{3}}{2}\right) + 2\sqrt{3}\ln{\sqrt{2-\sqrt{3}}}
\end{align}

調和数の母函数[編集]

調和数の列の母函数

\sum_{n=1}^\infty z^n H_n = 
\frac {-\ln(1-z)}{1-z}

で与えられる(ln(z) は自然対数)。また、冪指数型母函数は

\sum_{n=1}^\infty \frac {z^n}{n!} H_n =
 -e^z \sum_{k=1}^\infty \frac{1}{k} \frac {(-z)^k}{k!} = 
e^z \mbox {Ein}(z)

となる。ここで Ein(z) は整指数積分で、

\text{Ein}(z) = \text{E}_1(z) + \gamma + \ln z = 
\Gamma (0,z) + \gamma + \ln z

が成り立つものである(ただし、Γ(0, z) は不完全ガンマ函数)。

応用[編集]

調和数は、二重ガンマ函数に対する

 \psi(n) = H_{n-1} - \gamma

のような、いくつかの特殊函数に関する計算公式に現れる。このような関係式は、しばしば調和数のパラメータ n を整数以外に拡張するための定義式としても利用される。先の節で述べたような極限によって、調和数から定数 γ を定義することがよく行われるが、

 \gamma = \lim_{n \to \infty}{\left(H_n - \ln\left(n+{1 \over 2}\right)\right)}

としたほうが収斂が早い。

2002年にジェフリー・ラガリアス英語版は、リーマン予想が「不等式

 \sigma(n) \le H_n + \ln(H_n)e^{H_n}

が任意の自然数 n に対して成立し、かつ n > 1 のときは真の(等号無しの)不等式として成立する」という主張に等価であることを示した。ここで σ(n) は n約数和である。

一般化[編集]

一般化調和数[編集]

n-番目の m-次一般化調和数 (generalized harmonic number)

H_{n,m}=\sum_{k=1}^n \frac{1}{k^m}

で与えられる。n を無限大に飛ばした極限が存在するのは m > 1 の時に限られることに注意。一般化調和数を表す記号としては

H_{n,m}= H_n^{(m)} = H_m(n).

なども使われることがある。なお、m = 1 の場合が通常の調和数であり、添字 m を落として

H_n= \sum_{k=1}^n \frac{1}{k}

と書く。また、n → ∞ の極限で一般化調和数はリーマン・ゼータ函数に収斂する。つまり

\lim_{n\to \infty} H_{n,m} = \zeta(m)

が成り立つ。一般化調和数はベルヌーイ数を調べる際に現れ、またスターリング数を調べる際にも現れる。一般化調和数の母函数は

\sum_{n=1}^\infty z^n H_{n,m} = 
\frac {\mbox{Li}_m(z)}{1-z}

である。ここで Lim(z) は多重対数函数で |z| < 1 とする。この式で m = 1 としたものは、先に述べた調和数列の母函数に一致する。

複素平面への一般化[編集]

調和数についてのオイラーの積分公式は次の積分等式

\int_a^1 \frac {1-x^s}{1-x} dx = 
- \sum_{k=1}^\infty \frac {1}{k} {s \choose k} (a-1)^k

から従うが、この式は s を一般の複素数としても(二項係数を適切に拡張すれば)成り立つ。a = 0 とすれば、この公式から調和数を補間して複素平面へ拡張した函数の積分表示と級数表示が両方得られる。この積分等式自体はニュートン級数(ニュートンの一般二項定理

\sum_{k=0}^\infty {s \choose k} (-x)^k = (1-x)^s

から簡単な操作で得られる。調和数を補間する函数は、実は二重ガンマ函数 ψ(x) をつかって

\psi(s+1)+\gamma = \int_0^1 \frac {1-x^s}{1-x} dx

と書ける(γ はオイラー-マスケローニ定数)。この積分の過程を繰り返せば

H_{s,2}=-\sum_{k=1}^\infty \frac {(-1)^k}{k} {s \choose k} H_k

を得る。

関連項目[編集]

参考文献[編集]

この記事は、クリエイティブ・コモンズ・ライセンス 表示-継承 3.0 非移植のもと提供されているオンライン数学辞典『PlanetMath』の項目Harmonic numberの本文を含む