Математика / Интерполяция

Константа Лебега: усиление ошибки интерполяции

Константа Лебега измеряет наибольшее усиление ошибок значений при полиномиальной интерполяции на фиксированной сетке. Она зависит от расположения узлов и выбранного отрезка, но не от конкретных значений функции.

Опубликовано: 31 октября 2026 г.Обновлено: 31 октября 2026 г.

Формула

L_n(x)=\sum_{j=0}^n|\ell_j(x)|,\qquad\Lambda_n=\max_{x\in[a,b]}L_n(x)

Схема Функция Лебега трёх узлов

Для узлов −1,0,1 функция равна 1+|x|−x². Максимумы 1,25 достигаются при x=±0,5.

Обозначения

$ℓ_j(x)$: базисный полином Лагранжа, равный единице в своём узле и нулю в остальных, безразмерный
$L_n(x)$: локальный коэффициент усиления абсолютной погрешности данных, безразмерный
$Λ_n$: максимум функции Лебега на отрезке, безразмерный
$n$: максимальная степень интерполянта; узлов n+1, безразмерный

Условия применения

Узлы x₀,…,xₙ должны быть попарно различны и принадлежать фиксированному отрезку [a;b]. Рассматривается интерполяция полиномом степени не выше n с обычным базисом Лагранжа. Граница относится к максимальной абсолютной ошибке данных, одинаково ограниченной во всех узлах.

Ограничения

Константа характеризует чувствительность математической задачи к изменению узловых значений и не является ошибкой конкретного программного алгоритма. Малое Λₙ не гарантирует точного приближения произвольной негладкой функции. Перебор по конечной сетке даёт лишь нижнюю оценку максимума, если нет дополнительного доказательства.

Подробное объяснение

Полином Лагранжа имеет вид p(x)=Σfⱼℓⱼ(x). Если данные изменились на δfⱼ, то δp(x)=Σδfⱼℓⱼ(x). Из неравенства треугольника при |δfⱼ|≤ε сразу получается |δp(x)|≤εΣ|ℓⱼ(x)|. Сумма модулей и есть функция Лебега: она показывает возможное локальное усиление абсолютной ошибки.

Оценка в фиксированной точке достижима: знаки ошибок данных можно выбрать равными знакам базисных полиномов в этой точке. Поэтому Λₙ — норма оператора интерполяции из узловых значений в непрерывный полином с равномерной нормой. В каждом узле Lₙ=1, но между узлами сумма модулей может быть больше единицы.

Для базиса выполняется Σℓⱼ(x)=1. Без модулей сумма всегда единица и ничего не говорит об усилении; большие положительные и отрицательные слагаемые могут взаимно сокращаться. Именно наличие отрицательных базисных значений вне своего узла делает Lₙ(x)>1. Смена единиц координаты с одновременным линейным переносом всего отрезка не меняет значения константы.

Нужно различать две ошибки: изменение полинома от неточных данных и различие между точной функцией и её интерполянтом. Для непрерывной функции вторая ограничена (1+Λₙ)Eₙ, где Eₙ — ошибка лучшего равномерного приближения полиномами степени n. Такая связь объясняет, почему удачная сетка полезна, но не заменяет свойства самой функции и отдельный анализ округления.

Как пользоваться формулой

Зафиксируйте отрезок и все различные узлы; не сравнивайте константы, вычисленные на разных областях аргумента.
Постройте каждый базис Лагранжа или используйте устойчивое вычисление базисных значений через барицентрические веса.
В каждой точке суммируйте модули всех базисных значений, а не модуль их общей суммы, которая равна единице.
Найдите максимум на отрезке; для небольшого примера разберите знаки базиса и проверьте внутренние экстремумы.
Умножьте найденное Λₙ на абсолютную границу ошибки данных и отдельно оцените ошибку приближения функции.

Историческая справка

Название константы связано с математической традицией изучения операторов приближения и равномерной сходимости. В современной теории интерполяции это норма линейного оператора, который переводит набор узловых значений в полином. Такая постановка отделяет свойства сетки от конкретной функции и позволяет сравнивать разные способы выбора точек без заранее заданных наблюдений. Приведённая статья не приписывает всю полиномиальную интерполяцию Анри Лебегу.

