第5章 3.牛顿插值公式.ppt

牛顿插值公式 邹昌文 问题的提出 一般地构造以下基函数 问题 求作n次多项式 使满足 均差/ divided difference/的定义 牛顿基本插值公式 * (2) (1) 为了得到计算系数ci的一般方法,下面引进一般均差的概念. 1 1 1 0 1 0 1 0 1 1 1 0 1 0 ] , , ... , [ ] , , ... , [ ] , , ... , [ ] , ... , , [ ] , ... , [ + + - - + + + - - = - - = k k k k k k k k k k k x x x x x f x x x f x x x x x f x x x f x x f 均差表 …… ┊ ┊　 ┊　 ┊　 ┊　 ┊ ? f(x0,x1,x2,x3,x4 ) f(x1,x2,x3,x4 ) f(x2,x3,x4 ) f(x3,x4 ) f(x4) x4 ? ? f(x0,x1,x2,x3 ) f(x1,x2,x3 ) f(x2,x3 ) f(x3) x3 ? ? ? f(x0,x1,x2 ) f(x1,x2 ) f(x2) x2 ? ? ? ? f(x0,x1 ) f(x1) x1 ? ? ? ? ? f(x0) x0 ? 4阶 3阶 2阶 1阶 f(xk) xi 计算规律：任一个k(≥1) 阶均差的数值等于一个分式的值，其分子为所求均差左侧的数减去左上侧的数，分母为所求均差同一行最左边的基点值减去由它往上数第k个基点值。 注意：均差表中，对角线上的均差是构造牛顿型插值公式的重要数据。 粗线框出的部分在计算机上可存入二维数组 例 已知函数y=f(x)的观测数据如表，试构造差商表，并求f[2,4,5]及f[2,4,5,6]的值。 1 5 9 - 4 13 f(x) 0 2 4 5 6 x 解 n=4, 构造差商表 xi f(xi ) 1阶 2阶 3阶 4阶 0 2 4 5 6 2 1 1 5 9 -4 13 2 -13 17 0 -5 15 -1 5 f[2,4,5]= -5 f[2,4,5,6]=5 1 2 … … … … n?1 1 + (x ? x0) ? 2 + … … + (x ? x0)…(x ? xn?1) ? n?1 Nn(x) Rn(x) ai = f [ x0, …, xi ] *