newpoly.m

来自「一些matlab程序」· M 代码 · 共 19 行

function [c,d]=newpoly(x,y)
clear all
x=[1 2 3 4 5]';
y=[1 4 7 8 6];
n=length(x);
d=zeros(n,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)=d(k,k)+c(m);
end