Формул с постоянным шагом.

Может быть представлена в виде диагональной таблицы разностей:

x	y	y	2y	3y	4y	Примечание
x-2	y-2		2y-3		4y-4	2-я ИФН
		y–2		3y-3
x-1	y–1		2y-2		4y-3
		y-1		3y-2
x0	y0		2y-1		4y-2	ф. Стерлинга
		y0		3y-1		ф. Бесселя
x1	y1		2y0		4y-1
		y1		3y0
x2	y2		2y1		4y0	1-я ИФН

Мы рассмотрели интерполяционные формулы для равностоящих узлов интерполирования.

Рассмотрим формулы для произвольно заданных узлов интерполирования.

Наиболее часто используется формула Лагранжа.

Интерполяционная формула Лагранжа.

Пусть на отрезке [a;b] даны n+1 различных значений аргумента x: x₀, x₁,…, x_n и известны соответствующие их значению функции y=f(x) : f(x₀)=y₀, f(x₁)=y₁, f(x_n)=y_n. Требуется построить полином степени не выше , имеющий в заданных узлах те же значения, что и функция , т.е. L_n(x_i)=y_i при i=1,n

;

,

где L_i⁽ⁿ⁾- коэффициенты Лагранжа.

Следует отметить, если узлы равностоящие, то интерполяционный полином Лагранжа совпадает с интерполяционной формулой Ньютона.

Примечательно то, что формула Лагранжа зависит лишь от y_i, а не от разностей.

Частные случаи.

n=1

При n=1 имеем 2 точки: (x₀;y₀) и (x₁;y₁).

прямая, проходящая через эти точки-

n=2 (x₀;y₀), (x₁;y₁), (x₂;y₂)

Пример:

L₃(x)=x³+x²-x+2

Для вычисления лагранжевых подмножеств удобно составлять следующую таблицу разности:

Обозначим произведение элементов i-ой строки через D_i , а произведение главной диагонали П_n+1(x). Отсюда следует, что:

П_n+1(x)=(x-x₀)(x-x₁)…(x-x_n)

при i=1,n

Для упрощения вычислений можно использовать инвариантность (при равноотстоящих точках лагранжевых коэффициентов),если

x= at+b

x_j= at_j+b при j=0,n

то L_i⁽ⁿ⁾(x)= L_i⁽ⁿ⁾(t)

Схема Эйткена

Чаще всего требуется найти не общее выражение L_n(x) , а значение его при конкретных x , тогда будет удобно пользоваться интерполяционной схемой Эйткена:

Последовательно вычисляются многочлены:

и т.д.

Вычисления по схеме Эйткена удобно расположить в таблице:

Вычисления по схеме Эйткена обычно ведутся до тех пор, пока последовательные значения L_01…n(x) и L_01…n(n+1) не совпадут по заданной точности.

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

Пример: x=27, =0,1