Сплайн (математика) - Spline (mathematics)

Для инструмент для рисования, видеть Плоский шлиц.

Отдельные узлы в 1/3 и 2/3 образуют сплайн из трех кубических многочленов, пересекающихся с C² преемственность. Тройные узлы на обоих концах интервала гарантируют, что кривая интерполирует конечные точки.

В математика, а сплайн это особенный функция определенный кусочно к многочлены.В интерполирующий проблемы, сплайн-интерполяция часто предпочитают полиномиальная интерполяция потому что он дает аналогичные результаты, даже при использовании полиномов низкой степени, избегая при этом Феномен Рунге для более высоких степеней.

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

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

Вступление

Термин «сплайн» используется для обозначения широкого класса функций, которые используются в приложениях, требующих интерполяции и / или сглаживания данных. Данные могут быть одномерными или многомерными. Сплайн-функции для интерполяции обычно определяются как минимизаторы подходящих мер шероховатости (например, интегрального квадрата кривизны) с учетом ограничений интерполяции. Сглаживающие сплайны можно рассматривать как обобщения интерполяционных сплайнов, в которых функции определяются для минимизации взвешенной комбинации среднего квадрата ошибки аппроксимации по наблюдаемым данным и меры шероховатости. Для ряда содержательных определений меры шероховатости сплайн-функции оказываются конечномерными по своей природе, что является основной причиной их полезности в вычислениях и представлении. В оставшейся части этого раздела мы полностью сосредоточимся на одномерных полиномиальных сплайнах и используем термин «сплайн» в этом ограниченном смысле.

Определение

Начнем с того, что ограничимся обсуждением одномерный полиномиальный случай. В этом случае сплайн - это кусочно полиномиальная функция Эту функцию назовите S, принимает значения из интервала [а,б] и сопоставляет их с ${\displaystyle \mathbb {R} }$, набор действительные числа,

${\displaystyle S:[a,b]\to \mathbb {R} .}$

Мы хотим S быть кусочно определенным. Для этого пусть интервал [а,б] быть покрытым k упорядоченные подынтервалы с попарно непересекающимися интерьеры,

${\displaystyle [t_{i},t_{i+1}]{\mbox{ , }}i=0,\ldots ,k-1}$

${\displaystyle [a,b]=[t_{0},t_{1}]\cup [t_{1},t_{2}]\cup \cdots \cup [t_{k-2},t_{k-1}]\cup [t_{k-1},t_{k}]}$

${\displaystyle a=t_{0}\leq t_{1}\leq \cdots \leq t_{k-1}\leq t_{k}=b}$

По каждому из этих k "кусочки" [а,б], мы хотим определить многочлен, назовем его п_я.

${\displaystyle P_{i}:[t_{i},t_{i+1}]\to \mathbb {R} }$.

На я-й подынтервал [а,б], S определяется п_я,

${\displaystyle S(t)=P_{0}(t){\mbox{ , }}t_{0}\leq t<t_{1},}$

${\displaystyle S(t)=P_{1}(t){\mbox{ , }}t_{1}\leq t<t_{2},}$

${\displaystyle \vdots }$

${\displaystyle S(t)=P_{k-1}(t){\mbox{ , }}t_{k-1}\leq t\leq t_{k}.}$

Данный к + 1 точки т_я называются узлы. Вектор${\displaystyle {\mathbf {t}}=(t_{0},\dots ,t_{k})}$ называется узел вектор для сплайна. Если узлы равномерно распределены в интервале [а,б] мы говорим, что сплайн униформа, иначе мы говорим, что это неоднородный.

Если полиномиальные части п_я у каждого есть высшее образование п, то говорят, что сплайн имеет степень ${\displaystyle \leq n}$ (или изпорядок п + 1).

Если ${\displaystyle S\in C^{r_{i}}}$ в районе т_я, то говорят, что сплайн гладкость (по меньшей мере) ${\displaystyle C^{r_{i}}}$ в т_я. То есть при т_я две части п_я-1 и п_я разделяют значения общих производных от производной порядка 0 (значение функции) до производной порядка р_я (другими словами, две соседние части полинома соединяются с потеря гладкости не более п - р_я).

