1.4. Задача интерполяции как простейшая задача приближения функций

Пусть дана сетка узлов , где , и известны значения функции в узлах сетки: , причем (N+1) узловых точек – попарно различны.

Интерполяция обобщенными полиномами.

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

- «обобщенный полином N-го порядка», для которого выполняется условие совпадения значений полинома в узлах сетки со значениями . Это требование приводит к системе линейных уравнений:

.

(4)

Систему (4) называют Нормальной системой уравнений.

Для разрешимости системы (4) необходимо, чтобы определитель системы

.

Рассмотрим частный случай задачи интерполяции – Интерполяцию алгебраическими полиномами.

В этом случае в качестве базисной системы функции выступают степенные функции:

.

Обозначим - искомый интерполяционный полином N-ой степени

,

И запишем условия совпадения значений полинома с табличными значениями функции :

.

(5)

Определитель этой системы - определитель Вандермонда – отличен от нуля:

,

Поэтому система имеет единственное решение.

Система (5) Плохо обусловлена.(см. п. п 3.4.3.). Покажем, как можно построить искомый интерполяционный полином другим способом, не решая систему. Для этого построим так называемые Фундаментальные полиномы степени N: , удовлетворяющие условию:

.

(6)

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

Условие (6) непосредственно проверяется. Отсюда следует, что искомый интерполяционный полином можно записать в виде:

.

(7)

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

Приведем еще одно доказательство единственности полинома Лагранжа (7) независимо от логики решения нормальной системы (2).

Теорема 1.1. Полином Лагранжа (4), проходящий через все табличные (N+1) значения функции Y(X) – единственный.

От противного. Пусть Еще один полином степени N, решающий ту же задачу интерполяции. Рассмотрим разность - полином порядка . Очевидно, что этот полином имеет на отрезке ровно (N+1) корень, что противоречит основной теореме алгебры. Значит,

Теорема 1.2. (О погрешности интерполяции). Пусть функция , задана сетка узлов , – интерполяционный полином Лагранжа (4), построенный по значениям функции Y(X) в узлах сетки Xn. Тогда для погрешности интерполяции справедливы следующие оценки:

	(8)
– Теоретическая погрешность в точке ,
	(9)
– абсолютная погрешность в точке,
	(10)

– максимальная абсолютная погрешность на всем отрезке,

Где

–

– специальный полином (N+1)-ой степени, построенный по узлам как по нулям;

–

– существует в силу определения класса функций .

Запишем Y(X) в виде:

,

Где – погрешность интерполяции в точке . Очевидно по условию, что .

Отсюда следует, что погрешность (остаточный член) интерполяции можно искать в виде:

,

Где R– некоторая функция.

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

,

(11)

Где T – свободная переменная.

Положим в (11) T=X.

.

Т. е. функция обращается в 0 в точке T=X. Положим далее последовательно Получаем:

.

Т. о. мы получили, что на отрезке функция обращается в 0 ровно в (N+2) точках. Отсюда по теореме Ролля следует, что на интервале (A, b)

· существует, по крайней мере, (N+1) точка, в которой обращается в 0;

· существует, по крайней мере, N точек, в которых обращается в 0;

· …………………………

· существует, по крайней мере, 1 точка , в которой обращается в 0.

Продифференцируем формулу (11) по T n+1 раз и положим . Получим:

.

(12)

Учтем, что , т. к. степень полинома равна N. Далее получаем:

,

Т. к. - многочлен -Ой степени специального вида с коэффициентом при старшей степени, равным 1.

Подставляя эти результаты в (12), получаем:

,

Откуда следует

.

Подставляя в выражение для , получаем

,

Откуда следует формула (8). Оценки (9), (10) вытекают из (8) автоматически.

Пример 8. Пусть , . Построить интерполяционный полином Лагранжа L2(X) по узлам и оценить погрешность интерполяции в точке X=116 и на всем отрезке [A,B].

Построение полинома Лагранжа (соответствует квадратичной интерполяции) – отложим на семинарское занятие. Оценим погрешность интерполяции. Очевидно, что заданная функция дифференцируема на любое число раз. Нам достаточно оценить 3-ю производную.

Последовательно дифференцируя , получаем:

.

Подставляя в формулу (9), получим:

.

Оценим погрешность в точке X=116:

, откуда следует

.

Оценку на всем отрезке предлагается вычислить на семинаре.

Приведем еще одну форму записи интерполяционного полинома:

((10

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

((11 (13)

Система (13) легко решается, начиная с . Построенный в такой форме интерполяционный полином называется интерполяционным полиномом Ньютона. Далее будут рассмотрены частные виды полинома Ньютона для сетки равноотстоящих узлов.

< Предыдущая		Следующая >