clear all
clc
% generate the numbers by :
% [url]http://muchong.com/bbs/viewthread.php?tid=3918964&fpage=1[/url]
n1=100;
n2=n1*10;
x=rand(n1,1);
y=rand(n1,1);
z=rand(n2,1);


x_new=repmat(x,[n2/n1,1]);
y_new=repmat(y,[n2/n1,1]);

% transform the data
data=[x_new';y_new';z']';

% calculate the mean value of each point.
for i=1:n1
     ind=[data(:,1)==x(i)&data(:,2)==y(i)];
     z_temp=data(ind,:);
     z_mean(i)=mean(z_temp(:,3));
end

% output the data.
x_temp=[min(x):(max(x)-min(x))/n1:max(x)];
y_temp=[min(y):(max(y)-min(y))/n1:max(y)];

[X,Y]=meshgrid(x_temp,y_temp);
Z = griddata(x,y,z_mean,X,Y);

mesh(X,Y,Z)
hold on
plot3(x,y,z_mean,'o')
% if u feel the original point is too chaos, u can ctrl+r the following
% codes
plot3(x_new,y_new,z,'r.')
hold off

% typed by : mm
% typed on : 19:28 2011-12-10