log2kriging.m

用GIEWANK检验Kriging精度的算法

function[a,b,s]=log2(x)
a=0.5:0.5:10;
a=a';%输入数据S
b=log2(a);%输入数据Y
s=0;%精度
theta=10;lob=0.01;upb=20;
[dmodel, perf]=dacefit(a,b,@regpoly0,@corrgauss,theta,lob,upb);
c=5.51:0.01:6;
c=c';%检验点
d=log2(c);
e=predictor(c, dmodel);
figure(1);
plot(c,d,'r',c,e,'g')
figure(2);
plot(c,d,'r')
figure(3);
plot(c,e,'g')
f=d-e;
n=(6-5.5)/0.01;
for i=1:n
    if s<=f(i)*f(i)
       s=f(i)*f(i);
    end
end
s=sqrt(s)