griddata 调用方法：
ZI = griddata(x,y,z,XI,YI)
[XI,YI,ZI] = griddata(x,y,z,XI,YI)
[...] = griddata(...,method)
[...] = griddata(...,method,options)

method 的值 为
'linear' -- 则，以三角形为基础的线性内插
'cubic' -- 则，以三角形为基础的三次方程内插
'nearest' -- 则，用最邻近的点 内插
'v4'-- -- 则，MATLAB 4 格点样条函数内插
默认'linear' 线性内插

三角形为基础，就是按Delaunay方法先找出内插点四周的3个点，构成三角形，内插点在三角形内。然后线性内插，或三次方程内插。

'cubic' 和 'v4' 插值结果构成的曲面较光滑，'linear'和 'nearest' 插值结果构成的曲面不光滑不连续。

前3种方法具体算法见
[1] Barber, C. B., D.P. Dobkin, and H.T. Huhdanpaa, "The Quickhull Algorithm for Convex Hulls," ACM Transactions on Mathematical Software, Vol. 22, No. 4, Dec. 1996, p. 469-483. Available in PDF format at http://www.acm.org/pubs/citations/journals/toms/1996-22-4/p469-barber/.

第4种方法具体算法见
[2] Sandwell, David T., "Biharmonic Spline Interpolation of GEOS-3 and SEASAT Altimeter Data", Geophysical Research Letters, 14, 2, 139-142,1987.

(参考了 MathWorks 主站材料）