Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
3 семестр / Лекции / Лекция 11 - Приближение функций по методу интерполяции.doc
Скачиваний:
0
Добавлен:
16.05.2024
Размер:
577.54 Кб
Скачать

§11.4 Многочлен ньютона с разделенными разностями .

Пусть функция задана на множестве дискретных точек , , … значений аргумента с произвольным (не обязательно постоянным ) шагом.

ОПРЕДЕЛЕНИЕ . Величины

; = (11.8)

называются разделенными разностями первого порядка. Разделенные разности второго порядка определяются формулой:

; ; = (11.9)

Аналогично определяются разделенные разности третьего и более высоких порядков. Общее определение разделенной разности порядка таково:

(11.10)

Аналогично таблице конечных разностей строят таблицу разделенных разностей

Разделенные разности обладают следующими свойствами:

1.Разделенная разность является симметричной функцией своих аргументов, то есть ее значение не меняется при любой перестановке своих аргументов , , ,… .

2.Пусть функция имеет на отрезке, содержащем точки , , ,… , производную порядка k.Тогда справедливо равенство:

.

3. В случае, когда таблица значений аргумента имеет постоянный шаг h,разделенная и конечная разности связаны равенством:

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

+…

+ (11.11)

ПРИМЕР.11.4 Таблично заданную функцию запишем в виде вертикальной таблицы:

Теперь запишем многочлен:

Найдем значения : ,

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

, где

С учетом этого формулу для погрешности интерполяции в точке x, не являющейся узловой, можно записать следующим образом:

Когда величина мала, а функция достаточно гладкая, то справедливо приближенное равенство:

, из которого следует, что

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

Рассмотрим решение следующей задачи (РЗ 16 )

ПРИМЕР 11.5. Пусть функция задана таблицей. Вычислить значение в точке

и оценить погрешность интерполяции.

x

0.5

0.7

1.2

2

f(x)

0.87

1.23

1.41

0.5

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

Для приближения функции многочленом первой степени возьмем два ближайших к точке

, перенумеровав их так: , . Для построения многочлена 2-ой степени

добавим следующий узел . И последняя точка будет иметь номер

Перепишем таблицу с учетом новой нумерации точек. Составим таблицу разделенных разностей, обозначая их более компактным образом:

=0.257

=-0.9

=0.226

=

--0.606(6)

Оценка погрешности.

ОТВЕТ:

ЗАМЕЧАНИЕ. Для построения таблицы была использована функция .