Интерполяционный многочлен Ньютона для произвольной сетки узлов
Используя форму Лагранжа, представим интерполяционный многочлен в следующем виде:
Здесь ; ; (k = 1,2,…n) – интерполяционные многочлены в форме Лагранжа, построенные по узлам .
Рассмотрим разности
Таким образом, используя формулу (3.37), получим
(4.40)
а интерполяционный многочлен принимает форму
(4.41)
Эта форма называется интерполяционным многочленом Ньютона с раздельными разностями.
Выражение для погрешности имеет тот же вид, что и в случае многочлена Лагранжа [см.формулу(4.9)].
Отметим, что в формуле (4.41) интерполяционного многочлена на узлы накладывается единственное условие - их несовпадение. Поэтому нумерацию узлов можно произвести в произвольном порядке. Например, индексом «0» часто обозначают последующий узел таблицы, за принимают предпоследний узел и обозначают его и т.д. В этом случае многочлен (4.41) принимает форму
(4.42)
и ее называют многочленом Ньютона для интерполирования назад.
Сравнение форм Лагранжа и Ньютона интерполяционного многочлена позволяет рекомендовать использование представления в форме Лагранжа, во-первых, в теоретических исследованиях, например при изучении вопроса о сходимости к ; во-вторых, при интерполировании нескольких функций на одной и той же сетке узлов, поскольку в этом случае можно один раз вычислить множители Лагранжа и использовать их для интерполяции всех функций.
Представление в форме Ньютона оказывается более удобным в практических расчетах. Действительно, число используемых узлов и степень интерполяционного многочлена часто заранее не известно, а при переходе от n узлов к (n+1)-му узлу в форме Ньютона добавляется лишь один член, имеющий смысл поправки к уже вычисленному значению. В то же время в форме Лагранжа добавление еще одного слагаемого сопровождается полным пересчетом полученного ранее результата. Кроме того, в вычислительной практике интерполяция обычно осуществляется на не большом отрезке длиной h<1. При этом слагаемые формы Ньютона имеет порядок ,…, т.е. расположены в порядке убывания, что оказывается полезным при определении точности результата интерполирования.
Дата добавления: 2015-11-06; просмотров: 1903;