Интерполяционный многочлен Лагранжа

Пусть имеется Интерполяционный многочлен Лагранжа - student2.ru значение ф-ции Интерполяционный многочлен Лагранжа - student2.ru в Интерполяционный многочлен Лагранжа - student2.ru различных точках Интерполяционный многочлен Лагранжа - student2.ru

Интерполяционный многочлен Лагранжа - student2.ru

Надо найти функцию, которая смогла бы дать значение Интерполяционный многочлен Лагранжа - student2.ru в произвольной точке Интерполяционный многочлен Лагранжа - student2.ru . Для этого часто используют алгебраический многочлен Интерполяционный многочлен Лагранжа - student2.ru степени Интерполяционный многочлен Лагранжа - student2.ru , который в точках Интерполяционный многочлен Лагранжа - student2.ru принимает заданные значения, т.е. Интерполяционный многочлен Лагранжа - student2.ru (6.15)

Такой многочлен называетсяинтерполяционным. Точки Интерполяционный многочлен Лагранжа - student2.ru называются узлами интерполяции.

Теорема существования. Существует единственный интерполяционный многочлен Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru степени, удовлетворяющий условиям (6.15).

Доказательство. Существование подобного многочлена устанавливается

непосредственно путем его выписывания. Пусть Интерполяционный многочлен Лагранжа - student2.ru . Тогда

Интерполяционный многочлен Лагранжа - student2.ru . (6.16)

При Интерполяционный многочлен Лагранжа - student2.ru

Интерполяционный многочлен Лагранжа - student2.ru . (6.17)

В общем случае, при любом натуральном Интерполяционный многочлен Лагранжа - student2.ru

Интерполяционный многочлен Лагранжа - student2.ru . (6.18)

Непосредственной подстановкой в (6.16 - 6.18) убеждаемся, что Интерполяционный многочлен Лагранжа - student2.ru .

Нетрудно увидеть, что многочлен Интерполяционный многочлен Лагранжа - student2.ru в (6.18) имеет степень Интерполяционный многочлен Лагранжа - student2.ru . Видно, что

Интерполяционный многочлен Лагранжа - student2.ru , a Интерполяционный многочлен Лагранжа - student2.ru .

Единственность интерполяционного многочлена (6.18) доказывается методом от противного. Пусть кроме Интерполяционный многочлен Лагранжа - student2.ru имеется ещё некоторый алгебраический многочлен Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru степени, удовлетворяющий условиям

Интерполяционный многочлен Лагранжа - student2.ru , Интерполяционный многочлен Лагранжа - student2.ru (6.19)

Тогда, согласно (6.15) и (6.19)

Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru (6.20) Интерполяционный многочлен Лагранжа - student2.ru

Если в точках Интерполяционный многочлен Лагранжа - student2.ru разность (6.20) не равна нулю, то мы получим многочлен

Интерполяционный многочлен Лагранжа - student2.ru степени не выше Интерполяционный многочлен Лагранжа - student2.ru и в силу основной теоремы высшей алгебры он имеет не более Интерполяционный многочлен Лагранжа - student2.ru корней. Но из (6.20) следует, что мы уже имеем Интерполяционный многочлен Лагранжа - student2.ru корень. Чтобы избежать противоречия, мы должны согласиться, что Интерполяционный многочлен Лагранжа - student2.ru .

Это и есть доказательство единственности интерполяционного многочлена.

Интерполяционный многочлен (6.18) называется интерполяционным многочленом Лагранжа, а функции Интерполяционный многочлен Лагранжа - student2.ru -лагранжевыми коэффициентами.

Для вычисления Интерполяционный многочлен Лагранжа - student2.ru удобно применять следующее расположение разностей, подчеркнув разности, расположенные на главной диагонали:

Интерполяционный многочлен Лагранжа - student2.ru

Интерполяционный многочлен Лагранжа - student2.ru

Обозначим произведение элементов Интерполяционный многочлен Лагранжа - student2.ru строки через Di , а произведение элементов главной диагонали через Интерполяционный многочлен Лагранжа - student2.ru , т.е Интерполяционный многочлен Лагранжа - student2.ru . Тогда

Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru .

Иногда бывает полезным для упрощения вычислений использовать инвариантность коэффициентов Лагранжа относительно линейной подстановки: если

Интерполяционный многочлен Лагранжа - student2.ru , то Интерполяционный многочлен Лагранжа - student2.ru .

В случае равностоящих узлов имеются таблицы для лагранжевых коэффициентов и процесс вычисления значительно облегчается.

Схема Эйткена.

