Учебное пособие: Численные методы
Строится система вложенных отрезков
¦( n +1) -производная (n+1)-го порядка
Пусть
(6.6)
Если ¦(c)-полином n-ой степени, то производная (n+1)-го порядка равна 0, тогда Rn (x)≡0 и мы получаем точную аппроксимацию.
Теорема:
Многочлен Лагранжа вида (6.4) для таблично заданной функции единственен.
Доказательство:
Пусть Qn (x)- многочлен Лагранжа, построенный для этой же функции ¦(c) по тем же узлам интерполяции. Qn (x)¹Pn (x) Qn (xi )=yi =Pn (xi ),
Рассмотрим многочлен Ln (x)= Qn (x)-Rn (x)-это многочлен n-ой степени, для которого точки xi , i=0,n являются корнями. Это противоречит основной теореме алгебры, которая говорит о том, что полином n-ой степени имеет ровно n корней . А Ln (x) имеет n+1 корней . Противоречие доказывает теорему.
Интерполяционная схема Эйткина
Поскольку при большом числе узлов интерполяции вычисление значения полинома Лагранжа по формуле (6.4) громоздко, необходимо получить рекуррентную формулу.
Пусть ¦(c)- непрерывна, узлы выбраны на отрезке [a;b] таким образом, что:
Введем функцию
xi -узлы интерполяции;
yi= ¦(c)
Полином Лагранжа: Pn (x) см. (6.4)
Таким образом, функция Q0,1 (x) представляет собой полином Лагранжа l-ой степени, построенной по узлам x0 ,x1 введем функцию вида
Функция Q1,2 (x)- интерполяционный полином Лагранжа, построенный по узлам x1 ,x2 .
Введем теперь функцию
Аналогично:
Q0,1,2 (x2 )= у2
В силу единственности полинома Лагранжа, построенного по узлам x0 , x1 ,x2
функция Q0,1,2 (x) представляет собой интерполяционный полином Лагранжа 2-ой степени, построенный по узлам x0 , x1 ,x2 .
Введем функцию: