牛顿插值法 second version

最新推荐文章于 2022-11-11 16:59:04 发布

原创最新推荐文章于 2022-11-11 16:59:04 发布 · 165 阅读

0 ·

本内容遵循CC 4.0 BY-SA版权协议

收录于

这段代码展示了数值计算过程，包括生成数列、计算累积乘积、构造矩阵、求矩阵差商，并最终进行数值积分的近似计算，最后绘制了计算结果的图形。

for i=1:10
    x(i)=i+1993
end
y=[67.052,68.008,69.803,72.024,73.400,72.063,74.669,74.487,74.065,74.777]
m=2010
p(1)=1
for i=2:length(y)
    p(i)=p(i-1)*(m-x(i-1))
end
disp(p(2))
Q=zeros(length(y))
for i=1:length(y)
    Q(i,1)=y(i)
end
for i=2:length(y)
    for j=2:i
        Q(i,j)=(Q(i,j-1)-Q(i-1,j-1))/(x(i)-x(i-j+1))
    end
end
disp(Q)
t=0

for i=1:length(y)
    t=t+Q(i,i)*p(i)
end

disp(vpa(t,5))
x(length(y)+1)=2010
y(length(y)+1)=t
plot(x,y)