Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Расходящиеся ряды и асимптотические теории

.pdf
Скачиваний:
0
Добавлен:
11.08.2026
Размер:
632 Кб
Скачать

1. КОЕ-ЧТО ИЗ ИСТОРИИ РАСХОДЯЩИХСЯ РЯДОВ

11

и обозначим суммы наших двух рядов S0 и S1 соответственно. Имеем S0 = = 1 − S0, откуда S0 = 12 и

S1 = 1 2 + 3 4 + 5 6 + . . .

= 0 + 1 2 + 3 4 + 5 − . . . ,

откуда 2S1 = 1 1 + 1 1 + 1 1 + . . . и S1 = 14 . Итак,

1 1 + 1 1 + 1 1 + . . . = 12 , 1 2 + 3 4 + 5 6 + . . . = 14 .

По-видимому, Лейбниц был первым, кто приписал сумме 1 1 + 1 1 +

+ 1 1 + . . . значение 12 .

Все время рассматривая наши два примера, можно угадать два способа суммирования, кажущиеся разумными и позволяющие оправдать только что проведенное чисто формальное рассуждение. Эти два способа — дающие возможность больших обобщений — суммирование в среднем и суммирование Абеля.

Суммирование в среднем

Идея этого суммирования — в попытке построить положительную ме-

ру d полной массы 1 на подмножествах N и определить сумму ряда

an

как

 

 

 

=0

 

snd (n),

n

 

 

где sn обозначает частичную сумму

 

 

 

 

 

n

 

 

sn =

 

 

 

ap.

 

p=0

Результат не должен зависеть от первых членов ряда. Следовательно, мера должна обращаться в 0 на любом конечном подмножестве N, что приводит к противоречию, если предполагать меру σ-аддитивной! Итак, приходится ослабить требования и заменить понятие меры на понятие плотности [FGA]. Если — подобная плотность, обозначим (E) массу E подмножества N (если она существует). Проще всего построить такую плот-

12 РАСХОДЯЩИЕСЯ РЯДЫ И АСИМПТОТИЧЕСКИЕ ТЕОРИИ

ность, используя среднее арифметическое по Чезаро: распределить единичную массу между n + 1 точкой 0, . . . , n6 и перейти к пределу по n. Средние

tn = s0 + s1 + . . . + sn , n + 1

соответствующие ряду 1 1 + 1 1 + 1 1 + . . ., стремятся к 12 при n → +

(sn принимает значения 1 и 0 поочередно).

Если обозначить через E множество четных чисел, получим (E) =

= lim μn(E) = 12 .

Ряд 1 2 + 3 4 + 5 6 + . . . этим способом не суммируется. Но, применив этот способ дважды, получим в качестве суммы 14 .

В более общем случае для каждого целого n можно распределить единичную массу по некоторому конечному подмножеству N и предполагать, что если E конечно, то (E) = 0. Если, кроме того, предположить, что μn сосредоточена на отрезке [1, n], задача сведется к заданию правильной теплицевой матрицы, то есть бесконечной треугольной матрицы

 

a00

0

 

0 . . . . . . . . . . . . . . .

 

 

a10

a11

 

0 0 . . .

. . . . . . . . .

 

 

...

...

... ... ...

. . . . . . . . .

 

A =

.

.

.

..

.

..

.

..

 

,

 

..

..

 

 

 

. . . . . . . . .

 

 

 

an1

. . . . . . ann

 

 

 

an0

0 0 . . .

 

 

 

 

 

 

 

 

 

 

 

 

 

.

.

 

. . .

 

 

 

..

..

 

.. .. ..

. . . 0 . . .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

коэффициенты которой — положительные действительные числа, суммы по строкам тождественно равны 1, а столбцы представляют собой последовательности, сходящиеся к нулю. Тогда матрица-столбец, составленная из tn, получается умножением A на матрицу-столбец частичных сумм sn.

Композиция двух суммирований в среднем, очевидно, соответствует

произведению матриц. Например,

1 0 0 . . .

12 C = 13

...

1

0 . . .

3

3

 

2

.

 

.

1

1

. . .

 

 

 

 

 

 

..

..

. . .

6Соответствующая мера равна μn = (δ0 + . . . + δn)/(n + 1).

1. КОЕ-ЧТО ИЗ ИСТОРИИ РАСХОДЯЩИХСЯ РЯДОВ

13

соответствует методу Чезаро, а

 

 

 

1

0

0 . . .

 

 

 

3

1

0 . . .

 

 

 

4

4

.

.

 

C2

=

.

.

 

11

5

 

2

. . .

 

 

 

 

 

 

 

18

18

18

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..

..

..

 

..

— его использованной выше итерации.

Преобразование Эйлера соответствует теплицевой матрице:

 

 

C1

 

 

0

 

T =

1 C10

 

2 ...

 

 

 

 

 

 

 

 

1

 

0

 

 

 

Cn

 

n

 

 

2

 

 

 

.

 

 

.

 

 

.

 

 

 

 

 

 

 

0

 

 

0

. . .

1 C1

 

0

 

0

2

 

.

1

.

 

 

.

 

 

 

 

 

.

 

.

 

 

 

.

 

 

.

 

.

 

 

.

 

 

 

 

 

1

 

C1 . . .

1

Cn

n

n

 

n

 

 

 

 

n

2

 

.

 

 

.

 

2

.

 

 

 

 

 

 

 

 

 

 

.

 

 

.

 

 

.

 

 

 

.

 

 

.

 

 

.

 

. . .

. . .

... .

0

...

Это преобразование усредняет ряды путем преобразования

b

n

=

1

(C0 a

0

+ C1 a

1

+ . . . + Cna

).

 

 

 

2n+1

n

n

n n

 

 

 

 

 

 

 

 

 

 

Оно усредняет формальные ряды, связанные с гомографическим преобразованием (считая 1 неподвижной точкой):

 

 

 

 

 

 

 

 

1 y

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

f (x) = n=0 anxn,

gˆ(y) = n=0 bnyn,

gˆ(y) = f 1

2

1 y

Это преобразование было введено Эйлером как ускоряющее сходи-

мость. Например, для ряда

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

1

 

 

 

 

 

 

 

 

1 2

+

3 − . . .

 

 

 

 

 

оно дает

1

 

1

 

 

1

 

 

 

 

 

 

 

+

 

+

+ . . . ,

 

 

 

 

 

 

 

 

233

 

 

 

 

 

 

2

222

 

 

 

 

 

 

 

 

что ускоряет вычисление log 2.

14

РАСХОДЯЩИЕСЯ РЯДЫ И АСИМПТОТИЧЕСКИЕ ТЕОРИИ

Вместо того чтобы использовать меры μn с параметром n N и конечным носителем, можно использовать меры μt с параметром t > 0 и произвольным носителем. Основным примером служит плотность Бореля, свя-

занная с семейством мер Пуассона:

μt = e−t tn δn.

n=0 n!

Этот метод, возникший благодаря Борелю (1899), сопоставляет ряду

an в качестве суммы предел (если он существует)

n=0

lim e−tS(t),

где S(t) =

sn tn .

 

 

 

t→+

 

n=0 n!

В дальнейшем мы вернемся к этому фундаментальному методу. Применяя его к «ряду Лейбница» 1 1 + 1 1 + . . ., найдем

lim e−t ch t = 1 .

t→+2

Методы Абеля

Описав методы усреднения, переходим к методам Абеля. Самый про-

 

 

стой из них основан на теореме Абеля: если ряд

an сходится, то его

сумма задается равенством

 

=0

n

lim

 

 

anxn.

 

x→1

 

 

x<1 n=0

Отсюда идея — которую можно найти у Эйлера — для некоторых расходя-

щихся рядов принять за определение суммы ряда

an предел (если он

существует)

 

 

 

 

