Интерполяция многочленами. В общем виде задача интерполяции многочленами формулируется следующим образом

В общем виде задача интерполяции многочленами формулируется следующим образом. Пусть на отрезке [ x ₁, x_n ] таблично задана функция абсциссами x ₁, x ₂,..., x_n и ординатами y ₁, y ₂,..., y_n своих узловых точек. Требуется построить степенной полином вида

,

значения которого в точках x ₁, x ₂,..., x_n совпадает со значениями этой табличной функции

Такой полином всегда существует и оказывается единственным. Для вычисления коэффициентов этого полинома пользуются условиями равенства его значений и значений таблично заданной функции в узлах интерполяции. Это даёт систему из n линейного алгебраического уравнения относительно коэффициентов a ₀, a ₁,..., a_n вида

которая в матричной форме записывается следующим образом

и которая может быть разрешена любым из известных методов решения систем линейных алгебраических уравнений.

Использование такого подхода имеет ряд недостатков. Во-первых, с увеличением количества n узлов интерполяции пропорционально возрастает и степень системы линейных алгебраических уравнений, что приводит к усложнению ее решения. Во-вторых, с ростом n возрастает вероятность «биения» функции между узлами интерполяции за счёт членов полинома с высокими степенями.

Рассмотрим применение этого подхода на примере кусочно-линейной интерполяции таблично заданной функ­ции и поиска её значения при аргументе х = 1.6.

Для решения задачи строится полином 3-го порядка

,

коэффициенты которого есть решение следующей системы линейных алгебраических уравнений

В матричной форме эта система имеет вид

.

Её решение получается любым доступным методом, например методом Гаусса

.

Таким образом, табличная функция в случае интерполяции полиномом представляется в виде

.

Её значение в заданной точке x = 1.6 будет

.

Ниже на рис.2 представлен фрагмент рабочей книги Excel с реализацией интерполяции полиномом. При построении графика приближающей функции аргумент х изменяется с шагом 0.2, а значения функции вычисляются по общей формуле, адаптированной под конкретные значения из диапазона изменения аргумента.

Многочлен Лагранжа (J.L.Lagrange, 1795)

Представляет собой случай полиномиального представления приближающей функции, когда она ищется в виде линейной комбинации базисных функций j_k (x), которые должны быть определены для всего отрезка интерполяции [ x ₁, x_n ], линейно независимы, и их количество должно быть равно числу узлов таблично заданной функции

.

Коэффициенты с ₁, с ₂,..., с_n определяются исходя из условий ра-

Рис 2.

венства значений приближающей и исходной функций при табличных значениях аргумента, что сводит задачу к системе n линейных алгебраических уравнений относительно них, а в качестве функций j_k (x) используются полиномы (n –1) степени

,

которые для пяти узловых точек записываются в виде

,

,

,

,

.

Для каждого полинома характерно то, что для всех значений x_i в узловых точках он принимает нулевые значения, кроме k -ой, где его значение равно единице. Графики этих полиномов пред­ставлены на рис.3.

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

.

Процесс построения интерполирующего многочлена Лагранжа для пяти узловых точек показан на рис.4.

Рассмотрим работу метода на приведённом выше примере интерполяции таблично заданной по 4-м точкам функции. Сначала строятся четыре базовых полинома:

,

,

,

.

Они позволяют записать интерполирующий многочлен Лагранжа в виде

.

Для аргумента x = 1.6 многочлен Лагранжа даёт

Ниже на рис.5 представлен фрагмент рабочей книги Excel с реализацией интерполяции с помощью многочлена Лагранжа. Как видно из рисунка для аргумента x = 1.6 многочлен Лагранжа дал значение 0.528.

Для контроля правильности вычислений многочлена Лагранжа полезно строить графики базовых полиномов. Для рассматриваемого примера они приведены на рис.6.