Welcome to AE Resources
INTERPOLATION

INTERPOLATION

These examples also serve to illustrate that with more complex functions, the use of 1st order interpolation is not advisable. Instead, the polynomials can be used as interpolators. These polynomials can also be used to extrapolate data, but these can be problematic, as is shown in Fig. 4↓.
figure images/exrap_error.jpg
Figure 4 Extrapolation not valid using some polynomial interpolation

A more mathematically (computationally) efficient polynomial uses the nested forms of the polynomial. These forms include the Newton and Lagrangian forms of the equations.

NEWTON INTERPOLATING POLYNOMIALS
Consider the equation
p(x) = b0 + (x − c1){b1 + (x − c2){b2 + (x − c3)...{bn + (x − cn)}}}
Newton’s polynomial can be used for interpolating to fit a nth polynomial to n+1 points.
p(x) = bo + b1(x − xo) + ... + bn(x − xo)(x − x1)...(x − xn)
Here, x is the point of interest and at each point xo, x1, ..., xn, the function p is known. In order to compute the unknown value, x, the coefficients bo, b1, ..., bn must be computed. These are found using the following set of equations:
bo  =  p(xo) b1  =  p(x1, xo) = (p(x1) − p(xo))/(x1 − xo) . . . bn  =  pn[xn...xo] = (p(xn, xn − 1, ..., x1) − p(xn − 1, ...xo))/(xn − xo)
As an example to illustrate computation of the b coefficients, consider the situation when n = 2
p(x)  =  bo + b1(x − xo) + b2(x − xo)(x − x1) bo  =  p(xo) b1  =  p(x1, xo) = (p(x1) − p(xo))/(x1 − xo) b2  =  p(x2, x1, xo) = (p(x2, x1) − p(x1, xo))/(x2 − xo)
This method uses the coefficients coefficients for lower n values for each successive larger n value. A two-dimensional array can be set up to compute the b(0:n) coefficients to use the function divided differences already computed. Similar to a Taylor series, a remainder or truncation error can be computed:
Rn = p(x, xn, ..., xo)(x − xo)(x − x1)...(x − xn)
Notice that the remainder calculation requires that the function be known. An alternate way to estimate the error when the function is not known is to use an additional known point for the error estimate:
Rn = (pn + 1)/((n + 1)!)x)n + 1
Next Page →
Next Page →
← Previous Page
← Previous Page