Вектор${\displaystyle {\mathbf {r}}=(r_{1},\dots ,r_{k-1})}$ так, чтобы сплайн имел гладкость ${\displaystyle C^{r_{i}}}$ в т_я за ${\displaystyle i=0,\ldots ,k-1}$ называется вектор гладкости для шлица.

Учитывая вектор узла ${\displaystyle {\mathbf {t}}}$, градус п, а вектор гладкости ${\displaystyle {\mathbf {r}}}$ за ${\displaystyle {\mathbf {t}}}$, можно рассматривать множество всех сплайнов степени ${\displaystyle \leq n}$ имеющий узловой вектор ${\displaystyle {\mathbf {t}}}$ и вектор плавности ${\displaystyle {\mathbf {r}}}$. Оснащенный операцией сложения двух функций (точечное сложение) и взятия действительных кратных функций, этот набор становится реальным векторным пространством. Этот сплайновое пространство обычно обозначается как ${\displaystyle S_{n}^{\mathbf {r}}({\mathbf {t}})}$.

При математическом изучении полиномиальных сплайнов вопрос о том, что происходит, когда два узла, скажем, т_я и т_я+1, перемещаются вместе, есть простой ответ. Полиномиальный кусокп_я(т) исчезает, и кускип_я−1(т) и п_я+1(т) объединяются с суммой потерь непрерывности длят_я и т_я+1.То есть,

${\displaystyle S(t)\in C^{n-j_{i}-j_{i+1}}[t_{i}=t_{i+1}],}$ куда ${\displaystyle j_{i}=n-r_{i}}$

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

${\displaystyle (t_{0},t_{1},\cdots ,t_{1},t_{2},\cdots ,t_{2},t_{3},\cdots ,t_{k-2},t_{k-1},\cdots ,t_{k-1},t_{k})}$

куда т_я повторяется j_я раз для ${\displaystyle i=1,\dots ,k-1}$.

А параметрическая кривая на интервале [а,б]

${\displaystyle G(t)=(X(t),Y(t)){\mbox{ , }}t\in [a,b]}$

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

Примеры

Предположим, что интервал [а,б] равно [0,3], а подынтервалы - [0,1], [1,2] и [2,3]. Предположим, что полиномиальные части должны иметь степень 2, а части на [0,1] и [1,2] должны соединяться по значению и первой производной (на т= 1), в то время как части на [1,2] и [2,3] соединяются просто по стоимости (при т = 2). Это определило бы тип сплайна S(т) для которого

${\displaystyle S(t)=P_{0}(t)=-1+4t-t^{2}{\mbox{ , }}0\leq t<1}$

${\displaystyle S(t)=P_{1}(t)=2t{\mbox{ , }}1\leq t<2}$

${\displaystyle S(t)=P_{2}(t)=2-t+t^{2}{\mbox{ , }}2\leq t\leq 3}$

будет членом этого типа, а также

${\displaystyle S(t)=P_{0}(t)=-2-2t^{2}{\mbox{ , }}0\leq t<1}$

${\displaystyle S(t)=P_{1}(t)=1-6t+t^{2}{\mbox{ , }}1\leq t<2}$

${\displaystyle S(t)=P_{2}(t)=-1+t-2t^{2}{\mbox{ , }}2\leq t\leq 3}$

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

Самый простой сплайн имеет степень 0. Его также называют сплайном. ступенчатая функция. Следующий по простоте сплайн имеет степень 1. Его также называют сплайном. линейный шлиц. Замкнутый линейный сплайн (т. Е. Первый узел и последний одинаковые) на плоскости - это просто многоугольник.

Обычный шлиц - это естественный кубический шлиц степени 3 с непрерывностью C².Слово "естественный" означает, что вторые производные полиномов сплайна устанавливаются равными нулю в конечных точках интервала интерполяции.

