3.7.Интерполяционная формула Ньютона (для равностоящих узлов)

Используем свойства конечных разностей при выводе формулы Ньютона для равноотстоящих узлов:

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

а). для интерполирования вперед:

(8)

в). для интерполирования назад:

(9)

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

Обозначим через интерполяционный многочлен, построенный по узлам

Многочлен можем представить в виде: (8`)

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

.

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

.

По свойству 5 конечных разностей,

.

По свойству 4, .

В узлах интерполяции значения многочлена совпадают со значениями функции f(x). Значения разности определяются узловыми значениями функции, поэтому

и

Следовательно, = , к = 1,2,…,n, .

Подставляя эти выражения в формулу (8'), получим интерполяционную формулу Ньютона:

Величины , входящие в формулу для интерполирования вперёд, находятся на верхней диагонали таблицы. Аналогичные величины для второй формулы находятся на нижней диагонали (см. таблицу).

Пример 11. Составить интерполяционные формулы Ньютона по приведённой таблице.

x
1,00	2,7183	0,2859 0,3159	0,0300
1,10	3,0042
1,20	3,3201

Вычислить .

a.) По формуле (8) для интерполирования вперёд:

b.) По формуле (9) для интерполирования назад:

3.8. Погрешность формулы Ньютона.

Полная погрешность формулы Ньютона по общему правилу выражается формулой:

Вычислим устранимую и неустранимую погрешности метода. Оценим разность . Интерполяционный многочлен Ньютона тождественно совпадает с многочленом Лагранжа, построенным по тем же узлам, следовательно, остаётся неизменной оценка:

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

Используя формулу (8), видим, что

так как

Если многочлен Ньютона, построенный по приближённым значениям , то

Таким образом,

Пример 12. Оценить погрешность вычисления в примере 10.

Имеем, , т.е., . Следовательно,

Далее, и, следовательно,

Таким образом,

Если округлить, полученные в примере 10 значения до четырёх верных знаков после запятой, то .

Тогда .

Следовательно

3.9. Интерполяционные сплайны. Параболические интерполяционные сплайны

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

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

Для приближённого представления функции f(x) на отрезке выберем на нём две системы узлов:

а) узлы , называемые узлами сплайна.

в) узлы , называемые узлами интерполяции.

Узлы выбираем так, что для i=1,2,…,n-1 имеем , т.е. между двумя внутренними узлами сплайна лежит один узел интерполяции.

Функция называется параболическим интерполяционным сплайном для функции f(x), построенным по указанным узлам, если:

1.Функция и её производная непрерывны на отрезке .

2. На каждом отрезке между двумя последовательными узлами сплайна функция является многочленом второй степени:

для

3. В узлах интерполяции значения сплайна совпадают со значениями функции f(x).

, для n=2 .

Запишем систему условий, определяющих параболический интерполяционный сплайн :

a) условия интерполяции

= (1)

Так как узел

b) условия непрерывности сплайна в узлах сплайна.

В каждом внутреннем узле сплайна стыкуются два участка и . На первом из них сплайн , а на втором . Условие непрерывности в точке даёт

= (2)

с) условия непрерывности производной во внутренних узлах сплайна приводят к соотношениям , (3)

Узлы сплайна разбивают весь отрезок на (n+1) частей, на каждой из которых сплайн определяется тремя коэффициентами , поэтому необходимо определить 3(n+1)=3n+3 неизвестных коэффициентов.

Указанные выше условия (1),(2),(3) содержат соответственно n+1, n и n уравнений. Число неизвестных на 2 превосходит число уравнений. Для того, чтобы система имела лишь одно решение, нужно добавить ещё два условия.

1. Если функция f(x) периодическая с периодом Т=b-a ,то должно выполняться условие периодичности производных сплайна:

и (4)

В развёрнутом виде: .

3. В крайних узлах сплайна должны совпадать значения первых производных сплайна и функции:

,т.е.

(5)

Можно доказать, что система уравнений (1)-(3), дополненная условиями (4) или (5), имеет единственное решение.

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

При выполнении указанных дополнительных условий (4) или (5) справедливы оценки:

, для любого , (6)

где , -верхняя граница значений на отрезке .

Из оценок (6) следует, что при , значения сплайна и его производной стремятся соответственно, к значению функции и её производной:

и

Получена погрешность метода f(x)-

Неустранимую оценку погрешности вычислять не будем.

Пример 13.

На отрезке выбраны:

А) узлы сплайна:

В) узлы интерполяции:

Найти параболический сплайн , для функции

Составим систему уравнений, определяющих параболический сплайн :

1)Условия интерполяции:

2)Условия непрерывности во внутренних узлах сплайна:

;

3)Условия непрерывности производной во внутренних узлах сплайна:

4)Дополнительное условие: в крайних узлах сплайна значение его производной совпадает со значением производной функции:

Взяв из таблиц значение функции ,sinx для указанных значений аргумента, получим:

Решая систему получим:

Следовательно, сплайн определяется следующим образом:

=