Остаточный член многочлена Лагранжа

Интерполяционная схема Эйткена

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

На первом этапе строится последовательность многочленов первой степени по двум рядом стоящим узлам:

Очевидно, что все построенные многочлены являются многочленами Лагранжа (здесь ). На втором этапе строится последовательность многочленов Лагранжа 2-ой степени, при этом используются многочлены, вычисленные на предыдущем этапе. Расчеты выполняются по формулам:

Количество таких многочленов будет на 1 меньше чем на предыдущем этапе. Наконец на n- ом этапе строится многочлен Лагранжа степени n по формуле

Это и будет значение многочлена Лагранжа степени n в точке , построенного по узлам .

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

Расчеты по интерполяционной схеме Эйткена можно представить в виде следующей таблицы

Таблица 2.2.

Пример 2.2. Рассмотрим задачу вычисления многочлена Лагранжа в заданной точке для таблицы 2.1. Расчеты выполним для точки . Схема Эйткена дает следующие результаты:

; ;

.

Подставив в аналитическое выражение многочлена Лагранжа, полученное в примере 2.1 значение , имеем тот же результат:

.

Если точка не совпадает с узлом интерполирования , то погрешность интерполяционного многочлена (остаточный член)

не равна нулю.

Теорема 2.1. Если непрерывная и раз дифференцируемая функция на интервале , то существует некоторая точка , такая что остаточный член формулы Лагранжа равен:

. (2.12)

Доказательство. Рассмотрим функцию

, (2.13)

где – некоторая, пока неизвестная константа. Учитывая, что и для , функция имеет корень . Выберем константу так, чтобы и точка , не совпадающая ни с каким узлом , тоже была корнем функции . Тогда очевидно, что

, (2.14)

функция будет иметь корня . Пусть точка , тогда на концах каждого из интервала , ,…, , ,…, функция равна нулю. Сформулируем теорему Ролля. Если , то для непрерывной и дифференцируемой функции на интервале найдется не менее одной точки , в которой . Применяя эту теорему к каждому из перечисленных выше интервалов, на концах которых значения функции равны, можно утверждать, что внутри каждого интервала имеется хотя бы один корень. А это означает, что функция на интервале имеет не менее корня. Аналогично можно доказать, что на интервале имеет не менее корней и, наконец, имеет не менее одного корня. Обозначим корень уравнения переменной , тогда дифференцируя (2.13) раз и учитывая, что , получим формулу:

,

заменяя в которой на , получим уравнение для определения :

. (2.15)

Тогда из (2.15), выразив переменную и подставив ее в выражение (2.14), получается формула для остаточного члена (2.12). Теорема доказана.

Оценивая максимальное значение , получим формулу для погрешности метода при интерполировании в точке

, (2.16)

где

. (2.17)

Максимальная погрешность интерполирования на отрезке оценивается значением

.

Пример 2.3. Для функции по значениям в интерполяционных узлах , , построен многочлен Лагранжа 2-ой степени. Требуется оценить погрешность при вычислении многочлена в точке . Оценим величину , где интервал интерполирования зададим по крайним узлам: , . Вычислив производную, имеем . Так как на интервале убывает, то . Погрешности метода, вычисленная по формуле (2.16), равна

.

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

Дата добавления: 2014-12-29; Просмотров: 867; Нарушение авторских прав?; Мы поможем в написании вашей работы!