${\displaystyle S''(a)\,=S''(b)=0.}$

Это заставляет сплайн быть прямой линией за пределами интервала, не нарушая при этом его гладкости.

Алгоритм вычисления натуральных кубических сплайнов

Кубические шлицы имеют вид ${\displaystyle {S}_{j}\left(x\right)=a_{j}+b_{j}\left(x-x_{j}\right)+c_{j}{\left(x-x_{j}\right)}^{2}+d_{j}{\left(x-x_{j}\right)}^{3}}$.
Заданный набор координат ${\displaystyle C=\left[\left({x}_{0},{y}_{0}\right),\left({x}_{1},{y}_{1}\right),....,\left({x}_{n},{y}_{n}\right)\right]}$ мы хотим найти набор ${\displaystyle n\,}$ шлицы ${\displaystyle {S}_{i}\left(x\right)}$ за ${\displaystyle i=0,\ldots ,n-1.}$

Они должны удовлетворять:

• ${\displaystyle S_{i}\left(x_{i}\right)=y_{i}=S_{i-1}\left(x_{i}\right),i=1,\ldots ,n-1.}$
• ${\displaystyle S_{i}'\left(x_{i}\right)=S_{i-1}'\left(x_{i}\right),i=1,\ldots ,n-1.}$
• ${\displaystyle S_{i}''\left(x_{i}\right)=S_{i-1}''\left(x_{i}\right),i=1,\ldots ,n-1.}$
• ${\displaystyle S_{0}''\left(x_{0}\right)=S_{n-1}''\left(x_{n}\right)=0}$.

Определим один кубический сплайн ${\displaystyle S\,}$ как пятерка ${\displaystyle (a,b,c,d,x_{t})\,}$ куда ${\displaystyle a,b,c\,}$ и ${\displaystyle d\,}$ соответствуют коэффициентам в форме, показанной ранее, и ${\displaystyle x_{t}\,}$ равно ${\displaystyle x_{j}.\,}$

Алгоритм вычисления естественных кубических сплайнов:
Ввод: набор координат ${\displaystyle C\,}$, с ${\displaystyle \left|C\right|=n+1}$
Выход: установить сплайны, состоящие из п 5-кортежи.

1. Создать новый массив а размера п + 1 и для ${\displaystyle i=0,\ldots ,n}$ набор ${\displaystyle a_{i}=y_{i}\,}$
2. Создать новые массивы б и d каждый размер п.
3. Создать новый массив час размера п и для ${\displaystyle i=0,\ldots ,n-1}$ набор ${\displaystyle h_{i}=x_{i+1}-x_{i}\,}$
4. Создать новый массив α размера п и для ${\displaystyle i=1,\ldots ,n-1}$ набор ${\displaystyle {\alpha }_{i}={\frac {3}{{h}_{i}}}\left({a}_{i+1}-{a}_{i}\right)-{\frac {3}{{h}_{i-1}}}\left({a}_{i}-{a}_{i-1}\right)}$.
5. Создать новые массивы c, л, μ, и z каждый размер ${\displaystyle n+1\,}$.
6. Набор ${\displaystyle l_{0}=1,{\mu }_{0}=z_{0}=0\,}$
7. За ${\displaystyle i=1,\ldots ,n-1\,}$
   1. Набор ${\displaystyle {l}_{i}=2\left({x}_{i+1}-{x}_{i-1}\right)-{h}_{i-1}{\mu }_{i-1}}$.
   2. Набор ${\displaystyle {\mu }_{i}={\frac {{h}_{i}}{{l}_{i}}}}$.
   3. Набор ${\displaystyle {z}_{i}={\frac {{\alpha }_{i}-{h}_{i-1}{z}_{i-1}}{{l}_{i}}}}$.