n=0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

anxn.

 

 

 

 

 

 

lim

 

 

 

 

 

 

x→1

 

 

 

 

 

 

 

x<1 n=0

 

 

 

 

 

 

 

Применяя этот метод к ряду Лейбница, находим

 

 

 

 

 

 

 

 

 

 

 

 

1

 

lim

 

1

=

1

 

n n

 

и

 

 

,

(1) x =

1 + x

 

x→1 1 + x

2

 

n=0

 

 

x<1

 

 

 

 

 

 

 

 

 

 

 

 

как и ожидалось.

1. КОЕ-ЧТО ИЗ ИСТОРИИ РАСХОДЯЩИХСЯ РЯДОВ

15

Для ряда 1 2 + 3 4 + 5 6 + . . . находим

 

2

1

 

 

1 2x + 3x

+ . . . =

 

 

 

(1 + x)2

 

с пределом 14 .

Этот метод обобщается следующим образом. Полагаем x = e−t, тогда xn = e−tn. Пусть n} — данная строго возрастающая последователь-

ность строго положительных действительных чисел, стремящаяся к +.

Сопоставляем ряду an функцию

n=0

ft(x) = ane−tλn

n=0

(если она существует для всех достаточно малых t > 0) и — обязательно при условии существования — определяем сумму как

lim ane−tλn .

t→0

t>0 n=0

В дальнейшем мы остановимся на следующих примерах:

λn =

0

если n = 0

(Линделеф);

 

n log n

если n > 0

 

 

 

 

 

 

0

если n = 0, 1, 2

 

λn = n log n log(log n) если n > 2

(Харди [H2], 1941).

1.2. Функциональное уравнение для ζ-функции. Ряд Эйлера

Более детальное рассмотрение изложенных ниже вопросов можно найти в [H1] и [Bar]. Мемуар Л. Эйлера [Eu1] начинается так:

Соотношение, которое я собираюсь здесь рассмотреть, касается сумм следующих двух обобщенных бесконечных рядов:

1m 2m + 3m 4m + 5m 6m + 7m 8m + . . . , 11m 21m + 31m 41m + 51m 61m + 71m 81m + . . . .

16 РАСХОДЯЩИЕСЯ РЯДЫ И АСИМПТОТИЧЕСКИЕ ТЕОРИИ

Задачей Эйлера было изучение равенства

 

 

1s−1

2s−1 + 3s−1 . . .

 

(s − 1)!(2s 1)

 

1

 

 

 

 

 

 

 

.. .

 

=

 

cos

2 sπ .

 

1

 

1

+

1

(2s−1 1)πs

 

s

s

s

 

 

1 2

3

 

 

 

 

 

 

Точнее, он собирался «проверить» это равенство для всех целых s, для s = 12 и s = 32 . Для этого нужно суммировать расходящиеся ряды

1m − −2m + 3m 4m + 5m 6m + 7m 8m + . . . .

Эйлер использует метод Абеля. Для m = 0 и 1 он получил два изученных выше ряда с суммами 12 и 14 соответственно. В более общем случае он получает при k > 1

12k 22k + 32k − . . . = 0,

12k−1 22k−1 + 32k−1 − . . . = (1)k−1 22k 1 Bk . 2k

В случае Re s > 0 ζ-функция Римана определена равенством

ζ(s) = 1−s + 2−s + 3−s + . . . .

Сопоставим ей функцию

η(s) = 1−s 2−s + 3−s − . . . .

Имеем

η(s) = (1 21−S )ζ(s),

и знаменитое функциональное равенство Б. Римана для ζ-функции запишется в виде

(2s−1 1)η(1 − S) = (2s 1)π−s cos 12 Γ(s)η(s).

Мы узнаем тождество, угаданное Эйлером за сто лет до Римана.

Заметим, что теоретический подсчет, основанный на использовании сумм расходящихся рядов, приводит к строгому доказательству формулы в случае целых s. (Эйлер не намеревался приводить такое доказательство.)

