Monotone cubic interpolation

In the mathematical subfield of numerical analysis, monotone cubic interpolation is a variant of cubic interpolation that preserves monotonicity of the data set being interpolated.

Monotonicity is preserved by linear interpolation but not guaranteed by cubic interpolation.

Monotone cubic Hermite interpolation

Example showing cubic and monotone cubic interpolation of a monotone data set. Note that the red line is not monotone.

Monotone interpolation can be accomplished using cubic Hermite spline with the tangents m_i modified to ensure the monotonicity of the resulting Hermite spline.

Data preprocessing

Let the data points be (x_k,y_k) for k = 1,...,n

1. Compute the slopes of the secant lines between successive points:

   

   for .
2. Initialize the tangents at every data point as the average of the secants,

   

   for ; these may be updated in further steps. For the endpoints, use one-sided differences:

   

3. For , if Δ_k = 0 (if two successive y_k = y_{k + 1} are equal), then set m_k = m_{k + 1} = 0, as the spline connecting these points must be flat to preserve monotonicity. Ignore step 4 and 5 for those k.
4. Let α_k = m_k / Δ_k and β_k = m_{k + 1} / Δ_k.
5. Now, since to prevent overshoot or undershoot, we are trying to restrict the position vector (α_k,β_k) to a circle of radius 3, if , then set m_k = τ_kα_kΔ_k and m_{k + 1} = τ_kβ_kΔ_k where .

Note that only one pass of the algorithm is required.

Cubic interpolation

After the preprocessing, evaluation of the interpolated spline is equivalent to cubic Hermite spline, using the data x_k, y_k, and m_k for k = 1,...,n.

To evaluate at x, find the smallest value larger than x, x_upper, and the largest value smaller than x, x_lower, among x_k such that . Calculate

h = x_upper − x_lower and

then the interpolant is

f_interpolated(x) = y_lowerh₀₀(t) + hm_lowerh₁₀(t) + y_upperh₀₁(t) + hm_upperh₁₁(t)

where h_ii are the basis functions for the cubic Hermite spline.

References

• Fritsch, F. N.; Carlson, R. E. (1980). "Monotone Piecewise Cubic Interpolation". SIAM Journal on Numerical Analysis (SIAM) 17 (2): 238–246. doi:10.1137/0717021.