8. Набор ${\displaystyle l_{n}=1;z_{n}=c_{n}=0.\,}$
9. За ${\displaystyle j=n-1,n-2,\ldots ,0}$
   1. Набор ${\displaystyle c_{j}=z_{j}-{\mu }_{j}c_{j+1}\,}$
   2. Набор ${\displaystyle b_{j}={\frac {{a}_{j+1}-{a}_{j}}{{h}_{j}}}-{\frac {{h}_{j}\left({c}_{j+1}+2{c}_{j}\right)}{3}}}$
   3. Набор ${\displaystyle d_{j}={\frac {{c}_{j+1}-{c}_{j}}{3{h}_{j}}}}$
10. Создайте новый набор сплайнов и назовите его output_set. Заполните его п шлицы S.
11. За ${\displaystyle i=0,\ldots ,n-1}$
    1. Набор S_я,а = а_я
    2. Набор S_я,б = б_я
    3. Набор S_я,c = c_я
    4. Набор S_я,d = d_я
    5. Набор S_я,Икс = Икс_я
12. Выход output_set

Примечания

Можно спросить, что означает больше, чем п несколько узлов в векторе узла имеют, так как это приведет к непрерывности, подобной

${\displaystyle S(t)\in C^{-m}{\mbox{ , }}m>0}$

на месте этой высокой множественности. По соглашению, любая такая ситуация указывает на простой разрыв между двумя соседними полиномиальными частями. Это означает, что если узел т_я появляется больше чем п +1 раз в расширенном векторе узла, все его экземпляры превышают (п + 1) th можно убрать без изменения характера сплайна, так как все кратности п + 1, п + 2, п + 3 и т. Д. Имеют то же значение. Принято считать, что таким образом отбирался любой узловой вектор, определяющий любой тип сплайна.

Классический сплайн типа степени п используется в численном анализе, имеет преемственность

${\displaystyle S(t)\in \mathrm {C} ^{n-1}[a,b],\,}$

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

Другой тип сплайна, который часто используется в графике, например, в программах для рисования, таких как Adobe Illustrator из Adobe Systems, имеет кубические части, но имеет непрерывность не более

${\displaystyle S(t)\in \mathrm {C} ^{1}[a,b].}$

Этот тип сплайна также используется в PostScript а также в определении некоторых компьютерных типографских шрифтов.

Многие системы автоматизированного проектирования, разработанные для высококачественной графики и анимации, используют, например, расширенные узловые векторы. майя из Псевдоним. В системах автоматизированного проектирования часто используется расширенная концепция шлицев, известная как Неравномерный рациональный B-сплайн (NURBS).

Если доступны выборки данных из функции или физического объекта, сплайн-интерполяция - это подход к созданию сплайна, который аппроксимирует эти данные.

Общее выражение для C² Интерполяция кубического сплайна

Общее выражение для яth C² интерполяция кубического сплайна в точке Икс с естественным условием можно найти по формуле

${\displaystyle S_{i}(x)={\frac {z_{i}(x-t_{i-1})^{3}}{6h_{i}}}+{\frac {z_{i-1}(t_{i}-x)^{3}}{6h_{i}}}+\left[{\frac {f(t_{i})}{h_{i}}}-{\frac {z_{i}h_{i}}{6}}\right](x-t_{i-1})+\left[{\frac {f(t_{i-1})}{h_{i}}}-{\frac {z_{i-1}h_{i}}{6}}\right](t_{i}-x)}$

куда

• ${\displaystyle z_{i}=f^{\prime \prime }(t_{i})}$ - значения второй производной на яй узел.
• ${\displaystyle h_{i}^{}=t_{i}-t_{i-1}}$
• ${\displaystyle f(t_{i}^{})}$ - значения функции на яй узел.

Представления и имена

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

${\displaystyle a=t_{0}<\underbrace {t_{1}=\cdots =t_{1}} _{j_{1}}<\cdots <\underbrace {t_{k-2}=\cdots =t_{k-2}} _{j_{k-2}}<t_{k-1}=b}$

${\displaystyle j_{i}\leq n+1~,~~i=1,\ldots ,k-2.}$

Размерность равна сумме степени плюс кратности

${\displaystyle d=n+\sum _{i=1}^{k-2}j_{i}.}$

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

Литература по шлицам изобилует названиями специальных типов шлицев, которые связаны с:

• Варианты, сделанные для представления сплайна, например:
  • с помощью основа функции для всего сплайна (дав нам имя B-шлицы )
  • с помощью Многочлены Бернштейна как использовал Пьер Безье для представления каждого полиномиального фрагмента (дав нам имя Шлицы Безье )
• Варианты, сделанные при формировании вектора расширенного узла, например:
  • используя отдельные узлы для C^п-1 непрерывность и равномерное расположение этих узлов на [а,б] (давая нам равномерные шлицы)
  • использование узлов без ограничения расстояния (что дает нам неоднородные шлицы)
• Любые особые условия на сплайн, например:
  • обеспечение нулевых вторых производных на а и б (давая нам естественные шлицы)
  • требуя, чтобы заданные значения данных были на сплайне (что дает нам интерполирующие сплайны)

Часто специальное имя выбиралось для типа сплайна, удовлетворяющего двум или более основным пунктам, указанным выше. Например, Шпонка Эрмита представляет собой сплайн, который выражается с помощью полиномов Эрмита для представления каждой из отдельных частей полинома. Они чаще всего используются с п = 3; то есть как Кубические шлицы Эрмита. В этой степени они могут быть дополнительно выбраны только непрерывными по касательной (C¹); откуда следует, что все внутренние узлы двойные. Было изобретено несколько способов подгонки таких сплайнов к заданным точкам данных; то есть превратить их в интерполирующие сплайны и сделать это путем оценки правдоподобных значений касательной в местах пересечения каждых двух частей полинома (что дает нам Кардинальные шлицы, Шлицы Catmull-Rom, и Шлицы Кочанека-Бартельса, в зависимости от используемого метода).

Для каждого из представлений должны быть найдены некоторые средства оценки, чтобы значения сплайна могли быть получены по запросу. Для тех представлений, которые выражают каждый отдельный полиномиальный кусок п_я(т) с точки зрения некоторой основы для получения степени п полиномы, это концептуально просто:

• Для данного значения аргумента т, найдите интервал, в котором он лежит ${\displaystyle t\in [t_{i},t_{i+1}]}$
• Найдите полиномиальный базис, выбранный для этого интервала ${\displaystyle P_{0},\ldots ,P_{k-2}}$
• Найдите значение каждого базисного полинома в т: ${\displaystyle P_{0}(t),\ldots ,P_{k-2}(t)}$
• Найдите коэффициенты линейной комбинации тех базисных многочленов, которые дают сплайн на этом интервале. c₀, ..., c_k-2
• Сложите эту линейную комбинацию значений базисного полинома, чтобы получить значение сплайна в т:

${\displaystyle \sum _{j=0}^{k-2}c_{j}P_{j}(t).}$

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

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

История

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

Принято считать, что первая математическая ссылка на сплайны - это работа 1946 г. Шенберг, что, вероятно, первое место, где слово «сплайн» используется в связи с гладкой, кусочно-полиномиальной аппроксимацией. Однако эти идеи уходят корнями в авиастроение и судостроение. В предисловии к (Bartels et al., 1987), Робин Форрест описывает "подъем ", метод, используемый в британской авиационной промышленности во время Вторая Мировая Война для создания шаблонов для самолетов, пропуская тонкие деревянные полоски (так называемые "шлицы ") через точки, расположенные на полу большого дизайнерского чердака, техника, заимствованная из проектирования корпуса корабля. В течение многих лет практика проектирования судов использовала модели для проектирования в малом. Успешный проект затем был нанесен на миллиметровую бумагу и ключевые точки графика были перенесены на увеличенную миллиметровую бумагу в полный размер. Тонкие деревянные полоски обеспечивали интерполяцию ключевых точек в плавные кривые. Полосы удерживались на месте в дискретных точках (которые Форрест назвал «утками» ; Шенберг использовал "собак" или "крыс") и между этими точками принимал формы с минимальной энергией деформации. По словам Форреста, одним из возможных стимулов для математической модели этого процесса была потенциальная потеря критических компонентов конструкции для всего самолета. на случай попадания в чердак вражеской бомбы. Это привело к «коническому подъему», при котором конические сечения использовались для моделирования положения кривой между утками. Конический подъем был заменен тем, что мы назвали бы сплайнами в начале 1960-х годов. Эд на работе Дж. К. Фергюсон в Боинг и (несколько позже) М.А. Сабин в Британская авиастроительная корпорация.

