Суммирование расходящихся рядов методами Абеля, Бореля, Чезаро и Дирихле
Перевод поста Давендра Кападия (Devendra Kapadia) "The ABCD of Divergent Series."
Выражаю благодарность за помощь в переводе Андрею Дудину.
Какова сумма всех натуральных чисел? Интуиция подсказывает, что ответ — бесконечность. В математическом анализе сумма натуральных чисел является простым примером расходящегося ряда. Тем не менее, математики и физики сочли полезным придать дробные, отрицательные и даже нулевые значения суммам таких рядов. Цель моей статьи — желание отодвинуть завесу тайны, окружающую результаты суммирования расходящихся рядов. В частности, я буду использовать функцию Sum (функция поиска частичных сумм, рядов и т. п. в Mathematica), а так же другие функции в Wolfram Language для того, чтобы объяснить в каком смысле стоит рассматривать следующие утверждения:
Важность обозначений формул буквами A, B, C, и D вскоре станет вам понятна.
Начнем с того, что напомним понятие сходящегося ряда, используя следующую бесконечно убывающую геометрическую прогрессию.
Общий член ряда, начиная с n = 0, определяется по формуле:
In[1]:=
Теперь зададим сумму членов ряда от i = 0 до некоторого конечного значения i = n.
In[2]:=
Эта конечная сумма называется частичной суммой ряда.
График значений таких частичных сумм показывает, что их значения приближаются к числу 2 с ростом n:
In[3]:=
Out[3]=
Применяя функцию Limit (поиск предела последовательности или функции в точке) найдем предел значения частичных сумм этого ряда при стремлении n к бесконечности, что подтвердит наши наблюдения.
In[4]:=
Out[4]=
Функция Sum даёт такой же результат, когда мы производим суммирование членов ряда в пределах от 0 до бесконечности.
In[5]:=
Out[5]=
Мы говорим, что данный ряд (сумма данной бесконечно убывающей геометрической прогрессии) сходится и что его сумма равна 2.
Вообще, бесконечный ряд сходится, если последовательность его частичных сумм стремится к некоторому значению при неограниченном увеличении номера частичной суммы. В этом случае, предельное значение частичных сумм называется суммой ряда.
Бесконечный ряд который не сходится называется расходящимся. По определению, сумма расходящегося ряда не может быть найдена с помощью рассмотренного выше метода частичных сумм. Тем не менее, математики разработали различные способы присваивания конечных числовых значений суммам этих рядов. Такая сумма называется регуляризованной суммой расходящегося ряда. Процесс вычисления регуляризованных сумм называется регуляризацией.
Теперь мы рассмотрим пример A из вступления.
“A” обозначает Абеля, знаменитого норвежского математика, который предложил одну из техник регуляризации расходящихся рядов. В ходе своей короткой жизни, он умер всего в 26 лет, Абель достиг впечатляющих результатов в решении одних из самых трудных математических задач. В частности, он показал, что решение алгебраического уравнения пятой степени не может быть найдено в радикалах, поставив тем самым точку в проблеме, которая оставалась нерешенной на протяжении 250 лет до него.
Для того чтобы применить метод Абеля, заметим, что общий член данного ряда имеет вид:
In[6]:=
Это можно легко проверить, найдя несколько первых значений a[n].
In[7]:=
Out[7]=
Как можно увидеть на графике ниже, частичные суммы ряда принимают значения, равные 1 или 0 в зависимости от того, четное n или нечетное.
In[8]:=
Out[8]=
Естественно, что функция Sum выдает сообщение, о том что ряд расходится.
In[9]:=
Out[9]=
Регуляризация Абеля может быть применена к этому ряду в два шага. Сначала мы строим соответствующий степенной ряд.
In[10]:=
Out[10]=
Затем мы берем предел этой суммы при x стремящемся к 1, заметим при этом, что соответствующий ряд сходится для значений x меньших, но не равных 1.
In[11]:=
Out[11]=
Эти два шага можно объединить, сформировав, по сути, определение суммы расходящегося ряда по Абелю.
In[12]:=
Out[12]=
Мы можем получить тот же ответ используя опцию Regularization для функции Sum следующим образом.
In[13]:=
Out[13]=
Значение 1/2 представляется разумным, так как оно является средней величиной из двух значений, 1 и 0, принимаемых частичной суммой данного ряда. Кроме того, используемый в данном методе предельный переход интуитивно понятен, т. к. при x = 1 степенной ряд совпадает с нашим расходящимся рядом. Однако, Абель был сильно обеспокоен отсутствием строгости, которое было присуще математическому анализу того времени, и выражал свою обеспокоенность об этом:
«Расходящиеся ряды — изобретение дьявола, и это стыдно на них ссылаться при каких бы то ни было доказательствах. С их помощью, можно сделать любой вывод, какой ему будет угоден, и именно поэтому эти ряды производят столько ошибок и столько парадоксов.» (Н. Х. Абель в письме к своему бывшему учителю Берндту Хольмбою, Январь 1826)
Обратимся теперь к примеру B, в котором утверждается, что:
“B” обозначает Бореля, французского математика, который работал в таких областях как теория меры и теория вероятностей. В частности, Борель связан с так называемой “теоремой о бесконечных обезьянах”, которая утверждает, что если абстрактная обезьяна будет случайным образом ударять по клавиатуре пишущей машинки на протяжении бесконечного количества времени, то вероятность того, что она напечатает некоторый конкретный текст, например, полное собрание сочинений Уильяма Шекспира, отлична от нуля.
Для того чтобы применить метод Бореля заметим, что общий член данного ряда имеет вид:
In[14]:=
Регуляризация Бореля может быть применена к быстро расходящимся рядам в два шага. На первом шаге мы вычисляем экспоненциальную производящую функцию для последовательности членов данного ряда. Стоящий в знаменателе факториал обеспечивает сходимость данного ряд при всех значениях параметра t.
In[15]:=
Out[15]=
Затем мы производим преобразование Лапласа нашей экспоненциальной производящей функции и ищем его значение в точке s=1.
In[16]:=
Out[16]=
Out[17]=
Эти шаги можно объединить, в итоге мы получим, по сути, определение суммы расходящегося ряда по Борелю.
In[18]:=
Out[18]=
Также мы можем использовать специализированные функции Wolfram Language для поиска экспоненциальной производящей функции и преобразования Лапласа:
In[19]:=
Out[19]=
При этом, ответ можно получить непосредственно с помощью Sum следующим образом.
In[20]:=
Out[20]=
Определение суммы по Борелю разумно, т. к. оно даёт тот же самый результат, что и обычный метод частичных сумм, если его применить к сходящемуся ряду. В этом случае можно поменять местами суммирование и интегрирование, и затем определить Гамма-функцию, при этом мы получим, что соответствующий интеграл будет равен 1 и останется просто, по сути, исходная сумма ряда:
In[21]:=
Out[21]=
Однако в случае с расходящимися рядами поменять местами знаки суммы и интеграла нельзя, что приводит к интересным результатам, которые даёт данный метод регуляризации.
Суммирование по Борелю представляет собой универсальный метод суммирования расходящихся рядов, который применяется, скажем, в квантовой теории поля. О применении суммирования по Борелю существует огромная коллекция литературы.
Пример C утверждает что:
“C” обозначает Чезаро (на англ. языке его фамилия пишется как Cesaro), итальянского математика, который внес значительный вклад в дифференциальную геометрию, теорию чисел и математическую физику. Чезаро был очень продуктивным математиком и написал около 80 работ в период с 1884 по 1886 г., до того, как получил степень PhD в 1887!
Для начала заметим, что общий член ряда, начиная с n = 0, имеет вид:
In[22]:=
График показывает сильную осцилляцию частичных сумм данного ряда.
In[23]:=
Out[23]=
Метод Чезаро использует последовательность средних арифметических значений частичных сумм ряда для того, чтобы подавить осцилляции, что демонстрирует следующий график.
In[24]:=
Out[24]=
Формально говоря, суммирование по Чезаро определяется как предел последовательности средних арифметических значений частичных сумм ряда. Вычисляя данный предел для ряда из примера C, мы получим ожидаемый нами результат -1/2 (см. график выше).
In[25]:=
Out[25]=
Сумма по Чезаро может быть получена непосредственно, если мы в функции Sum используем данный тип регуляризации, указав соответствующее значение опции Regularization.
In[26]:=
Out[26]=
Метод суммирования по Чезаро играет важную роль в теории рядов Фурье, в которых ряды на основе тригонометрических функций используются для представления периодических функций. Ряд Фурье для непрерывной функции может и не сходится, но соответствующая сумма по Чезаро (или чезаровское среднее, как её обычно называют) всегда будет сходиться к функции. Этот красивый результат называется теоремой Фейера.
Наш последний пример утверждает, что сумма натурального ряда равна -1/12.
“D” означает Дирихле, немецкого математика, который совершил огромный вклад в теорию чисел и ряд других областей математики. О широте вкладов Дирихле можно судить, просто введя в Mathematica 10 следующий код.
In[27]:=
Out[27]//TableForm=
Регуляризация по Дирихле получила свое название от понятия “ряд Дирихле”, который определяется следующим образом:
Специальным случаем данного ряда является дзета-функция Римана, которую можно определить так:
In[28]:=
In[29]:=
Out[29]=
Функция SumConvergence говорит нам, что этот ряд сходится в том случае, если действительная часть параметра s будет больше 1.
In[30]:=
Out[30]=
Однако, сама по себе дзета-функция Римана может быть определена и для других значений параметра s с помощью процесса аналитического продолжения, известного из теории функций комплексного переменного. Например, при s = -1, мы получим:
In[31]:=
Out[31]=
Но при s = -1, ряд, задающий дзета-функцию Римана и есть натуральный ряд. Отсюда мы и получаем, что:
In[32]:=
Out[32]=
Еще один способ осознания этого результата заключается в том, чтобы ввести бесконечно малый параметр ε в выражение члена нашего расходящегося ряда, а затем найти разложение полученной функции в ряд Маклорена с помощью функции Series, как показано ниже.
In[33]:=
Out[33]=
Первое слагаемое
Аналогично можно получить безумно странное значение 0 для расходящейся суммы квадратов натуральных чисел.
In[34]:=
Out[34]=
В этом случае в соответствующем разложении отсутствуют члены, не зависящие от параметра ε, в результате мы получаем 0.
In[35]:=
Out[35]=
Регуляризация Дирихле тесно связана с процессом дзета регуляризации, который используется в современной теоретической физике. В своей знаменитой работе, выдающийся британский физик Стивен Хокинг применил данный метод к задаче вычисления Фейнмановых интегралов в искривленном пространстве-времени. Статья Хокинга описывает процесс дзета-регуляризации очень системно и она приобрела большую популярность после публикации.
Наши знания о расходящихся рядах основаны на глубочайших теориях, разработанных одними из лучших мыслителей последних нескольких столетий. Тем не менее, я соглашусь со многими читателям, которые как и я, чувствуют некоторое непонимание, когда они видят их в современных физических теориях. Великий Абель, вероятно, был прав, когда назвал данные ряды “изобретением дьявола”. Не исключено, что какой-то будущий Эйнштейн, обладающий умом, свободным от всяческих устоев и авторитетов, отбросит преобладающие научные убеждения и переформулирует фундаментальную физику так, что в ней не не будет места для расходящихся рядов. Но даже если такая теория станет реальностью, расходящиеся ряды все равно будут давать нам богатый источник математических идей, освещая дорогу к более глубокому пониманию нашей Вселенной.