Product Code Database
barcode-scavenger
» » Wiki: Interpolation
Tag

# Interpolation  ( Video )

C O N T E N T S
Rank: 100%

In the field of numerical analysis, interpolation is a method of constructing new within the range of a of known data points.

In and , one often has a number of data points, obtained by sampling or , which represent the values of a function for a limited number of values of the independent variable. It is often required to interpolate, i.e., estimate the value of that function for an intermediate value of the independent variable.

A closely related problem is the approximation of a complicated function by a simple function. Suppose the formula for some given function is known, but too complicated to evaluate efficiently. A few data points from the original function can be interpolated to produce a simpler function which is still fairly close to the original. The resulting gain in simplicity may outweigh the loss from interpolation error.

Example
This table gives some values of an unknown function $f\left(x\right)$.
 ! x|width="10px" !colspan=3 align=center f( x) 0 8415 9093 1411 7568 9589 2794
Interpolation provides a means of estimating the function at intermediate points, such as $x=2.5$.

We describe some of interpolation, differing in such properties as: accuracy, cost, number of data points needed, and of the resulting function.

Piecewise constant interpolation
The simplest interpolation method is to locate the nearest data value, and assign the same value. In simple problems, this method is unlikely to be used, as linear interpolation (see below) is almost as easy, but in higher-dimensional multivariate interpolation, this could be a favourable choice for its speed and simplicity.

Linear interpolation
One of the simplest methods is interpolation (sometimes known as lerp). Consider the above example of estimating f(2.5). Since 2.5 is midway between 2 and 3, it is reasonable to take f(2.5) midway between f(2) = 0.9093 and f(3) = 0.1411, which yields 0.5252.

Generally, linear interpolation takes two data points, say ( x a, y a) and ( x b, y b), and the interpolant is given by:

$y = y_a + \left\left( y_b-y_a \right\right) \frac\left\{x-x_a\right\}\left\{x_b-x_a\right\} \text\left\{ at the point \right\} \left\left( x,y \right\right)$

$\frac\left\{y-y_a\right\}\left\{y_b-y_a\right\} = \frac\left\{x-x_a\right\}\left\{x_b-x_a\right\}$

$\frac\left\{y-y_a\right\}\left\{x-x_a\right\} = \frac\left\{y_b-y_a\right\}\left\{x_b-x_a\right\}$

This previous equation states that the slope of the new line between $\left(x_a,y_a\right)$ and $\left(x,y\right)$ is the same as the slope of the line between $\left(x_a,y_a\right)$ and $\left(x_b,y_b\right)$

Linear interpolation is quick and easy, but it is not very precise. Another disadvantage is that the interpolant is not at the point x k.

The following error estimate shows that linear interpolation is not very precise. Denote the function which we want to interpolate by g, and suppose that x lies between x a and x b and that g is twice continuously differentiable. Then the linear interpolation error is

$|f\left(x\right)-g\left(x\right)| \le C\left(x_b-x_a\right)^2 \quad\text\left\{where\right\}\quad C = \frac18 \max_\left\{r\inx_a,x_b\right\} |g\text{'}\text{'}\left(r\right)|.$

In words, the error is proportional to the square of the distance between the data points. The error in some other methods, including polynomial interpolation and spline interpolation (described below), is proportional to higher powers of the distance between the data points. These methods also produce smoother interpolants.

Polynomial interpolation
Polynomial interpolation is a generalization of linear interpolation. Note that the linear interpolant is a . We now replace this interpolant with a of higher degree.

Consider again the problem given above. The following sixth degree polynomial goes through all the seven points:

$f\left(x\right) = -0.0001521 x^6 - 0.003130 x^5 + 0.07321 x^4 - 0.3577 x^3 + 0.2255 x^2 + 0.9038 x.$

Substituting x = 2.5, we find that f(2.5) = 0.5965.

Generally, if we have n data points, there is exactly one polynomial of degree at most n−1 going through all the data points. The interpolation error is proportional to the distance between the data points to the power n. Furthermore, the interpolant is a polynomial and thus infinitely differentiable. So, we see that polynomial interpolation overcomes most of the problems of linear interpolation.

