%?根據問題的主要特點：三維-散點-曲面光滑-擬合，下面一個例子基本滿足要求，以供參考：


clcclearclose?all

%%?矩陣數據曲面光滑
%?初始數據
n1=6;?x=linspace（-22n1）;
n2=11;?y=linspace（05n2）;
[Y0X0]=meshgrid（yx）;
Z0=X0.^2+sqrt（Y0）;
noisy=rand（n1n2）;
Z1=Z0+noisy;

%樣條擬合
p=0.9?%
sp=csaps（{xy}Z1p）;
v=fnval（sp{xy}）;

%?mesh圖
figure
subplot（221）mesh（X0Y0Z0）title（‘基準數據Z0‘）;?xlabel（‘x‘）ylabel（‘y‘）
subplot（222）mesh（X0Y0Z1）title（‘待光滑處理數據Z1‘）;
subplot（223）mesh（X0Y0v）title（‘光滑后數據v‘）;
subplot（224）mesh（X0Y0v-Z1）title（‘誤差v-Z1‘）;


%%?散點數據曲面光滑
n=20;
x=rand（1n）*2;?y=rand（1n）*2;?z=x.^2+（y）;
noisy=rand（1n）;?z1=z+noisy;

%?散點插值生產矩陣數據
n1=21;?n2=21;
xi=linspace（min（x）+0.1max（x）-0.1n1）;
yi=linspace（min（y）+0.1max（y）-0.1n2）;
[YiXi]=meshgrid（yixi）;
Z1?=?griddata（xyz1XiYi）;
Z0?=?griddata（xyz