Читать онлайн «Maple 9.5/10 в математике, физике и образовании». Страница 63

Используя понятия разделенных разностей для полинома Ньютона можно получить выражение:

Nn(x) = f(x0) + (x-x0)f(x1, x0) + (x-x0)(x-x1)f(x0, x1, x2) + … + (x-x0)(x-x1)…(x-xn)f(x, x0, x1, …, xn)  (5.6)

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

Итерационно-интерполяционный метод Эйткена позволяет свести вычисления коэффициентов интерполяционного полинома Лагранжа, с учетом его равенства в узлах интерполяции с исходными данными к вычислению функциональных определителей второго порядка. При этом эффективность метода повышается в тех случаях, когда нет необходимости в получении приближенного аналитического выражения функции f(х), заданной таблично, а требуется лишь определить значение в некоторой точке х*, отличной от узловых точек. Этот метод заключается в последовательной линейной интерполяции. Процесс вычисления f(x*) состоит в следующем: необходимо пронумеровать узлы интерполяции, например, в порядке убывания их от х*. Затем для каждой узловой точки интерполяции строятся соотношения:

которые является интерполяционными полиномами, построенными соответственно по узлам хi, хj, хk. Продолжая этот процесс, имеем следующий полином:

   (5.7)

Полученный полином является интерполяционным полиномом, построенный по узлам хi, xj, …, хk, хm. Это утверждение верное, так как Рn-1ij…k(х) и Рn-1j…km(x) являются интерполяционными полиномами. При его реализации предполагается, что функция гладкая, а также критерием оценки погрешности определяется некоторое значение, определяемое условиями конкретной задачи.

Метод Чебышева был создан для оптимального выбора узлов интерполяции, если это возможно при решении конкретной задачи, и для получения минимально возможной погрешности аппроксимации. Предполагается, что в выборе расположения узлов интерполяции ограничений нет, и предполагается, что узлы выбираются произвольно. Ставится задача о наилучшем выборе узлов. Наилучшими узлами х, следует признать те, для которых выражение max[a,b]|ωn(x)| минимально для рассматриваемого класса функций (алгебраических полиномов). Определение этих узлов сводится к нахождению корней полинома, наименее уклоняющихся от нуля на [a, b]. Такой полином порождается полиномом Чебышева первого рода Tn+1.

Полиномы Чебышева определены в интервале [-1,1]. Для перевода интерполяции в интервале [a, b], выполняется линейная замена переменной х:

В качестве узлов интерполяции берутся корни полинома Чебышева:

   (5.8)

Тогда погрешность Чебышевской интерполяции определяется выражением:

   (5.9)

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

Наиболее известным и широко применяемым является случай сплайновой интерполяции, когда между двумя точками строится полином n-й степени

   (5.10)

который в узлах интерполяции принимает значения интерполируемой функции и непрерывен вместе со своими (n-1)-ми производными. Такой кусочно-непрерывный интерполяционный полином называется сплайном. Его коэффициенты находят из условий в узлах интерполяции — равенства значений сплайна и приближаемой функции, а также равенства (n-1)-й производной соответствующих полиномов. Максимальная по всем частичным отрезкам степень полинома является степенью сплайна.

Одним из наиболее распространенных интерполяционных сплайнов является кубический интерполяционный сплайн. Для вывода уравнения кубического интерполяционного сплайна можно воспользоваться его представлением в виде гибкой линейки, изогнутой таким образом, что она проходит через значения функции в узлах, то есть, является упругой рейкой в состоянии равновесия. Это его состояние описывается уравнением S'''(х)=0, где S'''(х) — четвертая производная. Из этого следует, что между каждой парой соседних узлов интерполяционная формула записывается в виде полинома третьей степени. Этот полином удобно представить следующим образом:

S(x) = аi + bi(x-xi-1) + с(х-хi-1)² + di(x–хi-1)³, xi-1≤х≤xi, i = 1, 2, ..., n.

Система Maple позволяет легко вычислять коэффициенты кубических полиномов. Метод сплайновой интерполяции дает хорошие результаты при интерполяции непрерывных функций с гладкими производными 1-ой и 2-ой степени. При этом кубическая сплайновая интерполяция, построенная по узлам fi=f(хi), i=0,1,…,n, будет иметь минимум кривизны по сравнению с любой интерполяционной функцией, имеющей непрерывные первую и вторую производные. Выполнение сплайн-интерполяции функций с резким изменением производных дает, как правило, большие ошибки. Сплайны более высоких порядков, чем третий, используется редко, так как при вычислении большого числа коэффициентов может накапливаться ошибка, приводящая к значительным погрешностям.

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

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

Итак, при задании f(х1), …, f(хn) приближение к f(x) ищется в виде

   (5.11)

Коэффициенты аi, bi находятся из совокупности соотношений R(хj)=f(xj) (j=1,…,n), которые можно записать в виде

Данное уравнение образует систему n линейных уравнений относительно n+1 неизвестных. Такая система всегда имеет нетривиальное решение.

Функция R(x) может быть записана в явном виде в случае n нечетное, если р=q, и n четное, если р-q=1. Для записи функции R(x) в явном виде следует вычислять так называемые обратные разделенные разности, определяемые условиями

и рекуррентным соотношением

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

Использование рациональной интерполяции часто целесообразнее интерполяции полиномами в случае функций с резкими изменениями характера поведения или особенностями производных в точках.