Sum of the first n whole number reciprocals; 1/1 + 1/2 + 1/3 + ... + 1/n
Номер гармоники
H
n
{\displaystyle H_{n}}
с
n
=
⌊
x
⌋
{\displaystyle n=\lfloor x\rfloor }
(красная линия) со своим асимптотическим пределом
γ
+
ln
(
x
)
{\displaystyle \gamma +\ln(x)}
(синяя линия), где
γ
{\displaystyle \gamma }
– постоянная Эйлера–Машерони .
В математике гармоническое n -е число представляет собой сумму обратных величин первым n натуральных чисел :
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}}.}
Начиная с n = 1 , начинается последовательность номеров гармоник:
1
,
3
2
,
11
6
,
25
12
,
137
60
,
…
{\displaystyle 1,{\frac {3}{2}},{\frac {11}{6}},{\frac {25}{12}},{\frac {137}{60}},\dots }
Числа гармоник связаны со средним гармоническим числом тем, что n -й номер гармоники также в n раз превышает среднее гармоническое первых n положительных целых чисел.
Гармонические числа изучаются с древности и играют важную роль в различных разделах теории чисел . Их иногда называют гармоническими рядами , они тесно связаны с дзета-функцией Римана и появляются в выражениях различных специальных функций .
Числа гармоник примерно аппроксимируют функцию натурального логарифма. [1] : 143 и, таким образом, связанный с ним гармонический ряд растет безгранично, хотя и медленно. В 1737 году Леонард Эйлер использовал расхождение гармонического ряда, чтобы предоставить новое доказательство бесконечности простых чисел . Его работа была расширена на комплексную плоскость Бернхардом Риманом в 1859 году, что привело непосредственно к знаменитой гипотезе Римана о распределении простых чисел .
Когда стоимость большого количества предметов имеет распределение по закону Ципфа , общая стоимость n наиболее ценных предметов пропорциональна n -му номеру гармоники. Это приводит к множеству удивительных выводов относительно « длинного хвоста» и теории сетевой стоимости .
Теорема Бертрана -Чебышева подразумевает, что, за исключением случая n = 1 , числа гармоник никогда не являются целыми числами. [2]
Первые 40 гармонических чисел
show
n
Harmonic number, Hn
expressed as a fraction
decimal
relative size
1
1
1
1
2
3
/2
1.5
1.5
3
11
/6
~1.83333
1.83333
4
25
/12
~2.08333
2.08333
5
137
/60
~2.28333
2.28333
6
49
/20
2.45
2.45
7
363
/140
~2.59286
2.59286
8
761
/280
~2.71786
2.71786
9
7 129
/2 520
~2.82897
2.82897
10
7 381
/2 520
~2.92897
2.92897
11
83 711
/27 720
~3.01988
3.01988
12
86 021
/27 720
~3.10321
3.10321
13
1 145 993
/360 360
~3.18013
3.18013
14
1 171 733
/360 360
~3.25156
3.25156
15
1 195 757
/360 360
~3.31823
3.31823
16
2 436 559
/720 720
~3.38073
3.38073
17
42 142 223
/12 252 240
~3.43955
3.43955
18
14 274 301
/4 084 080
~3.49511
3.49511
19
275 295 799
/77 597 520
~3.54774
3.54774
20
55 835 135
/15 519 504
~3.59774
3.59774
21
18 858 053
/5 173 168
~3.64536
3.64536
22
19 093 197
/5 173 168
~3.69081
3.69081
23
444 316 699
/118 982 864
~3.73429
3.73429
24
1 347 822 955
/356 948 592
~3.77596
3.77596
25
34 052 522 467
/8 923 714 800
~3.81596
3.81596
26
34 395 742 267
/8 923 714 800
~3.85442
3.85442
27
312 536 252 003
/80 313 433 200
~3.89146
3.89146
28
315 404 588 903
/80 313 433 200
~3.92717
3.92717
29
9 227 046 511 387
/2 329 089 562 800
~3.96165
3.96165
30
9 304 682 830 147
/2 329 089 562 800
~3.99499
3.99499
31
290 774 257 297 357
/72 201 776 446 800
~4.02725
4.02725
32
586 061 125 622 639
/144 403 552 893 600
~4.05850
4.0585
33
53 676 090 078 349
/13 127 595 717 600
~4.08880
4.0888
34
54 062 195 834 749
/13 127 595 717 600
~4.11821
4.11821
35
54 437 269 998 109
/13 127 595 717 600
~4.14678
4.14678
36
54 801 925 434 709
/13 127 595 717 600
~4.17456
4.17456
37
2 040 798 836 801 833
/485 721 041 551 200
~4.20159
4.20159
38
2 053 580 969 474 233
/485 721 041 551 200
~4.22790
4.2279
39
2 066 035 355 155 033
/485 721 041 551 200
~4.25354
4.25354
40
2 078 178 381 193 813
/485 721 041 551 200
~4.27854
4.27854
гармоник , Тождества включающие номера
По определению числа гармоник удовлетворяют рекуррентному соотношению
H
n
+
1
=
H
n
+
1
n
+
1
.
{\displaystyle H_{n+1}=H_{n}+{\frac {1}{n+1}}.}
Гармонические числа связаны с числами Стирлинга первого рода соотношением
H
n
=
1
n
!
[
n
+
1
2
]
.
{\displaystyle H_{n}={\frac {1}{n!}}\left[{n+1 \atop 2}\right].}
Номера гармоник удовлетворяют тождествам рядов
∑
k
=
1
n
H
k
=
(
n
+
1
)
H
n
−
n
{\displaystyle \sum _{k=1}^{n}H_{k}=(n+1)H_{n}-n}
и
∑
k
=
1
n
H
k
2
=
(
n
+
1
)
H
n
2
−
(
2
n
+
1
)
H
n
+
2
n
.
{\displaystyle \sum _{k=1}^{n}H_{k}^{2}=(n+1)H_{n}^{2}-(2n+1)H_{n}+2n.}
Эти два результата очень аналогичны соответствующим интегральным результатам.
∫
0
x
log
y
d
y
=
x
log
x
−
x
{\displaystyle \int _{0}^{x}\log y\ dy=x\log x-x}
и
∫
0
x
(
log
y
)
2
d
y
=
x
(
log
x
)
2
−
2
x
log
x
+
2
x
.
{\displaystyle \int _{0}^{x}(\log y)^{2}\ dy=x(\log x)^{2}-2x\log x+2x.}
Есть несколько бесконечных суммаций, включающих числа гармоник и степени π : [3] [ нужен лучший источник ]
∑
n
=
1
∞
H
n
n
⋅
2
n
=
π
2
12
∑
n
=
1
∞
H
n
2
(
n
+
1
)
2
=
11
360
π
4
∑
n
=
1
∞
H
n
2
(
n
+
1
)
2
=
11
360
π
4
∑
n
=
1
∞
H
n
n
3
=
π
4
72
{\displaystyle {\begin{aligned}\sum _{n=1}^{\infty }{\frac {H_{n}}{n\cdot 2^{n}}}&={\frac {\pi ^{2}}{12}}\\\sum _{n=1}^{\infty }{\frac {H_{n}^{2}}{(n+1)^{2}}}&={\frac {11}{360}}\pi ^{4}\\\sum _{n=1}^{\infty }{\frac {H_{n}^{2}}{(n+1)^{2}}}&={\frac {11}{360}}\pi ^{4}\\\sum _{n=1}^{\infty }{\frac {H_{n}}{n^{3}}}&={\frac {\pi ^{4}}{72}}\end{aligned}}}
Интегральное представление Эйлера [4] является
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+\cdots +x^{n-1}.}
Используя замену x = 1 − u , другое выражение для H n :
H
n
=
∫
0
1
1
−
x
n
1
−
x
d
x
=
∫
0
1
1
−
(
1
−
u
)
n
u
d
u
=
∫
0
1
[
∑
k
=
1
n
(
n
k
)
(
−
u
)
k
−
1
]
d
u
=
∑
k
=
1
n
(
n
k
)
∫
0
1
(
−
u
)
k
−
1
d
u
=
∑
k
=
1
n
(
n
k
)
(
−
1
)
k
−
1
k
.
{\displaystyle {\begin{aligned}H_{n}&=\int _{0}^{1}{\frac {1-x^{n}}{1-x}}\,dx=\int _{0}^{1}{\frac {1-(1-u)^{n}}{u}}\,du\\[6pt]&=\int _{0}^{1}\left[\sum _{k=1}^{n}{\binom {n}{k}}(-u)^{k-1}\right]\,du=\sum _{k=1}^{n}{\binom {n}{k}}\int _{0}^{1}(-u)^{k-1}\,du\\[6pt]&=\sum _{k=1}^{n}{\binom {n}{k}}{\frac {(-1)^{k-1}}{k}}.\end{aligned}}}
График, демонстрирующий связь между гармоническими числами и натуральным логарифмом . Номер гармоники H n можно интерпретировать как сумму Римана интеграла:
∫
1
n
+
1
d
x
x
=
ln
(
n
+
1
)
.
{\displaystyle \int _{1}^{n+1}{\frac {dx}{x}}=\ln(n+1).}
Номер n- гармоники примерно равен натуральному логарифму n й . Причина в том, что сумма аппроксимируется интегралом
∫
1
n
1
x
d
x
,
{\displaystyle \int _{1}^{n}{\frac {1}{x}}\,dx,}
значение которого равно
ln n .
Значения последовательности H n − ln n монотонно убывают к пределу
lim
n
→
∞
(
H
n
−
ln
n
)
=
γ
,
{\displaystyle \lim _{n\to \infty }\left(H_{n}-\ln n\right)=\gamma ,}
где
γ ≈ 0,5772156649 —
постоянная Эйлера–Машерони . Соответствующее
асимптотическое разложение имеет вид
H
n
∼
ln
n
+
γ
+
1
2
n
−
∑
k
=
1
∞
B
2
k
2
k
n
2
k
=
ln
n
+
γ
+
1
2
n
−
1
12
n
2
+
1
120
n
4
−
⋯
,
{\displaystyle {\begin{aligned}H_{n}&\sim \ln {n}+\gamma +{\frac {1}{2n}}-\sum _{k=1}^{\infty }{\frac {B_{2k}}{2kn^{2k}}}\\&=\ln {n}+\gamma +{\frac {1}{2n}}-{\frac {1}{12n^{2}}}+{\frac {1}{120n^{4}}}-\cdots ,\end{aligned}}}
где
Bk — числа
Бернулли .
Генерирующие функции [ править ]
Производящая функция для чисел гармоник:
∑
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
−
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-1}}{k}}{\frac {z^{k}}{k!}}=e^{z}\operatorname {Ein} (z)}
где Ein(
z ) — полный
экспоненциальный интеграл . Экспоненциальный интеграл также можно выразить как
Ein
(
z
)
=
E
1
(
z
)
+
γ
+
ln
z
=
Γ
(
0
,
z
)
+
γ
+
ln
z
{\displaystyle \operatorname {Ein} (z)=\mathrm {E} _{1}(z)+\gamma +\ln z=\Gamma (0,z)+\gamma +\ln z}
где Γ(0,
z ) —
неполная гамма-функция .
Арифметические свойства [ править ]
Гармонические числа обладают несколькими интересными арифметическими свойствами. Хорошо известно, что
H
n
{\textstyle H_{n}}
является целым числом тогда и только тогда, когда
n
=
1
{\textstyle n=1}
, результат, который часто приписывают Тэзингеру. [5] Действительно, используя 2-адическое нормирование , нетрудно доказать, что для
n
≥
2
{\textstyle n\geq 2}
числитель
H
n
{\textstyle H_{n}}
является нечетным числом, а знаменатель
H
n
{\textstyle H_{n}}
это четное число. Точнее,
H
n
=
1
2
⌊
log
2
(
n
)
⌋
a
n
b
n
{\displaystyle H_{n}={\frac {1}{2^{\lfloor \log _{2}(n)\rfloor }}}{\frac {a_{n}}{b_{n}}}}
с некоторыми нечетными целыми числами
a
n
{\textstyle a_{n}}
и
b
n
{\textstyle b_{n}}
.
Как следствие теоремы Вольстенхолма , для любого простого числа
p
≥
5
{\displaystyle p\geq 5}
числитель
H
p
−
1
{\displaystyle H_{p-1}}
делится на
p
2
{\textstyle p^{2}}
. Более того, Эйзенштейн [6] доказал, что для всех нечетных простых чисел
p
{\textstyle p}
оно держится
H
(
p
−
1
)
/
2
≡
−
2
q
p
(
2
)
(
mod
p
)
{\displaystyle H_{(p-1)/2}\equiv -2q_{p}(2){\pmod {p}}}
где
q
p
(
2
)
=
(
2
p
−
1
−
1
)
/
p
{\textstyle q_{p}(2)=(2^{p-1}-1)/p}
является
фактором Ферма , и следствием этого является то, что
p
{\textstyle p}
делит числитель
H
(
p
−
1
)
/
2
{\displaystyle H_{(p-1)/2}}
если и только если
p
{\textstyle p}
является
простым числом Вифериха .
В 1991 году Ишваратасан и Левин [7] определенный
J
p
{\displaystyle J_{p}}
как набор всех положительных целых чисел
n
{\displaystyle n}
такой, что числитель
H
n
{\displaystyle H_{n}}
делится на простое число
p
.
{\displaystyle p.}
Они доказали, что
{
p
−
1
,
p
2
−
p
,
p
2
−
1
}
⊆
J
p
{\displaystyle \{p-1,p^{2}-p,p^{2}-1\}\subseteq J_{p}}
для всех простых чисел
p
≥
5
,
{\displaystyle p\geq 5,}
и они определили
гармонические простые числа как простые числа
p
{\textstyle p}
такой, что
J
p
{\displaystyle J_{p}}
имеет ровно 3 элемента.
Ишваратасан и Левин также предположили, что
J
p
{\displaystyle J_{p}}
является конечным множеством для всех простых чисел
p
,
{\displaystyle p,}
и что существует бесконечно много гармонических простых чисел. Бойд [8] проверил, что
J
p
{\displaystyle J_{p}}
конечно для всех простых чисел до
p
=
547
{\displaystyle p=547}
кроме 83, 127 и 397; и он предложил эвристику, предполагающую, что плотность гармонических простых чисел в множестве всех простых чисел должна быть равна
1
/
e
{\displaystyle 1/e}
. Истинный [9] показало, что
J
p
{\displaystyle J_{p}}
имеет нулевую асимптотическую плотность , а Бин-Лин Ву и Юн-Гао Чен [10] доказал, что количество элементов
J
p
{\displaystyle J_{p}}
не превышающий
x
{\displaystyle x}
самое большее
3
x
2
3
+
1
25
log
p
{\displaystyle 3x^{{\frac {2}{3}}+{\frac {1}{25\log p}}}}
, для всех
x
≥
1
{\displaystyle x\geq 1}
.
Приложения [ править ]
Номера гармоник появляются в нескольких формулах расчета, таких как дигамма-функция.
ψ
(
n
)
=
H
n
−
1
−
γ
.
{\displaystyle \psi (n)=H_{n-1}-\gamma .}
Это соотношение также часто используется для определения расширения чисел гармоник до нецелых чисел
n . Числа гармоник также часто используются для определения
γ с использованием введенного ранее предела:
γ
=
lim
n
→
∞
(
H
n
−
ln
(
n
)
)
,
{\displaystyle \gamma =\lim _{n\rightarrow \infty }{\left(H_{n}-\ln(n)\right)},}
хотя
γ
=
lim
n
→
∞
(
H
n
−
ln
(
n
+
1
2
)
)
{\displaystyle \gamma =\lim _{n\to \infty }{\left(H_{n}-\ln \left(n+{\frac {1}{2}}\right)\right)}}
сходится быстрее.
В 2002 году Джеффри Лагариас доказал [11] что гипотеза Римана эквивалентна утверждению, что
σ
(
n
)
≤
H
n
+
(
log
H
n
)
e
H
n
,
{\displaystyle \sigma (n)\leq H_{n}+(\log H_{n})e^{H_{n}},}
верно для любого
целого числа n ≥ 1 со строгим неравенством, если
n > 1 ; здесь
σ ( n ) обозначает
делителей n
. сумму
Собственные значения нелокальной задачи
λ
φ
(
x
)
=
∫
−
1
1
φ
(
x
)
−
φ
(
y
)
|
x
−
y
|
d
y
{\displaystyle \lambda \varphi (x)=\int _{-1}^{1}{\frac {\varphi (x)-\varphi (y)}{|x-y|}}\,dy}
даны
λ
=
2
H
n
{\displaystyle \lambda =2H_{n}}
, где по соглашению
H
0
=
0
{\displaystyle H_{0}=0}
, а соответствующие собственные функции задаются
полиномами Лежандра
φ
(
x
)
=
P
n
(
x
)
{\displaystyle \varphi (x)=P_{n}(x)}
.
[12]
числа Обобщенные гармоник
Номер n- й обобщенной гармоники порядка m определяется выражением
H
n
,
m
=
∑
k
=
1
n
1
k
m
.
{\displaystyle H_{n,m}=\sum _{k=1}^{n}{\frac {1}{k^{m}}}.}
(В некоторых источниках это также может обозначаться как
H
n
(
m
)
{\textstyle H_{n}^{(m)}}
или
H
m
(
n
)
.
{\textstyle H_{m}(n).}
)
Частный случай m = 0 дает
H
n
,
0
=
n
.
{\displaystyle H_{n,0}=n.}
Частный случай m = 1 сводится к обычному номеру гармоники:
H
n
,
1
=
H
n
=
∑
k
=
1
n
1
k
.
{\displaystyle H_{n,1}=H_{n}=\sum _{k=1}^{n}{\frac {1}{k}}.}
Предел
H
n
,
m
{\textstyle H_{n,m}}
при n → ∞ конечно, если m > 1 , с обобщенным гармоническим числом, ограниченным и сходящимся к дзета-функции Римана
lim
n
→
∞
H
n
,
m
=
ζ
(
m
)
.
{\displaystyle \lim _{n\rightarrow \infty }H_{n,m}=\zeta (m).}
Наименьшее натуральное число k такое, что k н не делит знаменатель номера обобщенной гармоники H ( k , n ), а знаменатель чередующегося номера обобщенной гармоники H ' ( k , n ) не равен для n = 1, 2, ... :
77, 20, 94556602, 42, 444, 20, 104, 42, 76, 20, 77, 110, 3504, 20, 903, 42, 1107, 20, 104, 42, 77, 20, 2948, 110, 136, 20, 76, 42, 903, 20, 77, 42, 268, 20, 7004, 110, 1752, 20, 19203, 42, 77, 20, 104, 42, 76, 20, 370, 110, 1107, 20, ... (последовательность A128670 в OEIS )
Соответствующая сумма
∑
k
=
1
n
k
m
{\displaystyle \sum _{k=1}^{n}k^{m}}
встречается при изучении чисел Бернулли ; гармонические числа появляются также при изучении чисел Стирлинга .
Некоторые интегралы от чисел обобщенных гармоник:
∫
0
a
H
x
,
2
d
x
=
a
π
2
6
−
H
a
{\displaystyle \int _{0}^{a}H_{x,2}\,dx=a{\frac {\pi ^{2}}{6}}-H_{a}}
и
∫
0
a
H
x
,
3
d
x
=
a
A
−
1
2
H
a
,
2
,
{\displaystyle \int _{0}^{a}H_{x,3}\,dx=aA-{\frac {1}{2}}H_{a,2},}
где
A —
постоянная Апери ζ (3),
и
∑
k
=
1
n
H
k
,
m
=
(
n
+
1
)
H
n
,
m
−
H
n
,
m
−
1
for
m
≥
0.
{\displaystyle \sum _{k=1}^{n}H_{k,m}=(n+1)H_{n,m}-H_{n,m-1}{\text{ for }}m\geq 0.}
Каждое обобщенное гармоническое число порядка m можно записать как функцию чисел гармоник порядка
m
−
1
{\displaystyle m-1}
с использованием
H
n
,
m
=
∑
k
=
1
n
−
1
H
k
,
m
−
1
k
(
k
+
1
)
+
H
n
,
m
−
1
n
{\displaystyle H_{n,m}=\sum _{k=1}^{n-1}{\frac {H_{k,m-1}}{k(k+1)}}+{\frac {H_{n,m-1}}{n}}}
например:
H
4
,
3
=
H
1
,
2
1
⋅
2
+
H
2
,
2
2
⋅
3
+
H
3
,
2
3
⋅
4
+
H
4
,
2
4
{\displaystyle H_{4,3}={\frac {H_{1,2}}{1\cdot 2}}+{\frac {H_{2,2}}{2\cdot 3}}+{\frac {H_{3,2}}{3\cdot 4}}+{\frac {H_{4,2}}{4}}}
для Производящая функция чисел обобщенных гармоник:
∑
n
=
1
∞
z
n
H
n
,
m
=
Li
m
(
z
)
1
−
z
,
{\displaystyle \sum _{n=1}^{\infty }z^{n}H_{n,m}={\frac {\operatorname {Li} _{m}(z)}{1-z}},}
где
Li
m
(
z
)
{\displaystyle \operatorname {Li} _{m}(z)}
—
полилогарифм , а
| г | < 1 . Приведенная выше производящая функция для
m = 1 является частным случаем этой формулы.
Дробный аргумент для обобщенных чисел гармоник можно ввести следующим образом:
Для каждого
p
,
q
>
0
{\displaystyle p,q>0}
целое число и
m
>
1
{\displaystyle m>1}
целое или нет, мы имеем из полигамма-функций:
H
q
/
p
,
m
=
ζ
(
m
)
−
p
m
∑
k
=
1
∞
1
(
q
+
p
k
)
m
{\displaystyle H_{q/p,m}=\zeta (m)-p^{m}\sum _{k=1}^{\infty }{\frac {1}{(q+pk)^{m}}}}
где
ζ
(
m
)
{\displaystyle \zeta (m)}
—
дзета-функция Римана . Соответствующее рекуррентное соотношение
H
a
,
m
=
H
a
−
1
,
m
+
1
a
m
.
{\displaystyle H_{a,m}=H_{a-1,m}+{\frac {1}{a^{m}}}.}
Некоторые особые значения
H
1
4
,
2
=
16
−
5
6
π
2
−
8
G
H
1
2
,
2
=
4
−
π
2
3
H
3
4
,
2
=
16
9
−
5
6
π
2
+
8
G
H
1
4
,
3
=
64
−
π
3
−
27
ζ
(
3
)
H
1
2
,
3
=
8
−
6
ζ
(
3
)
H
3
4
,
3
=
(
4
3
)
3
+
π
3
−
27
ζ
(
3
)
{\displaystyle {\begin{aligned}H_{{\frac {1}{4}},2}&=16-{\tfrac {5}{6}}\pi ^{2}-8G\\H_{{\frac {1}{2}},2}&=4-{\frac {\pi ^{2}}{3}}\\H_{{\frac {3}{4}},2}&={\frac {16}{9}}-{\frac {5}{6}}\pi ^{2}+8G\\H_{{\frac {1}{4}},3}&=64-\pi ^{3}-27\zeta (3)\\H_{{\frac {1}{2}},3}&=8-6\zeta (3)\\H_{{\frac {3}{4}},3}&=\left({\frac {4}{3}}\right)^{3}+\pi ^{3}-27\zeta (3)\end{aligned}}}
где
G —
константа Каталана . В частном случае, когда
p
=
1
{\displaystyle p=1}
, мы получаем
H
n
,
m
=
ζ
(
m
,
1
)
−
ζ
(
m
,
n
+
1
)
,
{\displaystyle H_{n,m}=\zeta (m,1)-\zeta (m,n+1),}
где
ζ
(
m
,
n
)
{\displaystyle \zeta (m,n)}
– дзета-функция Гурвица . Это соотношение используется для численного расчета номеров гармоник.
Формулы умножения [ править ]
Теорема умножения применима к гармоническим числам. Используя полигамма- функции, получаем
H
2
x
=
1
2
(
H
x
+
H
x
−
1
2
)
+
ln
2
H
3
x
=
1
3
(
H
x
+
H
x
−
1
3
+
H
x
−
2
3
)
+
ln
3
,
{\displaystyle {\begin{aligned}H_{2x}&={\frac {1}{2}}\left(H_{x}+H_{x-{\frac {1}{2}}}\right)+\ln 2\\H_{3x}&={\frac {1}{3}}\left(H_{x}+H_{x-{\frac {1}{3}}}+H_{x-{\frac {2}{3}}}\right)+\ln 3,\end{aligned}}}
или, в более общем смысле,
H
n
x
=
1
n
(
H
x
+
H
x
−
1
n
+
H
x
−
2
n
+
⋯
+
H
x
−
n
−
1
n
)
+
ln
n
.
{\displaystyle H_{nx}={\frac {1}{n}}\left(H_{x}+H_{x-{\frac {1}{n}}}+H_{x-{\frac {2}{n}}}+\cdots +H_{x-{\frac {n-1}{n}}}\right)+\ln n.}
Для обобщенных чисел гармоник имеем
H
2
x
,
2
=
1
2
(
ζ
(
2
)
+
1
2
(
H
x
,
2
+
H
x
−
1
2
,
2
)
)
H
3
x
,
2
=
1
9
(
6
ζ
(
2
)
+
H
x
,
2
+
H
x
−
1
3
,
2
+
H
x
−
2
3
,
2
)
,
{\displaystyle {\begin{aligned}H_{2x,2}&={\frac {1}{2}}\left(\zeta (2)+{\frac {1}{2}}\left(H_{x,2}+H_{x-{\frac {1}{2}},2}\right)\right)\\H_{3x,2}&={\frac {1}{9}}\left(6\zeta (2)+H_{x,2}+H_{x-{\frac {1}{3}},2}+H_{x-{\frac {2}{3}},2}\right),\end{aligned}}}
где
ζ
(
n
)
{\displaystyle \zeta (n)}
—
дзета-функция Римана .
Гипергармонические числа [ править ]
Следующее обобщение обсуждалось Дж. Х. Конвеем и Р. К. Гаем в их книге «Книга чисел» 1995 года . [1] : 258 Позволять
H
n
(
0
)
=
1
n
.
{\displaystyle H_{n}^{(0)}={\frac {1}{n}}.}
Тогда n-е
гипергармоническое число порядка
r (
r>0 ) определяется рекурсивно как
H
n
(
r
)
=
∑
k
=
1
n
H
k
(
r
−
1
)
.
{\displaystyle H_{n}^{(r)}=\sum _{k=1}^{n}H_{k}^{(r-1)}.}
В частности,
H
n
(
1
)
{\displaystyle H_{n}^{(1)}}
это обычный номер гармоники
H
n
{\displaystyle H_{n}}
.
гармонические числа Римские
Римские гармонические числа , [13] названные в честь Стивена Романа , были представлены Даниэлем Лебом и Джан-Карло Ротой в контексте обобщения умбрального исчисления с логарифмами. [14] Существует множество возможных определений, но одно из них —
n
,
k
≥
0
{\displaystyle n,k\geq 0}
, является
c
n
(
0
)
=
1
,
{\displaystyle c_{n}^{(0)}=1,}
и
c
n
(
k
+
1
)
=
∑
i
=
1
n
c
i
(
k
)
i
.
{\displaystyle c_{n}^{(k+1)}=\sum _{i=1}^{n}{\frac {c_{i}^{(k)}}{i}}.}
Конечно,
c
n
(
1
)
=
H
n
.
{\displaystyle c_{n}^{(1)}=H_{n}.}
Если
n
≠
0
{\displaystyle n\neq 0}
, они удовлетворяют
c
n
(
k
+
1
)
−
c
n
(
k
)
n
=
c
n
−
1
(
k
+
1
)
.
{\displaystyle c_{n}^{(k+1)}-{\frac {c_{n}^{(k)}}{n}}=c_{n-1}^{(k+1)}.}
Формулы закрытой формы
c
n
(
k
)
=
n
!
(
−
1
)
k
s
(
−
n
,
k
)
,
{\displaystyle c_{n}^{(k)}=n!(-1)^{k}s(-n,k),}
где
s
(
−
n
,
k
)
{\displaystyle s(-n,k)}
-
числа Стирлинга первого рода, обобщенные на отрицательный первый аргумент, и
c
n
(
k
)
=
∑
j
=
1
n
(
n
j
)
(
−
1
)
j
−
1
j
k
,
{\displaystyle c_{n}^{(k)}=\sum _{j=1}^{n}{\binom {n}{j}}{\frac {(-1)^{j-1}}{j^{k}}},}
который нашел
Дональд Кнут .
Фактически, эти числа были определены в более общем виде с использованием римских чисел и римских факториалов , которые включают отрицательные значения для
n
{\displaystyle n}
. Это обобщение было полезно в их исследовании по определению гармонических логарифмов .
Гармонические числа для действительных и комплексных значений [ править ]
Формулы, приведенные выше,
H
x
=
∫
0
1
1
−
t
x
1
−
t
d
t
=
∑
k
=
1
∞
(
x
k
)
(
−
1
)
k
−
1
k
{\displaystyle H_{x}=\int _{0}^{1}{\frac {1-t^{x}}{1-t}}\,dt=\sum _{k=1}^{\infty }{x \choose k}{\frac {(-1)^{k-1}}{k}}}
являются интегралом и последовательным представлением функции, которая интерполирует гармонические числа и посредством
аналитического продолжения расширяет определение на комплексную плоскость, отличную от отрицательных целых чисел
x . Интерполирующая функция на самом деле тесно связана с
дигамм-функцией.
H
x
=
ψ
(
x
+
1
)
+
γ
,
{\displaystyle H_{x}=\psi (x+1)+\gamma ,}
где
ψ ( x ) — дигамма-функция, а
γ —
постоянная Эйлера–Машерони . Процесс интегрирования можно повторить, чтобы получить
H
x
,
2
=
∑
k
=
1
∞
(
−
1
)
k
−
1
k
(
x
k
)
H
k
.
{\displaystyle H_{x,2}=\sum _{k=1}^{\infty }{\frac {(-1)^{k-1}}{k}}{x \choose k}H_{k}.}
для Ряд Тейлора чисел гармоник равен
H
x
=
∑
k
=
2
∞
(
−
1
)
k
ζ
(
k
)
x
k
−
1
for
|
x
|
<
1
{\displaystyle H_{x}=\sum _{k=2}^{\infty }(-1)^{k}\zeta (k)\;x^{k-1}\quad {\text{ for }}|x|<1}
который получается из ряда Тейлора для дигамма-функции (
ζ
{\displaystyle \zeta }
–
дзета-функция Римана ).
При попытке аппроксимировать H x для комплексного числа x эффективно сначала вычислить H m для некоторого большого целого числа m . Используйте это как приближение значения H m + x . Затем используйте рекурсивное соотношение H n = H n −1 + 1/ n назад m раз, чтобы развернуть его до приближения для H x . Более того, это приближение является точным в пределе, когда m стремится к бесконечности.
В частности, для фиксированного целого числа n это тот случай, когда
lim
m
→
∞
[
H
m
+
n
−
H
m
]
=
0.
{\displaystyle \lim _{m\rightarrow \infty }\left[H_{m+n}-H_{m}\right]=0.}
Если n не является целым числом, невозможно сказать, верно ли это уравнение, поскольку мы еще (в этом разделе) не определили числа гармоник для нецелых чисел. Однако мы получаем уникальное расширение гармонических чисел до нецелых чисел, настаивая на том, что это уравнение продолжает выполняться, когда произвольное целое число n заменяется произвольным комплексным числом x ,
lim
m
→
∞
[
H
m
+
x
−
H
m
]
=
0
.
{\displaystyle \lim _{m\rightarrow \infty }\left[H_{m+x}-H_{m}\right]=0\,.}
Поменяв местами порядок двух частей этого уравнения и затем вычитая их из
H x , получим
H
x
=
lim
m
→
∞
[
H
m
−
(
H
m
+
x
−
H
x
)
]
=
lim
m
→
∞
[
(
∑
k
=
1
m
1
k
)
−
(
∑
k
=
1
m
1
x
+
k
)
]
=
lim
m
→
∞
∑
k
=
1
m
(
1
k
−
1
x
+
k
)
=
x
∑
k
=
1
∞
1
k
(
x
+
k
)
.
{\displaystyle {\begin{aligned}H_{x}&=\lim _{m\rightarrow \infty }\left[H_{m}-(H_{m+x}-H_{x})\right]\\[6pt]&=\lim _{m\rightarrow \infty }\left[\left(\sum _{k=1}^{m}{\frac {1}{k}}\right)-\left(\sum _{k=1}^{m}{\frac {1}{x+k}}\right)\right]\\[6pt]&=\lim _{m\rightarrow \infty }\sum _{k=1}^{m}\left({\frac {1}{k}}-{\frac {1}{x+k}}\right)=x\sum _{k=1}^{\infty }{\frac {1}{k(x+k)}}\,.\end{aligned}}}
Этот бесконечный ряд сходится для всех комплексных чисел x, за исключением целых отрицательных чисел, что не удается, поскольку попытка использовать рекурсивное соотношение H n = H n −1 + 1/ n в обратном направлении через значение n = 0 включает деление на ноль. Согласно этой конструкции, функция, определяющая номер гармоники для комплексных значений, является единственной функцией, которая одновременно удовлетворяет условиям (1) H 0 = 0 , (2) H x = H x −1 + 1/ x для всех комплексных чисел x , кроме неположительные целые числа и (3) lim m →+∞ ( H m + x − H m ) = 0 для всех комплексных значений x .
Эту последнюю формулу можно использовать, чтобы показать, что
∫
0
1
H
x
d
x
=
γ
,
{\displaystyle \int _{0}^{1}H_{x}\,dx=\gamma ,}
где
γ —
константа Эйлера–Машерони или, в более общем смысле, для каждого
n имеем:
∫
0
n
H
x
d
x
=
n
γ
+
ln
(
n
!
)
.
{\displaystyle \int _{0}^{n}H_{x}\,dx=n\gamma +\ln(n!).}
Специальные значения для дробных аргументов [ править ]
Существуют следующие специальные аналитические значения для дробных аргументов от 0 до 1, определяемые интегралом
H
α
=
∫
0
1
1
−
x
α
1
−
x
d
x
.
{\displaystyle H_{\alpha }=\int _{0}^{1}{\frac {1-x^{\alpha }}{1-x}}\,dx\,.}
Дополнительные значения могут быть сгенерированы из рекуррентного отношения.
H
α
=
H
α
−
1
+
1
α
,
{\displaystyle H_{\alpha }=H_{\alpha -1}+{\frac {1}{\alpha }}\,,}
или из отношения отражения
H
1
−
α
−
H
α
=
π
cot
(
π
α
)
−
1
α
+
1
1
−
α
.
{\displaystyle H_{1-\alpha }-H_{\alpha }=\pi \cot {(\pi \alpha )}-{\frac {1}{\alpha }}+{\frac {1}{1-\alpha }}\,.}
Например:
H
1
2
=
2
−
2
ln
2
H
1
3
=
3
−
π
2
3
−
3
2
ln
3
H
2
3
=
3
2
+
π
2
3
−
3
2
ln
3
H
1
4
=
4
−
π
2
−
3
ln
2
H
3
4
=
4
3
+
π
2
−
3
ln
2
H
1
6
=
6
−
3
2
π
−
2
ln
2
−
3
2
ln
3
H
1
8
=
8
−
1
+
2
2
π
−
4
ln
2
−
1
2
(
ln
(
2
+
2
)
−
ln
(
2
−
2
)
)
H
1
12
=
12
−
(
1
+
3
2
)
π
−
3
ln
2
−
3
2
ln
3
+
3
ln
(
2
−
3
)
{\displaystyle {\begin{aligned}H_{\frac {1}{2}}&=2-2\ln 2\\H_{\frac {1}{3}}&=3-{\frac {\pi }{2{\sqrt {3}}}}-{\frac {3}{2}}\ln 3\\H_{\frac {2}{3}}&={\frac {3}{2}}+{\frac {\pi }{2{\sqrt {3}}}}-{\frac {3}{2}}\ln 3\\H_{\frac {1}{4}}&=4-{\frac {\pi }{2}}-3\ln 2\\H_{\frac {3}{4}}&={\frac {4}{3}}+{\frac {\pi }{2}}-3\ln 2\\H_{\frac {1}{6}}&=6-{\frac {\sqrt {3}}{2}}\pi -2\ln 2-{\frac {3}{2}}\ln 3\\H_{\frac {1}{8}}&=8-{\frac {1+{\sqrt {2}}}{2}}\pi -4\ln {2}-{\frac {1}{\sqrt {2}}}\left(\ln \left(2+{\sqrt {2}}\right)-\ln \left(2-{\sqrt {2}}\right)\right)\\H_{\frac {1}{12}}&=12-\left(1+{\frac {\sqrt {3}}{2}}\right)\pi -3\ln {2}-{\frac {3}{2}}\ln {3}+{\sqrt {3}}\ln \left(2-{\sqrt {3}}\right)\end{aligned}}}
Они вычисляются с помощью дигамм-теоремы Гаусса , которая, по сути, утверждает, что для натуральных чисел p и q с p < q
H
p
q
=
q
p
+
2
∑
k
=
1
⌊
q
−
1
2
⌋
cos
(
2
π
p
k
q
)
ln
(
sin
(
π
k
q
)
)
−
π
2
cot
(
π
p
q
)
−
ln
(
2
q
)
{\displaystyle H_{\frac {p}{q}}={\frac {q}{p}}+2\sum _{k=1}^{\lfloor {\frac {q-1}{2}}\rfloor }\cos \left({\frac {2\pi pk}{q}}\right)\ln \left({\sin \left({\frac {\pi k}{q}}\right)}\right)-{\frac {\pi }{2}}\cot \left({\frac {\pi p}{q}}\right)-\ln \left(2q\right)}
Римана с дзета - Связь функцией
Некоторые производные дробных чисел гармоник имеют вид
d
n
H
x
d
x
n
=
(
−
1
)
n
+
1
n
!
[
ζ
(
n
+
1
)
−
H
x
,
n
+
1
]
d
n
H
x
,
2
d
x
n
=
(
−
1
)
n
+
1
(
n
+
1
)
!
[
ζ
(
n
+
2
)
−
H
x
,
n
+
2
]
d
n
H
x
,
3
d
x
n
=
(
−
1
)
n
+
1
1
2
(
n
+
2
)
!
[
ζ
(
n
+
3
)
−
H
x
,
n
+
3
]
.
{\displaystyle {\begin{aligned}{\frac {d^{n}H_{x}}{dx^{n}}}&=(-1)^{n+1}n!\left[\zeta (n+1)-H_{x,n+1}\right]\\[6pt]{\frac {d^{n}H_{x,2}}{dx^{n}}}&=(-1)^{n+1}(n+1)!\left[\zeta (n+2)-H_{x,n+2}\right]\\[6pt]{\frac {d^{n}H_{x,3}}{dx^{n}}}&=(-1)^{n+1}{\frac {1}{2}}(n+2)!\left[\zeta (n+3)-H_{x,n+3}\right].\end{aligned}}}
А используя ряд Маклорена , мы имеем для x < 1, что
H
x
=
∑
n
=
1
∞
(
−
1
)
n
+
1
x
n
ζ
(
n
+
1
)
H
x
,
2
=
∑
n
=
1
∞
(
−
1
)
n
+
1
(
n
+
1
)
x
n
ζ
(
n
+
2
)
H
x
,
3
=
1
2
∑
n
=
1
∞
(
−
1
)
n
+
1
(
n
+
1
)
(
n
+
2
)
x
n
ζ
(
n
+
3
)
.
{\displaystyle {\begin{aligned}H_{x}&=\sum _{n=1}^{\infty }(-1)^{n+1}x^{n}\zeta (n+1)\\[5pt]H_{x,2}&=\sum _{n=1}^{\infty }(-1)^{n+1}(n+1)x^{n}\zeta (n+2)\\[5pt]H_{x,3}&={\frac {1}{2}}\sum _{n=1}^{\infty }(-1)^{n+1}(n+1)(n+2)x^{n}\zeta (n+3).\end{aligned}}}
Для дробных аргументов от 0 до 1 и для a > 1,
H
1
/
a
=
1
a
(
ζ
(
2
)
−
1
a
ζ
(
3
)
+
1
a
2
ζ
(
4
)
−
1
a
3
ζ
(
5
)
+
⋯
)
H
1
/
a
,
2
=
1
a
(
2
ζ
(
3
)
−
3
a
ζ
(
4
)
+
4
a
2
ζ
(
5
)
−
5
a
3
ζ
(
6
)
+
⋯
)
H
1
/
a
,
3
=
1
2
a
(
2
⋅
3
ζ
(
4
)
−
3
⋅
4
a
ζ
(
5
)
+
4
⋅
5
a
2
ζ
(
6
)
−
5
⋅
6
a
3
ζ
(
7
)
+
⋯
)
.
{\displaystyle {\begin{aligned}H_{1/a}&={\frac {1}{a}}\left(\zeta (2)-{\frac {1}{a}}\zeta (3)+{\frac {1}{a^{2}}}\zeta (4)-{\frac {1}{a^{3}}}\zeta (5)+\cdots \right)\\[6pt]H_{1/a,\,2}&={\frac {1}{a}}\left(2\zeta (3)-{\frac {3}{a}}\zeta (4)+{\frac {4}{a^{2}}}\zeta (5)-{\frac {5}{a^{3}}}\zeta (6)+\cdots \right)\\[6pt]H_{1/a,\,3}&={\frac {1}{2a}}\left(2\cdot 3\zeta (4)-{\frac {3\cdot 4}{a}}\zeta (5)+{\frac {4\cdot 5}{a^{2}}}\zeta (6)-{\frac {5\cdot 6}{a^{3}}}\zeta (7)+\cdots \right).\end{aligned}}}
Примечания [ править ]
^ Перейти обратно: а б
Джон Х., Конвей; Ричард К., Гай (1995). Книга чисел . Коперник.
^ Грэм, Рональд Л.; Кнут, Дональд Э.; Паташник, Орен (1994). Конкретная математика . Аддисон-Уэсли.
^ Сондоу, Джонатан и Вайсштейн, Эрик В. «Гармоническое число». Из MathWorld — веб-ресурса Wolfram. http://mathworld.wolfram.com/HarmonicNumber.html
^ Сандифер, К. Эдвард (2007), Как Эйлер сделал это , MAA Spectrum, Математическая ассоциация Америки, стр. 206, ISBN 9780883855638 .
^ Вайсштейн, Эрик В. (2003). CRC Краткая математическая энциклопедия . Бока-Ратон, Флорида: Chapman & Hall/CRC. п. 3115. ИСБН 978-1-58488-347-0 .
^ Эйзенштейн, Фердинанд Готхольд Макс (1850). «Новый род теоретико-числовых функций, которые зависят от двух элементов и определяются некоторыми линейными функциональными уравнениями». Об этом сообщает Royal. Пруссия. Академическая наука Берлин . 15 :36-42.
^ Ишваратасан, Арулаппа; Левин, Юджин (1991). «p-целые гармонические суммы» . Дискретная математика . 91 (3): 249–257. дои : 10.1016/0012-365X(90)90234-9 .
^ Бойд, Дэвид В. (1994). «Р-адическое исследование частичных сумм гармонического ряда» . Экспериментальная математика . 3 (4): 287–302. CiteSeerX 10.1.1.56.7026 . дои : 10.1080/10586458.1994.10504298 .
^ Санна, Карло (2016). «О p-адической оценке гармонических чисел» (PDF) . Журнал теории чисел . 166 : 41–46. дои : 10.1016/j.jnt.2016.02.020 . hdl : 2318/1622121 .
^ Чен, Юн-Гао; Ву, Бин-Лин (2017). «О некоторых свойствах гармонических чисел». Журнал теории чисел . 175 : 66–86. дои : 10.1016/j.jnt.2016.11.027 .
^ Джеффри Лагариас (2002). «Элементарная проблема, эквивалентная гипотезе Римана». амер. Математика. Ежемесячно . 109 (6): 534–543. arXiv : math.NT/0008177 . дои : 10.2307/2695443 . JSTOR 2695443 .
^ Э.О. Так (1964). «Некоторые методы обтекания тупых тонких тел». Дж. Гидромеханика . 18 (4): 619–635. Бибкод : 1964JFM....18..619T . дои : 10.1017/S0022112064000453 . S2CID 123120978 .
^ Сесма, Дж. (2017). «Возвращение к римским гармоническим числам» . Журнал теории чисел . 180 : 544–565. arXiv : 1702.03718 . дои : 10.1016/j.jnt.2017.05.009 . ISSN 0022-314X .
^ Леб, Дэниел Э; Рота, Джан-Карло (1989). «Формальный степенной ряд логарифмического типа» . Достижения в математике . 75 (1): 1–118. дои : 10.1016/0001-8708(89)90079-0 . ISSN 0001-8708 .
Внешние ссылки [ править ]
Эта статья включает в себя материал из числа Harmonic на сайте PlanetMath , который доступен под лицензией Creative Commons Attribution/Share-Alike License .