Первая интерполяционная формула Ньютона
практически неудобна для интерполирования функции вблизи узлов таблицы. В этом
случае обычно применяется вторая интерполяционная
формула Ньютона.
Описание задачи. Пусть
имеем последовательность значений функции
,
для равноотстоящих значений
аргумента , где - шаг интерполяции. Построим
полином следующего вида:
,
или, используя обобщённую
степень, получаем:
. (1)
Тогда, при выполнении равенства , , получим
, .
Подставим эти значения в формулу
(1). Тогда, окончательно, вторая
интерполяционная формула Ньютона имеет вид:
. (2)
Введём более удобную запись формулы (2).
Пусть , тогда , и т. д.
Подставив эти значения в формулу (2), получим:
. (3)
Это и есть обычный вид второй интерполяционной формулы Ньютона.
Для приближённого вычисления значений функции полагают:
.
Как первая, так и вторая интерполяционные
формулы Ньютона могут быть использованы для экстраполирования функции, т. е.
для нахождения значений функции для значений
аргументов , лежащих вне пределов таблицы. Если и близко к , то выгодно применять первую интерполяционную формулу
Ньютона, причём тогда . Если же и близко к , то удобнее пользоваться второй интерполяционной формулой
Ньютона, причём . Таким
образом, первая интерполяционная формула Ньютона обычно используется для интерполирования вперёд и экстраполирования назад, а вторая интерполяционная
формула Ньютона, наоборот, - для интерполирования
назад и экстраполирования вперёд.
Заметим, что операция экстраполирования,
вообще говоря, менее точна, чем операция интерполирования в узком смысле
слова.
Пример. Приняв шаг , построить интерполяционный полином Ньютона для функции , заданной таблицей
|
0,5
|
0,55
|
0,6
|
0,65
|
0,7
|
0,75
|
0,8
|
|
0,875
|
0,7088
|
0,5361
|
0,3572
|
0,173
|
-0,0156
|
-0,2081
|
Решение. Составляем таблицу разностей (таблица 1). Так как разности
третьего порядка практически постоянны, то в формуле (3) полагаем . Приняв , , будем иметь:
,
или
,
где .
Это и есть искомый интерполяционный полином Ньютона.
Таблица 1
|
|
|
|
|
0,5
0,55
0,6
0,65
0,7
0,75
0,8
|
0,875
0,7088
0,5361
0,3572
0,173
-0,0156
-0,20
|
-0,1662
-0,1727
-0,1789
-0,1842
-0,1886
-0,1925
|
-0,0065
-0,0062
-0,0053
-0,0044
-0,0039
|
0,0003
0,0009
0,0009
0,0005
|
|