Оценка погрешностей интерполяционных формул




 

Согласно условиям задачи интерполирования, построенные для аппроксимируемой функции многочлены , и принимают значения , ,…, , совпадающие со значениями функции только в точках x 0, x 1, x 2 ,…, xn интервала [ a, b ]. Для других значений значения интерполяционного многочлена, очевидно, будут отличаться от значений на некоторую величину

.

Возникает вопрос, какова величина остаточного члена интерполяционной формулы. Величина остаточного члена зависит от ряда факторов - свойств функции , числа и расположения узлов xi, шага h, положения точки x * интерполирования. Поэтому в самом общем случае этот вопрос является открытым и подлежит исследованию для каждого конкретного аппроксимационного процесса.

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

,

т. е. отклонение функции от ее интерполяционного многочлена в точке x = x *.

Если точка x * интерполирования заранее неизвестна, а интерполирование осуществляется на всем отрезке [ a, b ], то за меру погрешности принимают максимальное значение абсолютной величины разности между функцией и интерполяционным многочленом на отрезке [ a, b ] (рисунок 2)

.

В этом случае считается, что многочлен равномерно аппроксимирует функцию с точностью ∆ на отрезке [ a, b ]. При этом величина ∆ называется равномерной оценкой отклонения функции от ее интерполяционного многочлена.

Заметим, что в условиях задачи интерполирования не накладываются какие-либо ограничения на функцию . Однако, если предполагать, что аппроксимируемая функция на отрезке [ a, b ] непрерывна и имеет все производные , ,…, до (n +1) – го порядка включительно, то равномерная оценка ее погрешности интерполирования может быть найдена, исходя из вида используемого интерполяционного многочлена [2]. Так, при интерполировании функции по формуле Лагранжа для заданных произвольных значений аргумента x: x 0, x 1, x 2 ,…, xn получена следующая равномерная оценка ее абсолютной погрешности:

(26)

где

; (27)

Для фиксированной точки x = x * оценка погрешности производится по формуле

.

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

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

, , .

По формуле (27) получим

.

По формуле (26) находим:

.

При использовании формул Ньютона (20), (25) оценка погрешности найденных значений интерполируемой функции несколько упрощается. Эти формулы получены для равноотстоящих узлов интерполирования x 0, x 1, x 2 ,…, xn, образующих арифметическую прогрессию: , ,…, , где h – шаг изменения аргумента. Выполнив постановку в формуле (26) переменной t

,

получим

. (28)

Эта формула позволяет определить равномерную оценку первой интерполяционной формулы Ньютона.

Аналогично, полагая в формуле (26) , получим равномерную оценку второй интерполяционной формулы Ньютона

.

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

,

приближенно можно принять

.

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

.

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

.

 

 

Литература.

1. Данилина Н.И. и др. Вычислительная математика: Учеб. пособие для техникумов. – М.: Высш. шк., 1985, с. 256-287, 294-298, 306-308.

2. Демидович Б.П. Основы вычислительной математики: Учеб. пособие для высш. техн. учеб. заведений. М.: ФИЗМАТЛИТ, 1960, с. 497-519, 527-539.

3. Копченова Н.Н., Марон И.А. Вычислительная математика в примерах и задачах. Учеб. пособие. – СПб.: Изд-во «Лань», 2008, с. 100-107, 113-117.

4. Турчак Л.И., Плотников П.В. Основы численных методов: Учеб. пособие. – М.: ФИЗМАТЛИТ, 2005, с. 31-36,47-58.

5. Гутер Р.С., Овчинский Б.В. Элементы численного анализа и математической обработки результатов опыта. М.: Изд-во «Наука» 1970, с. 93-128, 133-135.

 

Контрольные вопросы и упражнения для приобретения



Поделиться:




Поиск по сайту

©2015-2024 poisk-ru.ru
Все права принадлежать их авторам. Данный сайт не претендует на авторства, а предоставляет бесплатное использование.
Дата создания страницы: 2019-01-11 Нарушение авторских прав и Нарушение персональных данных


Поиск по сайту: