1. Introduction

2. Definition
    Simple splines  Composing simple splines 

3. Examples
    Algorithm for computing natural cubic splines 

4. Continuity levels
     Natural continuity    Reduced continuity    Locally negative continuity  

5. General expression for a C² interpolating cubic spline

6. Representations and names

7. History

8. References

9. External links

In mathematics, a spline is a function defined piecewise by polynomials. In interpolating problems, spline interpolation is often preferred to polynomial interpolation because it yields similar results, even when using low degree polynomials, while avoiding Runge's phenomenon for higher degrees.

In the computer science subfields of computer-aided design and computer graphics, the term spline more frequently refers to a piecewise polynomial parametric curve{{citation needed|date=January 2018}}. Splines are popular curves in these subfields because of the simplicity of their construction, their ease and accuracy of evaluation, and their capacity to approximate complex shapes through curve fitting and interactive curve design.{{citation needed|date=January 2018}}

The term spline comes from the flexible spline devices used by shipbuilders and draftsmen to draw smooth shapes.^[1] A backronym of "Smooth Polynomial Lines Interpolating Numerical Estimates" has also been proposed{{Citation needed|date=March 2019}}.

Introduction

The term "spline" is used to refer to a wide class of functions that are used in applications requiring data interpolation and/or smoothing. The data may be either one-dimensional or multi-dimensional. Spline functions for interpolation are normally determined as the minimizers of suitable measures of roughness (for example integral squared curvature) subject to the interpolation constraints. Smoothing splines may be viewed as generalizations of interpolation splines where the functions are determined to minimize a weighted combination of the average squared approximation error over observed data and the roughness measure. For a number of meaningful definitions of the roughness measure, the spline functions are found to be finite dimensional in nature, which is the primary reason for their utility in computations and representation. For the rest of this section, we focus entirely on one-dimensional, polynomial splines and use the term "spline" in this restricted sense.

Definition

Simple splines

A simple spline is a piecewise polynomial function in a single variable. Such a spline S maps values from an interval [a,b] to the set of real numbers

For the piecewise definition of S the interval [a,b] is partitioned in n subintervals that together make up the whole interval. This is done by selecting additional n - 1 numbers t_k (1 ≤ k ≤ n - 1), representing points within this interval, such that

yielding n subintervals

such that

For each of these subintervals a polynomial P_i will be constructed

Excluding the right endpoints of each subinterval (with exception of the rightmost) allows the piecewise definition of S on [a,b] by the n polynomials P_i on the ith subinterval

If the n polynomials P_i each have degree at most m, then the spline is said to be of degree m (or is of order m+1).

The n+1 points t_j (0 ≤ j ≤ n) are called knots. The vector

is called a knot vector for the spline. If the knots are equidistantly distributed in the interval [a,b] the spline is called uniform, otherwise it is non-uniform.

Polynomial functions are continuous everywhere, they are said to belong to the class Polynomials are also r times continuously differentiable, and so they belong also to the class (for arbitrary r). Since for a polynomial function of grade m all derivatives of order m + 1 and higher vanish to the zero function, one can sensibly state that m is the highest (non-trivial) degree of smoothness for this polynomial function.

At each point common to two intervals, the two polynomials and are glued together for defining S. If the values of the two polynomials agree there, then S is continuous at If the two polynomials also have the same tangent there, that is, when also their first derivatives agree at then S is said to be in in a neighborhood of If the first derivatives of the two polynomials agree at then S is said to be in in a neighborhood of The values describe the smoothness of S at the gluing points Assuming a spline of degree m, the loss of smoothness (with respect to the maximal nontrivial smoothness) at the point is given by A loss of m at , i.e. establishes continuity there by imposing the condition a smaller loss, yields additional equations for involving the derivatives up to the desired smoothness. Collecting the n - 1 values yields the smoothness vector of the spline S.

An inner knot t_k can be "deleted" by making it equal to a neighboring knot t_{k ± 1}. By identifying t_k with t_{k + 1} the polynomial piece

disappears, and the pieces

join with the sum of the continuity losses for

t_k and t_k+1.

That is,

where

This leads to a more general understanding of a knot vector. The continuity loss at any point can be considered to be the result of multiple knots located at that point, and a spline type can be completely characterized by its degree m and its extended knot vector

where t_k is repeated j_k times

for . Obviously, n does not refer anymore to the number of components in this vector.

Given a knot vector , a degree m, and a smoothness vector for this knot vector, one can consider the set of all splines of degree having this knot vector and this smoothness vector . Equipped with the operation of adding two piecewise defined polynomial functions (pointwise addition) and taking their real multiples, together with the zero function, this set becomes a real vector space. This spline space is commonly denoted by .

Composing simple splines

A parametric curve in can be defined as a continuous map from an interval, the domain of the parameter, to the space For example, a plane curve is a mapping from the interval [a,b] to the Euclidean plane

with and being real functions on

If both functions and are simple spline functions of the same degree with the same extended knot vectors on that interval, then

specifies a spline curve in .

Examples

Suppose the interval [a,b] is [0,3] and the subintervals are [0,1], [1,2], and [2,3]. Suppose the polynomial pieces are to be of degree 2, and the pieces on [0,1] and [1,2] must join in value and first derivative (at t=1) while the pieces on [1,2] and [2,3] join simply in value (at t = 2).

This would define a type of spline S(t) for which

would be a member of that type, and also

would be a member of that type.

(Note: while the polynomial piece 2t is not quadratic, the result is still called a quadratic spline. This demonstrates that the degree of a spline is the maximum degree of its polynomial parts.) The extended knot vector for this type of spline would be (0, 1, 2, 2, 3).

The simplest spline has degree 0. It is also called a step function. The next most simple spline has degree 1. It is also called a linear spline. A closed linear spline (i.e, the first knot and the last are the same) in the plane is just a polygon.

A common spline is the natural cubic spline of degree 3 with continuity C². The word "natural" means that the second derivatives of the spline polynomials are set equal to zero at the endpoints of the interval of interpolation

Algorithm for computing natural cubic splines

Cubic splines have polynomial pieces of the form {{nowrap|Given }} coordinates we find polynomials, and which satisfy {{nowrap|for :}}

and

One such polynomial is given by a 4-tuple where and correspond to the coefficients as used above.

Input: a set of coordinates

Output: a spline as a set of polynomial pieces, each represented by a 5-tuple.

1. Create a new array a of size n + 1, and for set
2. Create new arrays b, d and μ each of size n
3. Create a new array h of size k and for set
4. Create a new array β of size n-1 and for set
5. Create new arrays c, l, and z each of size .
6. Set
7. For
   1. Set
   2. Set
   3. Set
8. Set
9. For
   1. Set
   2. Set
   3. Set
10. Create the spline as a new set of polynomials and call it output_set. Populate it with n 4-tuples for the polynomials P_i.
11. For
    1. Set P_i,a = a_i
    2. Set P_i,b = b_i
    3. Set P_i,c = c_i
    4. Set P_i,d = d_i
12. Output output_set

Continuity levels

If sampled data from a function or a physical object are available, spline interpolation is an approach to creating a spline that approximates those data.

Natural continuity

The classical spline type of degree n used in numerical analysis has continuity

which means that every two adjacent polynomial pieces meet in their value and first n - 1 derivatives at each knot. The mathematical spline that most closely models the flat spline is a cubic (n = 3), twice continuously differentiable (C²), natural spline, which is a spline of this classical type with additional conditions imposed at endpoints a and b.

Reduced continuity

Another type of spline that is much used in graphics, for example in drawing programs such as Adobe Illustrator from Adobe Systems, has pieces that are cubic but has continuity only at most

This spline type is also used in PostScript as well as in the definition of some computer typographic fonts.

Many computer-aided design systems that are designed for high-end graphics and animation use extended knot vectors, for example Maya from Alias. Computer-aided design systems often use an extended concept of a spline known as a non-uniform rational B-spline (NURBS).

Locally negative continuity

It might be asked what meaning more than n multiple knots in a knot vector have, since this would lead to continuities like

at the location of this high multiplicity. By convention, any such situation indicates a simple discontinuity between the two adjacent polynomial pieces. This means that if a knot t_i appears more than n + 1 times in an extended knot vector, all instances of it in excess of the (n + 1)th can be removed without changing the character of the spline, since all multiplicities n + 1, n + 2, n + 3, etc. have the same meaning. It is commonly assumed that any knot vector defining any type of spline has been culled in this fashion.

General expression for a C² interpolating cubic spline

The general expression for the ith C² interpolating cubic spline at a point x with the natural condition can be found using the formula

where

• are the values of the second derivative at the ith knot.
• 
• are the values of the function at the ith knot.

Representations and names

For a given interval [a,b] and a given extended knot vector on that interval, the splines of degree n form a vector space. Briefly this means that adding any two splines of a given type produces spline of that given type, and multiplying a spline of a given type by any constant produces a spline of that given type. The dimension of the space containing all splines of a certain type can be counted from the extended knot vector:

The dimension is equal to the sum of the degree plus the multiplicities

If a type of spline has additional linear conditions imposed upon it, then the resulting spline will lie in a subspace. The space of all natural cubic splines, for instance, is a subspace of the space of all cubic C² splines.

The literature of splines is replete with names for special types of splines. These names have been associated with:

• The choices made for representing the spline, for example:
  • using basis functions for the entire spline (giving us the name B-splines)
  • using Bernstein polynomials as employed by Pierre Bézier to represent each polynomial piece (giving us the name Bézier splines)
• The choices made in forming the extended knot vector, for example:
  • using single knots for C^n-1 continuity and spacing these knots evenly on [a,b] (giving us uniform splines)
  • using knots with no restriction on spacing (giving us nonuniform splines)
• Any special conditions imposed on the spline, for example:
  • enforcing zero second derivatives at a and b (giving us natural splines)
  • requiring that given data values be on the spline (giving us interpolating splines)

Often a special name was chosen for a type of spline satisfying two or more of the main items above. For example, the Hermite spline is a spline that is expressed using Hermite polynomials to represent each of the individual polynomial pieces. These are most often used with n = 3; that is, as cubic Hermite splines. In this degree they may additionally be chosen to be only tangent-continuous (C¹); which implies that all interior knots are double. Several methods have been invented to fit such splines to given data points; that is, to make them into interpolating splines, and to do so by estimating plausible tangent values where each two polynomial pieces meet (giving us cardinal splines, Catmull–Rom splines, and Kochanek–Bartels splines, depending on the method used).

For each of the representations, some means of evaluation must be found so that values of the spline can be produced on demand. For those representations that express each individual polynomial piece P_i(t) in terms of some basis for the degree n polynomials, this is conceptually straightforward:

• For a given value of the argument t, find the interval in which it lies
• Look up the polynomial basis chosen for that interval
• Find the value of each basis polynomial at t:
• Look up the coefficients of the linear combination of those basis polynomials that give the spline on that interval c₀, ..., c_k-2
• Add up that linear combination of basis polynomial values to get the value of the spline at t:

However, the evaluation and summation steps are often combined in clever ways. For example, Bernstein polynomials are a basis for polynomials that can be evaluated in linear combinations efficiently using special recurrence relations. This is the essence of De Casteljau's algorithm, which features in Bézier curves and Bézier splines.

For a representation that defines a spline as a linear combination of basis splines, however, something more sophisticated is needed. The de Boor algorithm is an efficient method for evaluating B-splines.

History

Before computers were used, numerical calculations were done by hand. Although piecewise-defined functions like the sign function or step function were used, polynomials were generally preferred because they were easier to work with. Through the advent of computers splines have gained importance. They were first used as a replacement for polynomials in interpolation, then as a tool to construct smooth and flexible shapes in computer graphics.

It is commonly accepted that the first mathematical reference to splines is the 1946 paper by Schoenberg, which is probably the first place that the word "spline" is used in connection with smooth, piecewise polynomial approximation. However, the ideas have their roots in the aircraft and shipbuilding industries. In the foreword to (Bartels et al., 1987), Robin Forrest describes "lofting", a technique used in the British aircraft industry during World War II to construct templates for airplanes by passing thin wooden strips (called "splines") through points laid out on the floor of a large design loft, a technique borrowed from ship-hull design. For years the practice of ship design had employed models to design in the small. The successful design was then plotted on graph paper and the key points of the plot were re-plotted on larger graph paper to full size. The thin wooden strips provided an interpolation of the key points into smooth curves. The strips would be held in place at discrete points (called "ducks" by Forrest; Schoenberg used "dogs" or "rats") and between these points would assume shapes of minimum strain energy. According to Forrest, one possible impetus for a mathematical model for this process was the potential loss of the critical design components for an entire aircraft should the loft be hit by an enemy bomb. This gave rise to "conic lofting", which used conic sections to model the position of the curve between the ducks. Conic lofting was replaced by what we would call splines in the early 1960s based on work by J. C. Ferguson at Boeing and (somewhat later) by M. A. Sabin at British Aircraft Corporation.

The word "spline" was originally an East Anglian dialect word.{{citation needed|date=September 2017}}

The use of splines for modeling automobile bodies seems to have several independent beginnings. Credit is claimed on behalf of de Casteljau at Citroën, Pierre Bézier at Renault, and Birkhoff, Garabedian, and de Boor at General Motors (see Birkhoff and de Boor, 1965), all for work occurring in the very early 1960s or late 1950s. At least one of de Casteljau's papers was published, but not widely, in 1959. De Boor's work at General Motors resulted in a number of papers being published in the early 1960s, including some of the fundamental work on B-splines.

Work was also being done at Pratt & Whitney Aircraft, where two of the authors of (Ahlberg et al., 1967)—the first book-length treatment of splines—were employed, and the David Taylor Model Basin, by Feodor Theilheimer. The work at General Motors is detailed nicely in (Birkhoff, 1990) and (Young, 1997). Davis (1997) summarizes some of this material.

References

• {{cite journal |last=Schoenberg |first=Isaac J. |url=http://www.ams.org/journals/qam/1946-04-01/S0033-569X-1946-15914-5/S0033-569X-1946-15914-5.pdf |title=Contributions to the Problem of Approximation of Equidistant Data by Analytic Functions: Part A.—On the Problem of Smoothing or Graduation. A First Class of Analytic Approximation Formulae |journal=Quarterly of Applied Mathematics |volume=4 |issue=2 |pages=45–99 |year=1946 |ref=harv }}
• {{cite journal |last=Schoenberg |first=Isaac J. |url=http://www.ams.org/journals/qam/1946-04-02/S0033-569X-1946-16705-2/S0033-569X-1946-16705-2.pdf |title=Contributions to the Problem of Approximation of Equidistant Data by Analytic Functions: Part B.—On the Problem of Osculatory Interpolation. A Second Class of Analytic Approximation Formulae |journal=Quarterly of Applied Mathematics |volume=4 |issue=2 |pages=112–141 |year=1946 |ref=harv }}
• {{cite journal |last=Ferguson |first=James C. |title=Multivariable Curve Interpolation |journal=Journal of the ACM |volume=11 |issue=2 |pages=221–228 |year=1964 |doi=10.1145/321217.321225 |ref=harv }}
• {{cite book |last=Ahlberg |first=J. Harold |last2=Nielson |first2=Edwin N. |last3=Walsh |first3=Joseph L. |title=The Theory of Splines and Their Applications |location=New York |publisher=Academic Press |year=1967 |isbn=0-12-044750-9 |ref=harv }}
• {{cite book |last=Birkhoff |first= |chapter=Fluid dynamics, reactor computations, and surface representation |editor-first=Steve |editor-last=Nash |title=A History of Scientific Computation |location=Reading |publisher=Addison-Wesley |year=1990 |isbn=0-201-50814-1 |ref=harv }}
• {{cite book |last=Bartels |first=Richard H. |last2=Beatty |first2=Brian A. |last3=Barsky |first3=John C. |title=An Introduction to Splines for Use in Computer Graphics and Geometric Modeling |location=San Mateo |publisher=Morgan Kaufmann |year=1987 |isbn=0-934613-27-3 |ref=harv }}
• {{cite book |last=Birkhoff |last2=de Boor |chapter=Piecewise polynomial interpolation and approximation |editor-first=H. L. |editor-last=Garabedian |title=Proc. General Motors Symposium of 1964 |pages=164–190 |publisher=Elsevier |location=New York and Amsterdam |year=1965 |ref=harv }}
• {{cite journal |last=Davis |url=http://www.wpi.edu/~pwdavis/sinews/spline17.htm |title=B-splines and Geometric design |journal=SIAM News |volume=29 |issue=5 |year=1997 |ref=harv }}
• {{cite journal |last=Epperson |url=http://www.netlib.org/na-digest-html/98/v98n26.html#1 |title=History of Splines |journal=NA Digest |volume=98 |issue=26 |pages= |year=1998 |ref=harv }}
• {{cite book |last=Stoer |last2=Bulirsch |title=Introduction to Numerical Analysis |location=New York |publisher=Springer-Verlag |pages=93–106 |isbn=0-387-90420-4 |ref=harv }}
• {{cite journal |last=Schoenberg |title=Contributions to the problem of approximation of equidistant data by analytic functions |journal=Quart. Appl. Math. |volume=4 |pages=45–99 and 112–141 |year=1946 |ref=harv }}
• {{cite journal |last=Young |title=Garrett Birkhoff and applied mathematics |journal=Notices of the AMS |volume=44 |issue=11 |pages=1446–1449 |year=1997 }}
• {{cite book |last=Chapra |last2=Canale |title=Numerical Methods for Engineers |location=Boston |publisher=McGraw-Hill |edition=5th |year=2006 |isbn=0-07-291873-X |ref=harv }}

Specific

External links

• Cubic Splines Module Prof. John H. Mathews California State University, Fullerton
• An Interactive Introduction to Splines, ibiblio.org

• XLL Excel Addin Function Implementation of cubic spline

• Online Cubic Spline Interpolation Utility
• Learning by Simulations Interactive simulation of various cubic splines
• Symmetrical Spline Curves, an animation by Theodore Gray, The Wolfram Demonstrations Project, 2007.
• [https://richardfuhr.neocities.org/BusyBCurves.html Exploring Bezier and Spline Curves], an interactive web app

• Notes, PPT, Mathcad, Maple, Mathematica, Matlab, Holistic Numerical Methods Institute
• various routines, NTCC
• Sisl: Opensource C-library for NURBS, SINTEF
• VBA Spline Interpolation, vbnumericalmethods.com

• Kot Pir Abdul Aziz Halt railway station
• Kot pri Damlju
• Kotpuli Nayanar
• Kotputli (Rajasthan Assembly constituency)
• Kotra, Berasia tehsil, Bhopal
• Kot Radha Kishan railway station
• Kot Radha Kishan Tehsil
• Kot Radha Kishn railway station
• Kotra, Huzur tehsil, Bhopal
• Kotramide
• Kotramide I
• Kotra River
• Kotravai
• Kotraza (Stragari)
• Kotraža (Stragari)
• Kotri Barrage
• Kotri barrage
• Kotri Bridge
• Kotrigur
• Kotrigurs
• Kotri Jn railway station
• Kotri Junction railway station
• Kotringo
• Kotromanic dynasty
• Kotromanović