О многочленах Берштейна для любителей математики (+ немного формул)

 МНОГОЧЛЕНЫ БЕРНШТЕЙНА В ГЕОМЕТРИИ ГАУССА–ЧЕБЫШЕВА
Популярный обзор: идеи, геометрия и ключевые соотношения

8 ноября 2026 года
Подготовлено в математическом обсуждении с Ю. В. Крякиным.

СТРАНИЦА 1. ПАРАДОКС ПЛОСКОГО ОТРЕЗКА И МОНЕТЫ БЕРНШТЕЙНА

В тысяча девятьсот двенадцатом году Сергей Натанович Бернштейн предложил удивительно наглядный способ приближать любую непрерывную функцию \(f\in C[0,1]\) на единичном отрезке алгебраическими многочленами. Его идея выросла из схемы независимых испытаний Бернулли. Представим себе несимметричную монету, у которой вероятность выпадения герба равна координате точки \(x\in[0,1]\). Если подбросить такую монету \(n\) раз, то вероятность получить ровно \(k\) гербов задаётся биномиальным весом. Бернштейн разбил отрезок на \(n\) равных частей и усреднил значения функции в равномерных узлах \(k/n\):
\[
B_n(f;x)=\sum_{k=0}^n f\left(\frac{k}{n}\right)\binom{n}{k}x^k(1-x)^{n-k}.
\] (1)

Казалось бы, раз узлы делят отрезок совершенно поровну, то и приближать функцию такой многочлен должен равномерно во всех точках. Однако в действительности всё происходит иначе. В самой середине отрезка случайный разброс числа гербов максимален, поэтому биномиальное усреднение захватывает широкую окрестность точки. Но по мере приближения к концам \(0\) и \(1\) монета становится почти детерминированной. Дисперсия распределения и классический предел Вороновской устроены неоднородно:
\[
B_n\bigl((t-x)^2;x\bigr)=\frac{x(1-x)}{n},\qquad
\lim_{n\to\infty}n\bigl(B_n(g;x)-g(x)\bigr)=\frac{1}{2}x(1-x)g''(x).
\] (2)
Внутренний шаг усреднения пропорционален \(\sqrt{x(1-x)/n}\), а на самих концах погрешность обращается в нуль.

Из-за этого вырождения обычный модуль непрерывности с постоянным шагом вдоль отрезка оказывается слишком грубым. Чтобы согласовать меру гладкости с весом \(x(1-x)\), в работах З. Дицана и В. Тотика шаг разности вдоль отрезка стали искусственно умножать на \(\sqrt{x(1-x)}\). На этом пути исследователей ждала серьёзная преграда: прямолинейный шаг вблизи границы неизбежно упирается в конец интервала и выводит точку за пределы \([0,1]\). Это заставляло делить отрезок на зоны и громоздить тяжёлые технические конструкции там, где задача по своей природе должна быть простой и прозрачной.

\newpage

СТРАНИЦА 2. ПОДЪЁМ НА ОКРУЖНОСТЬ: ГЕОМЕТРИЯ ГАУССА–ЧЕБЫШЕВА

Причина описанных трудностей носит чисто геометрический характер. Мы пытаемся описать поведение многочлена Бернштейна, двигаясь вдоль плоского диаметра, тогда как сам этот многочлен живёт в геометрии полуокружности, построенной над отрезком \([0,1]\).

Свяжем каждую точку отрезка с углом \(\theta\in[0,\pi]\) на дуге подстановкой \(x=\sin^2(\theta/2)\). Если поворачивать точку по дуге на один и тот же малый угол \(u\), то её проекция на диаметр смещается на величину \(\Delta x\approx u\sqrt{x(1-x)}\) — в точности по закону весового шага Бернштейна! Главное преимущество окружности в том, что у неё нет краёв: при повороте через конец диаметра точка плавно отражается обратно внутрь отрезка, никогда не покидая его пределов.

Поднимем функцию на дугу равенством \(F(\theta)=f(\sin^2(\theta/2))\). Симметричный поворот в обе стороны задаёт оператор чебышевского сдвига \(\mathcal{T}_u\) и второй модуль непрерывности Гаусса–Чебышева на единичном отрезке:
\[
(\mathcal{T}_u f)\left(\sin^2\frac{\theta}{2}\right)=\frac{F(\theta+u)+F(\theta-u)}{2},\qquad
\omega^*_2(f,t)=\sup_{0<u\le t}\|(I-\mathcal{T}_u)f\|_{C[0,1]}.
\] (3)

