
- •Глава 4. Многочленная интерполяция
- •4.1.Интерполяционный многочлен Лагранжа Постановка задачи интерполирования
- •Задача многочленной интерполяции. Ищется алгебраический многочлен n-й степени , удовлетворяющий условиям интерполяции
- •Интерполяционный многочлен в форме Лагранжа
- •Оценка погрешности интерполяции
- •Теорема 1. Пусть функция имеет на отрезкепроизводные до-го порядка включительно. Тогда для любогох из отрезка найдется точкатакая, что
- •Доказательство. Рассмотрим вспомогательную функцию
- •Первое свойство разделенных разностей. Разделенные разности – симметричные функции своих аргументов, то есть они не меняются при перестановке аргументов.
- •Интерполяционная формула Ньютона
- •Доказательство. Пусть представляет собой некоторый многочленn-й степени . Учитывая, что, перепишем равенство (4.2.4) для случая, когда:
- •Алгоритм построения интерполяционного многочлена Ньютона
- •4.3. Кратные узлы интерполяции. Интерполяционный многочлен Эрмита Кратные узлы интерполяции. Формулировка задачи интерполяции с кратными узлами
- •Предел интерполяционного многочлена Ньютона при слиянии узлов
- •Интерполяционная формула Эрмита
- •Доказательство. Требуется доказать, что для полученного предела интерполяционного многочлена выполняются условия интерполяции (4.3.3). Для этого мы зафиксируем произвольный номер и докажем, что
- •Оценка погрешности интерполяционного многочлена Эрмита
- •Теорема 2. Пусть функция имеет на отрезкепроизводные до-го порядка включительно. Тогда для любогох из отрезка найдется точкатакая, что
- •4.4. Обратное интерполирование. Кусочно - многочленная интерполяция. Интерполирование сплайнами Обратное интерполирование
- •Кусочно-многочленная интерполяция
- •Интерполирование сплайнами Сплайны. Пусть на отрезке задана сетка точек, которая делитнаn отрезков ,. Сплайномn-го порядка называется функция
- •4.5. Многочлены Чебышева. Чебышевские узлы интерполяции Многочлены Чебышева
- •Контрольные вопросы и задания
Первое свойство разделенных разностей. Разделенные разности – симметричные функции своих аргументов, то есть они не меняются при перестановке аргументов.
Для разностей первого порядка это проверяется элементарно:
.
Для разделенных разностей более высоких порядков это свойство доказывается методом математической индукции.
Второе свойство
разделенных разностей.
Пусть
–
произвольный многочлен n-го
порядка,
заданные попарно различные точки, а
точка х
не совпадает ни с одной из них. Тогда
–многочлен
-й
степени,
–многочлен
-й
степени,
…
–многочлен нулевой
степени, то есть постоянная, равная
,
.
Доказательство.
По определению
.
Числитель этой дроби,
,
представляет собой многочленn-го
порядка, который обращается в нуль в
точке
.
Следовательно,
– многочлен
-й
степени, что и требовалось доказать.
Рассмотрим следующую
разделенную разность. По определению
.
Числитель этой дроби
представляет собой многочлен
-го
порядка, который обращается в нуль в
точке
.
Следовательно,
– многочлен
-й
степени, что и требовалось доказать.
Продолжим этот процесс. Каждая последующая разделенная разность, очевидно, представляет собой многочлен, степень которого на единицу меньше, чем у предыдущего. Поэтому
представляет собой многочлен нулевой степени, то есть постоянную, значение которой не зависит от х. Поэтому
.
Отсюда
.
Второе свойство полностью доказано.
Интерполяционная формула Ньютона
Пусть
–
некоторая функция с известными значениями
в узлах
.
Рассмотрим разделенную разность
.
Выразим из этого равенства
:
.
(4.2.1)
Рассмотрим
разделенную разность
.
Выразим из этого равенства
:
Подставим это выражение в (4.2.1):
.
(4.2.2)
Рассмотрим
разделенную разность
.
Выразим из этого равенства
:
.
Подставим это выражение в равенство (4.2.2):
.
(4.2.3)
Продолжим этот процесс пока возможно и получим
…………………………………………………
.
(4.2.4)
Сумма всех членов
правой части равенства (4.2.4), кроме
последнего, представляет собой
алгебраический многочлен n-й
степени. Обозначим его через
.
…
….
(4.2.5)
Следовательно, формулу (4.2.4) можно записать в виде
.
(4.2.6)
Заданный формулой
(4.2.5) алгебраический многочлен
представляет собой интерполяционный
многочлен для функции
с узлами интерполяции
.
Доказательство. Пусть представляет собой некоторый многочленn-й степени . Учитывая, что, перепишем равенство (4.2.4) для случая, когда:
…
…...
(4.2.7)
Это представление
справедливо для любого многочлена n-й
степени, в том
числе и для интерполяционного многочлена.
Пусть теперь
- интерполяционный многочлен. Тогда, по
условиям интерполяции,
.
Отсюда следует, что
,
,
,
…………………………………….
.
(4.2.8)
Заменив в представлении (4.2.7) разделенные разности по формулам (4.2.8), получим представление интерполяционного многочлена по формуле (4.2.5). Что и требовалось доказать.
Интерполяционный многочлен, представленный в виде (4.2.5), называют интерполяционным многочленом Ньютона.
Заметим, что из равенства (4.2.6) получается еще одно выражение для погрешности интерполяции
.
(4.2.9)
Сравнивая это
выражение с выражением (4.1.10), легко
видеть, что должно быть справедливо
следующее утверждение. Для
любого
найдется
такое, что будет справедливо равенство
.
(4.2.10)
Интерполяционную формулу Ньютона (4.2.5) можно записать и более компактно
.
(4.2.11)