sininterpl.m

用GIEWANK检验Kriging精度的算法

function[a,b,s]=sin(x)
a=0.5:0.5:10;
a=a';%输入数据S
b=sin(a);%输入数据Y
s=0;%精度
c=4.01:0.01:6;
c=c';%检验点
d=sin(c);
e=interp1(a,b,c,'spline');
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)