
- •Глава 4. Многочленная интерполяция
- •4.1.Интерполяционный многочлен Лагранжа Постановка задачи интерполирования
- •Задача многочленной интерполяции. Ищется алгебраический многочлен n-й степени , удовлетворяющий условиям интерполяции
- •Интерполяционный многочлен в форме Лагранжа
- •Оценка погрешности интерполяции
- •Теорема 1. Пусть функция имеет на отрезкепроизводные до-го порядка включительно. Тогда для любогох из отрезка найдется точкатакая, что
- •Доказательство. Рассмотрим вспомогательную функцию
- •Первое свойство разделенных разностей. Разделенные разности – симметричные функции своих аргументов, то есть они не меняются при перестановке аргументов.
- •Интерполяционная формула Ньютона
- •Доказательство. Пусть представляет собой некоторый многочленn-й степени . Учитывая, что, перепишем равенство (4.2.4) для случая, когда:
- •Алгоритм построения интерполяционного многочлена Ньютона
- •4.3. Кратные узлы интерполяции. Интерполяционный многочлен Эрмита Кратные узлы интерполяции. Формулировка задачи интерполяции с кратными узлами
- •Предел интерполяционного многочлена Ньютона при слиянии узлов
- •Интерполяционная формула Эрмита
- •Доказательство. Требуется доказать, что для полученного предела интерполяционного многочлена выполняются условия интерполяции (4.3.3). Для этого мы зафиксируем произвольный номер и докажем, что
- •Оценка погрешности интерполяционного многочлена Эрмита
- •Теорема 2. Пусть функция имеет на отрезкепроизводные до-го порядка включительно. Тогда для любогох из отрезка найдется точкатакая, что
- •4.4. Обратное интерполирование. Кусочно - многочленная интерполяция. Интерполирование сплайнами Обратное интерполирование
- •Кусочно-многочленная интерполяция
- •Интерполирование сплайнами Сплайны. Пусть на отрезке задана сетка точек, которая делитнаn отрезков ,. Сплайномn-го порядка называется функция
- •4.5. Многочлены Чебышева. Чебышевские узлы интерполяции Многочлены Чебышева
- •Контрольные вопросы и задания
4.5. Многочлены Чебышева. Чебышевские узлы интерполяции Многочлены Чебышева
Рассмотрим множество функций
,
,
(4.5.1)
определенных
на отрезке
и имеющих множества значений, принадлежащих
отрезку
.
По формуле (4.5.1) непосредственно получим:
,
.
(4.5.2)
Покажем,
что
совпадает с алгебраическим многочленомn-й
степени при любом натуральном значении
n.
Обозначим
.
Тогда
.
(4.5.3)
По известной тригонометрической формуле
.
Учитывая равенство (4.5.3), последнюю формулу можно записать в виде
.
Выразим
из последней формулы
и, учитывая равенства (4.5.2), получим
рекуррентную формулу
.
(4.5.4)
Вычисляя
по этой рекуррентной формуле функции
последовательно при
получим
,
,
,
,
…
Легко
видеть, что, продолжая эти вычисления,
мы на каждом шаге будем получать
алгебраические многочлены, причем
степени их каждый раз будет увеличиваться
на 1. Таким образом, функция
совпадает в области определения с
алгебраическим многочленомn-й
степени при любом натуральном значении
n.
Функции
получили названиемногочленов
Чебышева.
Рассмотрим некоторые свойства многочленов
Чебышева.
Найдем корни многочленов Чебышева. Для этого решим уравнение
при
.
Из него непосредственно получим
.
(4.5.5)
Здесь
величина k
может принимать любые целые значения.
Но уравнение (4.5.5) будет иметь решения
относительно x
только при тех значениях k,
при которых значение величины
будет заключено между 0 и
.
Поэтому уравнение (4.5.5) будет иметь
решения относительноx
только при тех значениях k,
при которых
.
(4.5.6)
Условия
(4.5.6) будут выполняться при
(
).
Например, при
величинаk
может принять только одно значение
.
При
.
При
и так далее. Зафиксируем некоторое
натуральное значениеn.
Для каждого из значений
уравнение (4.5.5) будет иметь единственное
решение:
,
.
(4.5.7)
Таким
образом, функция
(
)
будет иметьn
корней на отрезке
и все эти корни получаются по формулам
(4.5.7).
Найдем
коэффициенты при старшей степени
многочленов Чебышева. Рассматривая
формулы для
при
легко заметить, что при
единственный коэффициент многочлена
равен 1, а при
коэффициенты при старшей степени
многочленов
равны
.
Методом математической индукции с
помощью формулы (4.5.4) легко доказать,
что коэффициенты при старшей степени
многочленов
равны
при любых натуральных значенияхn.
Наряду
с введенными многочленами Чебышева,
часто используются и нормированные
многочлены Чебышева, получаемые из
путем деления их на коэффициенты при
старшей степени:
,
=
,
(4.5.8)
многочлены с коэффициентами при старшей степени, равными 1.
Функции
(
),
так же как и
,
будут иметьn
корней на отрезке
,
и все эти корни получаются по формулам
(4.5.7). Кроме того, они обладают одним
замечательным свойством, которое мы
приведем без доказательства.
Теорема 1 (теорема
Чебышева).
Из всех многочленов степени n+1
с коэффициентами при старшей степени,
равными 1 (многочленов вида
),
нормированный многочлен Чебышева
наименее уклоняется от нуля на отрезке
.
То есть
принимает наименьшее значение, если
.
Наименьшее значение
равно 1 при
,
а при всех остальных натуральных
значенияхn
наименьшее значение
будет равно
.
Чебышевские узлы интерполяции
В
параграфе 4.1 была доказана оценка
погрешности многочленной интерполяции.
Если интерполируемая функция
имеет на отрезке
ограниченную производнуюn+1-го
порядка и существует положительная
постоянная
такая, что
на отрезке
для любогоx,
принадлежащего отрезку
будет справедлива оценка погрешности
интерполяционного многочлена
(4.1.15)
.
В
этой оценке только величина
зависит отx.
Величина
неотрицательна. Поэтому наименьшее
значение погрешность интерполяции на
отрезке
будет принимать при тех же значениях
,
при которых принимает наименьшее
значение величина
.
Вычислим эти значения. Введем для этого
линейное преобразование
,
(4.5.9)
которое
отображает отрезок
на отрезок
.
Обратное преобразование, отображающее
на
,
очевидно, имеет вид
.
(4.5.10)
Введем точки
,
(4.5.11)
соответствующие
точкам
.
Поскольку точки
,
соответствующие им точки
.
Зная точки
,
можно получить соответствующие точки
по формуле
,
.
(4.5.12)
Сделаем
в выражении
замену переменных по формулам (4.5.9),
(4.5.12):
.
(4.5.13)
Согласно
теореме Чебышева,
достигается, если
совпадают с корнями многочленов Чебышева
:
,
.
(4.5.14)
Чтобы
получить эту формулу, необходимо в
формуле (4.5.7) заменить x
на t,
а n
на n+1.
Причем этот минимум равен
.Таким
образом, если
выбрать узлы интерполяции
на
отрезке
по
формулам
(4.5.14), (4.5.12),
то оценка
погрешности интерполяции
примет
наименьшее значение, равное
.
Эти узлы получили названиечебышевских
узлов интерполяции.
Если
выбраны чебышевские узлы интерполяции
,
то будет справедлива оценка погрешности
интерполяции
:
.
(4.5.15)
Эту
оценку называют наилучшей
равномерной оценкой погрешности
интерполяции.
Ее можно использовать для того, чтобы
ответить на вопрос об условиях сходимости
интерполяционных многочленов Лагранжак
при
.
Пусть функция
имеет производные любого порядка на
отрезке
и существует постояннаяМ
такая, что для
:
(здесь
мажорантные оценки модуля производной
на отрезке
).
Тогда
.
А
поскольку длявыполняется неравенство
,
то согласно теореме о промежуточной последовательности,
.
Таким
образом, если функция
имеет производные любого порядка на
отрезке
и существует постояннаяМ
такая, что
:
,
а для интерполяции на отрезке
выбираются чебышевские узлы, то
погрешность интерполяции будет стремиться
к 0 при
.