newpoly.m

各种插值拟合的子程序

function [C,D]=newpoly(x,y) %牛顿插值
n=length(x);
D=zeros(n);
D(:,1)=y';
for j=2:n  %计算差商矩阵
    for k=j:n
        D(k,j)=(D(k,j-1)-D(k-1,j-1))/(x(k)-x(k-j+1)); 
    end
end
C=D(n,n);
for k=(n-1):-1:1  %构造插值多项式
    C=conv(C,poly(x(k)));
    m=length(C);
    C(m)=C(m)+D(k,k);
end