Численные методы математического анализа
1. Задача интерполяции. Многочлен Лагранжа
1.1. Постановка задачи
Пусть задана функция
. Часто нахождение
значений этой функции может оказаться трудоемкой задачей. Например,
- параметр в
некоторой сложной задаче, после решения которой определяется значение f(x), или
f(x) измеряется в дорогостоящем эксперименте. В этих случаях можно получить
небольшую таблицу значений функции, но прямое нахождение ее значений при большом
количестве значений аргумента нереально. В такой ситуации f(x) заменяется
приближенной формулой
, которая в
определенном смысле близка к функции f(x). Близость обеспечивается введением в
функцию
свободных параметров и их
соответствующим выбором.
Итак, известны значения функции f(x) в точках ![]()
,
. Потребуем, чтобы для некоторой функции ![]()
выполнялись равенства:
(1)
Если (1) рассматривать как систему для определения
, то этот
способ называется интерполяцией (Лагранжевой).
Если
зависит от
нелинейно, то
интерполяция нелинейная, иначе интерполяция линейная. В случае линейной
интерполяции можно записать
(2)
(где
- система
линейно-независимых функций)
Подставим (2) в (1). Относительно
получаем линейную
систему уравнений:
,
(3)
Для однозначной разрешимости системы должно быть
.
Для того, чтобы задача нтерполирования имела единственное
решение, система функций ![]()
должна для любых
несовпадающих
удовлетворять условию:
(4)
Система функций, удовлетворяющая условию (4), называется чебышевской.
Наиболее простой и (для многих случаев) удобной является
система функций
,
. Функция
при этом представляет
собой многочлен степени
(интерполяционный
многочлен) с коэффициентами
.
Система уравнений (3) в этом случае имеет вид:
,
(5)
Определителем этой системы является отличный от нуля
определитель Вандермонда: 
Отсюда следует, что интерполяционный многочлен существует и единственен.
Непосредственное решение системы (5) для нахождения
aj уже при небольших n приводит к сильному
искажению значений aj. Получим явный вид
нтерполяционного многочлена
, не решая систему
(5).
Если
многочлен степени n
такой что
, то
- искомый интерполяционный многочлен степени
, т.к.
.
Так как
при
, то
делится на
для любых
, то
есть
. Так как
, то
.
Таким образом,
(6)
Такая форма записи интерполяционного многочлена называется
многочленом Лагранжа и обозначается, как правило,
.
Существуют и другие формы записи того же самого интерполяционного многочлена.
Если обозначить
, то
.
Тогда (6) можно записать в виде 
В узлах
многочлен Лагранжа
совпадает с заданной функцией, в остальных точках в общем случае
не совпадает с
(кроме
случая, когда
многочлен степени не выше
). Разность
-
- остаточный член. Запишем ее в виде
.
При
. Найдем постоянную
такую, чтобы
в некоторой фиксированной точке
, в которой мы
рассматриваем погрешность.
Относительно
будем предполагать
кратную
дифференцируемость.
Значение c, при котором
существует и равно
. Тогда
функция
равна нулю по крайней мере в
точках
.
По теореме Ролля производная
равна нулю по крайней
мере в
точках
. Далее,
равна нулю по крайней мере в n точках и т.д.
Для
- производной получаем,
что существует по крайней мере одна точка
такая, что
.
![]()
Отсюда получаем при
. В этом случае в
точке ![]()
, т.е. остаточный член в точке
имеет вид:
![]()
Значение
зависит от
- точки, в
которой рассматривается погрешность.
,
где ![]()