Интерполяционный полином Лагранжа. Наиболее общей формулой полиномиального интерполирования является интерполяционная формула Лагранжа

Наиболее общей формулой полиномиального интерполирования является интерполяционная формула Лагранжа. Задача параболиче­ского интерполирования в этом случае формулируется следующим об­разом: на отрезке [a, b] в узлах интерполяции х₀ x₁..., х_n задается функция f(x) cвоими n + 1 значениями y₀=f(x₀), y₁=f(x₁),...y_n=f(x_n) требуется построить многочлен L(x) так, чтобы в узлах интерполяции x_0, x_1, .... х_n его значения совпадали со значениями заданной функции, т. е.

l(x_o)=y_o, l(x₁)=y₁,..., l(x_n)=y_n.

Следует отметить, что в такой постановке задачи узлы интерполя­ции х₀, x₁..., х_n могут произвольно отстоять друг от друга на отрез­ке [а, b], иными словами, узлы интерполяции неравноотстоящие, т. е. h=x_i+1-x_i≠ const (i=0, 1,.., n-1); величина h называ­ется шагом интерполяции.

Задача имеет решение, если степень многочлена L(x), которым мы заменяем неизвестную функцию, не выше n.

Пусть нам даны значения функции f(x) в n+1 точке:

y_i = f(x_i) при i = 0, 1, … n.

Тогда существует полином L(x) = a₀ + a₁х + … + a_nхⁿ, значения которого совпадают со значениями f(х) во всех точках х_i. Этот полином называется полиномом Лагранжа и служит примером глобальной интерполяции.

Коэффициенты a_j можно найти, решив следующую систему:

(4)

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

Неизвестные найдем по формулам Крамера:

гдеD— определитель системы (4).

Найдя коэффициенты, можно представить интерполяционный многочлен в виде

Этот многочлен можно представить в другой форме:

где L_j(x) - полиномы n-ой степени, называемые лагранжевыми коэффициентами и имеющие вид:

. (5)

Легко убедиться в справедливости формулы (5). Действительно, степень каждого полинома L_j в точности равна n, сам он равен 1 в точке х = x_j и обращается в нуль в остальных узлах интерполяции. Формула (5) позволяет вычислять значения полинома Лагранжа и без нахождения его коэффициентов.

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

Пример 1. Положим n = 1. Ясно, что мы имеем в этом случае две точки и интерполяционная формула Лагранжа дает уравнение прямой, проходящей через две заданные точки.

.

Примем n = 2. Тогда получим уравнение параболы, проходящей через три точки.

.

n= 3. Пусть заданы значения x₀= 1, x ₁ =3,  x ₂= 7, x ₃ = 12 и y₀= 5.6, y ₁ =6.7,   y ₂= 8.1, y ₃ = 10.3. Определить значение неизвестной функции для х =6, 5.

Для данного случая, когда мы имеем четыре значения функции, интерполяционная формула Лагранжа представляется так:

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

Определим значение функции при х = 6, 5:

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

Для удобства вычислений составим вспомогательную таблицу где х_о, x₁, …, x_n — узлы интерполяции, а х — значение аргумент, для которого определяется приближенное значение по интерполяционной формуле Лагранжа.

,

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

			…
			…
			…
…	…	…	…	…	…
			…

Пример 2. Функция f(х) задана таблично:

Пользуясь интерполяционным многочленом Лагранжа, найти ее значение в точке x = 1.

Решение. Составим таблицу и найдем П₅ (1) = - 48. Приближенное значение функции в точке х =1, т. е f(1) L₄(1), определим по формуле.