数值积分上机报告牛顿差值.pdfVIP

中国矿业大学 （北京）理学院 数值分析实验报告 实验名称 Newton 插值法计算函数值 实验时间 2014 年4 月3 组 组 长 员 一、实验目的，内容 二、相关背景知识介绍 三、代码 四、数值结果 五、计算结果的分析 六、计算中出现的问题，解决方法及体会 一、实验目的，内容  通过编程实现Newton 插值法，构造线性、二次、三次、四次Newton 插值多项式，加深 对Newton 插值法的理解；  应用所编程序解决实际算例；  通过与Lagrange 插值结果作比较，比较两种算法的结果精度。 二、相关背景知识介绍 给定插值点序列(x , f (x )), i 0,1, , n 。构造Newton 插值多项式N (u) 。输入要计算的 i i n x, 函数点 并计算N (x) 的值，利用Newton 插值公式，当增加一个节点时，只需在后面多计 n 算一项，而前面的计算仍有用；另一方面N (x) 的各项系数恰好又是各阶均差，而各阶均差 n 可用均差公式来计算。 f (x ) f (x )   1 0 为f (x) 的一阶均差。 f x ,x  0 1 x x 1 0     f x ,L,x ,x f x , x ,L,x f x ,x ,L,x  0 k2 k 0 1 k1 为f (x) 的k 阶均差。  0 1 k  x x k k1 图一：均差表 Newton 插值法计算步骤： n , (i 1,2,3) x 1.输入 值及(x ,f (x )) ，要计算的函数点 。 i i 2.对给定的 ，由 x N (x)  f (x ) (xx ) f x , x (xx )(xx ) f x , x , x (xx ) n 0 0  0 1  0 1  0 1 2  0 (xx ) (xx ) f x , x , x 1 n1  0 1 n  计算N (x) 的值。 n 3.输出N (x) 。 n Newton 图二： 插值法流程图 三、代码 问题一： , 1 2 3 4 Newton