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)
0 ~0 B" Z' @, G8 X* [format long
. B6 c: A0 u, B  C' ?n = length(x);
for j = 1:n     %fill in y colum of Newton triangle
' Q- X2 b5 b1 ?, [0 }2 ~) \    v(j,1) = y(j);
3 `" r' Q: {7 P) v# w0 t5 O& K4 cend
- z6 _8 L( r' w7 f* I7 k5 J7 ffor 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
( M2 L: b, ?  z7 Q  L8 Xend
  V0 K- D! k  _) ~! K* Mfor i = 1:n
5 U# }- Y# \" U9 d    v(i) = v(1,i);     %read along top of triangle      
8 I* Y- p1 Y6 r" s7 h0 lend     %output coefficients
, ^" s' d3 N1 W, J% _y0 = v(n)*(x0 - x(n-1)) + v(n-1);  %initial nest
, ~9 Q( {) g  e: N6 B+ vfor i = 1:n-2                          %use nest
    y0 = y0*(x0 - x(n-1-i)) + v(n-1-i);
4 ]1 U" a; T6 e5 X1 v8 P1 ]end

youOK

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