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

haidaowang

% page 151  3
% 写出一个插值多项式
% input :vector(x,y),interpolation point x0
% output:the result at x = x0
function y0 = page_157_1(x ,y,x0)
6 z% k) w2 \1 c5 S. |) |2 zformat long
n = length(x);
/ ~2 A% R7 M' r  l7 @% `  g1 ^for j = 1:n     %fill in y colum of Newton triangle
6 w; X+ ?9 z( X6 B- p  E    v(j,1) = y(j);
$ |5 Y5 K+ y6 C$ k2 a- g6 `+ A& l4 @end
for i = 2:n             %for colum i
    for j = 1:n+1-i             %fill in colum from top to bottom
$ y& N, @3 D! k  Q5 W) R0 ~        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      
6 A/ I( W7 n7 [end     %output coefficients
8 A5 p+ D# F4 K0 W' b0 |: Zy0 = v(n)*(x0 - x(n-1)) + v(n-1);  %initial nest
6 i' B) @) a0 o4 X# Nfor i = 1:n-2                          %use nest
: b! Q% h! x6 z3 ]( J    y0 = y0*(x0 - x(n-1-i)) + v(n-1-i);
end
& n( \+ h9 L# @$ p

youOK

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