函数griddata的超牛改进版，值得一下哦

% Gridfit demo script

% This script file is designed to be used in cell mode
% from the matlab editor, or best of all, use the publish
% to HTML feature from the matlab editor. Older versions
% of matlab can copy and paste entire blocks of code into
% the Matlab command window.

Topographic dataload bluff_data;
x=bluff_data(:,1);
y=bluff_data(:,2);
z=bluff_data(:,3);
% Two ravines on a hillside. Scanned from a
% topographic map of an area in upstate New York.
plot3(x,y,z,'.')


% Turn the scanned point data into a surface
gx=0:4:264;
gy=0:4:400;
g=gridfit(x,y,z,gx,gy);
figure
colormap(hot(256));
surf(gx,gy,g);
camlight right;
lighting phong;
shading interp
line(x,y,z,'marker','.','markersize',4,'linestyle','none');
title 'Use topographic contours to recreate a surface'


Fitting a trigonometric surfaceclear

n1 = 15;
n2 = 15;
theta = rand(n1,1)*pi/2;
r = rand(1,n2);

x = cos(theta)*r;
y = sin(theta)*r;
x=x(:);
y=y(:);

x = [[0 0 1 1]';x;x;1-x;1-x];
y = [[0 1 0 1]';y;1-y;y;1-y];
figure
plot(x,y,'.')
title 'Data locations in the x-y plane'


z = sin(4*x+5*y).*cos(7*(x-y))+exp(x+y);

xi = linspace(0,1,51);
[xg,yg]=meshgrid(xi,xi);
zgd = griddata(x,y,z,xg,yg);

figure
surf(xi,xi,zgd)
colormap(hot(256))
camlight right
lighting phong
title 'Griddata on trig surface'
% Note the wing-like artifacts along the edges, due
% to the use of a Delaunay triangulation in griddata.

zgrid = gridfit(x,y,z,xi,xi);

figure
surf(xi,xi,zgrid)
colormap(hot(256))
camlight right
lighting phong
title('Gridfit to trig surface')


The trig surface with highly different scalings on the x and y axesxs = x/100;
xis = xi/100;

ys = y*100;
yis = xi*100;

% griddata has problems with badly scaled data
[xg,yg]=meshgrid(xis,yis);
zgd = griddata(xs,ys,z,xg,yg);

figure
surf(xg,yg,zgd)
colormap(hot(256))
camlight right
lighting phong
title 'Serious problems for griddata on badly scaled trig surface'

% autoscaling on (the default)
zgrids = gridfit(xs,ys,z,xis,yis,'autoscale','on');

% plot the autoscaled result
figure
surf(xis,yis,zgrids)
colormap(hot(256))
camlight right
lighting phong
title 'Gridfit (automatically scaled) on trig surface'


Fitting the "peaks" surfaceclear

n = 100;
x = (rand(n,1)-.5)*6;
y = (rand(n,1)-.5)*6;

z = peaks(x,y);

xi = linspace(-3,3,101);
zpgf = gridfit(x,y,z,xi,xi);

[xg,yg]=meshgrid(xi,xi);
zpgd = griddata(x,y,z,xg,yg,'cubic');

figure
surf(xi,xi,zpgd)
colormap(jet(256))
camlight right
lighting phong
title 'Griddata (method == cubic) on peaks surface'

figure
surf(xi,xi,zpgf)
colormap(hsv(256))
camlight right
lighting phong
title('Gridfit to peaks surface')

我的理解，gridata是把一块有限大的布铺在数据点上，而这个改进的函数是把一块无限大的木板钉在数据点上，这样看起来不仅曲面很光滑而且可以外推作简单的预报。来自www.。

学习一下，谢谢楼主分享成果。

非常好，谢谢你的分享！

这么用啊这个？哪位大神说一声啊？

谢谢楼主分享···

下载顶起！

x=0:400:5600;
y=0:400:4800;
z=[370 470 550 600 670 690 670 620 580 450 400 300 100 150 250
   510 620 730 800 850 870 850 780 720 650 500 200 300 350 320
   650 760 880 970 1020 1050 1020 830 800 700 300 500 550 480 350
   740 880 1080 1130 1250 1280 1230 1040 900 500 700 780 750 650 550
   830 980 1180 1320 1450 1420 1400 1300 700 900 850 840 380 780 750
   880 1060 1230 1390 1500 1500 1400 900 1100 1060 950 870 900 930 950
   910 1090 1270 1500 1200 1100 1350 1450 1200 1150 1010 880 1000 1050 1100
   950 1190 1370 1500 1200 1100 1550 1600 1550 1380 1070 900 1050 1150 1200
   1430 1450 1460 1500 1550 1600 1550 1600 1600 1600 1550 1500 1500 1550 1500
   1420 1430 1450 1480 1500 1550 1510 1430 1300 1200 980 850 750 550 500
   1380 1410 1430 1450 1470 1320 1280 1200 1080 940 780 620 460 370 350
   1370 1390 1410 1430 1440 1140 1110 1050 950 820 690 540 380 300 210
   1350 1370 1390 1400 1410 960 940 880 800 690 570 430 290 210 150
];
xi=0:4:5600;
yi=0:4:4800;
[xii,yii]=meshgrid(xi,yi);
zii=gridfit(x,y,z,xii,yii);
楼主能帮忙看看这么写报错了