1. КОЕ-ЧТО ИЗ ИСТОРИИ РАСХОДЯЩИХСЯ РЯДОВ

17

Метод, которым Эйлер обосновывал возможность вычислений с помощью расходящихся рядов, во многом взят из его исследования ряда Эйлера:

(1)nn!.

n=0

Это обоснование, по большей части «экспериментальное», очень любопытно. Много деталей, связанных с этим примером, можно найти в [Bar]. В современных терминах общую идею Эйлера можно изложить следующим образом: пусть дан ряд с общим членом an; следует найти программу (или программы), порождающие an. Процитируем Эйлера (1749):

Но я уже отмечал по другому поводу, что слову «сумма» надо придать более широкое значение и понимать под ним дробь или другое аналитическое выражение, которое, будучи разложено [в ряд] согласно принципам анализа, приводит к тому же самому ряду, чью сумму мы ищем.

В другом месте он пишет:

...summa cujusque seriei est valor expressionis illius finitae, ex cujus evolutione illa series oritur7.

В [Eu2] Эйлер описал четыре (новых) способа для суммирования ряда

(1)nn!. Путем прямого подсчета можно проверить, как согласуются

n=0

эти разные приближения. Мы не будем описывать эти способы (что очень хорошо сделано в [Bar]). В связи с ними нас больше всего интересует тот факт, что соответствующий формальный ряд

f (x) = (1)nn!xn+1

n=0

является формальным решением дифференциального уравнения Эйлера

x2y + y = x.

Среди других методов рассмотрим преобразование в непрерывную сходящуюся дробь

11 + 11 + 11 + 21 + 21 + 31 + 31 + . . . ,

7. . . сумма некоторого ряда — это значение конечного выражения, чье разложение приводит

кданному ряду.

18

РАСХОДЯЩИЕСЯ РЯДЫ И АСИМПТОТИЧЕСКИЕ ТЕОРИИ

и итерация преобразования Эйлера следует из суммирования до наименьшего члена, описанного в следующем параграфе. (Этот метод очень хорошо объяснен в трактате по анализу, принадлежащем Лакруа.) Метод, основанный на дифференциальном уравнении Эйлера, связан с суммированием Бо- реля–Лапласа и будет играть основополагающую роль в дальнейшем. (Само же суммирование Бореля–Лапласа связано с борелевской плотностью, описанной выше [H1].)

Позднее Харди [H2] суммировал ряд Эйлера, используя описанный выше метод Абеля:

 

 

 

 

 

+

 

 

lim

 

 

1!x2

+ 2!x3

 

(

1)nn!xn+1e−tn log n log(log n) .

x

+

t→0

 

 

 

 

 

t>0

 

 

 

n=3

1.3.Эйлер, Коши, Пуанкаре и суммирование до наименьшего члена

Вначале главы VIII второго тома «Новых методов небесной механики» [P] А. Пуанкаре пишет:

Между геометрами и астрономами есть определенного рода недоразумения в связи со значением слова сходимость. Геометры, озабоченные

совершенной строгостью и зачастую слишком безразличные к громоздкости запутанных вычислений, с помощью которых они обосновывают возможность, не задумываясь о ее эффективном использовании, говорят, что ряд сходится, когда сумма членов стремится к определенному пределу, даже если первые члены убывают слишком медленно. Астрономы, напротив, имеют привычку говорить, что ряд сходится, когда, к примеру, двадцать первых членов убывают очень быстро, даже если следующие члены должны неопределенно возрастать. Так, в качестве простого примера рассмотрим два ряда с общими членами 1000n/n! и n!/1000n.

Геометры скажут, что первый ряд сходится, причем быстро,. . .; но второй ряд они сочтут расходящимся. . .

Астрономы, напротив, сочтут первый ряд расходящимся,. . . а второй — сходящимся.

Оба эти правила корректны: первое — для теоретических исследований, второе — для численных приложений.

