系数.txt

来自「数值计算方法--拉格朗日插值法和牛顿插值法」· 文本 代码 · 共 30 行

% 拉格朗日差值多项式系数

x=[0 0.1 0.2 0.3 0.4];
y=[0.5 0.5398 0.5793 0.6179 0.6554];
n=5;
for i=1:n
 l(i)=y(i);
  for j=1:n
   if i~=j
   l(i)=l(i)/(x(i)-x(j));
   end;   
  end;
end;


%牛顿差值多项式系数

x=[0 0.1 0.2 0.3 0.4];
y=[0.5 0.5398 0.5793 0.6179 0.6554];
n=5;
nn=y;
for j=2:n
for i=n:-1:j
nn(i)=(nn(i)-nn(i-1))/(x(i)-x(i-j+1));
end
end