オアの調和数については「調和数 」を、各項の逆数が等差数列であるような数列については「調和数列 」をご覧ください。
数学 において、n -番目の調和数 （ちょうわすう、英 : harmonic number n までの自然数 の逆数和
n = ⌊x ⌋ に対する調和数 H n ,1x )（青）に漸近収斂する。
H
n
=
1
+
1
2
+
1
3
+
⋯
+
1
n
=
∑
k
=
1
n
1
k
{\displaystyle 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
=
∫
0
1
1
−
x
n
1
−
x
d
x
{\displaystyle H_{n}=\int _{0}^{1}{\frac {1-x^{n}}{1-x}}\,dx}
はオイラー による。この等式は簡単な代数的等式
1
−
x
n
1
−
x
=
1
+
x
+
…
+
x
n
−
1
{\displaystyle {\frac {1-x^{n}}{1-x}}=1+x+\ldots +x^{n-1}}
を使えば明らかである。また、積分の変数を単純に x = 1 − u と変換すれば、H n
H
n
=
∫
0
1
1
−
(
1
−
u
)
n
u
d
u
=
∫
0
1
[
∑
k
=
1
n
(
−
1
)
k
−
1
(
n
k
)
u
k
−
1
]
d
u
=
∑
k
=
1
n
(
−
1
)
k
−
1
(
n
k
)
∫
0
1
u
k
−
1
d
u
=
∑
k
=
1
n
(
−
1
)
k
−
1
1
k
(
n
k
)
{\displaystyle {\begin{aligned}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{aligned}}}
が得られる。同じ表現は、第三レトケシュ恒等式 （英語版 ） x 1 = 1, ..., x n
Π
k
(
1
,
…
,
n
)
=
(
−
1
)
n
−
k
(
k
−
1
)
!
(
n
−
k
)
!
{\displaystyle \Pi _{k}(1,\ldots ,n)=(-1)^{n-k}(k-1)!(n-k)!}
なる事実を用いることでも得られる。すなわち
H
n
=
H
n
,
1
=
∑
k
=
1
n
1
k
=
(
−
1
)
n
−
1
n
!
∑
k
=
1
n
1
k
2
Π
k
(
1
,
…
,
n
)
=
∑
k
=
1
n
(
−
1
)
k
−
1
1
k
(
n
k
)
{\displaystyle 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 {n}{k}}}
が成り立つ。また、レトケシュ恒等式を x 1 = 12 , ..., x n n 2 に対して用いれば、この場合
Π
k
(
1
2
,
2
2
,
…
,
n
2
)
=
(
−
1
)
n
−
k
(
n
−
k
)
!
(
n
+
k
)
!
2
k
2
{\displaystyle \Pi _{k}(1^{2},2^{2},\ldots ,n^{2})=(-1)^{n-k}{\frac {(n-k)!(n+k)!}{2k^{2}}}}
となるので、ζ(2) の第 n -部分和についての類似の公式
H
n
,
2
=
∑
k
=
1
n
1
k
2
=
2
∑
k
=
1
n
(
−
1
)
k
−
1
1
k
2
(
n
k
)
(
n
+
k
k
)
{\displaystyle 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 {n}{k}}{\binom {n+k}{k}}}}
を得る。H n n の自然対数 ln(n ) と同程度の速さである。このことは、H n
∫
1
n
d
x
x
(
=
ln
(
n
)
)
{\displaystyle \int _{1}^{n}{dx \over x}\quad (=\ln(n))}
で近似することによって確認できる。数列 (H n n )) は単調に減少して、
lim
n
→
∞
H
n
−
ln
(
n
)
=
γ
{\displaystyle \lim _{n\to \infty }H_{n}-\ln(n)=\gamma }
なる定数（この定数 γ はオイラー・マスケローニ定数 と呼ばれ、その値は 0.5772156649... である）を極限 にもち、これに対応する漸近展開 は
H
n
=
ln
n
+
γ
+
1
2
n
−
1
−
1
12
n
−
2
+
1
120
n
−
4
+
O
(
n
−
6
)
{\displaystyle H_{n}=\ln {n}+\gamma +{\frac {1}{2}}n^{-1}-{\frac {1}{12}}n^{-2}+{\frac {1}{120}}n^{-4}+{\mathcal {O}}(n^{-6})}
で与えられる。
分数パラメータに対する特殊値
編集
調和数 H n n を積分
H
α
=
∫
0
1
1
−
x
α
1
−
x
d
x
{\displaystyle H_{\alpha }=\int _{0}^{1}{\frac {1-x^{\alpha }}{1-x}}dx}
によって拡張すれば、0 と 1 の間の分数値をもつパラメータ α に対する解析的な特殊値を定めることができる。あるいはさらに漸化式
H
α
=
H
α
−
1
+
1
α
{\displaystyle H_{\alpha }=H_{\alpha -1}+{\frac {1}{\alpha }}}
によって拡張することもでき、結局は任意の x > 0 に対して（x が整数のときもそうでないときも）
H
x
=
x
∑
k
=
1
∞
1
k
(
x
+
k
)
{\displaystyle H_{x}=x\sum _{k=1}^{\infty }{\frac {1}{k(x+k)}}}
が成立する。いくつかの特殊値について計算すれば、以下のようになる。
H
3
/
4
=
4
3
−
3
ln
2
+
π
2
,
H
2
/
3
=
3
2
(
1
−
ln
3
)
+
π
6
3
,
H
1
/
2
=
2
−
2
ln
2
,
H
1
/
3
=
3
−
π
2
3
−
3
2
ln
3
,
H
1
/
4
=
4
−
π
2
−
3
ln
2
,
H
1
/
6
=
6
−
π
2
3
−
2
ln
2
−
3
2
ln
(
3
)
,
H
1
/
8
=
8
−
π
2
−
4
ln
2
−
1
2
{
π
+
ln
(
2
+
2
)
−
ln
(
2
−
2
)
}
,
H
1
/
12
=
12
−
3
(
ln
2
+
ln
3
2
)
−
π
(
1
+
3
2
)
+
2
3
ln
2
−
3
{\displaystyle {\begin{aligned}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{aligned}}}
調和数の母函数
編集
調和数の列の母函数 は
∑
n
=
1
∞
z
n
H
n
=
−
ln
(
1
−
z
)
1
−
z
{\displaystyle \sum _{n=1}^{\infty }z^{n}H_{n}={\frac {-\ln(1-z)}{1-z}}}
で与えられる（ln(z ) は自然対数 ）。また、冪指数型母函数は
∑
n
=
1
∞
z
n
n
!
H
n
=
−
e
z
∑
k
=
1
∞
1
k
(
−
z
)
k
k
!
=
e
z
Ein
(
z
)
{\displaystyle \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 ) は整指数積分 で、
Ein
(
z
)
=
E
1
(
z
)
+
γ
+
ln
z
=
Γ
(
0
,
z
)
+
γ
+
ln
z
{\displaystyle {\text{Ein}}(z)={\text{E}}_{1}(z)+\gamma +\ln z=\Gamma (0,z)+\gamma +\ln z}
が成り立つものである（ただし、Γ(0, z ) は不完全ガンマ函数 ）。
調和数は、ディガンマ関数 に対する
ψ
(
n
)
=
H
n
−
1
−
γ
{\displaystyle \psi (n)=H_{n-1}-\gamma }
のような、いくつかの特殊函数に関する計算公式に現れる。このような関係式は、しばしば調和数のパラメータ n を整数以外に拡張するための定義式としても利用される。先の節で述べたような極限によって、調和数から定数 γ を定義することがよく行われるが、
γ
=
lim
n
→
∞
(
H
n
−
ln
(
n
+
1
2
)
)
{\displaystyle \gamma =\lim _{n\to \infty }{\left(H_{n}-\ln \left(n+{1 \over 2}\right)\right)}}
としたほうが収斂が早い。
2002年にジェフリー・ラガリアス は、リーマン予想 が「不等式
σ
(
n
)
≤
H
n
+
ln
(
H
n
)
e
H
n
{\displaystyle \sigma (n)\leq H_{n}+\ln(H_{n})e^{H_{n}}}
が任意の自然数 n に対して成立し、かつ n > 1 のときは真の（等号無しの）不等式として成立する」という主張に等価であることを示した。ここで σ(n ) は n の約数和 である。
一般化調和数
編集
n -番目の m -次一般化調和数 (generalized harmonic number ) は
H
n
,
m
=
∑
k
=
1
n
1
k
m
{\displaystyle H_{n,m}=\sum _{k=1}^{n}{\frac {1}{k^{m}}}}
で与えられる。n を無限大に飛ばした極限が存在するのは m > 1 の時に限られることに注意。一般化調和数を表す記号としては
H
n
,
m
=
H
n
(
m
)
=
H
m
(
n
)
.
{\displaystyle H_{n,m}=H_{n}^{(m)}=H_{m}(n).}
なども使われることがある。なお、m = 1 の場合が通常の調和数であり、添字 m を落として
H
n
=
∑
k
=
1
n
1
k
{\displaystyle H_{n}=\sum _{k=1}^{n}{\frac {1}{k}}}
と書く。また、n → ∞ の極限で一般化調和数はリーマンゼータ函数 に収斂する。つまり
lim
n
→
∞
H
n
,
m
=
ζ
(
m
)
{\displaystyle \lim _{n\to \infty }H_{n,m}=\zeta (m)}
が成り立つ。一般化調和数はベルヌーイ数 を調べる際に現れ、またスターリング数 を調べる際にも現れる。一般化調和数の母函数は
∑
n
=
1
∞
z
n
H
n
,
m
=
Li
m
(
z
)
1
−
z
{\displaystyle \sum _{n=1}^{\infty }z^{n}H_{n,m}={\frac {{\mbox{Li}}_{m}(z)}{1-z}}}
である。ここで Lim z ) は多重対数函数 で |z | < 1 とする。この式で m = 1 としたものは、先に述べた調和数列の母函数に一致する。
複素平面への一般化
編集
調和数についてのオイラーの積分公式は次の積分等式
∫
a
1
1
−
x
s
1
−
x
d
x
=
−
∑
k
=
1
∞
1
k
(
s
k
)
(
a
−
1
)
k
{\displaystyle \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 とすれば、この公式から調和数を補間して複素平面へ拡張した函数の積分表示と級数表示が両方得られる。この積分等式自体はニュートン級数 （ニュートンの一般二項定理 ）
∑
k
=
0
∞
(
s
k
)
(
−
x
)
k
=
(
1
−
x
)
s
{\displaystyle \sum _{k=0}^{\infty }{s \choose k}(-x)^{k}=(1-x)^{s}}
から簡単な操作で得られる。調和数を補間する函数は、実はディガンマ関数 ψ(x ) をつかって
ψ
(
s
+
1
)
+
γ
=
∫
0
1
1
−
x
s
1
−
x
d
x
{\displaystyle \psi (s+1)+\gamma =\int _{0}^{1}{\frac {1-x^{s}}{1-x}}dx}
と書ける（γ はオイラー-マスケローニ定数）。この積分の過程を繰り返せば
H
s
,
2
=
−
∑
k
=
1
∞
(
−
1
)
k
k
(
s
k
)
H
k
{\displaystyle H_{s,2}=-\sum _{k=1}^{\infty }{\frac {(-1)^{k}}{k}}{s \choose k}H_{k}}
を得る。
関連項目
編集
参考文献
編集
Arthur T. Benjamin, Gregory O. Preston, Jennifer J. Quinn, A Stirling Encounter with Harmonic Numbers 75 (2) pp 95-103.
Donald Knuth . The Art of Computer Programming , Volume 1: Fundamental Algorithms , Third Edition. Addison-Wesley, 1997. ISBN 0-201-89683-4 . Section 1.2.7: Harmonic Numbers, pp.75–79.Ed Sandifer, How Euler Did It -- Estimating the Basel problem
Weisstein, Eric W. "Harmonic Number" . mathworld.wolfram.com (英語). Peter Paule and Carsten Schneider, Computer Proofs of a New Family of Harmonic Number Identities
Wenchang CHU, A Binomial Coefficient Identity Associated with Beukers' Conjecture on Apery Numbers The Electronic Journal of Combinatorics , 11 , #N15.
Ayhan Dil and Istvan Mezo, A Symmetric Algorithm for Hyperharmonic and Fibonacci Numbers 206 , 942--951.
Zoltán Retkes, "An extension of the Hermite–Hadamard Inequality ", Acta Sci. Math. (Szeged) この記事は、クリエイティブ・コモンズ・ライセンス 表示-継承 3.0 非移植 のもと提供されているオンライン数学辞典『PlanetMath 』の項目Harmonic number の本文を含む
![{\displaystyle {\begin{aligned}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{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8ce56fa34d90b72e4da5efb26629fc1eccbf23b4)