Итак, можно говорить о сходящихся рядах «в смысле геометров» или «в смысле астрономов». Заметим, что на практике, в приложениях, можно констатировать, что прежде всего сходящиеся «в смысле астрономов» ряды имеют очень быстро возрастающий общий член, хотя изначально он убывал. Итак, то, что Пуанкаре рассматривал как возможность, на самом

1. КОЕ-ЧТО ИЗ ИСТОРИИ РАСХОДЯЩИХСЯ РЯДОВ

19

деле — правило. Он не отдавал себе отчет в этом явлении — хотя был очень хорошо с ним знаком, — приводя определение асимптотических разложений. Наибольшие недостатки теории Пуанкаре связаны с этим замечанием. Чтобы преодолеть их, надо обратиться к работам других математиков. Уже

упомянутый нами ряд Эйлера (1)nn!xn+1 (смежный с рядом, на кото-

n=0

рый ссылается Пуанкаре) позволит нам начать понимать, что надо делать. В конце концов мы получим теорию, точную асимптотическую теорию, в которой противоречие между двумя отмеченными Пуанкаре понятиями сходимости исчезнет.

Следующий анализ ряда Эйлера считается классическим (ср., например, Олвера [O]). Пусть

+

e−t/x

f (x) = 0

1 + t dt.

Заметим, что

 

fn(x) = x − 1!x2 + 2!x3 + . . . + (1)n−1(n − 1)!xn,

 

+

tne−t/x

 

n

Rn(x) = (1)

0

1 + t dt.

Легко проверить, что

 

 

 

 

f (x) = fn(x) + Rn(x),

 

 

+

|Rn(x)| < 0

tne−t/xdt = n!xn+1.

Итак, остаток Rn(x) имеет тот же знак, что и первое из опущенных слагаемых (1)nn!xn, и его модуль мажорируется модулем этого слагаемого. Рассуждая как в случае знакочередующегося ряда при проверке «особого критерия», получаем оценки

f2p(x) < f (x) < f2p+1(x)

для всех целых p > 1.

20

РАСХОДЯЩИЕСЯ РЯДЫ И АСИМПТОТИЧЕСКИЕ ТЕОРИИ

Если x > 0 фиксировано, то это не позволяет (как в случае проверки особого критерия для ряда, общий член которого стремится к нулю) получить сколь угодно точную оценку для f (x) (здесь общий член быстро стремится к +). Но мы получаем оценку, качество которой зависит от x. Точнее, мы увидим, что она экспоненциально хороша, когда x «мало».

Ясно, что при фиксированном x > 0 наилучшее приближение f (x) последовательностью fn(x) получается (этим методом), когда скачок

 

 

 

 

|f2p+1(x) − f2p(x)| = (2p)!x2p

наименьший

из возможных, то есть когда общий член ряда Эйлера

n

n!x

n+1

наименьший из возможных по абсолютной величине.

n=0(1)

 

 

В таком случае есть смысл остановить суммирование на соответствующем индексе: это и есть суммирование до наименьшего члена.

Здесь отношение двух последовательных членов равно nx по абсолютной величине. Если фиксированное x > 0 «мало», а n меняется, то это соотношение меньше 1 при n < x1 1 и больше 1 при n > x1 1. Итак,

общий член убывает до n = N x1, а затем неопределенно возрастает. При малых x получаем ряд, сходящийся в смысле астрономов. Качество аппроксимации равно N !xN N !N −N . По формуле Стирлинга,

N ! 2πN N N e−N ,

и мы видим, что порядок аппроксимации равен

2πN e−N 2π/xe1/x,

то есть экспоненциален.

Следуя концепции, выдвинутой Пуанкаре в 1886 году (раздел 1.7), скажем, что ряд Эйлера

f (x) =

 

(1)nn!xn+1

 

n=0

 

 

 

это асимптотическое разложение функции

 

+

e−t/x

f (x) = 0

1 + t dt

в начале координат.

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]