4.3.8 Интерполяционный многочлен Ньютона для произвольной сетки узлов

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

Здесь ;; (k = 1,2,…n) – интерполяционные многочлены в форме Лагранжа, построенные по узлам .

Рассмотрим разности

Таким образом, используя формулу (3.37), получим

(4.40)

а интерполяционный многочлен принимает форму

(4.41)

Эта форма называется интерполяционным многочленом Ньютона с раздельными разностями.

Выражение для погрешности имеет тот же вид, что и в случае многочлена Лагранжа [см.формулу(4.9)].

Отметим, что в формуле (4.41) интерполяционного многочлена на узлы накладывается единственное условие - их несовпадение. Поэтому нумерацию узлов можно произвести в произвольном порядке. Например, индексом «0» часто обозначают последующий узел таблицы, за принимают предпоследний узел и обозначают егои т.д. В этом случае многочлен (4.41) принимает форму

(4.42)

и ее называют многочленом Ньютона для интерполирования назад.

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

Представление в форме Ньютона оказывается более удобным в практических расчетах. Действительно, число используемых узлов и степень интерполяционного многочлена часто заранее не известно, а при переходе от n узлов к (n+1)-му узлу в форме Ньютона добавляется лишь один член, имеющий смысл поправки к уже вычисленному значению. В то же время в форме Лагранжа добавление еще одного слагаемого сопровождается полным пересчетом полученного ранее результата. Кроме того, в вычислительной практике интерполяция обычно осуществляется на не большом отрезке длиной h<1. При этом слагаемые формы Ньютона имеет порядок ,…, т.е. расположены в порядке убывания, что оказывается полезным при определении точности результата интерполирования.