- •Глава 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 получим . Абсолютная погрешность этого приближенного значения функции равна. Ее называют такжеабсолютной погрешностью интерполяции с кратными узлами. Будем искать оценку абсолютной погрешности интерполяции. Она основывается на теореме, которую мы приведем без доказательства.