World Library  
Flag as Inappropriate
Email this Article

Spline (mathematics)

Article Id: WHEBN0000457680
Reproduction Date:

Title: Spline (mathematics)  
Author: World Heritage Encyclopedia
Language: English
Subject: Smoothing spline, Interpolation, PDIFF, Carl R. de Boor, List of mathematical shapes
Collection: Interpolation, Splines (Mathematics)
Publisher: World Heritage Encyclopedia
Publication
Date:
 

Spline (mathematics)

A quadratic spline composed of six polynomial segments. Between point 0 and point 1 a straight line. Between point 1 and point 2 a parabola with second derivative = 4. Between point 2 and point 3 a parabola with second derivative = -2. Between point 3 and point 4 a straight line. Between point 4 and point 5 a parabola with second derivative = 6. Between point 5 and point 6 a straight line.
A cubic spline composed of seven polynomial segments. This shape is used as a pulse in the article Pulse (physics)
The second derivative of the cubic spline above.

In mathematics, a spline is a numeric function that is piecewise-defined by polynomial functions, and which possesses a sufficiently high degree of smoothness at the places where the polynomial pieces connect (which are known as knots).[1][2]

In interpolating problems, spline interpolation is often preferred to polynomial interpolation because it yields similar results to interpolating with higher degree polynomials while avoiding instability due to Runge's phenomenon. In computer graphics, parametric curves whose coordinates are given by splines are popular 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.

The most commonly used splines are cubic spline, i.e., of order 3—in particular, cubic B-spline, which is equivalent to C2 continuous composite Bézier curves.[3] They are common, in particular, in spline interpolation simulating the function of flat splines.

The term spline is adopted from the name of a flexible strip of metal commonly used by drafters to assist in drawing curved lines.[4]

Contents

  • Examples 1
  • History 2
  • Definition 3
  • Derivation of a cubic spline interpolating between points 4
  • See also 5
  • References 6
  • Further reading 7
  • External links 8
    • Theory 8.1
    • Excel functions 8.2
    • Online utilities 8.3
    • Computer code 8.4

Examples

The bell-shaped Irwin-Hall spline
The second derivative of the spline above

A simple example of a quadratic spline (a spline of degree 2) is

S(t) = \begin{cases} (t+1)^2-1 & -2 \le t < 0\\ 1-(t-1)^2 & 0 \le t \le 2 \end{cases}

for which S'(0)=2.

A simple example of a cubic spline is

S(t) = \left|t\right|^3

as

S(t) = \begin{cases} t^3 & t \ge 0\\ -t^3 & t < 0 \end{cases}

and

S'(0) =\ 0
S''(0) =\ 0

An example of using a cubic spline to create a bell shaped curve is the Irwin-Hall distribution polynomials:

f_X(x)= \begin{cases} \frac{1}{4}(x+2)^3 & -2\le x \le -1\\ \frac{1}{4}\left(3|x|^3 - 6x^2 +4 \right)& -1\le x \le 1\\ \frac{1}{4}(2-x)^3 & 1\le x \le 2 \end{cases}

History

Before computers were used, numerical calculations were done by hand. Functions such as the step function were used but polynomials were generally preferred. With the advent of computers, splines first replaced polynomials in interpolation, and then served in construction of smooth and flexible shapes in computer graphics.[5]

The word "spline" originally meant a thin wood or metal slat in East Anglian dialect. By 1895 it had come to mean a flexible ruler used to draw curves.[6] These splines were used in the aircraft and shipbuilding industries. For years ship designers had employed models to design hulls 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 the key points (using lead weights called "ducks"[7] (see Spline Ducks for illustration) or "dogs" or "rats"[8]) and between these points would assume shapes of minimum strain energy. In the foreword to Bartels et al. (1987),[7] 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.

It is commonly accepted that the first mathematical reference to splines is the 1946 paper by Schoenberg,[8] which is probably the first place that the word "spline" is used in connection with smooth, piecewise polynomial approximation. 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[9] at Boeing and (somewhat later) by M.A. Sabin at British Aircraft Corporation.

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,[10] Garabedian, and de Boor at General Motors (see Birkhoff and de Boor, 1965),[11] 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.[12]

Work was also being done at Pratt & Whitney Aircraft, where two of the authors of the first book-length treatment of splines (Ahlberg et al., 1967)[13] 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).[14] Davis (1997) summarizes some of this material.

Definition

A spline is a piecewise-polynomial real function

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

on an interval [a,b] composed of k subintervals [t_{i-1}, t_i] with

a = t_0 < t_1 < \cdots < t_{k-1} < t_k = b.

