标签：rand dmodel 模型 40 YX matlab theta X1 kriging

kriging工具箱：https://orbit.dtu.dk/en/publications/dace-a-matlab-kriging-toolbox

x=rand(1,100)*5;
y=rand(1,100)*5;
z=x./(y+1)+0.01*rand(1,100);
data=[x',y',z'];
scatter(x,y,25,z);
colorbar;

%模型参数设置
theta = [5 5]; lob = [1e-1 1e-1]; upb = [20 20];
%变异函数模型为高斯模型
[dmodel, perf] = dacefit(data(:,1:2), data(:,3), @regpoly0, @corrgauss, theta, lob, upb);
%创建一个40*40的网格，范围为0-5
X = gridsamp([0 0;5 5], 40);
%格网点的预测值返回在矩阵YX中，预测点的均方根误差返回在矩阵MSE中
[YX,MSE] = predictor(X, dmodel);
X1 = reshape(X(:,1),40,40); X2 = reshape(X(:,2),40,40);
YX = reshape(YX, size(X1));         %size(X1)=40*40
mesh(X1, X2, YX);         %绘制预测表面

theta是初始值，lob和upb是参数范围，计算结果存在dmodel.theta中。

>> dmodel.theta
ans =
   0.903390452322443   4.718047509652705

当对模型不确定或者变量波动剧烈时，所给范围宜较宽。

若初始值选取不合理，则可能造成曲面波动。

如以下情况：

标签：rand,dmodel,模型,40,YX,matlab,theta,X1,kriging
来源： https://www.cnblogs.com/dingdangsunny/p/15658914.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。