MATLAB Mathematics - SERC - Index of

Interpolation
The next step is to use interpolation to compute function values over a grid. Use
meshgrid to create the grid, and griddatan to do the interpolation:
delta = 0.05;
d = -1:delta:1;
[x0,y0,z0] = meshgrid(d,d,d);
X0 = [x0(:), y0(:), z0(:)];
v0 = griddatan(X,v,X0);
v0 = reshape(v0, size(x0));
Then use isosurface and related functions to visualize the surface that
consists of the (x,y,z) values for which the function takes a constant value. You
could pick any value, but the example uses the value 0.6. Since the function is
the squared distance from the origin, the surface at a constant value is a
sphere:
p = patch(isosurface(x0,y0,z0,v0,0.6));
isonormals(x0,y0,z0,v0,p);
set(p,'FaceColor','red','EdgeColor','none');
view(3);
camlight;
lighting phong
axis equal
title('Interpolated sphere from scattered data')
Note A smaller delta produces a smoother sphere, but increases the
compute time.
2-35