zxecexample.m

来自「数值线性代数的Matlab应用程序包 共13个程序函数」· M 代码 · 共 11 行

x=1:10;
y=[2,3,5,7,10,2,9,13,4,5];
N=9;
xp=x(1):0.01:x(length(x));
yp=zeros(1,length(xp));
p=ZXEC(x,y,N-1);
for i=1:N
  yp=yp+p(i)*xp.^(N-i);
end
plot(x,y,'*',xp,yp);