Syntax: |
w = BIVINTERP(x,y,z,u,v)
|
This function interpolates, from values of the function given at input grid points in an x-y plane and for a given set of points in the plane, the values of a single-valued bivariate function z = z(x,y). The method is based on a piece-wise function composed of a set of bicubic polynomials in x and y. Each polynomial is applicable to a rectangle of the input grid in the x-y plane. Each polynomial is determined locally.
The first two input parameters are vectors. Vector x
contains the x-coordinates
of the input grid points, in ascending order. Vector y
contains the
y-coordinates of the input grid points, in ascending order. Both x
and
y
must be monotonically increasing. The third parameter is a matrix, z
,
which contains the values of the function at the input grid points, z[i][j]
is the
data value at (x[i],y[j])
. The last two parameters are vectors. Vector
u
contains the x-coordinates of the desired
points, and vector v
contains the y-coordinates of the desired points. Vectors
u
and v
must have the same number of elements. The output is a vector,
w
, containing the interpolated values, w[i]
is the interpolated value at the
location (u[i],v[i])
.
Algorithm derived from an article by Hiroshi Akima, Communications of the ACM, volume 17, number 1, January 1974, pp. 26-31.