Интерполяционные полиномы Гаусса, Стирлинга и Бесселя.

При построении интерполяционных формул Ньютона используются лишь значения функции, которые расположенные по одну сторону от выбранного начального значения, то есть эти формулы носят односторонний характер. Во многих случаях оказываются полезные интерполяционные формулы, которые содержат в себе как следующие, так и предшествующие значения функции по отношению к ее начальному значению. Чаще используются разности, расположенные в горизонтальной строке диагональной таблицы разностей данной функции, которые отвечают начальным значениям x₀ и y₀, или в строках, которые непосредственно примыкают к ней. Эти разности Δ f_-1, Δ f_0, Δ f²_-1... называют центральными разностями (табл.4), где

Таблица 4.

Соответствующие интерполяционные формулы носят название интерполяционных формул с центральными разностями. К ним относятся формулы Гаусса, Стирлинга, Бесселя.

Пусть есть 2n+1 равноотстоящих узлов интерполирования х_-1, х _-(n-1),..., х _-1, х₀, х₁,... х_n-1, х_n, где х_i=х_i+1–x_i = h = const (i= -n, -(n-1),..., n –1)) и для функций f(x) известные ее значения в этих узлах f(x_i) (i = 0, ± 1,..., ± n). Нужно построить многочлен степени не выше 2n такой, что P(x_i) = f(x_i), при i = 0, ± 1,..., ± n.

Приведем интерполяционные формулы, опуская их вывод. Первая интерполяционная формула Гаусса:

, (12)

где а – центральные разности, которые образуют в табл.4 нижнюю ломаную строку.

Вторая интерполяционная формула Гаусса:

, (13)

где – центральные разности, которые образуют в табл.4 верхнюю ломаную строку.

Формулы Гаусса используются для интерполяции в середине таблицы вблизи x₀: первая при х > х₀, вторая — для х < x₀.

Интерполяционная формула Стирлинга получается, если взять среднее арифметическое интерполяционных формул Гаусса (12) и (13):

Из второй интерполяционной формулы Гаусса можно получить интерполяционную формулу Бесселя в виде

При более детальном рассмотрении интерполяционных формул оказывается, что при | t |≤ 0, 25 целесообразно использовать формулу Стирлинга, а при 0, 25 ≤ t ≤ 0, 75 — формулу Бесселя.

Остаточный член интерполяционных формул Гаусса и Стирлинга для порядка 2n максимальных использованных разностей из Таблицы 4 и при х Î [x₀ – nh, x₀ +nh ] имеет вид:

Если же аналитическое выражение функции неизвестно, то при малом h:

Если (2n+1) -порядок максимально использованной разности из таблицы и х Î [x₀–nh, x₀ + (n+1)h], то остаточный член интерполяционной формулы Бесселя запишется в виде:

Если же функция f(x) заданна таблично и шаг h мал, то

В частности, при t= 0, 5 получим погрешность интерполяции на середину

или

.

Пример