牛顿插值先求出差商表,再去利用差商表的对角线的值去求多项式。
function [A,C,L,wcgs,Cw]= newploy_final(X,Y)
n=length(X);%向量X中元素xi的个数,此处n=5
A=zeros(n,n); %把变量A定义为一个n行n列的零矩阵用来存储各阶差商
A(:,1)=Y'; %将Y转置赋值给零矩阵A的第一列
s=0.0; p=1.0; q=1.0; c1=1.0;
for j=2:n
for i=j:n
A(i,j)=(A(i,j-1)-A(i-1,j-1))/(X(i)-X(i-j+1));
end
b=poly(X(