Интегралы Борвейна: паттерн, который ломается — и почему это не баг

Интегралы Борвейна: паттерн, который ломается — и почему это не баг

Спикер: 3Blue1Brown | Длительность: 17:25

Транскрипт

Загадочный паттерн

Перед нами последовательность вычислений. Каждое из них — интеграл, и каждое оказывается равным числу π. Паттерн выглядит так, будто он будет продолжаться вечно. Но в определённый момент он ломается: вместо π мы получаем значение, которое лишь едва-едва меньше π.

Главный герой — функция sinc(x) = sin(x)/x. Это осциллирующая синусоида, затухающая к краям. В точке x = 0 предел равен 1, и мы определяем sinc(0) = 1 для непрерывности.

Интеграл sinc от минус до плюс бесконечности равен ровно π. Если мы умножим sinc(x) на sinc(x/3) — растянутую версию, — интеграл по-прежнему π. Добавим sinc(x/5) — снова π. И так далее через нечётные числа.

Но на 15-м множителе паттерн ломается. Точное значение — дробь от π, где и числитель, и знаменатель — числа порядка 400 миллиардов миллиардов. Отец и сын Борвейны описали это в статье. Коллега-исследователь, обнаруживший это на компьютере, решил, что это баг системы. Но это не баг — это реальное явление.

Аналогия со скользящим средним

Рассмотрим прямоугольную функцию rect(x): она равна 1 на отрезке [-1/2, 1/2] и 0 вне его. Применим к ней скользящее среднее с окном шириной 1/3. Получим сглаженную функцию с плато, шириной 1 - 1/3 = 2/3.

Применим ещё одно скользящее среднее с окном 1/5 — плато сузится ещё на 1/5. Затем 1/7, 1/9, 1/11, 1/13... Каждый раз плато уменьшается.

Критический момент: когда сумма 1/3 + 1/5 + 1/7 + ... + 1/13 всё ещё меньше 1, плато существует и значение функции в нуле равно ровно 1. Но при добавлении 1/15 сумма превышает 1, плато исчезает, и значение в нуле становится чуть меньше 1.

Этот коэффициент, чуть меньший единицы, — в точности тот же множитель, который стоит перед π в «сломанном» интеграле. Два явления количественно идентичны.

Случай с 2cos(x): начальное плато имеет ширину 2, и сумма 1/3 + 1/5 + ... должна превысить 2 — это происходит при 1/113.

Связь через преобразование Фурье

Подставим x → πx: задача «доказать интеграл = π» превращается в «доказать интеграл = 1». Нормализованная sinc-функция с π внутри аргумента имеет площадь 1.

Три ключевых факта:

1. Sinc и rect связаны преобразованием Фурье — одна является Фурье-образом другой.
2. Интеграл функции от -∞ до +∞ равен значению её Фурье-образа в нуле.
3. Теорема о свёртке: Фурье-образ произведения функций = свёртка их Фурье-образов.

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

Практические задания

Задание 1: Численная проверка интегралов Борвейна

Используйте Python с библиотекой scipy для вычисления последовательности интегралов. Начните с ∫sinc(x)dx, затем ∫sinc(x)·sinc(x/3)dx, и так далее. На каждом шаге сравнивайте результат с π. Зафиксируйте шаг, на котором равенство нарушается, и величину отклонения.

Задание 2: Визуализация скользящих средних

Реализуйте на Python последовательное применение скользящего среднего к прямоугольной функции. Постройте 8 графиков на одном полотне, показывая сужение плато на каждом шаге. Отметьте значение f(0) на каждом графике.

Задание 3: Кумулятивные суммы и пороги

Вычислите кумулятивные суммы ряда 1/3 + 1/5 + 1/7 + ... до 30 членов. Найдите момент превышения порогов 1 и 2. Сопоставьте с порогами ломки паттерна (15 и 113).

Задание 4: Преобразование Фурье функции sinc

Вычислите FFT дискретизированной sinc-функции и покажите, что результат — прямоугольная функция. Поэкспериментируйте с sinc(x/k) для разных k и покажите, как меняется ширина «прямоугольника».

Задание 5: Эксперимент с произвольными последовательностями

Замените обратные нечётные числа на другую последовательность (например, 1/2, 1/4, 1/8 или 1/π, 1/e, 1/√2). Предскажите порог ломки по кумулятивной сумме и проверьте численно через интегралы.

Лучшие цитаты

«Иногда кажется, что вселенная просто издевается над тобой» — 3Blue1Brown

«Это не баг, это реальное явление» — 3Blue1Brown

«Когда коллега-исследователь вычислял эти интегралы с помощью системы компьютерной алгебры, он решил, что это какой-то баг» — 3Blue1Brown

«Вы не можете просто случайным образом подставлять новые множители в интеграл и ожидать, что ничего не изменится» — 3Blue1Brown

«Два явления не просто качественно похожи — они количественно совпадают» — 3Blue1Brown

«Это систематический способ сменить перспективу, при которой трудные задачи иногда выглядят проще» — 3Blue1Brown

«Вы получаете огромную отдачу от одного вычисления» — 3Blue1Brown

«Когда паттерн ломается на значении 15, он делает это на ничтожно малую величину» — 3Blue1Brown