Matlab之用牛顿均差求插值多项式

haidaowang · 发表于 2020-8-27 13:55

EDA365欢迎您登录！

您需要登录才可以下载或查看，没有帐号？注册

x

% page 151  3
% 写出一个插值多项式
% input :vector(x,y),interpolation point x0
% output:the result at x = x0
function y0 = page_157_1(x ,y,x0)
format long
n = length(x);
for j = 1:n    %fill in y colum of Newton triangle
v(j,1) = y(j);
end
for i = 2:n          %for colum i
for j = 1:n+1-i          %fill in colum from top to bottom
      v(j,i) = (v(j+1,i-1) - v(j,i-1))/(x(j+i-1)-x(j));
end
end
for i = 1:n
v(i) = v(1,i);    %read along top of triangle
end    %output coefficients
y0 = v(n)*(x0 - x(n-1)) + v(n-1);  %initial nest
for i = 1:n-2                         %use nest
y0 = y0*(x0 - x(n-1-i)) + v(n-1-i);
end

youOK · 发表于 2020-8-27 14:43

Matlab之用牛顿均差求插值多项式

帐号		自动登录	找回密码
密码			注册

Matlab之用牛顿均差求插值多项式

EDA365欢迎您登录！

推荐内容 /1