计算方法牛顿插值法

牛顿插值法

一:实验目的：

1.matlab中多项式的表示及多项式运算

2.用matlab实现牛顿插值法

二:实验代码

(一)function [p]=Newton_Ployfit(X,Y)

if size(X) ~= size(Y)

error;

end

format long g

r=size(X);n=r(2);

M=ones(n,n);

M(:,1)=Y';

for i=2:n

for j=i:n

M(j,i)=(M(j,i-1)-M(j-1,i-1))/(X(j)-X(j-i+1)); end

end

M

p0=[zeros(1,n-1) M(1,1)];p=p0;

for i=1:n-1

p1=M(i+1,i+1).*poly(X(1:i));

p0=[zeros(1,n-i-1) p1];

p=p+p0;

end

(二)x0=linspace(0,2*pi,10);

y0=sin(x0);

p=Newton_Ployfit(x0,y0);

x=0:0.2:2*pi;

y1=sin(x);

y2=polyval(p,x);

plot(x,y1,'co',x,y2,'r');

三：流程图

四:实验结果(一)

(二)