sinkriging.m

来自「用GIEWANK检验Kriging精度的算法」· M 代码 · 共 25 行

function[a,b,s]=sin(x)
a=0.5:0.5:10;
a=a';%输入数据S
b=sin(a);%输入数据Y
s=0;%精度
theta=10;lob=0.01;upb=20;
[dmodel, perf]=dacefit(a,b,@regpoly0,@corrgauss,theta,lob,upb);
c=4.01:0.01:6;
c=c';%检验点
d=sin(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-4)/0.01;
for i=1:n
    if s<=f(i)*f(i)
       s=f(i)*f(i);
    end
end
s=sqrt(s)