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为均方误差。