[Kriging](二)Matlab Kriging工具箱dace简单使用

工具箱

工具箱下载来自于百度网盘，由于是别人分享的这里就不做引用了。

配置

将工具箱解压缩，然后放置在matlab的toolbox路径下，然后在matlab内设置路径，如下图所示

简单测试

load data1.mat

其中S相当于xy的坐标点，Y相当于每个坐标点对应的值。换言之，如果要使用自己的测试样例数据，加入我们有一个data由XYZ组成，那么我们可以先把XY存储为Y变量，然后，再把Z单独存为S变量，再把S和Y共同保存在一个data.mat文件中以供dace工具箱使用，达到自定义数据的效果。

读进来的就是工具箱内的一个样例数据data1.mat，本博文先简单测试下如何使用工具箱进行kriging，选择POLY0回归方法以及GAUSS相关方法，假设点为各向异性。然后使用dacefit方法来进行建模
使用dmodel内的参数进行新数据点的值估计，使用了predictor方法。

代码

load data1.mat;
% 初始化参数
theta = [10 10]; lob = [1E-1 1E-1]; upb = [20 20];
[dmodel, perf] = dacefit(S, Y, @regpoly0, @corrgauss, theta, lob, upb);

X = gridsamp([0 0;100 100], 40);
[YX, MSE] = predictor(X, dmodel);

X1 = reshape(X(:,1),40,40); X2 = reshape(X(:,2),40,40);
YX = reshape(YX, size(X1));

figure(1); mesh(X1, X2, YX);
hold on;
plot3(S(:,1),S(:,2),Y,'.k', 'MarkerSize',10);
hold off;

运行结果

前文

(一)网上下载的kriging克里金的C++程序的初步调试