lq.m

Lagrange插值、Newton插值函数等插值函数

function scyty=lq(x,sx,ys,h,m)
c1=clock;
k=1;
for i=1:size(x,2)
    for j=k:size(sx,2)-1
        if x(i)<=sx(j+1)
            break;
        end
    end
    scyty(i)=m(j)*(sx(j+1)-x(i))^3/6/h(j)+m(j+1)*(x(i)-sx(j))^3/6/h(j)+(ys(j)-m(j)*h(j)^2/6)*(sx(j+1)-x(i))/h(j)+(ys(j+1)-m(j+1)*h(j)^2/6)*(x(i)-sx(j))/h(j);
    k=j;
end
c2=clock;
c=c2-c1;
time=c(6)