关于#MATLAB#克里金插值的问题，如何解决？

MATLAB克里金插值出现问题

clc,clear
data=load('E:\pcl\代码\测试1\zhitongdaochu6.txt');
S=data(1:100,1:2);
Z=data(1:100,3);
shu=size(Z);

[mS,mZ]=dsmerge(S,Z,1e-14);
theta = [-100 100]; lob = [1e-1 1e-1]; upb = [20 20];%参数
%调用克里金插值算法工具箱
%进行拟合操作
[dmodel, perf] = dacefit(mS, mZ, @regpoly0, @corrspherical, theta, lob, upb);

%创建一个40*40的格网，标注范围为0-100，即格网间距为2.5
%S存储了点位坐标值，Y为观测值
p=100;
X = gridsamp([0 0;20000 20000], p);  

% X=[83.731    32.36];     %单点预测的实现
%格网点的预测值返回在矩阵YX中，预测点的均方根误差返回在矩阵MSE中
[YX,MSE] = predictor(X, dmodel); 
%转换坐标 n行n列
X1 = reshape(X(:,1),p,p); X2 = reshape(X(:,2),p,p);
YX = reshape(YX, size(X1));         %size(X1)=40*40
figure(1), mesh(X1, X2, YX)         %绘制预测表面

z轴一直是一个值，没有变化，但读入的z轴数据是正确变化的

希望平台的大家能指点一下