package scipy

One-dimensional interpolating spline for a given set of data points.

Fits a spline y = spl(x) of degree `k` to the provided `x`, `y` data. Spline function passes through all provided points. Equivalent to `UnivariateSpline` with s=0.

Parameters ---------- x : (N,) array_like Input dimension of data points -- must be strictly increasing y : (N,) array_like input dimension of data points w : (N,) array_like, optional Weights for spline fitting. Must be positive. If None (default), weights are all equal. bbox : (2,) array_like, optional 2-sequence specifying the boundary of the approximation interval. If None (default), ``bbox=x[0], x[-1]``. k : int, optional Degree of the smoothing spline. Must be 1 <= `k` <= 5. ext : int or str, optional Controls the extrapolation mode for elements not in the interval defined by the knot sequence.

* if ext=0 or 'extrapolate', return the extrapolated value. * if ext=1 or 'zeros', return 0 * if ext=2 or 'raise', raise a ValueError * if ext=3 of 'const', return the boundary value.

The default value is 0.

check_finite : bool, optional Whether to check that the input arrays contain only finite numbers. Disabling may give a performance gain, but may result in problems (crashes, non-termination or non-sensical results) if the inputs do contain infinities or NaNs. Default is False.

See Also -------- UnivariateSpline : Superclass -- allows knots to be selected by a smoothing condition LSQUnivariateSpline : spline for which knots are user-selected splrep : An older, non object-oriented wrapping of FITPACK splev, sproot, splint, spalde BivariateSpline : A similar class for two-dimensional spline interpolation

Notes ----- The number of data points must be larger than the spline degree `k`.

Examples -------- >>> import matplotlib.pyplot as plt >>> from scipy.interpolate import InterpolatedUnivariateSpline >>> x = np.linspace(-3, 3, 50) >>> y = np.exp(-x**2) + 0.1 * np.random.randn(50) >>> spl = InterpolatedUnivariateSpline(x, y) >>> plt.plot(x, y, 'ro', ms=5) >>> xs = np.linspace(-3, 3, 1000) >>> plt.plot(xs, spl(xs), 'g', lw=3, alpha=0.7) >>> plt.show()

Notice that the ``spl(x)`` interpolates `y`:

>>> spl.get_residual() 0.0

Construct a new spline representing the antiderivative of this spline.

Parameters ---------- n : int, optional Order of antiderivative to evaluate. Default: 1

Returns ------- spline : UnivariateSpline Spline of order k2=k+n representing the antiderivative of this spline.

Notes -----

.. versionadded:: 0.13.0

See Also -------- splantider, derivative

Examples -------- >>> from scipy.interpolate import UnivariateSpline >>> x = np.linspace(0, np.pi/2, 70) >>> y = 1 / np.sqrt(1 - 0.8*np.sin(x)**2) >>> spl = UnivariateSpline(x, y, s=0)

The derivative is the inverse operation of the antiderivative, although some floating point error accumulates:

>>> spl(1.7), spl.antiderivative().derivative()(1.7) (array(2.1565429877197317), array(2.1565429877201865))

Antiderivative can be used to evaluate definite integrals:

>>> ispl = spl.antiderivative() >>> ispl(np.pi/2) - ispl(0) 2.2572053588768486

This is indeed an approximation to the complete elliptic integral :math:`K(m) = \int_0^\pi/2 1 - m\sin^2 x^

1/2

}

dx`:

>>> from scipy.special import ellipk >>> ellipk(0.8) 2.2572053268208538

Construct a new spline representing the derivative of this spline.

Parameters ---------- n : int, optional Order of derivative to evaluate. Default: 1

Returns ------- spline : UnivariateSpline Spline of order k2=k-n representing the derivative of this spline.

See Also -------- splder, antiderivative

Notes -----

.. versionadded:: 0.13.0

Examples -------- This can be used for finding maxima of a curve:

>>> from scipy.interpolate import UnivariateSpline >>> x = np.linspace(0, 10, 70) >>> y = np.sin(x) >>> spl = UnivariateSpline(x, y, k=4, s=0)

Now, differentiate the spline and find the zeros of the derivative. (NB: `sproot` only works for order 3 splines, so we fit an order 4 spline):

>>> spl.derivative().roots() / np.pi array( 0.50000001, 1.5 , 2.49999998)

This agrees well with roots :math:`\pi/2 + n\pi` of :math:`\cos(x) = \sin'(x)`.

Return all derivatives of the spline at the point x.

Parameters ---------- x : float The point to evaluate the derivatives at.

Returns ------- der : ndarray, shape(k+1,) Derivatives of the orders 0 to k.

Examples -------- >>> from scipy.interpolate import UnivariateSpline >>> x = np.linspace(0, 3, 11) >>> y = x**2 >>> spl = UnivariateSpline(x, y) >>> spl.derivatives(1.5) array(2.25, 3.0, 2.0, 0)

Return spline coefficients.

Return positions of interior knots of the spline.

Internally, the knot vector contains ``2*k`` additional boundary knots.

Return weighted sum of squared residuals of the spline approximation.

This is equivalent to::

sum((wi * (yi-spl(xi)))**2, axis=0)

Return definite integral of the spline between two given points.

Parameters ---------- a : float Lower limit of integration. b : float Upper limit of integration.

Returns ------- integral : float The value of the definite integral of the spline between limits.

Examples -------- >>> from scipy.interpolate import UnivariateSpline >>> x = np.linspace(0, 3, 11) >>> y = x**2 >>> spl = UnivariateSpline(x, y) >>> spl.integral(0, 3) 9.0

which agrees with :math:`\int x^2 dx = x^3 / 3` between the limits of 0 and 3.

A caveat is that this routine assumes the spline to be zero outside of the data limits:

>>> spl.integral(-1, 4) 9.0 >>> spl.integral(-1, 0) 0.0

Return the zeros of the spline.

Restriction: only cubic splines are supported by fitpack.

Continue spline computation with the given smoothing factor s and with the knots found at the last call.

This routine modifies the spline in place.

Print the object to a human-readable representation.

Print the object to a human-readable representation.

Pretty-print the object to a formatter.