Рядом с этим разностным модулем действует функционал Петре \(K_1(f,s)\), измеряющий приближение функции гладкими функциями с контролем второй производной по углу \(\theta\). Поскольку периодическое сглаживание Стеклова на окружности сохраняет чётность и не требует никаких граничных сшивок, прямая оценка приближения многочленами Бернштейна доказывается в несколько строк сразу для всех \(n\ge 1\) с явной малой константой:
\[
\|B_nf-f\|_{C[0,1]}\le 3\,\omega^*_2\left(f,\sqrt{\frac{2}{n}}\right).
\] (4)

\newpage

СТРАНИЦА 3. КРАЕВОЙ НАКЛОН, СЖАТИЕ И ВЕРОЯТНОСТНЫЙ СМЫСЛ

Можно ли обратить оценку (4) снизу непосредственно для классического многочлена Бернштейна? Здесь обнаруживается тонкое различие на самой границе. Вторая производная по углу \(-D_\theta^2\) при пересчёте к координате \(x\) порождает родной дифференциальный оператор Гаусса–Чебышева:
\[
\mathcal{A}g(x)=-x(1-x)g''(x)+\left(x-\frac{1}{2}\right)g'(x).
\] (5)
В концах \(0\) и \(1\) старший член \(x(1-x)g''(x)\) исчезает, но остаётся первая производная \(\mp\frac{1}{2}g'\). Угловой поворот смещает конец отрезка внутрь на величину \(\sin^2(u/2)\approx u^2/4\), поэтому модуль \(\omega^*_2\) чувствует краевой наклон графика. Классический же оператор \(B_n\) прибит к концам намертво: его предел Вороновской (2) не содержит члена с \(g'(x)\) и слеп к наклону у границы.

Чтобы освободить многочлен Бернштейна от зажатия на концах и вернуть ему угловую геометрию, достаточно слегка сжать его аргумент к центру отрезка \(1/2\) на множитель \(\cos(1/\sqrt{n})\approx 1-\frac{1}{2n}\):
\[
B_n^*(f;x)=B_n(f;x^*),\qquad
x^*=\left(1-\frac{1}{2n}\right)\left(x-\frac{1}{2}\right)+\frac{1}{2}.
\] (6)
Сумма базисных весов при этом остаётся тождественно равной единице: \((x^*+(1-x^*))^n=1\), но предел Вороновской для \(B_n^*\) в точности совпадает с половиной родного оператора \(\mathcal{A}\). В результате оценка (4) становится полностью обратимой в строгой форме Тотика:
\[
\sup_{k\ge n}\|B_k^*f-f\|_{C[0,1]}\asymp\omega^*_2\left(f,\frac{1}{\sqrt{n}}\right).
\] (7)

У сжатия \(x^*\) есть глубокий вероятностный смысл. В статистике оно отвечает байесовской регуляризации Лапласа, которая уводит оценку вероятности от наивных крайних значений \(0\) и \(1\) к априорному центру \(1/2\). А в физике броуновского движения по окружности проекция случайного шага \(\cos(\theta\pm u)\) из-за кривизны дуги в среднем умножается на \(\cos u\approx 1-u^2/2\), создавая естественный центростремительный дрейф Ито к экватору.

\newpage

СТРАНИЦА 4. ОТ ПРОГРЕССИЙ ФИЛЛИПСА К ЧЕБЫШЕВСКОЙ СЕТКЕ НА ДУГЕ

В 1996–1997 годах Джордж Филлипс в Сент-Эндрюсе построил обобщённые многочлены Бернштейна по геометрической прогрессии, заменив целые числа \(k\) на гауссовы \(q\)-числа \([k]_q=(1-q^k)/(1-q)\). При вещественном \(0<q<1\) узлы Филлипса \([k]_q/[n]_q\) сгущаются только к правому концу отрезка, разрушая симметрию задачи.

