01,,n
1,,n
1,,)n x
及数值分析各算法流程图
一、插值
流程图转换为ns图1、 拉格朗日插值流程图:( 相应程序:lagrintp(x,y,xx))
2,,n ,,j n 1,2,,n 1,,)n 2、 牛顿插值流程图
(1)产生差商表的算法流程图(相应程序:divdiff(x,y))
注:1、另一程序divdiff1(x,y),输出的矩阵包含了节点向量。而divdiff(x,y)不含节点向量。    2、另一程序tableofdd(x,y,m),输出的是表格形式,添加了表头。
1,,),,n m 及1,,m (2)非等距节点的牛顿插值流程图(相应程序:newtint11(x,y,xx,m)) 、
注:1、虽然程序newtint11(x,y,xx,m)考虑了多种情形,看上去很复杂,但基本流程结构还是如上图所示。
2、程序中调用的子程序是divdiff 。
若调用的子程序是divdiff1的话,流程图中的第三,第四,第五步要相应的改一下数字。
2,3,,1m +1,,j
1,2,,n
=1,2,,)n m 及(3)求差分表的流程图(相应程序:difference(y,m))
注:1、difference 输出的是矩阵D 。而另一程序tableofd(y,m),输出的是带有表头的差分表。
(4)牛顿向前插值的流程图(相应程序:newtforward1(x,y,xx,m))
n x m
1,,),,
1,,m
注:1、程序newtforward1(x,y,xx,m))的结构与上述流程图一致,xx可以是数组。
2、另一程序newtforward(x,y,xx,m))先求出插值多项式,再求插值多项式在插值点的函数值。基本结构还是和上面的流程图一样。

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198,我们将在24小时内删除。