B-сплайн

B-сплайн (от англ. basis spline — базисный сплайн) — это математическая функция, представляющая собой кусочно-полиномиальную кривую минимальной степени гладкости, обладающую свойством локальности и образующая базис в пространстве сплайнов. B-сплайны являются обобщением кривых Безье и широко применяются в компьютерной графике, системах автоматизированного проектирования (САПР), численном анализе и обработке сигналов для моделирования и интерполяции сложных кривых и поверхностей.

Определение и математическая основа

B-сплайн степени n определяется как параметрическая кривая, заданная набором контрольных точек и вектором узлов. Формально, B-сплайн C(t) вычисляется как взвешенная сумма контрольных точек P_i с весами, определяемыми базисными функциями N_i,n(t):

C(t) = Σ_i=0^m P_i N_i,n(t*)

где t — параметр, изменяющийся в диапазоне от 0 до 1 или в пределах заданного интервала узлов.

Базисные функции

Базисные функции B-сплайна N_i,n(t) рекурсивно определяются по формуле Кокса — де Бура. Для степени n = 0 (кусочно-постоянная функция) они имеют вид:

N_i,0(t) = 1, если t_i ≤ t < t_i+1, и 0 в противном случае.

Для степени n > 0 используется рекурсия:

N_i,n(t) = ((t — t_i) / (t_i+n — t_i)) N_i,n-1(t) + ((t_i+n+1 — t) / (t_i+n+1 — t_i+1)) N_i+1,n-1(t)

где t_i — элементы вектора узлов.

Вектор узлов

Вектор узлов — это неубывающая последовательность значений параметра t, которая определяет, где и как соединяются полиномиальные сегменты сплайна. Вектор может быть:

• Равномерным: узлы расположены на одинаковом расстоянии друг от друга (например, [0, 1, 2, 3, 4]).
• Неравномерным: расстояния между узлами различны, что позволяет локально изменять форму кривой.
• Открытым (разомкнутым): первый и последний узлы повторяются n+1 раз (например, для кубического сплайна: [0, 0, 0, 0, 1, 2, 3, 4, 4, 4, 4]). Это гарантирует, что кривая начинается в первой контрольной точке и заканчивается в последней.
• Замкнутым (периодическим): узлы образуют циклическую последовательность, что позволяет создавать замкнутые кривые.

Классификация B-сплайнов

B-сплайны классифицируются по нескольким признакам:

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

• Линейные (степень 1): кусочно-линейные кривые, состоящие из отрезков прямых.
• Квадратичные (степень 2): кривые второго порядка, обеспечивающие непрерывность первой производной.
• Кубические (степень 3): наиболее распространённый тип, обеспечивающий непрерывность второй производной и оптимальный баланс между гладкостью и вычислительной сложностью.
• Сплайны высших степеней: используются редко из-за склонности к осцилляциям (феномен Рунге) и высокой вычислительной стоимости.

По типу контрольных точек

• Интерполирующие B-сплайны: кривая проходит через все заданные точки.
• Аппроксимирующие B-сплайны: кривая приближается к контрольным точкам, но не обязательно проходит через них (наиболее распространённый тип).

По структуре

• Рациональные B-сплайны (NURBS): кривые, заданные с весовыми коэффициентами для каждой контрольной точки. Позволяют точно представлять конические сечения (окружности, эллипсы) и другие кривые, не представимые обычными B-сплайнами.
• Однородные B-сплайны: все веса равны единице.

Свойства B-сплайнов

B-сплайны обладают рядом важных свойств, определяющих их применение:

• Локальность: изменение положения одной контрольной точки влияет только на часть кривой, ограниченную n+1 сегментами. Это позволяет выполнять локальную правку формы без перестроения всей кривой.
• Выпуклая оболочка: кривая целиком лежит внутри выпуклой оболочки своих контрольных точек. Это обеспечивает предсказуемость и стабильность формы.
• Непрерывность: B-сплайн степени n имеет непрерывные производные до порядка n-1 в узлах (C^n-1). Например, кубический сплайн имеет непрерывную вторую производную (C²).
• Инвариантность к аффинным преобразованиям: при параллельном переносе, повороте или масштабировании всех контрольных точек кривая преобразуется аналогично.
• Базисность: набор B-сплайнов заданной степени и с заданным вектором узлов образует базис в пространстве всех кусочно-полиномиальных функций с заданными условиями гладкости.

Применение

B-сплайны нашли широкое применение в различных областях науки и техники:

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

• Моделирование кривых и поверхностей: B-сплайны используются для создания гладких контуров объектов, траекторий движения и анимации персонажей.
• Шрифты: форматы шрифтов, такие как TrueType и PostScript, используют квадратичные и кубические B-сплайны для описания контуров глифов.
• Компьютерные игры: для создания плавных анимаций и траекторий движения объектов.

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

• Проектирование автомобилей, самолётов и кораблей: B-сплайны и NURBS являются стандартом для описания сложных аэродинамических поверхностей.
• Архитектура: создание нестандартных архитектурных форм и фасадов.
• Промышленный дизайн: моделирование корпусов бытовой техники, мебели и других изделий.

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

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

Медицина

• Компьютерная томография и МРТ: B-сплайны применяются для реконструкции трёхмерных моделей органов и тканей.
• Хирургическое планирование: моделирование разрезов и имплантатов.

История

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

Сравнение с кривыми Безье

B-сплайны являются обобщением кривых Безье. Основные отличия:

• Локальность: кривая Безье степени n определяется n+1 контрольными точками, и изменение одной точки влияет на всю кривую. B-сплайн обеспечивает локальность.
• Количество контрольных точек: для кривой Безье число точек строго фиксировано степенью, для B-сплайна — произвольно.
• Гладкость: B-сплайн может иметь более высокую степень гладкости в точках соединения сегментов, чем составная кривая Безье.

Примеры реализации

В компьютерной графике B-сплайны реализованы в большинстве графических API, включая OpenGL, DirectX и Vulkan. В математических пакетах, таких как MATLAB и GNU Octave, существуют встроенные функции для работы со сплайнами (например, spapi и spmak). В языках программирования Python библиотеки SciPy и NumPy предоставляют инструменты для создания и оценки B-сплайнов.

Критика и ограничения

Несмотря на широкое распространение, B-сплайны имеют некоторые ограничения:

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

Источники

• Де Бур, К. «Практическое руководство по сплайнам». — М.: Радио и связь, 1985.
• Роджерс, Д., Адамс, Дж. «Математические основы машинной графики». — М.: Мир, 2001.
• Piegl, L., Tiller, W. «The NURBS Book». — Springer, 1997.
• Farin, G. «Curves and Surfaces for CAGD: A Practical Guide». — Academic Press, 2002.
• Шёнберг, И. «Contributions to the problem of approximation of equidistant data by analytic functions». — Quarterly of Applied Mathematics, 1946.