Перенесём конструкцию Филлипса на окружность, положив параметр деформации равным повороту \(q=e^{i\pi/n}\). Тогда гауссовы числа превращаются в ядра Дирихле, а квадраты модулей отношений \([k]_q/[n]_q\) задают на отрезке \([0,1]\) симметричную сетку узлов Гаусса–Чебышева:
\[
t_{n,k}=\left|\frac{[k]_q}{[n]_q}\right|^2=\sin^2\frac{k\pi}{2n},\qquad k=0,1,\dots,n.
\] (8)
Эта сетка есть проекция равномерного деления дуги \(\theta_k=k\pi/n\). Она квадратично сгущает точки измерения функции одновременно к обоим концам \(0\) и \(1\) по закону \(\Delta t_{n,k}\approx\frac{\pi}{n}\sqrt{t_{n,k}(1-t_{n,k})}\).

Подстановка чебышевских узлов (8) и угловой доли дуги \(u(x)=\frac{2}{\pi}\arcsin\sqrt{x}\) рождает дуговой оператор Бернштейна–Гаусса–Чебышева:
\[
BC_n(f;x)=\sum_{k=0}^n f\left(\sin^2\frac{k\pi}{2n}\right)\binom{n}{k}x^k(1-x)^{n-k}.
\] (9)
На функциях с корневыми особенностями у края переход к узлам \(t_{n,k}\) даёт колоссальный выигрыш. В узлах (8) сингулярная функция \(f(x)=\sqrt{x}\) превращается в гладкую синусоиду \(\sin\frac{k\pi}{2n}\), и скорость её приближения на дуге возрастает с медленного порядка \(O(1/\sqrt{n})\) до быстрого порядка \(O(1/n)\).

\newpage

СТРАНИЦА 5. КУЗОВА «РЕНО», ДВОЯКАЯ ВЫПУКЛОСТЬ И МНОГОМЕРНЫЙ МИР

В конце 1950-х годов Поль де Кастельжо на «Ситроене» и Пьер Безье на «Рено» положили базис Бернштейна в основу компьютерного проектирования автомобильных кузовов. Им требовались два свойства: вычисление кривой устойчивыми выпуклыми усреднениями (алгоритм де Кастельжо) и строгое сохранение выпуклости каркаса, чтобы блики света на металле не давали ложных волн.

И классический оператор \(B_n\), и сжатый оператор \(B_n^*\), и чебышевский оператор \(BC_n\) полностью сохраняют монотонность функции и не увеличивают числа перемен знака. Характер же сохраняемой выпуклости определяется второй производной:
\[
(B_nf)''(x)=n(n-1)\sum_{k=0}^{n-2}\Delta^2f\left(\frac{k}{n}\right)p_{n-2,k}(x),\qquad
(BC_nf)''(x)=n(n-1)\sum_{k=0}^{n-2}\Delta^2f(t_{n,k})p_{n-2,k}(x).
\] (10)
У операторов \(B_n\) и \(B_n^*\) вторые разности берутся по равномерным узлам диаметра, поэтому они сохраняют обычную плоскую выпуклость \(f''(x)\ge 0\), нужную на пологих панелях дверей и крыши. 

У оператора \(BC_n\) узлы \(t_{n,k}\) равномерны по дуге, поэтому он сохраняет дуговую выпуклость Гаусса–Чебышева \(F''(\theta)=-\mathcal{A}f(x)\ge 0\). Именно такая выпуклость присуща реальным обводам капотов, крыльев и колёсных арок, которые подходят к габаритным краям по вертикальной касательной. Например, для полуокружности (эллиптического кузова) плоская вторая производная взрывается на краях, тогда как дуговой оператор Чебышева идеален и регулярен:
\[
f(x)=\sqrt{x(1-x)}=\frac{1}{2}\sin\theta\quad\Longrightarrow\quad
F''(\theta)=-(\mathcal{A}f)(x)=-\sqrt{x(1-x)}\le 0.
\] (11)

Наконец, вся эта одномерная геометрия непосредственно переносится в многомерный мир: над квадратной панелью кузова \([0,1]^2\) поднимается двумерный тор, а над треугольным симплексом и шаром возвышается многомерная сфера, где роль нашего оператора \(\mathcal{A}\) и сжатия \(x^*\) берут на себя сферический оператор Лапласа–Бельтрами и сжатие к барицентру треугольника.