Глобальная интерполяция.
В случае глобальной интерполяции отыскивается единый полином на всем интервале , то есть строится полином, который используется для интерполяции функции ) на всем интервале изменения аргумента . Будем искать интерполирующую функцию в виде полинома (многочлена) – ой степени:
Какова должна быть степень многочлена, чтобы удовлетворить всем условиям интерполяции? Допустим, что заданы две точки: и , то есть . Через эти точки можно провести единственную прямую, т.е. интерполирующей функцией будет полином первой степени
Через три точки можно провести параболу
и т.д. Рассуждая таким способом, можно предположить, что искомый полином должен иметь степень N .
Для того, чтобы доказать это, выпишем систему уравнений на коэффициенты. Уравнения системы представляют собой условия интерполяции при каждом :
Данная система является линейной относительно искомых коэффициентов . Известно, что СЛАУ имеет решение, если ее определитель отличен от нуля. Определитель данной системы
носит имя определителя Вандермонда. Из курса математического анализа известно, что он отличен от нуля, если (т.е. все узлы интерполяции различные). Таким образом, доказано, что система имеет решение. Мы показали, что для нахождения коэффициентов надо решить СЛАУ, что является сложной задачей. Но есть другой способ построения полинома N–й степени, который не требует решения такой системы.
Полином Лагранжа.
Решение ищем в виде:
где – базисные полиномы –й степени, для которых выполняется условие:
.
Убедимся в том, что если такие полиномы построены, то будет удовлетворять условиям интерполяции:
Для построения базисных полиномов определим
Легко понять, что
Функция является полиномом –й степени от и для нее выполняются условия "базисности":
Таким образом, нам удалось решить задачу о построении интерполирующего полинома N– й степени, и для этого не нужно решать СЛАУ. Полином Лагранжа можно записать в виде компактной формулы:
Погрешность этой формулы можно оценить, если исходная функция имеет производные до N+1 порядка:
Из этой формулы следует, что погрешность метода зависит от свойств функции , а также от расположения узлов интерполяции и точки z. Как показывают расчетные эксперименты, полином Лагранжа имеет малую погрешность при небольших значениях N<20. При бόльших N погрешность начинает расти, что свидетельствует о том, что метод Лагранжа не сходится (т.е. его погрешность не убывает с ростом N). Рассмотрим частные случаи. Пусть N=1, т.е. заданы значения функции только в двух точках. Тогда базовые полиномы имеют вид:
получаем формулы кусочно–линейной интерполяции. Пусть N=2.
В результате мы получили формулы так называемой квадратичной или параболической интерполяции.
Пример. Заданы значений некоторой функции:
x | 3,5 | |||
f | -1 | 0,2 | 0,5 | 0,8 |
Требуется найти значение функции при , используя интерполяционный полином Лагранжа. Для этого случая , т.е. полином Лагранжа имеет третий порядок. Вычислим значения базисных полиномов при :
Дата добавления: 2015-11-06; просмотров: 3278;