The restriction of S to an interval i is a polynomial

P_i: [t_{i-1}, t_i] \to \mathbb{R},

so that

S(t) = P_1 (t) \mbox{ , } t_0 \le t < t_1,
S(t) = P_2 (t) \mbox{ , } t_1 \le t < t_2,
\vdots
S(t) = P_k (t) \mbox{ , } t_{k-1} \le t \le t_{k}.

The highest order of the polynomials P_i (t) is said to be the order of the spline S. The spline is said to be uniform if all subintervals are of the same length, and non-uniform otherwise.[15]

The idea is to choose the polynomials in a way that guarantees sufficient smoothness of S. Specifically, for a spline of order n, S is required to be both continuous and continuously differentiable to order n-1 at the interior points t_i: for i=1, \dots, k-1 and j=0, \dots, n-1

P_i^{(j)} (t_i) = P_{i+1}^{(j)} (t_i).

Derivation of a cubic spline interpolating between points

Spline interpolation is one of the most common uses of splines.

See also

References

  1. ^ Judd, Kenneth L. (1998). Numerical Methods in Economics. MIT Press. p. 225.  
  2. ^ Chen, Wai-Kai (2009). Feedback, Nonlinear, and Distributed Circuits. CRC Press. pp. 9–20.  
  3. ^ http://books.google.ro/books?id=9bQ0f8sYqaAC&pg=PA215
  4. ^ Katz, Mitchell H. (2011). Multivariable Analysis: A Practical Guide for Clinicians and Public Health Researchers. Cambridge University Press. p. 82.  
  5. ^ Epperson, History of Splines, NA Digest, vol. 98, no. 26, 1998.
  6. ^ "spline".  
  7. ^ a b Bartels, Beatty, and Barsky, An Introduction to Splines for Use in Computer Graphics and Geometric Modeling, 1987.
  8. ^ a b Schoenberg, Contributions to the problem of approximation of equidistant data by analytic functions, Quart. Appl. Math., vol. 4, pp. 45–99 and 112–141, 1946.
  9. ^ Ferguson, James C, Multi-variable curve interpolation, J. ACM, vol. 11, no. 2, pp. 221-228, Apr. 1964.
  10. ^ Birkhoff, Fluid dynamics, reactor computations, and surface representation, in: Steve Nash (ed.), A History of Scientific Computation, 1990.
  11. ^ Birkhoff and de Boor, Piecewise polynomial interpolation and approximation, in: H. L. Garabedian (ed.), Proc. General Motors Symposium of 1964, pp. 164–190. Elsevier, New York and Amsterdam, 1965.
  12. ^ Davis, B-splines and Geometric design, SIAM News, vol. 29, no. 5, 1997.
  13. ^ Ahlberg, Nilson, and Walsh, The Theory of Splines and Their Applications (1967)
  14. ^ Young, Garrett Birkhoff and applied mathematics, Notices of the AMS, vol. 44, no. 11, pp. 1446–1449, 1997.
  15. ^ Fan, Jianqing & Yao, Qiwei (2005). "Spline Methods". Nonlinear time series: nonparametric and parametric methods. Springer. p. 247.  

Further reading

  • Smith, P. L. (1979). "Splines as a Useful and Convenient Statistical Tool". The American Statistician 33 (2): 57–62.  
  • Stoer; Bulirsch. Introduction to Numerical Analysis. Springer Science+Business Media. pp. 93–106.  
  • Chapra, Canale. Numerical Methods for Engineers (5th ed.). 

External links

Theory

Excel functions

  • Open source Excel cubic spline User Defined Function
  • SRS1 Cubic Spline for Excel - Free Excel cubic spline function (with utility to embed spline function code into any workbook)

Online utilities

  • Visual, interactive comparison of various types of splines
  • 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.

Computer code

  • Notes, PPT, Mathcad, Maple, Mathematica, Matlab, Holistic Numerical Methods Institute
  • various routines, NTCC
  • Sisl: Opensource C-library for NURBS, SINTEF
  • Closed Bezier Spline, C#, WPF, Oleg V. Polikarpotchkin
  • Bezier Spline from 2D Points, C#, WPF, Oleg V. Polikarpotchkin
This article was sourced from Creative Commons Attribution-ShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and USA.gov, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for USA.gov and content contributors is made possible from the U.S. Congress, E-Government Act of 2002.
 
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
 
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a non-profit organization.
 


Copyright © World Library Foundation. All rights reserved. eBooks from Project Gutenberg are sponsored by the World Library Foundation,
a 501c(4) Member's Support Non-Profit Organization, and is NOT affiliated with any governmental agency or department.