Whitney and Carbon Асимптотика

 Как приблизить сложную кривую простым многочленом? Эстафета идей от Уитни до искусственного интеллекта

Представьте, что перед вами на экране нарисован сложный график некоторой непрерывной функции \( f(x) \) на отрезке \( [0, 1] \). Это может быть профиль крыла самолета, запись звуковой волны или колебание температуры. Мы хотим заменить эту сложную функцию простым алгебраическим многочленом \( P_{n-1}(x) \) степени не выше \( n-1 \):

\[ P_{n-1}(x) = c_0 + c_1 x + c_2 x^2 + \dots + c_{n-1} x^{n-1}. \]

Но как заранее узнать, насколько точно многочлен способен повторить форму функции \( f(x) \)?

Если функция прямая линия, ее идеально приблизит многочлен первой степени. Чтобы проверить, является ли функция прямой, достаточно взять три точки на равном расстоянии \( h \) друг от друга — скажем, \( x \), \( x+h \) и \( x+2h \) — и проверить, лежит ли средняя точка на отрезке между крайними. Величина

\[ \Delta_h^2 f(x) = f(x) - 2f(x+h) + f(x+2h) \]

называется второй конечной разностью. На любой прямой она равна нулю.

Точно так же для любого натурального числа \( n \) можно составить конечную разность \( n \)-го порядка с шагом \( h \), коэффициенты которой берутся из треугольника Паскаля с чередующимися знаками:

\[ \Delta_h^n f(x) = \sum_{j=0}^n (-1)^j \binom{n}{j} f(x + jh). \]

Замечательное свойство этой суммы состоит в том, что для любого многочлена степени меньше \( n \) она тождественно равна нулю при любых \( x \) и \( h \). А для произвольной функции \( f \) максимальное по модулю значение такой разности, когда все ее точки \( x, x+h, \dots, x+nh \) помещаются внутри отрезка \( [0, 1] \), называется модулем непрерывности \( n \)-го порядка и обозначается \( \omega_n(f) \):

\[ \omega_n(f) = \sup_{x, x+nh \in [0, 1]} |\Delta_h^n f(x)|. \]

Величина \( \omega_n(f) \) — это своего рода «детектор кривизны» \( n \)-го порядка. Если она мала, функция на любом участке ведет себя почти как многочлен степени \( n-1 \). Возникает естественный вопрос: верно ли, что тогда и на всем отрезке \( [0, 1] \) сразу функцию \( f(x) \) можно приблизить одним многочленом с ошибкой, не превышающей \( W_n \cdot \omega_n(f) \), где \( W_n \) — некоторое число, зависящее только от порядка \( n \)?

1957 год: Хасслер Уитни и Арне Берлинг. Задача обретает форму

Для \( n=1 \) (приближение константой) и \( n=2 \) (приближение прямой) ответ прост и был найден английским математиком Х. Буркиллом в 1952 году: \( W_1 = W_2 = 1/2 \). Буркилл предположил, что конечные константы \( W_n \) существуют для всех \( n \).

В 1957 году выдающийся американский математик Хасслер Уитни доказал эту гипотезу. В его работе проявилось удивительное противоречие между бесконечной полупрямой и конечным отрезком. Уитни привел изящное интегральное тождество, подсказанное ему шведским математиком Арне Берлингом. Если функция задана на всей бесконечной полуоси \( [0, +\infty) \), то можно взять разность \( \Delta_h^n f(x) \) и усреднить ее по шагу \( h \) на огромном интервале от \( 0 \) до \( A \). При \( A \to \infty \) почти все промежуточные значения функции взаимно уничтожаются, и в две строчки получается оценка колебания функции с константой 2 для всех \( n \) сразу!

