sinpolyfit.m

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

function[a,b,s]=sin(x)
a=0.5:0.5:10;
a=a';%输入数据S
b=sin(a);%输入数据Y
s=0;%精度
[p,m]=polyfit(a,b,3);
[f,deta]= polyval(p,a,m);
plot(a,b,'o',a,f,'-')
d=f-b;
n=10/0.5;
for i=1:n
    if s<=d(i)*d(i)
       s=d(i)*d(i);
    end
end
s=sqrt(s)