DACE-Kriging 模型用法

1. 关于theta取值

    Dace工具箱内dacefit函数自带theta 优化，只需在初始设置 theta0、lob 和upb即可。

   猜测初始值theta0的选取，不影响最终结果。但是影响计算效率和加点数。 

    Dace工具箱手册内给的算例为：

    theta=[10 10]; lob=[1e-1 1e-1]; upb=[20 20];

    [dmodel, perf]=dacefit(S,Y,@repoly0,@corrgauss,theta,lob,upb)

2. 关于multiple design site 的问题：

    错误原因为：Kriging 建模时 不允许有重复设计点。

    解决办法为：dsmerge

    [mS,mY]=dsmerge(S,Y)

    [mS,mY]=dsmerge(S,Y,ds)

    [mS,mY]=dsmerge(S,Y,ds,nms)

3. 关于Kriging建模梯度的问题：

    Dace 工具箱可以直接输出梯度。

    [y,dy,mse]=predictor(x,dmodel)

    其中： y为响应预测值；dy为梯度；mse为均方误差。