Но на конечном отрезке \( [0, 1] \) шаг \( h \) нельзя устремить в бесконечность: как только точка \( x + nh \) упирается в правый конец отрезка \( 1 \), шагать дальше некуда. Уитни обошел это препятствие дискретным перебором по мелкой сетке, но платой стал гигантский рост константы с ростом \( n \). Даже для кубических разностей (\( n=3 \)) точное значение не поддавалось, и сам Уитни написал фразу, ставшую вызовом на десятилетия: «Задача нахождения констант \( W_n \) и \( W'_n \), по-видимому, чрезвычайно трудна...».

Здесь \( W'_n \) — константа для самого естественного способа построения многочлена: интерполяции Лагранжа \( L_{n-1}f \), когда многочлен просто проводят через \( n \) равноотстоящих точек отрезка \( x_m = m/(n-1) \). Наряду с ней важнейшую роль играет константа \( \widetilde{W}_n \) для проектора по средним \( Q_{n-1}f \), когда многочлен подбирают так, чтобы его средние значения (интегралы) на \( n \) равных ячейках отрезка совпадали со средними значениями самой функции. Очевидно, что наилучшее приближение всегда не хуже любого конкретного способа, поэтому \( W_n \le \widetilde{W}_n \) и \( W_n \le W'_n \).

1980-е годы: Гипотезы Благовеста Сендова, операторы Иванова — Такева и прорыв к тройке

Четверть века спустя болгарский математик академик Благовест Сендов выдвинул смелую идею: константы Уитни вообще не растут с ростом степени \( n \), даже когда \( n \) стремится к бесконечности! Более того, Сендов сформулировал две знаменитые гипотезы:

Первая гипотеза Сендова: \( W_n \le 1 \) для всех \( n \).
Вторая гипотеза Сендова: \( W'_n \le 2 \) для всех \( n \).

Доказать ограниченность \( W_n \) пытались многие. Первые оценки давали рост порядка \( n 2^n \) (Ю. А. Брудный) и \( (n+1)n^n \). Решающий сдвиг произошел в 1985 году в Софии. К. Иванов и М. Такев придумали, как перенести непрерывное усреднение Уитни — Берлинга на конечный отрезок: нужно усреднять разность \( \Delta_y^n f(x - iy) \) по шагу \( y \), выбирая индекс \( i \) в зависимости от того, в какой по счету ячейке отрезка находится точка \( x \). Опираясь на эти операторы, Сендов доказал, что \( W_n \le 6 \).

Вскоре выяснилось, что конструкцию Сендова можно радикально улучшить. Если перейти к проектору по средним и вычесть многочлен с теми же интегралами на ячейках \( [0, 1/n], [1/n, 2/n], \dots \), то первообразная остатка \( \Phi_n(y) = \int_0^y g_n(s) ds \) обращается в нуль во всех узлах сетки \( j/n \). Анализ леммы Сендова привел Ю. В. Крякина в 1989 году к модифицированному тождеству (независимо оценку объявили также Ю. А. Брудный и Б. Сендов), которое сразу снизило общую границу до \( W_n \le 3 \).

Однако на этом пути обнаружилась главная геометрическая ловушка задачи. Во всех внутренних ячейках отрезка, от \( 1/n \) до \( 1 - 1/n \), модифицированная лемма Сендова легко давала идеальную оценку \( |g_n(x)| \le \omega_n(f) \) с константой 1. Вся трудность сосредоточилась в узких крайних ячейках у концов отрезка — на интервалах \( [0, 1/n) \) и \( (1-1/n, 1] \).

1990-е — 2002 год: Тождество с первообразной и барьер Гилевича — Крякина — Шевчука

Почему именно крайняя ячейка \( [0, 1/n) \) так упорно сопротивляется? Когда точка \( x \) очень близка к левому концу \( 0 \), мы не можем шагнуть влево — там граница отрезка. Мы вынуждены пускать лучи конечной разности далеко вправо, с большим шагом \( h = (1-x)/n \).

В работах Ю. В. Крякина 1990-х годов это длинное усреднение \( A_0 g_n(x) = \int_0^1 \Delta_{rh}^n g_n(x) dr \) было точно раскрыто через первообразную \( \Phi_n \):

\[ A_0 g_n(x) = g_n(x) + \frac{H_n}{h} \Phi_n(x) + \frac{1}{h} \sum_{j=1}^n \frac{(-1)^j}{j} \binom{n}{j} \Phi_n(x + jh), \]

где \( H_n = 1 + 1/2 + 1/3 + \dots + 1/n \approx \ln n \) — гармоническое число. В то же время для самой первообразной \( \Phi_n(y) \), благодаря нулям в узлах \( j/n \), работает красивое тождество Жука — Натансона, дающее точную огибающую:

\[ |\Phi_n(y)| \le |q_n(y)| \omega_n(f), \quad q_n(y) = y \prod_{j=1}^n \left(1 - \frac{ny}{j}\right). \]

Беда в том, что если в знаковой сумме из раскрытия \( A_0 g_n(x) \) заменить каждое слагаемое \( \Phi_n(x+jh) \) на его верхнюю оценку \( |q_n(x+jh)| \) по модулю (без учета знаков), то при малых \( x \) набегает лишняя единица, а вместе с дополнительным слагаемым от леммы И. А. Шевчука при больших \( n \) возникает добавка \( e^{-2} \). В итоге в фундаментальной работе Я. Гилевича, Ю. В. Крякина и И. А. Шевчука (2002), объединившей этот метод с идеей промежуточного приближения М. Такева, были установлены рекорды, продержавшиеся почти четверть века:

\[ \widetilde{W}_n \le 2 + e^{-2} \approx 2{,}135, \quad W'_n \le 3. \]

Чтобы пробить стену двойки для средних и тройки для Лагранжа при больших \( n \), нужно было перестать оценивать знаковую сумму по модулю в каждой точке по отдельности.

Осень 2026 года: Пограничный слой, второе усреднение (Argon / GPT) и ключевая нормировка Carbon

В октябре — ноябре 2026 года в математическом диалоге Ю. В. Крякина с моделями искусственного интеллекта (линия Argon — GPT 6.1 — Carbon) произошел каскад открытий, который полностью прояснил устройство задачи на краю отрезка при \( n \to \infty \).

Шаг первый: открытие масштаба \( \Lambda_n = n H_n \).
Оказалось, что в первой ячейке ширины \( 1/n \) самое опасное поведение функции разворачивается в еще более узком пограничном слое у самого нуля, на расстоянии порядка \( x = t / (n H_n) \), где \( t > 0 \) — безразмерная координата. В этом масштабе при \( n \to \infty \) множитель при первообразной \( \frac{H_n}{h} \Phi_n(x) \) превращается в предельную функцию \( U(t) = \lim \Lambda_n \Phi_n(t/\Lambda_n) \), а огибающая Жука — Натансона принимает универсальный вид:

\[ |U(t)| \le t e^{-t}. \]

Шаг второй: второе положительное усреднение.
Как поймать знаковую сумму \( T_n(x) = h^{-1} \sum_{j=1}^n \frac{(-1)^j}{j} \binom{n}{j} \Phi_n(x+jh) \), не теряя сокращения знаков? Было введено семейство взвешенных усреднений

\[ A_k g_n(x) = \int_0^1 w_k(r) \Delta_{rh}^n g_n(x) dr, \quad w_k(r) = \frac{k}{k-1}(1 - r^{k-1}), \quad k \ge 2. \]

Поскольку вес \( w_k(r) \ge 0 \) и его среднее равно 1, оператор \( A_k \) тоже не превосходит единицы по норме модуля непрерывности. Но при повторном интегрировании по частям знаменатели \( j \) в знаковой сумме превращаются в степени \( j^k \), и вклад всех дальних ячеек в пределе при фиксированном \( k \) и \( n \to \infty \) сворачивается в аккуратный интеграл Лапласа — гамма-свертку от одной и той же предельной первообразной \( U \) в первой ячейке:

\[ A_k g_n(t/\Lambda_n) = g_n(t/\Lambda_n) + \frac{k}{k-1} U_n(t) - I_k(t) + o(1), \quad I_k(t) = \frac{1}{(k-1)!} \int_0^\infty w^{k-2} U(t+w) dw. \]

Оценка интеграла \( I_k(t) \) через огибающую \( w e^{-w} \) сначала оставляла неизвестный экспоненциальный хвост \( L e^{-t} \), что позволило спустить оценку для средних с \( 2 + e^{-2} \) сначала до \( (3 + \ln 2)/2 \approx 1{,}846 \), а затем до \( 2 - \frac{1}{2}\ln 2 \approx 1{,}653 \).

Шаг третий: идея Carbon — нормировка в конце отрезка и обнуление знаковой суммы!
Решающее прояснение внесла идея Carbon. Сначала обсуждалась попытка алгебраически сократить общие значения первообразной в соседних точках за счет гладкого выбора шага, но это оставляло слишком большой знакопеременный хвост третьего порядка. Тогда внимание переключили на прямое согласование двух операторов — обычного усреднения \( A_0 \) и взвешенного \( A_k \) — непосредственно в конце отрезка \( x = 0 \)!

Посмотрите, как просто и красиво работает этот аргумент:
1. Вес \( w_k(r) = \frac{k}{k-1}(1 - r^{k-1}) \) при больших \( k \) почти на всем интервале \( [0, 1] \) очень близок к единице: интеграл от модуля разности \( |w_k(r) - 1| \) не превосходит \( 2/k \). Поэтому для любой функции разница между операторами мала:

\[ |A_k g_n(x) - A_0 g_n(x)| \le \frac{2}{k}. \]

2. В самой крайней точке \( x = 0 \) шаг равен \( h = 1/n \), поэтому все узлы \( jh = j/n \) в формуле для \( A_0 g_n(0) \) попадают точно в концы ячеек, где первообразная равна нулю: \( \Phi_n(j/n) = 0 \)! Значит, в точке \( x = 0 \) выполняется точное равенство без всяких остатков:

\[ A_0 g_n(0) = g_n(0). \]

3. Сравнивая это с формулой для \( A_k \) при \( t = 0 \), получаем, что интеграл \( I_k(0) \) обязан стремиться к нулю при \( k \to \infty \), а вместе с ним по оценке гамма-плотности равномерно стремится к нулю и весь остаток \( I_k(t) \to 0 \)!

Следствие оказалось ошеломляющим: в масштабе пограничного слоя \( x = t / (n H_n) \) вся сложная знаковая сумма \( T_n(x) \) из длинного тождества асимптотически исчезает (стремится к нулю со скоростью \( O(H_n^{-2/5}) \))! Вся формула длинного усреднения схлопывается в элементарное соотношение:

\[ A_0 g_n(t/\Lambda_n) = g_n(t/\Lambda_n) + U_n(t) + o(1). \]

Итоговые формулы: вплотную к гипотезам Сендова

Как только знаковая сумма исчезла, асимптотические оценки посыпались как спелые яблоки.

1. Для интерполяции по средним \( \widetilde{W}_n \) и наилучшего приближения \( W_n \):
Поскольку \( A_0 g_n(x_n) \le 1 \), в точке максимума \( t = c \) получаем

\[ g_n(x_n) \le 1 - U(c) + o(1) \le 1 + c e^{-c} + o(1). \]

Максимум функции \( c e^{-c} \) достигается при \( c = 1 \) и равен \( e^{-1} = 1/e \approx 0{,}367879 \). Отсюда немедленно следует новая асимптотическая граница:

\[ \limsup_{n \to \infty} W_n \le \limsup_{n \to \infty} \widetilde{W}_n \le 1 + e^{-1} \approx 1{,}367879. \]

Более того, количественная оценка интегралов Лапласа дает эту границу для каждого \( n \) с явным остатком \( \widetilde{W}_n \le 1 + e^{-1} + O(H_n^{-2/5}) \).

2. Для интерполяции Лагранжа \( W'_n \) (Вторая гипотеза Сендова в пределе!):
В задаче Лагранжа ошибка в точке \( x_n = c/\Lambda_n \) складывается из самого значения \( g_n(x_n) \), вклада крайнего узла \( g_n(0) e^{-c} \) и суммы по остальным узлам, не превышающей \( c e^{-c} \). Сравнение двух положительных усреднений \( A_0 \) и \( A_2 \) в точках \( 0 \) и \( x_n \) связывает дефициты обоих концов через интеграл \( \int_0^c U(t) dt \), а сравнение разностей вдоль двух близких аффинных сеток дает дифференциальное ограничение \( |(e^t U(t))'| \le e^t \). После подстановки в формулу ошибки все неизвестные члены сокращаются, оставляя при каждом \( c > 0 \) оценку

\[ E \le 1 + \left(1 + c + \frac{c^2}{2}\right) e^{-c} < 2. \]

Это означает, что в пределе при \( n \to \infty \) Вторая гипотеза Сендова верна асимптотически:

\[ \limsup_{n \to \infty} W'_n \le 2, \quad W'_n \le 2 + O(H_n^{-2/5}), \]

причем приближение к двойке возможно лишь в сверхтонком подслое \( x_n = o((n \ln n)^{-1}) \), где \( g_n(0) \to -1 \), \( g_n(x_n) \to 1 \), а предельная первообразная \( U(t) \) тождественно равна нулю!

Так семидесятилетняя цепочка идей — от бесконечного интеграла Уитни — Берлинга и локальных усреднений Иванова — Такева — Сендова через граничные тождества Крякина — Шевчука к анализу пограничного слоя и краевой нормировке в тандеме математика и ИИ — позволила укротить хаос знакопеременных сумм и вплотную приблизиться к полному решению проблемы констант Уитни.