
- •Интерполяция
- •1. Постановка задачи и основные определения
- •2. Интерполяция многочленами
- •2.1 Интерполяционная формула Лагранжа
- •2.2 Погрешность интерполяционного многочлена
- •2.3 О сходимости процесса интерполяции
- •3. Интерполяция сплайнами
- •3.1 Что такое сплайн
- •3.2 Построение кубического сплайна
- •3.3 Сплайн-квадратура
- •Приложение. Текст библиотечного модуля splinunt на языке Паскаль для интерполяции кубическими сплайнами.
2.2 Погрешность интерполяционного многочлена
2.3 О сходимости процесса интерполяции
3. Интерполяция сплайнами
Сплайны – сравнительно недавнее изобретение, получившее широкое распространение во многих областях вычислительной математики, начиная с 60-х годов [1]. Они обладают целым рядом замечательных свойств. В частности, если говорить об интерполяции, сплайны (в отличие от многочленов) гарантируют сходимость при равномерном распределении узлов. Кроме того, естественный кубический сплайн является самой гладкой среди всех функций, интерполирующих заданный набор точек.
Более подробно с теорией и применениями сплайнов, помимо классической монографии [Error: Reference source not found], можно ознакомиться по книгам [2, 3].
3.1 Что такое сплайн
Сплайн (от англ. spline) или сплайн-функция – это непрерывная кусочно-полиномиальная функция, определенная на заданном отрезке и имеющая на этом отрезке несколько непрерывных производных. Простейшим примером сплайна может служить ломаная, проведенная через заданные точки (понятно, что в этом случае можно говорить только о непрерывности функции, но не производных). Однако, наибольший интерес представляют сплайны, построенные на основе многочленов более высоких степеней.
Вероятно, наиболее известными являются кубические сплайны. Собственно говоря, возникновение самого термина «сплайн» связано именно с кубическими сплайнами. Дело в том, что английское слово spline первоначально обозначало чертежный инструмент для проведения плавных кривых через заданные точки. Этот инструмент представляет собой упругую линейку (рейку) с приспособлением для закрепления в нужных точках. Решение соответствующей механической задачи показывает, что форма, которую приобретает такая рейка, в промежутках между точками фиксации описывается кубическими многочленами, причем на всем протяжении кривой первая и вторая производные сохраняют непрерывность. В крайних точках вторая производная, характеризующая кривизну, принимает нулевое значение, поскольку концы рейки, закрепленные только с одной стороны, распрямляются. Математическую функцию, воспроизводящую форму рейки, называют естественным кубическим сплайном.
Рассмотренная выше механическая аналогия позволяет легко объяснить гладкость сплайн-функции. Очевидно, что рейка, закрепленная в отдельных точках, примет ту форму, которой отвечает минимум потенциальной энергии. Из механики известно, что потенциальная энергия упругой деформации элемента рейки, связанной с его изгибом, пропорциональна квадрату кривизны. Для рейки в целом потенциальная энергия равна интегралу от квадрата кривизны. Таким образом, ее форма отвечает минимальной кривизне (в интегральном смысле).