关于应用Isight做代理模型及EI采点+优化

1.代理模型应用MATLAB的DACE工具箱制作

制作流程代码如下：

load('shuju1.txt');
X(:,1)=shuju1(:,1);
X(:,2)=shuju1(:,2);
YY(:,1)=shuju1(:,3);
theta=[10 10];
lob=[1e-1 1e-1];
upb=[20 20];
[dmodel,perf]=dacefit(X,YY,@regpoly2,@corrgauss,theta,lob,upb);
[y MSE]=predictor([x1;x2],dmodel);

一：加载txt数据并赋值

二：初始化k模型,theta初始值和上下限，并选择回归模型（高斯）和模型（二阶）

三：倒数第二句是求得代理模型dmodel

四：倒数第一句可以由代理模型，输入任意X，获得对应的输出Y.

z=(x1^2+x2-11)^2+(x1+x2^2-7)^2;
dzz=abs(z-y);
contz=1;
if(dzz>1||dzz<-1)
load('shuju1.txt');
fout=fopen('shuju1.txt');
[m,n]=size(shuju1);
for k=1:m
fprintf(fout,'%d %d %d\r\n',shuju1(k,1),shuju1(k,2),shuju1(k,3));
end
fprintf(fout,'%d %d %d\r\n',x1,x2,z