Interpolating Cubic Polynomial Splines

Top  Previous  Next

 

Cubic Polynomial Splines

 

Assume that image002ap values image004ap of a function image006ap are known at the nodes image008ap, image010ap. We desire to construct a smooth curve passing through the given points image012ap with polynomial splines of third degree, also called  cubic splines. If we assume that the nodes are ordered image014ap, then we can represent the desired curve by a spline function image016ap with image018ap, which is composed of cubic polynomials image020ap for image022ap and image024ap. These image020ap are differently for each subinterval and must satisfy certain connecting conditions at the nodes.

 

Definition of Interpolating Cubic Spline Functions

 

Depending on the conditions imposed at image026ap and image028ap, different types of spline functions image016ap will result. For a set of ordered nodes image008ap with image014ap, the spline functions mentioned previously can be defined by the following conditions:

 

(1)image030ap.
(2)image016ap is a cubic polynomial image020ap in each subinterval image032ap, image024ap.
(3)image016ap interpolates image034ap, that is, image036ap for image010ap.

 

Types of Interpolating Cubic Spline Functions

 

(a)Natural Cubic Splines: For image038ap or image040ap, image016ap is represented by the tangent of image016ap at image026ap or image042ap, thus image044ap.
(b)Complete Cubic Splines: Cubic spline with specified first end point derivatives — For the end point conditions image046ap, image048ap.
(c)Generalized Natural Cubic Spline: Cubic spline with specified second end point derivatives — For the end point conditions image050ap, image052ap.
(d)Cubic Spline with Specified Third End-Point Derivatives: For the end point conditions image054ap, image056ap.
(e)Cubic Spline with Not-a-Knot Condition: If image016ap satisfies the endpoint conditions image058ap and image060ap. This condition means that the third derivative of the spline function is continuous at the nodes image062ap and image064ap. Thus image062ap and image064ap are not “genuine” nodes for the spline function.

 

To construct a function image016ap that satisfies (1) and (2) above, we set

 

image066ap, for image022ap, image024ap

 

A cubic polynomial spline is twice continuously differentiable, depending on 4 parameters. Hence it is a spline function of odd degree image070ap for image072ap [image074ap-times continuously differentiable, image076ap-parametric].

 

Properties (1) and (3) for image016ap lead to image078ap conditions for the image020ap:

 

(A)image080ap,                image010ap
(B)image082ap,                image010ap
(C)image084ap,        image024ap
(D)image087ap,        image010ap

 

Here we have set image089ap and image091ap. Two additional end point conditions can be specified at will. Then we have to satisfy image093ap conditions for the  image093ap coefficients image095ap for image024ap and image097ap. Property (1) is the single most rigorous requirement of a spline function image016ap. It gives rise to conditions (C) and (D) and effects a smooth transition between the graphs of the polynomials image020ap and image099ap at the node image101ap for each image103ap: The graphs of adjacent polynomials image020ap and  image099ap have the same curvature at image101ap. This property of spline functions is especially useful for approximating a function image034ap which is known empirically and which one can draw using a flexible ruler, commonly known as a shipbuilders spline.