Источник: http://bspu.secna.ru/~pvv/mathpage/calcmath/complex/lection_a.htm - WWW - версия учебно - методического комплекса по численным методам разработана на кафедре Информатики и Математического Обеспечения Петрозаводского государственного университета

Численные методы математического анализа

1. Задача интерполяции. Многочлен Лагранжа

1.1. Постановка задачи

Пусть задана функция . Часто нахождение значений этой функции может оказаться трудоемкой задачей. Например, - параметр в некоторой сложной задаче, после решения которой определяется значение f(x), или f(x) измеряется в дорогостоящем эксперименте. В этих случаях можно получить небольшую таблицу значений функции, но прямое нахождение ее значений при большом количестве значений аргумента нереально. В такой ситуации f(x) заменяется приближенной формулой , которая в определенном смысле близка к функции f(x). Близость обеспечивается введением в функцию свободных параметров и их соответствующим выбором.

Итак, известны значения функции f(x) в точках

, . Потребуем, чтобы для некоторой функции

выполнялись равенства:

(1)

Если (1) рассматривать как систему для определения , то этот способ называется интерполяцией (Лагранжевой).

Если зависит от нелинейно, то интерполяция нелинейная, иначе интерполяция линейная. В случае линейной интерполяции можно записать

(2)

(где - система линейно-независимых функций)

Подставим (2) в (1). Относительно получаем линейную систему уравнений:

, (3)

Для однозначной разрешимости системы должно быть .

Для того, чтобы задача нтерполирования имела единственное решение, система функций должна для любых несовпадающих удовлетворять условию:

(4)

Система функций, удовлетворяющая условию (4), называется чебышевской.

1.2. Построение интерполяционного многочлена Лагранжа

Наиболее простой и (для многих случаев) удобной является система функций , . Функция при этом представляет собой многочлен степени (интерполяционный многочлен) с коэффициентами .

Система уравнений (3) в этом случае имеет вид:

, (5)

Определителем этой системы является отличный от нуля определитель Вандермонда:

Отсюда следует, что интерполяционный многочлен существует и единственен.

Непосредственное решение системы (5) для нахождения aj уже при небольших n приводит к сильному искажению значений aj. Получим явный вид нтерполяционного многочлена , не решая систему (5).

Если многочлен степени n такой что , то - искомый интерполяционный многочлен степени , т.к. .

Так как при , то делится на для любых , то есть

. Так как , то .

Таким образом,
(6)

Такая форма записи интерполяционного многочлена называется многочленом Лагранжа и обозначается, как правило, .

Существуют и другие формы записи того же самого интерполяционного многочлена.

Если обозначить , то

.

Тогда (6) можно записать в виде

1.3. Остаточный член

В узлах многочлен Лагранжа совпадает с заданной функцией, в остальных точках в общем случае не совпадает с (кроме случая, когда многочлен степени не выше ). Разность - - остаточный член. Запишем ее в виде .

При . Найдем постоянную такую, чтобы в некоторой фиксированной точке , в которой мы рассматриваем погрешность.

Относительно будем предполагать кратную дифференцируемость.

Значение c, при котором существует и равно . Тогда функция равна нулю по крайней мере в точках .

По теореме Ролля производная равна нулю по крайней мере в точках . Далее, равна нулю по крайней мере в n точках и т.д.

Для - производной получаем, что существует по крайней мере одна точка такая, что .

Отсюда получаем при . В этом случае в точке , т.е. остаточный член в точке имеет вид:

Значение зависит от - точки, в которой рассматривается погрешность.

Оценка остаточного члена:

,

где