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

haidaowang

% page 151  3
% 写出一个插值多项式
% input :vector(x,y),interpolation point x0
) S% ?4 g: J& b- A: T$ p# s' \% output:the result at x = x0
function y0 = page_157_1(x ,y,x0)
format long
/ }) L6 G" G: w" T% }n = length(x);
; T; u0 Q( h. v/ N  ~9 P; ]0 K9 sfor j = 1:n     %fill in y colum of Newton triangle
2 E' v! V3 D0 x    v(j,1) = y(j);
: J3 _. \$ F' D/ }! @end
6 I* B7 u/ W1 Z9 a7 kfor i = 2:n             %for colum i
    for j = 1:n+1-i             %fill in colum from top to bottom
0 x$ O% r0 a9 A4 ^& o        v(j,i) = (v(j+1,i-1) - v(j,i-1))/(x(j+i-1)-x(j));
    end
end
4 U" h. X3 q0 I( I4 X# J+ S" bfor i = 1:n
    v(i) = v(1,i);     %read along top of triangle      
end     %output coefficients
+ W5 @3 _# r0 g: a, s* ty0 = 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
, N% U& h3 D6 D2 Z" h

youOK

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