数值分析r研.ppt

牛顿插值多项式 解 xk f (xk) 一阶均差 二阶均差 2.0 1.414214 2.1 1.449138 0.34924 2.2 1.483240 0.34102 -0.04110 Newton插值公式： 取x=2.15，得N2(2.15)=1.466292. 利用误差估计公式，注意到 * 牛顿插值多项式 x x0 x1 x2 x3 f(x) f(x0) f(x1) f(x2) f(x3) 求各阶差商 差商表 k xk 0 1 2 3 x0 x1 x2 x3 f[xk] f[xk , xk+1] f[xk , xk+1 , xk+2] f[xk , xk+1 , xk+2 , xk+3] f[x0] f[x1] f[x2] f[x3] f[x0, x1] f[x1 , x2] f[x2 , x3] f[x0 , x1 , x2] f[x1 , x2 , x3] f[x0 , x1 , x2 , x3] 牛顿插值多项式 由各阶差商的定义 一阶差商 二阶差商 同理，得 …… 将以上各式分别乘以 然后将各式相加，并消去左右两边相等的部分，得 牛顿插值多项式 其中，Nn(x)是至多n次的多项式。 牛顿插值多项式 牛顿插值多项式 由插值多项式的唯一性知，n次牛顿插值多项式与拉格朗日多项式是相等的，即 牛顿插值的优点：每增加一个节点，插值多项式只增加一项，即 便于递推运算。 以上推导过程只强调节点位置的不同，各节点按任意大小排列 1）列表计算各阶均差，如下表 k xk 0 1 2 3 x0 x1 x2 x3 f[xk] f[xk , xk+1] f[xk , xk+1 , xk+2] f[xk , xk+1 , xk+2 , xk+3] f[x0] f[x1] f[x2] f[x3] f[x0, x1] f[x1 , x2] f[x2 , x3] f[x0 , x1 , x2] f[x1 , x2 , x3] f[x0 , x1 , x2 , x3] 2）将表中带有下划线的三角形的顶边项与该行最后一列的对应项相乘后，再相加，即得Newton插值多项式。 作出Newton插值多项式的步骤 输入数据 计算各阶差商 对 ，计算 (3) (4) 输出 ，停机。 Newton插值多项式的算法 插值多项式为 例5.5：利用差商求通过点(0, 1), (2, 2), (3, 4)的插值多项式。 解：构造差商表 一阶差商 二阶差商 0 1 1 2 2 0.5 3 4 2 0.5 读出多项式的系数1, 0.5, 0.5。插值多项式可以写出： 牛顿插值多项式 例5.6：在例5.5中，增加第4个数据点(1,0) ，构造三次插值多项式。 解：在差商表中，增加一个数据点，可得 差商均差 二阶差商 三阶差商 0 1 1 2 2 0.5 3 4 2 0.5 1 0 2 0 -0.5 牛顿插值多项式 牛顿插值多项式 x0 f[x0] x1 f[x1] f[x0,x1] x2 f[x2] f[x1,x2] f[x0,x1,x2] x3 f[x3] f[x2,x3] f[x1,x2,x3] f[x0,x1,x2,x3] 例5.7 由已知节点x0, x1, x2 , …, xm求得m次插值多项式Nm(x)，现增加一个任意位置的节点xm＋1 ，求m＋1次插值多项式Nm＋1(x)。 解 差商表（m=3） x4 f[x4] f[x3,x4] f[x2,x3,x4] f[x1,x2,x3 ,x4] f[x0,x1,x2,x3 ,x4] 牛顿插值多项式 例5.8 x 30o 45o 60o sin(x) 1/2 求二次牛顿插值多项式，并计算sin(50) 。 解： k xk f[xk] f[xk, xk+1] f[xk, xk+1 , xk+2] 0 30 1 45 2 60 一阶差商: 二阶差商: 牛顿插值多项式 二次牛顿插值多项式为 N2(50)=0.76543 经过n+1个点的d次插值多项式 经过n+1个点的d次插值项式到底有多少个？ 时，插值多项式的存在唯性定理已经回答了这一问题，即存在一个且仅有一个次数不超过n的插值多项式。 当dn时，有几个插值多项式存在？ 插值多项式 考虑例5.6，在例5.5中给出了3个点(0, 1), (2, 2), (3, 4)，求出了一个唯一的2次插值多项式。 增加了一个点(1,0)之后，可以得到一个三次插值多项式，该三次插值多项式显然经过前3个点。 改变第4个点，三次多项式不同 （无穷多个