xn=[0,1,2,4]';
yn=[3,6,11,51]';
x=0.5;
n=length(xn);
A(:,1)=yn;
for j=2:n
for i=j:n
A(i,j)=(A(i,j-1)-A(i-1,j-1))/(xn(i)-xn(i-j+1));
end
end
y=0;
for i=1:n
t=1;
for j=1:i-1
t=t*(x-xn(j));
end
y=y+A(i,i)*t;
end
y
数值计算方法牛顿插值matlab实现
最新推荐文章于 2023-11-07 21:27:46 发布