2. Аппроксимация функций Интерполяционный многочлен Лагранжа

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

Задача приближения функций называется задачей аппроксимации, функция – аппроксимирующей. Аппроксимация, при которой приближение строится на заданном дискретном множестве точек , называется точечной.

Процесс вычисления значений функции в точках , отличных от значений , называется интерполяцией. Если аргумент находится за пределами отрезка , то задача определения значения функции в точке называется экстраполированием.

Пусть задана система точек , всего значений некоторой функции . Требуется построить функцию такую, что ее значения в точках , называемых узлами интерполяции, совпадали с заданными значениями функции:

, , т.е.

требуется строгое совпадение значений функций в узлах интерполяции. Геометрически это означает, что график интерполирующей функции проходит через заданные точки.

Будем искать интерполирующую функцию в виде канонического многочлена степени

. (2.1)

Полагая, что , получим:

. (2.2)

Таким образом, чтобы многочлен (2.1) был интерполяционным, нужно, чтобы его коэффициенты удовлетворяли системе уравнений:

(2.3)

Получили систему линейных уравнений порядка (n + 1), её определитель

так называемый определитель Вандермонда, не равен нулю, если все различны, следовательно, существует единственный интерполирующий многочлен канонического вида (2.1). Однако практическое построение интерполяционного многочлена таким путем малоэффективно.

Построим интерполирующий полином Лагранжа в виде линейной комбинации:

,

где – базисные полиномы такие, что:

.

Значение базисного полинома в узлах интерполяции будет:

тогда интерполяционный многочлен Лагранжа имеет вид:

.

При по двум точкам построим линейный многочлен:

.

При по трем точкам построим многочлен Лагранжа второй степени:

Погрешность интерполяции полинома Лагранжа

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

. (2.3)

В точках , будет ,следовательно, , значит функция имеет по меньшей мере (n + 1) нулей на рассматриваемом отрезке.

Пусть , не совпадает ни с одной точкой. Подберем значение так, чтобы функция в этой точке имела еще один корень:

.

Следовательно, можно подобрать единственное значение коэффициента так, чтобы , , тогда имеет по крайней мере (n + 2) корней. По теореме Ролля между нулями функции существует, по крайней мере, один нуль производной, тогда производная имеет (n + 1) нуль. Применим теорему Ролля к производной, тогда имеет нулей и т.д. Производная имеет один корень, следовательно, существует точка такая, что .

Продифференцируем равенство (2.3) (n + 1) раз:

Следовательно, учитывая, что , получим . Поскольку – любая точка из промежутка, не совпадающая с узловой, то заменим ее произвольной точкой , тогда:

. .

Если производная порядка ограничена , то оценка погрешности полиномиальной интерполяции:

.

Рассмотрим линейную интерполяцию

,

если .

Оценим погрешность линейной интерполяции, для этого найдем :

,

Преобразуем .

Найдем точку, в которой производная этого выражения имеет нуль:

,

тогда и оценка погрешности имеет вид:

.