Если требуется найти не общее выражение Ln(x), а лишь его значения при конкретных х, и при этом значения функции даны в достаточно большом количестве узлов, то удобно пользоваться интерполяционной схемой Эйткена.

Согласно этой схеме последовательно вычисляются многочлены

Интерполяционный многочлен Лагранжа - student2.ru ,

Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru ,

----------------------------------------------------------------------------------------

­­­­­­­­­­­­­­­­­­­­­­­­­­­

Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru и т.д.

Интерполяционный многочлен Интерполяционный многочлен Лагранжа - student2.ru степени (принимающий в точках Интерполяционный многочлен Лагранжа - student2.ru значения Интерполяционный многочлен Лагранжа - student2.ru ) записывается следующим образом:

Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru .

Вычисления по схеме Эйткена удобно располагать в такой таблице

Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru
Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru        
Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru      
Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru    
Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru  
Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru Интерполяционный многочлен Лагранжа - student2.ru

Вычисления по схеме Эйткена обычно ведут до тех пор, пока последовательные значения Интерполяционный многочлен Лагранжа - student2.ru и Интерполяционный многочлен Лагранжа - student2.ru не совпадут в пределах заданной точности. Схема Эйткена легко реализуется на ЭВМ и обеспечивает возможность автоматического контроля точности вычислений.

Замечания.

1. Непосредственной проверкой можно убедиться, что Интерполяционный многочлен Лагранжа - student2.ru .

Это тождество может служить контролем при вычислении лагранжевых коэффициентов Интерполяционный многочлен Лагранжа - student2.ru .

2. Интерполяционный многочлен (6.18) линейно зависит от значений функции Интерполяционный многочлен Лагранжа - student2.ru .

Поэтому интерполяционный многочлен для суммы двух функций равен сумме интерполяционных многочленов для слагаемых.

Пример. Построить интерполяционный многочлен Лагранжа по табличным данным:

Интерполяционный многочлен Лагранжа - student2.ru
Интерполяционный многочлен Лагранжа - student2.ru
Интерполяционный многочлен Лагранжа - student2.ru

Решение. Согласно (6.18) при Интерполяционный многочлен Лагранжа - student2.ru имеем

Интерполяционный многочлен Лагранжа - student2.ru

Интерполяционный многочлен Лагранжа - student2.ru .

Погрешность интерполяции.

Можно написать следующее равенство.

Интерполяционный многочлен Лагранжа - student2.ru ,

где Интерполяционный многочлен Лагранжа - student2.ru - остаточный член или погрешность интерполяции.

Для оценки остаточного члена предположим, что

Интерполяционный многочлен Лагранжа - student2.ru ,

где Интерполяционный многочлен Лагранжа - student2.ru - отрезок, содержащий все узлы интерполяции Интерполяционный многочлен Лагранжа - student2.ru .

Будем искать Интерполяционный многочлен Лагранжа - student2.ru в виде:

Интерполяционный многочлен Лагранжа - student2.ru , (6.22)

где Интерполяционный многочлен Лагранжа - student2.ru , (6.23)

rn(x) - некоторая функция от Интерполяционный многочлен Лагранжа - student2.ru .

Воспользуемся следующим приемом:

Введем в рассмотрение фун-цию Интерполяционный многочлен Лагранжа - student2.ru .

Интерполяционный многочлен Лагранжа - student2.ru (6.24)

оставив в составе Интерполяционный многочлен Лагранжа - student2.ru функцию rn(x) при некотором произвольном фиксированном Интерполяционный многочлен Лагранжа - student2.ru но Интерполяционный многочлен Лагранжа - student2.ru .

Ф-ция Интерполяционный многочлен Лагранжа - student2.ru обращается в нуль при Интерполяционный многочлен Лагранжа - student2.ru , Интерполяционный многочлен Лагранжа - student2.ru , и Интерполяционный многочлен Лагранжа - student2.ru (на основании формул (6.15),(6.21), (6.22), (6.23) ).

Т.е. она принимает нулевые значения не менее чем в Интерполяционный многочлен Лагранжа - student2.ru точках отрезка Интерполяционный многочлен Лагранжа - student2.ru , на котором изменяется Интерполяционный многочлен Лагранжа - student2.ru . По теореме Ролля первая производная по t Интерполяционный многочлен Лагранжа - student2.ru обращается в нуль, по крайней мере в Интерполяционный многочлен Лагранжа - student2.ru точке интервала Интерполяционный многочлен Лагранжа - student2.ru , Интерполяционный многочлен Лагранжа - student2.ru равна нулю минимум в Интерполяционный многочлен Лагранжа - student2.ru точках этого интервала и т.д. Продолжая эти рассуждения и дойдя до Интерполяционный многочлен Лагранжа - student2.ru производной можно сказать, что найдется хотя бы одна точка Интерполяционный многочлен Лагранжа - student2.ru , в которой Интерполяционный многочлен Лагранжа - student2.ru .

