Интерполирование функции алгебраическими многочленами

Заданы различные точки (числа ) x₀,x₁, x₂,…, x_n лежащие в области определения функции f(x), Точки x_i, i=0,1,…n, называются узлами интерполяции. Требуется построить многочлен L_n(x), степени не выше n, значения которого в узлах интерполяции совпадают со значениями функции f(x), то есть

L_n(x_i)=f(x_i), i=0,1,…,n (1)

Многочлен L_n(x_i), удовлетворяющий условиям (1), называется интерполяционным многочленом, простроенным для функции f(x) по заданным узлам.

Множеством F служит множество функций, определенных в узлах интерполяции; множеством H – множество многочленов, степени не выше n; мерой близости служит мера . Для интерполяционного многочлена : , поэтому L_n(x) наилучшее приближение по этой мере.

Если искомый многочлен L_n(x) представить в виде:

L_n(x)= a_nxⁿ+a_n-1x^n-1+..+a₁x+a_0,

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

a_nx₀ⁿ+a_n-1x₀^n-1+..+a₁x₀+a₀=f(x₀)

a_nx₁ⁿ+a_n-1x₁^n-1+..+a₁x₁+a₀=f(x₁) (2)

…………..

a_nx_nⁿ+a_n-1x_n^n-1+..+a₁x_n+a₀=f(x_n)

Система (2) содержит (n+1) уравнений и (n+1) неизвестных a₀,a₁,…, a_{n .}

Разрешимость поставленной задачи сводится к разрешимости системы линейных уравнений (2). Система разрешима при любых правых частях, если ее определитель отличен от нуля.

Определитель системы является определителем Вандермонда, который не равен нулю, поэтому система уравнений имеет единственное решение.

Рассмотрим пример составления интерполяционного многочлена непосредственно из его определения.

Пример 8.

В узлах x₀=0, x₁=1,x₂=3 функция f(x) принимает значения y₀=-1 y₁=2, y₂= 5. Построить интерполяционный многочлен, отвечающий этим данным.

В нашем случае n=2. Следовательно, требуется построить интерполяционный многочлен L₂(x)=a₂x²+a₁x+a_0.

Условия задачи приводят к системе уравнений:

a₂0²+a₁0+a₀=-1

a₂1²+a₁1+ a₀ =2

a₂3²+a₁3+a₀=5

Ее решение: a₂ = -1/2, a₁ = 7/2, a₀ = -1.

Следовательно, L₂(x)= .