matlab源码--散点数据生成格网、等高线并三维显示（通过克吕金插值）

kriging程序，用于散点数据生成格网，且可以绘制等高线图，二维或者三维显示。内有matlab的程序代码以及原理讲解的PDF

function F = barnes （xi yi zi c g）

%	F = barnes （xi yi zi c g）
% 
% Filter kriged data.
% 
% 	F（ij） = f0（ij） + sum （wn * Dn） / sum （wn）
%
% o?f0（ij）= barnes first estimation 
%	wn     = exp（-h^2 / （4*c*g）） = weigth function for x and y data
%	Dn     = zi-f0 （of coordinates （xiyi））
%	cg    = filter constants
%	F（ij） = filtered data 	
%
% refs:	Tintor?et al. （1991） Mesoscale Dynamics and Vertical Motion in
%		the Alboran Sea Journal of Physical Oceanigraphy 21811:823.
%
%	Maddox （1980） An objective Technique for Separating Macroscale and 
%		Mesoscale Features in Meteorological Data Montly Weather 
%		Review 108 1108:1121.

% Weigth function 

  wf0 = ‘exp（-h.^2 ./ （4*c））‘;
  ww  = ‘exp（-h.^2 ./ （4*c*g））‘;

% Input format: position
  
  nx = length（xi）;
  ny = length（yi）;
  p  = kregrid（min（xi）xi（2）-xi（1）max（xi）...
		min（yi）yi（2）-yi（1）max（yi））;
  x = p（:1）;
  y = p（:2）;

% Input format: data
  
  [mn] = size（zi）;
  z = []; 
  for i = 1:n
	z = [z;zi（i:）‘];  
  end

% Sum over N: all data

  N = length（z）;

% Filter to get f0（ij）

 for j = 1:N
 
    dx = x - x（j）;
    dy = y - y（j）;
    h  = （dx.^2 + dy.^2） .^（0.5）;
    wn = eval（[wf0]）;

    wf   = sum（wn .* z）;
    w    = sum（wn）;
    f0（j） = wf / w;

  end	

% Filter to get F（ij）

  for j = 1:N
 
    dx = x - x（j）;
    dy = y - y（j）;
    h  = （dx.^2 + dy.^2） .^（0.5）;
    wn = eval（[ww]）;

    wf   = sum（wn .* （z-f0‘））;
    w    = sum（wn）;
    f（j） = wf / w;

  end	

% Output

  F = deplie（f0+fnxny）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         438  1998-05-15 15:38  chiaux.m
     文件        2071  1998-05-15 15:38  chitable.m
     文件        1270  1998-05-15 15:38  chiprob.m
     文件        1450  1998-05-15 15:39  contents_chi.m
     文件       18616  1998-05-15 15:38  a2.mat
     文件        2396  1998-05-15 15:39  variogr.m
     文件        3238  1998-05-15 15:38  Vario3dr.m
     文件        3435  1998-05-15 15:38  Vario3di.m
     文件        3219  1998-05-15 15:38  Vario2dr.m
     文件        3027  1998-05-15 15:38  Vario2di.m
     文件        3144  1998-05-15 15:39  var2diuv.m
     文件        1613  1998-05-15 15:39  trans.m
     文件        7144  2001-07-25 14:14  TMP.mat
     文件         692  1998-05-15 15:39  tintore.m
     文件         655  1998-05-15 15:39  stn2utm.m
     文件        1154  1998-05-15 15:39  sal15.mat
     文件        1324  1998-05-15 15:39  readme
     文件        3034  1998-05-15 15:39  plt90stn.m
     文件        1299  1998-05-15 15:39  outvario.m
     文件         576  1998-05-15 15:39  mrqminv.mat
     文件        2827  1998-05-15 15:39  mrqmin.m
     文件         985  1998-05-15 15:39  mrqcof.m
     文件         481  1998-05-15 15:39  means.m
     文件         297  1998-05-15 15:39  mat4dp.m
     文件         259  1998-05-15 15:39  mat3dp.m
     文件        3638  1998-05-15 15:39  ll2mtu.m
     文件        1332  2001-07-25 13:37  KStest.m
     文件      154112  2001-07-25 15:00  ksone.dll
     文件        6339  2001-07-25 14:51  kridemo.m
     文件       13967  1998-05-15 15:39  kri.mat
     文件         536  1998-05-15 15:39  kregrid3.m
............此处省略27个文件信息