Перепишем теперь последнее уравнение, предварительно продифференцировав Интерполяционный многочлен Лагранжа - student2.ru по Интерполяционный многочлен Лагранжа - student2.ru его левую и правую части:

Интерполяционный многочлен Лагранжа - student2.ru .

Так как Интерполяционный многочлен Лагранжа - student2.ru , а Интерполяционный многочлен Лагранжа - student2.ru в любой точке Интерполяционный многочлен Лагранжа - student2.ru , как производная многочлена Интерполяционный многочлен Лагранжа - student2.ru степени, то для точки Интерполяционный многочлен Лагранжа - student2.ru можно написать равенство

Интерполяционный многочлен Лагранжа - student2.ru ,

где Интерполяционный многочлен Лагранжа - student2.ru - производная Интерполяционный многочлен Лагранжа - student2.ru - го порядка от Интерполяционный многочлен Лагранжа - student2.ru .

Следовательно,

Интерполяционный многочлен Лагранжа - student2.ru (6.25)

или Интерполяционный многочлен Лагранжа - student2.ru , (6.26)

и Интерполяционный многочлен Лагранжа - student2.ru , (6.27)

где Интерполяционный многочлен Лагранжа - student2.ru - некоторая точка интервала Интерполяционный многочлен Лагранжа - student2.ru , положение которой зависит от рассматриваемого значения Интерполяционный многочлен Лагранжа - student2.ru

Из последнего равенства (6.27) можно оценить погрешность интерполяции в текущей точке Интерполяционный многочлен Лагранжа - student2.ru :

Интерполяционный многочлен Лагранжа - student2.ru . (6.28)

Здесь Интерполяционный многочлен Лагранжа - student2.ru .

Оценка максимальной погрешности интерполяции на всем отрезке Интерполяционный многочлен Лагранжа - student2.ru :

Интерполяционный многочлен Лагранжа - student2.ru . (6.29)

Пример. Оценить погрешность приближения функции Интерполяционный многочлен Лагранжа - student2.ru в точке Интерполяционный многочлен Лагранжа - student2.ru и на всем отрезке Интерполяционный многочлен Лагранжа - student2.ru , где Интерполяционный многочлен Лагранжа - student2.ru , с помощью интерполяционного многочлена Лагранжа второй степени Интерполяционный многочлен Лагранжа - student2.ru , построенного с узлами Интерполяционный многочлен Лагранжа - student2.ru .

Решение. Здесь Интерполяционный многочлен Лагранжа - student2.ru . Следовательно, нам потребуются производные до 3- го порядка включительно.

Интерполяционный многочлен Лагранжа - student2.ru .

Интерполяционный многочлен Лагранжа - student2.ru .

На основании (6.28)

Интерполяционный многочлен Лагранжа - student2.ru .

А в силу оценки (6.29):

Интерполяционный многочлен Лагранжа - student2.ru

Здесь опущено нахождение максимума Интерполяционный многочлен Лагранжа - student2.ru .

6.4. Линейная интерполяция.

Перепишем многочлен Лагранжа для случая Интерполяционный многочлен Лагранжа - student2.ru .

Интерполяционный многочлен Лагранжа - student2.ru (6.30)

Интерполяция в виде такого многочлена называется линейной.

 
  Интерполяционный многочлен Лагранжа - student2.ru

Рис.6.3. К вопросу о линейной интерполяции.

Введем обозначения Интерполяционный многочлен Лагранжа - student2.ru . Тогда Интерполяционный многочлен Лагранжа - student2.ru .

Интерполяционный многочлен Лагранжа - student2.ru .

Итак: Интерполяционный многочлен Лагранжа - student2.ru . (6.31).

Величина Интерполяционный многочлен Лагранжа - student2.ru называется фазой интерполяции, которая изменяется в пределах от 0 до 1, когда Интерполяционный многочлен Лагранжа - student2.ru пробегает значения от Интерполяционный многочлен Лагранжа - student2.ru до Интерполяционный многочлен Лагранжа - student2.ru .

Геометрически линейная интерполяция означает (рис. 6.3) замену графика функции на отрезке Интерполяционный многочлен Лагранжа - student2.ru хордой, соединяющей точки Интерполяционный многочлен Лагранжа - student2.ru .

Наши рекомендации