现代科学与工程计算-孟大志-第一篇1 2 2.pptVIP

★ Newton插值的算法 用牛顿插值公式，首先要计算各阶差商值，然后再计算插值。牛顿插值由承袭性质容易实现，关键是实现差商。 1. 由插值点序列{x, f(xi)}, i =0,1,2,…, n, 形成差商表 f [x0,…,xk], 求出待定系数ai 2. 形成插值基函数， 写出牛顿插值 Nn(x) ★ 差商与导数的关系 例 对 f (x) = x7?x4+3x+1， 求 f [20,21], f [x,20,21,…,26]和 f [x,20,21,…,27]。 解： 显然， f (7) (x) = 7!, f (8) (x) = 0, 由性质4得 性质1(差分与函数值的关系) 各阶差分均可表示为函值 fj = f(xj)的线性组合: 性质4(差分与差商的关系): * * * 回顾：Lagrange 插值的优缺点： 优点：具有严格的规律性,便于记忆。 缺点：计算量大、不具有承袭性。 Lagrange 插值虽然易算，但若要增加一个节点时， 全部基函数 lk(x) 都需重新算过。 为了克服这一缺点， 需建立具有承袭性的插值公式 —— Newton插值公式. 鲁东大学数 学 与 信 息 学 院 1.2.3 Newton插值 已知条件：设某个函数y = f (x) 在区间 [a, b] 上是连续的，且已知区间 [a, b] 上n +1个互异节点处的函数值，即 目的：构造一个不超过n次的Newton插值多项式 使得在插值节点处满足条件 (1.2.11) 而且在增加一个节点时插值多项式只需要增加一项。 鲁东大学数 学 与 信 息 学 院 1．插值多项式的构造 构造思想: 将插值多项式看成插值基函数的线性组合，定义符合Newton插值 多项式目的的插值基函数，利用待定系数法确定组合系数来构造 插值多项式。 为了表示组合系数的计算，需要给出差商的表示。 定义1.2.7 给定[a, b]中互不相同的点 x0, x1, x2, …, 以及 f (x) 在这些点处相应的函数值f (x0), f (x1), f (x2), …, 1阶差商 零阶差商 鲁东大学数 学 与 信 息 学 院 2阶差商 k阶差商： 依此递推法生成 可见差商的定义具有递推性，即高阶差商是两个低一阶 差商的差商。 计算过程如表1.2.1。 表1.2.1差商表 鲁东大学数 学 与 信 息 学 院 规定函数值为零阶差商 差商表的“金字塔”结构 a0 a1 a2 a3 a4 鲁东大学数 学 与 信 息 学 院 差商的基本性质： 性质2 各阶差商具有对称性,即任意调换节点的次序,差商的值不变。 鲁东大学数 学 与 信 息 学 院 性质3：（差商与导数的关系）若 在 上存在n 阶导数，且节点 ，则n阶差商与导数 的关系如下 (2) 基函数构造方法： 定义1.2.8 为了使构造的插值多项式具有承袭性，定义 为关于节点 的Newton插值基函数，它们适合节点 的增加的要求，并且显然它们是互相是线性无关的。 鲁东大学数 学 与 信 息 学 院 基本思想：利用待定系数法构造插值多项式 定义1.2.9 插值基函数固定以后，可将插值多项式表示为 插值基函数的线性组合，即 只需要确定组合系数 由插值条件(1.2.11)可得三角形方程组 ，就可以完成多项式的构造。 用从上向下依次代入的方法求解方程组。引入差商 的记法，可得组合系数为： 鲁东大学数 学 与 信 息 学 院 则 (1.2.12) 称为关于节点 的Newton插值多项式。 组合系数可以通过建立如表1.2.1所示的差商表来计算。 鲁东大学数 学 与 信 息 学 院 引入差商后我们就可以重新改写Newton插值公式 改进：增加一个节点 xn+1时，只需加一项 2．Newton插值余项 由插值多项式的存在唯一性可知 插值多项式 Ln(x) ，即 本质上就是n次 两者只是具体表现形式不同而已，故二者插值余项必定相同，即 由差商与导数的关系可得Newton插值余项公式的形式 鲁东大学数 学 与 信 息 学 院 （1.2.13） Newton插值余项不要求函数的高阶导数存在，因此更具有一般性，它对于f (x)是由离散点给出的情形或f (x)导数不存在的情形均适用。 例3 已知 的函数表 0.52110 0.30452 0.20134 0 f(xi) 0.50 0.30 0.20 0 xi 求出三次Newton插值多项式，计算f(0.23)的近似值并用均差的余