However, polynomial interpolation also has some disadvantages. Calculating the interpolating polynomial is computationally expensive (see computational complexity) compared to linear interpolation. Furthermore, polynomial interpolation may exhibit oscillatory artifacts, especially at the end points (see Runge's phenomenon).

Polynomial interpolation can estimate local maxima and minima that are outside the range of the samples, unlike linear interpolation. For example, the interpolant above has a local maximum at x ≈ 1.566, f( x) ≈ 1.003 and a local minimum at x ≈ 4.708, f( x) ≈ −1.003. However, these maxima and minima may exceed the theoretical range of the function—for example, a function that is always positive may have an interpolant with negative values, and whose inverse therefore contains false vertical asymptotes.

More generally, the shape of the resulting curve, especially for very high or low values of the independent variable, may be contrary to commonsense, i.e. to what is known about the experimental system which has generated the data points. These disadvantages can be reduced by using spline interpolation or restricting attention to Chebyshev polynomials.

Spline interpolation
Remember that linear interpolation uses a linear function for each of intervals x k, x k+1. Spline interpolation uses low-degree polynomials in each of the intervals, and chooses the polynomial pieces such that they fit smoothly together. The resulting function is called a spline.

For instance, the natural cubic spline is cubic and twice continuously differentiable. Furthermore, its second derivative is zero at the end points. The natural cubic spline interpolating the points in the table above is given by

$f\left(x\right) = \begin\left\{cases\right\}$
-0.1522 x^3 + 0.9937 x, & \text{if } x \in 0,1, \\ -0.01258 x^3 - 0.4189 x^2 + 1.4126 x - 0.1396, & \text{if } x \in 1,2, \\ 0.1403 x^3 - 1.3359 x^2 + 3.2467 x - 1.3623, & \text{if } x \in 2,3, \\ 0.1579 x^3 - 1.4945 x^2 + 3.7225 x - 1.8381, & \text{if } x \in 3,4, \\ 0.05375 x^3 -0.2450 x^2 - 1.2756 x + 4.8259, & \text{if } x \in 4,5, \\ -0.1871 x^3 + 3.3673 x^2 - 19.3370 x + 34.9282, & \text{if } x \in 5,6. \end{cases}

In this case we get f(2.5) = 0.5972.

Like polynomial interpolation, spline interpolation incurs a smaller error than linear interpolation and the interpolant is smoother. However, the interpolant is easier to evaluate than the high-degree polynomials used in polynomial interpolation. However, the global nature of the basis functions leads to ill-conditioning. This is completely mitigated by using splines of compact support, such as are implemented in Boost.Math and discussed in Kress.

Via Gaussian processes
is a powerful non-linear interpolation tool. Many popular interpolation tools are actually equivalent to particular Gaussian processes. Gaussian processes can be used not only for fitting an interpolant that passes exactly through the given data points but also for regression, i.e., for fitting a curve through noisy data. In the geostatistics community Gaussian process regression is also known as .

Other forms
Other forms of interpolation can be constructed by picking a different class of interpolants. For instance, rational interpolation is interpolation by rational functions using Padé approximant, and trigonometric interpolation is interpolation by trigonometric polynomials using . Another possibility is to use .

The Whittaker–Shannon interpolation formula can be used if the number of data points is infinite.

Sometimes, we know not only the value of the function that we want to interpolate, at some points, but also its derivative. This leads to Hermite interpolation problems.

When each data point is itself a function, it can be useful to see the interpolation problem as a partial problem between each data point. This idea leads to the displacement interpolation problem used in transportation theory.

In higher dimensions
Multivariate interpolation is the interpolation of functions of more than one variable. Methods include bilinear interpolation and bicubic interpolation in two dimensions, and trilinear interpolation in three dimensions. They can be applied to gridded or scattered data.

In digital signal processing
In the domain of digital signal processing, the term interpolation refers to the process of converting a sampled digital signal (such as a sampled audio signal) to that of a higher sampling rate () using various digital filtering techniques (e.g., convolution with a frequency-limited impulse signal). In this application there is a specific requirement that the harmonic content of the original signal be preserved without creating aliased harmonic content of the original signal above the original Nyquist limit of the signal (i.e., above fs/2 of the original signal sample rate). An early and fairly elementary discussion on this subject can be found in Rabiner and Crochiere's book Multirate Digital Signal Processing. R.E. Crochiere and L.R. Rabiner. (1983). Multirate Digital Signal Processing. Englewood Cliffs, NJ: Prentice–Hall.

Related concepts
The term is used to find data points outside the range of known data points.

In problems, the constraint that the interpolant has to go exactly through the data points is relaxed. It is only required to approach the data points as closely as possible (within some other constraints). This requires parameterizing the potential interpolants and having some way of measuring the error. In the simplest case this leads to approximation.

Approximation theory studies how to find the best approximation to a given function by another function from some predetermined class, and how good this approximation is. This clearly yields a bound on how well the interpolant can approximate the unknown function.

Generalization
If we consider $x$ as a variable in a topological space, with and the function $f\left(x\right)$ mapping to a , then the problem is treated as "interpolation of operators".Colin Bennett, Robert C. Sharpley, Interpolation of Operators, Academic Press 1988 The classical results about interpolation of operators are the Riesz–Thorin theorem and the Marcinkiewicz theorem. There are also many other subsequent results.

• Barycentric coordinates – for interpolating within on a triangle or tetrahedron
• Bilinear interpolation
• Brahmagupta's interpolation formula
• Imputation (statistics)
• Lagrange interpolation
• Multivariate interpolation
• Newton–Cotes formulas
• Polynomial interpolation
• Simple rational approximation

 Page 1 of 1 1
 Page 1 of 1 1

### Account

Social:
Pages:  ..   ..
Items:  ..

General: Atom Feed  ..
Help:  ..   ..
Category:  ..   ..
Media:  ..   ..
Posts:  ..   ..   ..

Page:  ..
Summary:  ..
1 Tags
10/10 Page Rank
5 Page Refs
4s Time