Сплайн

Сплайн — это функция, область определения которой разбита на конечное число отрезков, на каждом из которых функция совпадает с некоторым алгебраическим многочленом. Сплайны относятся к классу кусочно-полиномиальных функций и используются для аппроксимации, интерполяции и сглаживания данных. Основное свойство сплайна — обеспечение высокой степени гладкости в точках соединения (узлах) полиномиальных сегментов, что позволяет получать кривые и поверхности, лишённые резких изломов и осцилляций, характерных для интерполяции одним полиномом высокой степени. Термин происходит от английского spline — гибкое лекало, использовавшееся в черчении для проведения плавных кривых через заданные точки.

История

Прообразом математического сплайна послужило физическое устройство — гибкая рейка (сплайн), которую чертёжники фиксировали в опорных точках для вычерчивания плавных линий. В механике форма такой рейки описывается дифференциальным уравнением изгиба балки, что даёт кубический полином между точками опоры.

Математическая теория сплайнов начала активно развиваться в середине XX века. В 1946 году американский математик Исаак Шёнберг опубликовал работу, в которой впервые ввёл понятие сплайна как математического объекта и заложил основы теории. В 1960-е годы, с развитием вычислительной техники, сплайны стали широко применяться в задачах численного анализа, компьютерной графики и автоматизированного проектирования (САПР). Ключевой вклад в развитие теории и практического применения сплайнов внесли Карл де Бур (алгоритмы вычисления B-сплайнов), Пьер Безье (кривые Безье) и другие исследователи.

Основные свойства и классификация

Сплайны классифицируются по нескольким признакам: степени полиномов, условиям гладкости в узлах и типу базиса.

По степени полинома

• Линейный сплайн (степень 1): на каждом отрезке функция является линейной. Обеспечивает непрерывность самой функции, но не её первой производной (график — ломаная линия).
• Квадратичный сплайн (степень 2): на каждом отрезке — квадратичный полином. Обеспечивает непрерывность функции и её первой производной.
• Кубический сплайн (степень 3): наиболее распространённый тип. На каждом отрезке — кубический полином. Обеспечивает непрерывность функции, первой и второй производных. Для построения кубического сплайна обычно требуются дополнительные граничные условия на концах интервала (например, «естественный сплайн» — вторая производная на концах равна нулю, или «зажатый сплайн» — заданы значения первой производной на концах).
• Сплайны высших степеней (4, 5 и выше): используются реже, так как при увеличении степени растёт вычислительная сложность и могут возникать нежелательные осцилляции.

По типу базиса

• Интерполяционные сплайны: проходят точно через все заданные точки данных (узлы интерполяции).
• Сглаживающие сплайны: не обязательно проходят через узлы, а минимизируют некоторый функционал, балансирующий между близостью к данным и гладкостью кривой. Используются для обработки зашумлённых данных.
• B-сплайны (базисные сплайны): каждая сегментная функция строится на основе локального базиса, что обеспечивает высокую вычислительную эффективность и устойчивость. B-сплайны лежат в основе многих алгоритмов компьютерной графики и САПР.
• Кривые Безье: частный случай сплайнов, задаваемый контрольными точками. Кривая не проходит через все контрольные точки, а лишь аппроксимирует их, при этом первая и последняя точки являются началом и концом кривой. Кривые Безье широко используются в векторной графике (шрифты TrueType, PostScript, Adobe Illustrator) и анимации.
• NURBS (Non-Uniform Rational B-Splines) — неоднородные рациональные B-сплайны. Являются обобщением B-сплайнов, позволяющим точно представлять конические сечения (окружности, эллипсы) и другие кривые второго порядка. NURBS — стандарт де-факто в современной САПР и промышленном дизайне.

Применение

Сплайны находят применение в самых разных областях науки, техники и искусства.

Компьютерная графика и анимация

• Моделирование кривых и поверхностей: сплайны используются для создания плавных контуров трёхмерных объектов, автомобильных кузовов, корпусов кораблей и самолётов, ландшафта.
• Анимация: кривые Безье и B-сплайны применяются для задания траекторий движения объектов, интерполяции ключевых кадров (скелетная анимация), управления параметрами эффектов.
• Векторные шрифты: контуры букв в шрифтах TrueType и OpenType описываются квадратичными и кубическими кривыми Безье.

Автоматизированное проектирование (САПР)

• Дизайн промышленных изделий: сплайны являются основным инструментом для построения гладких обводов в автомобилестроении, авиастроении, судостроении, дизайне бытовой техники.
• Обратный инжиниринг: по облаку точек, полученному с 3D-сканера, строится сплайновая поверхность, повторяющая форму реального объекта.

Численный анализ и обработка данных

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

Обработка изображений и сигналов

• Интерполяция изображений: бикубические сплайны используются для масштабирования изображений (увеличения или уменьшения разрешения) с сохранением плавности переходов.
• Восстановление сигналов: сплайны применяются для восстановления аналогового сигнала по его дискретным отсчётам.

Статистика и машинное обучение

• Регрессионные сплайны: используются для построения нелинейных регрессионных моделей, позволяя гибко аппроксимировать сложные зависимости между переменными.
• MARS (Multivariate Adaptive Regression Splines): метод множественной адаптивной регрессии, основанный на сплайнах, для построения прогностических моделей.

Примеры

Интерполяция кубическим сплайном

Пусть заданы точки (0,0), (1,1), (2,0), (3,1). Для построения интерполяционного кубического сплайна необходимо найти четыре кубических полинома (по одному на каждый отрезок между точками), которые будут непрерывны вместе со своими первыми и вторыми производными в узлах (1,2,3). Решение системы линейных уравнений даёт коэффициенты полиномов. Полученная кривая будет плавно проходить через все четыре точки.

Кривая Безье в векторной графике

В программе Adobe Illustrator или Inkscape инструмент «Перо» позволяет создавать кривые Безье. Пользователь задаёт опорные точки и управляющие рычаги (касательные), которые определяют форму кривой. Изменяя длину и направление рычагов, можно добиться практически любой формы линии.

Интересные факты

• Математические сплайны, используемые в компьютерной графике, часто называют «кривыми Безье» в честь французского инженера Пьера Безье, который разработал их для проектирования кузовов автомобилей Renault.
• В САПР и системах трёхмерного моделирования (например, Blender, Autodesk Maya, 3ds Max) сплайны являются одним из основных типов геометрических примитивов, наряду с полигональными сетками и NURBS-поверхностями.
• Способность сплайнов минимизировать кривизну (для кубических сплайнов — минимизация интеграла квадрата второй производной) делает их аналогом физической гибкой рейки, которая принимает форму с минимальной энергией изгиба.

Источники

• Шёнберг И. Дж. Contributions to the problem of approximation of equidistant data by analytic functions. — 1946.
• де Бур К. Практическое руководство по сплайнам. — М.: Радио и связь, 1985.
• Завьялов Ю. С., Квасов Б. И., Мирошниченко В. Л. Методы сплайн-функций. — М.: Наука, 1980.
• Роджерс Д., Адамс Дж. Математические основы машинной графики. — М.: Мир, 2001.
• Фарен Г. Кривые и поверхности для компьютерной графики. — М.: Бином, 2005.