В 2010 году Ллойд Трефетен опубликовал в библиотеке примеров Chebfun материал Lebesgue functions and Lebesgue constants с вычислениями для равномерных и чебышёвских сеток. Там сравниваются рост констант и чувствительность интерполянтов. В книге Approximation Theory and Approximation Practice 2013 года эти идеи включены в алгоритмический подход к приближению функций. Современные программы визуализируют Lₙ(x), однако графический поиск максимума по-прежнему следует отличать от аналитически доказанного результата.

Историческая линия формулы

Имя Лебега закреплено за константой в теории приближений. Здесь используется современное определение через норму оператора Лагранжа; пример Chebfun Трефетена служит конкретным источником вычислительной интерпретации, а не свидетельством первоначального авторства определения.

Пример

Возьмём узлы −1, 0 и 1 на [−1;1], то есть n=2. Базисы: ℓ₀=x(x−1)/2, ℓ₁=1−x², ℓ₂=x(x+1)/2. На правой половине 0≤x≤1 первый неположителен, остальные неотрицательны. Поэтому L₂=−ℓ₀+ℓ₁+ℓ₂=1+x−x². Максимум достигается при производной 1−2x=0, то есть x=0,5, и равен 1,25. По симметрии слева такой же максимум при −0,5. Значит Λ₂=1,25. Если ошибка каждого узлового значения не больше 0,02 В, изменение интерполированного напряжения на всём отрезке не превышает 1,25·0,02=0,025 В. Это оценка влияния ошибок данных, а не отклонения от неизвестного физического сигнала.

Частая ошибка

Не вычисляйте |Σℓⱼ| вместо Σ|ℓⱼ|: первое выражение тождественно единице. Не выдавайте наибольшее значение на редкой пробной сетке за доказанный максимум. Не умножайте Λₙ на относительную погрешность без учёта масштаба данных: здесь вывод относится к абсолютной границе. Не называйте обусловленность сетки численной устойчивостью конкретной реализации.

Практика

Задачи с решением

Два концевых узла

Условие. Найдите константу Лебега линейной интерполяции в узлах 0 и 1 на [0;1].

Решение. Базисы равны 1−x и x. Оба неотрицательны на отрезке, поэтому L₁(x)=(1−x)+x=1 во всех точках.

Ответ. Λ₁=1.

Достижимость границы

Условие. Для узлов −1,0,1 в точке x=0,5 возьмите ошибки данных −0,02;0,02;0,02 В. Найдите ошибку полинома.

Решение. Базисные значения равны −0,125;0,75;0,375. Их скалярное произведение с ошибками даёт 0,0025+0,015+0,0075=0,025 В.

Ответ. δp(0,5)=0,025 В — граница достигается.

Дополнительные источники

https://www.chebfun.org/examples/approx/LebesgueConst.html
https://dlmf.nist.gov/3.3
https://people.math.sc.edu/burkardt/c_src/lebesgue/lebesgue.html

Связанные формулы

Математика

Узлы Чебышёва: выбор точек для интерполяции

$x_k=\frac{a+b}{2}+\frac{b-a}{2}\cos\frac{(2k+1)\pi}{2n},\quad k=0,\ldots,n-1$

Узлы Чебышёва — неравномерные точки, сгущённые около концов отрезка. Корни полинома Tₙ дают удобную сетку для интерполяционного полинома степени не выше n−1 и уменьшают наибольший модуль узлового произведения.

Математика

Интерполяция Лагранжа: полином через заданные точки

$P_n(x)=\sum_{k=0}^{n}y_kL_k(x),\qquad L_k(x)=\prod_{\substack{0\leq j\leq n\\j\ne k}}\frac{x-x_j}{x_k-x_j}$

Формула Лагранжа строит полином степени не выше n по n+1 различным узлам. Каждый базисный множитель выбирает значение своего узла и обращает вклады остальных узлов в ноль.

Математика

Барицентрическая интерполяция: веса и вычисление полинома

$P(x)=\frac{\sum_{k=0}^{n}\dfrac{w_ky_k}{x-x_k}}{\sum_{k=0}^{n}\dfrac{w_k}{x-x_k}},\qquad w_k=\frac{1}{\prod_{j\ne k}(x_k-x_j)},\quad x\ne x_k$

Барицентрическая форма вычисляет полином Лагранжа через заранее рассчитанные веса узлов. Она удобна для множества запросов по одной таблице; при точном совпадении аргумента с узлом нужно сразу вернуть узловое значение.