Слово «сплайн» изначально было Восточноанглийский диалектное слово.

Использование шлицев для моделирования автомобильных кузовов имеет несколько независимых начал. Кредит запрашивается от имени de Casteljau в Citroën, Пьер Безье в Renault, и Биркофф, Карабедян, и де Бур в Дженерал Моторс (см. Birkhoff and de Boor, 1965), все для работ, происходивших в самом начале 1960-х или в конце 1950-х годов. По крайней мере, одна из статей де Кастельжау была опубликована, но не широко, в 1959 году. Дженерал Моторс в результате в начале 1960-х годов был опубликован ряд статей, в том числе некоторые фундаментальные работы по B-шлицы.

Работа также велась в компании Pratt & Whitney Aircraft, где были задействованы два автора (Ahlberg et al., 1967) - первой трактовки шлицев в целую книгу, а также Модель бассейна Дэвида Тейлора, Федора Тейльхаймера. Работа на Дженерал Моторс подробно описано в (Birkhoff, 1990) и (Young, 1997). Дэвис (1997) резюмирует некоторые из этих материалов.

Рекомендации

• Фергюсон, Джеймс К., Многопараметрическая интерполяция кривой, J. ACM, т. 11, вып. 2, стр. 221-228, апрель 1964 г.
• Альберг, Нильсон и Уолш, Теория сплайнов и их приложения, 1967.
• Биркгоф, Гидродинамика, расчеты реакторов и представление поверхностей, в: Стив Нэш (ред.), История научных вычислений, 1990.
• Бартельс, Битти и Барский, Введение в сплайны для использования в компьютерной графике и геометрическом моделировании, 1987.
• Биркгоф и де Бур, Кусочно-полиномиальная интерполяция и аппроксимация, в: Х.Л. Гарабедян (ред.), Proc. Симпозиум General Motors 1964 г., С. 164–190. Эльзевир, Нью-Йорк и Амстердам, 1965 год.
• Дэвис, B-шлицы и геометрический дизайн, Новости SIAM, т. 29, нет. 5, 1997.
• Эпперсон, История сплайнов, Дайджест НС, т. 98, нет. 26, 1998.
• Стоер и Булирш, Введение в численный анализ. Springer-Verlag. п. 93-106. ISBN  0387904204
• Шенберг, Вклад в проблему аппроксимации эквидистантных данных аналитическими функциями, Кварта. Appl. Математика, т. 4. С. 45–99 и 112–141, 1946.
• Янг, Гаррет Биркгоф и прикладная математика, Уведомления AMS, т. 44, нет. 11. С. 1446–1449, 1997.
• Чапра, Канале, "Численные методы для инженеров", 5-е издание.

внешняя ссылка

Теория

• Интерактивное введение в сплайны, ibiblio.org

Функция Excel

• Дополнительная функция XLL Excel Реализация кубического сплайна

Онлайн-утилиты

• Онлайн-утилита интерполяции кубических сплайнов
• Обучение с помощью моделирования Интерактивное моделирование различных кубических шлицев
• Симметричные кривые сплайна, анимация Теодор Грей, Демонстрационный проект Wolfram, 2007.

Компьютерный код

• Заметки, PPT, Mathcad, Maple, Mathematica, Matlab, Институт целостных численных методов
• различные процедуры, NTCC
• Sisl: C-библиотека с открытым исходным кодом для NURBS, SINTEF
• Сплайн-интерполяция VBA, vbnumericalmethods.com