Интерполяционная формула Ньютона

Интерполяционная формула Ньютона позволяет выразить интерполяционный многочлен P_n(x) через значение f(x) в одном из узлов и через разделенные разности функции f(x), построенные по узлам x₀, x₁,…, x_n. Эта формула является разностным аналогом формулы Тейлора:

(7.4)

Прежде чем приводить формулу Ньютона, рассмотрим сведения о разделенных разностях. Пусть в узлах известны значения функцииf(x). Предполагаем, что среди точек x_k, k = 0, 1,…, n нет совпадающих. Тогда разделенными разностями первого порядка называются отношения

(7.5)

Будем рассматривать разделенные разности, составленные по соседним узлам, то есть выражения . По этим разделенным разностям первого порядка можно построить разделенные разности второго порядка:

(7.6)

Аналогично определяются разности более высокого порядка. То есть пусть известны разделенные разности k-го порядка тогда разделенная разность k+1-го порядка определяется как

(7.7)

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

(7.8)

Показано, что интерполяционный многочлен Лагранжа (7.3) совпадает с интерполяционным многочленом Ньютона (7.8).

Замечания

• В формуле (7.8) не предполагалось, что узлы x₀, x₁,…, x_n расположены в каком-то определенном порядке. Поэтому роль точки x₀ в формуле (7.8) может играть любая из точек x₀, x₁,…, x_n. Соответствующее множество интерполяционных формул можно получить из (7.8), перенумеровав узлы. Например, тот же самый многочлен P_n(x) можно представить в виде

(7.9)

• Если то (7.8) называется формулой интерполирования вперед, а (7.9) - формулой интерполирования назад.
• Интерполяционную формулу Ньютона удобнее применять в том случае, когда интерполируется одна и та же функция f(x), но число узлов интерполяции постепенно увеличивается. Если узлы интерполяции фиксированы и интерполируется не одна, а несколько функций, то удобнее пользоваться формулой Лагранжа

Сходимость интерполяционного процесса

Обсудим следующий вопрос: будет ли стремиться к нулю погрешность интерполирования f(x) – L_n(x), если число узлов n неограниченно увеличивать:

1. Свойства сходимости или расходимости интерполяционного процесса зависят как от выбора последовательности сеток, так и от гладкости функции f(x).

2. Известны примеры несложных функций, для которых интерполяционный процесс расходится.

Так последовательность интерполяционных многочленов, построенных для непрерывной функции по равноотстоящим узлам на отрезке
[-1, 1], не сходится к функции ни в одной точке отрезка [-1, 1], кроме точек –1, 0, 1. На рис. 7.2 в качестве иллюстрации изображен график многочленаL₉(x) при , построенного для функции по равноотстоящим узлам на отрезке [-1,1].

Рис. 7.2. Сходимость интерполяционных многочленов

3. Чтобы избежать этих некорректностей, в практике вычислений обычно избегают пользоваться интерполяционными многочленами высокой степени.

Сжимающие отображения.

Понятие сжимающего отбражения. Неподвижные точки. Метод простой итерации для операторного уравнения с сжимающим оператором. Оценка погрешности. Примеры: решение систем линейных алгебраических уравнений, решение нелинейных уравнений и систем.