
- •Глава 4. Многочленная интерполяция
- •4.1.Интерполяционный многочлен Лагранжа Постановка задачи интерполирования
- •Задача многочленной интерполяции. Ищется алгебраический многочлен n-й степени , удовлетворяющий условиям интерполяции
- •Интерполяционный многочлен в форме Лагранжа
- •Оценка погрешности интерполяции
- •Теорема 1. Пусть функция имеет на отрезкепроизводные до-го порядка включительно. Тогда для любогох из отрезка найдется точкатакая, что
- •Доказательство. Рассмотрим вспомогательную функцию
- •Первое свойство разделенных разностей. Разделенные разности – симметричные функции своих аргументов, то есть они не меняются при перестановке аргументов.
- •Интерполяционная формула Ньютона
- •Доказательство. Пусть представляет собой некоторый многочленn-й степени . Учитывая, что, перепишем равенство (4.2.4) для случая, когда:
- •Алгоритм построения интерполяционного многочлена Ньютона
- •4.3. Кратные узлы интерполяции. Интерполяционный многочлен Эрмита Кратные узлы интерполяции. Формулировка задачи интерполяции с кратными узлами
- •Предел интерполяционного многочлена Ньютона при слиянии узлов
- •Интерполяционная формула Эрмита
- •Доказательство. Требуется доказать, что для полученного предела интерполяционного многочлена выполняются условия интерполяции (4.3.3). Для этого мы зафиксируем произвольный номер и докажем, что
- •Оценка погрешности интерполяционного многочлена Эрмита
- •Теорема 2. Пусть функция имеет на отрезкепроизводные до-го порядка включительно. Тогда для любогох из отрезка найдется точкатакая, что
- •4.4. Обратное интерполирование. Кусочно - многочленная интерполяция. Интерполирование сплайнами Обратное интерполирование
- •Кусочно-многочленная интерполяция
- •Интерполирование сплайнами Сплайны. Пусть на отрезке задана сетка точек, которая делитнаn отрезков ,. Сплайномn-го порядка называется функция
- •4.5. Многочлены Чебышева. Чебышевские узлы интерполяции Многочлены Чебышева
- •Контрольные вопросы и задания
Интерполяционная формула Эрмита
Докажем, что полученный в предыдущем пункте предел интерполяционного многочлена при слиянии узлов, и есть искомый интерполяционный многочлен Эрмита
.
(4.3.14)
Доказательство. Требуется доказать, что для полученного предела интерполяционного многочлена выполняются условия интерполяции (4.3.3). Для этого мы зафиксируем произвольный номер и докажем, что
.
(4.3.15)
Интерполяционный многочлен не зависит от нумерации узлов. Поэтому и его предел также не будет зависеть от нумерации узлов и
…
+
…
+
.
(4.3.16)
Здесь
использовано стандартное представление
многочлена Ньютона и представление
разделенных разностей с кратными узлами.
Но разложение записано не до конца.
Сумма всех остальных членов записана
в виде
.
Функция
представляет собой алгебраический
многочлен степени
,
не обращающийся в 0 в точке
.
Произведение
и все его производные до
-го
порядка включительно, очевидно, обращаются
в 0 в точке
.
Поэтому из равенства (4.3.16) следуют
доказываемые равенства (4.3.15). Что и
требовалось доказать.
Итак, интерполяционный многочлен Эрмита построен. Алгоритм его построения аналогичен алгоритму построения многочлена Ньютона. Единственное отличие заключается в том, что у многочлена Эрмита имеются табличные разности с кратными узлами.
Пример 2
Для
функции
построить многочлен Эрмита с узлами
,
и значениями
,
,
,
.
Запишем многочлен Эрмита в общем виде
.
Вычислим табличные разности и запишем их в табл. 4.4
,
,
,
,
,
,
,
,
.
Таблица 4.4
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Осталось составить интерполяционный многочлен
.
Построенный
интерполяционный многочлен в рассмотренном
примере совпал с исходной функцией
.
Это не случайно и непосредственно
следует из единственности интерполяционного
многочлена Эрмита.
Пример 3
Даны значения
функции
и ее производных в узле
:
.
Требуется построить по этим данным
интерполяционный многочлен Эрмита
.
Запишем его в виде интерполяционного многочлена Ньютона с кратными узлами и, учитывая выражение (4.3.12), получим
.
Построенный многочлен Эрмита совпал с многочленом Тейлора. Таким образом, многочлен Эрмита является обобщением и объединяет в себе многочлен Тейлора и интерполяционный многочлен.
Оценка погрешности интерполяционного многочлена Эрмита
Пусть
функция
определена на заданном отрезке
и известны табличные значения функции
и ее производных в заданных узлах
,
содержащихся на отрезке
,
причем известными считаются значения
функции и ее производных
…….;
Здесь
,
.
.
По этим данным строится интерполяционный
многочлен Эрмита
,
который вместе со своими производными
имеет те же табличные значения в узлах.
Используем его в качестве приближения
для функции
на отрезке
.
Тогда точное значение функцииy
в точке x
будет равно
,
а в качестве приближенного значения
функцииy
в точке x
получим
.
Абсолютная погрешность этого приближенного
значения функции равна
.
Ее называют такжеабсолютной
погрешностью интерполяции с кратными
узлами. Будем
искать оценку абсолютной погрешности
интерполяции. Она основывается на
теореме, которую мы приведем без
доказательства.