Численные методы математического анализа
1. Задача интерполяции. Многочлен Лагранжа
1.1. Постановка задачи
Пусть задана функция . Часто нахождение
значений этой функции может оказаться трудоемкой задачей. Например,
- параметр в
некоторой сложной задаче, после решения которой определяется значение f(x), или
f(x) измеряется в дорогостоящем эксперименте. В этих случаях можно получить
небольшую таблицу значений функции, но прямое нахождение ее значений при большом
количестве значений аргумента нереально. В такой ситуации f(x) заменяется
приближенной формулой
, которая в
определенном смысле близка к функции f(x). Близость обеспечивается введением в
функцию
свободных параметров и их
соответствующим выбором.
Итак, известны значения функции f(x) в точках
,
. Потребуем, чтобы для некоторой функции
выполнялись равенства:
(1)
Если (1) рассматривать как систему для определения , то этот
способ называется интерполяцией (Лагранжевой).
Если зависит от
нелинейно, то
интерполяция нелинейная, иначе интерполяция линейная. В случае линейной
интерполяции можно записать
(2)
(где - система
линейно-независимых функций)
Подставим (2) в (1). Относительно получаем линейную
систему уравнений:
,
(3)
Для однозначной разрешимости системы должно быть .
Для того, чтобы задача нтерполирования имела единственное
решение, система функций должна для любых
несовпадающих
удовлетворять условию:
(4)
Система функций, удовлетворяющая условию (4), называется чебышевской.
Наиболее простой и (для многих случаев) удобной является
система функций ,
. Функция
при этом представляет
собой многочлен степени
(интерполяционный
многочлен) с коэффициентами
.
Система уравнений (3) в этом случае имеет вид:
,
(5)
Определителем этой системы является отличный от нуля
определитель Вандермонда:
Отсюда следует, что интерполяционный многочлен существует и единственен.
Непосредственное решение системы (5) для нахождения
aj уже при небольших n приводит к сильному
искажению значений aj. Получим явный вид
нтерполяционного многочлена , не решая систему
(5).
Если многочлен степени n
такой что
, то
- искомый интерполяционный многочлен степени
, т.к.
.
Так как при
, то
делится на
для любых
, то
есть
. Так как
, то
.
Таким образом, (6)
Такая форма записи интерполяционного многочлена называется
многочленом Лагранжа и обозначается, как правило, .
Существуют и другие формы записи того же самого интерполяционного многочлена.
Если обозначить , то
.
Тогда (6) можно записать в виде
В узлах многочлен Лагранжа
совпадает с заданной функцией, в остальных точках в общем случае
не совпадает с
(кроме
случая, когда
многочлен степени не выше
). Разность
-
- остаточный член. Запишем ее в виде
.
При . Найдем постоянную
такую, чтобы
в некоторой фиксированной точке
, в которой мы
рассматриваем погрешность.
Относительно будем предполагать
кратную
дифференцируемость.
Значение c, при котором существует и равно
. Тогда
функция
равна нулю по крайней мере в
точках
.
По теореме Ролля производная равна нулю по крайней
мере в
точках
. Далее,
равна нулю по крайней мере в n точках и т.д.
Для - производной получаем,
что существует по крайней мере одна точка
такая, что
.
Отсюда получаем при . В этом случае в
точке
, т.е. остаточный член в точке
имеет вид:
Значение зависит от
- точки, в
которой рассматривается погрешность.
,
где