数值计算CH3插值法与最小二乘法—3.4Newton插值法.ppt

华长生制作 在等距节点的前提下,差商与差分有如下关系: 差商与差分的关系 依此类推： 4-4 差分表示的 Newton 插值公式 1. Newton向前(差分)插值公式及余项公式 记插值点： 如果节点 是等距的, 即 ------(7) 由均差与向前差分的关系： -----(6) 以及 式(7)、(6)代入差商表示的插值公式： -----(8) Newton向前(差分)插值公式 得差分表示的 Newton 插值公式化为： Newton 插值多项式的余项公式 节点等距时, 化为 ------(9) 若由差分近似代替导数,上式可化为： ------(10) 实用的估计式： 2. Newton向后(差分)插值公式及余项公式 则： Newton 向后插值基本公式为： 当 x 位于表末时, 设： Newton向后（差分）插值公式： 同时可得 Newton向后（差分）插值余项公式： 导数形式余项公式： 差分形式余项公式： 实用的估计式： 例3: 给定 f (x) = cos x 的函数表如下： 用四次Newton插值多项式计算cos0.048及cos0.566的近似值，并估计误差. 0.82534 0.87758 0.92106 0.95534 0.98007 0.99500 1.0000 f(xk) 0.6 0.5 0.4 0.3 0.2 0.1 0.0 xk 6 5 4 3 2 1 0 k 先构造差分表, 解 再根据插值点在表中的位置选择插值节点. 0.0 0.1 0.2 0.3 0.4 0.5 0.6 1.0000 0.99500 0.98007 0.95534 0.92106 0.87758 0.82534 -0.00500 -0.01493 -0.02473 -0.03428 -0.04348 -0.05224 -0.00993 -0.00980 -0.00955 -0.00920 -0.00876 0.00013 0.00025 0.00035 0.00044 0.00012 0.00010 0.00009 -0.00002 -0.00001 0.00001 0.048 分别需要几个插值节点？如何选择？ 0.566 显然, h = 0.1, 当 x = 0.048 时, 用前向公式, * 第三章 插值法与 最小二乘法 张红梅 自动化学院 2010年3月 3.4 Newton 插值法 Lagrange 插值多项式的基函数： 优点：形式对称，有很强的规律性，便于记忆。 缺点： (1) 重复计算多, 导致计算量大； (2) 插值基函数 lj(x) 依赖于所有节点，当增加插值节点时，原来已算出的所有 lj(x) 都需要重新计算，使计算量加大。 ？ 是否存在其它的插值法，可以解决上述问题呢 ? Newton插值基函数 差商的定义和性质 差商表示的Newton插值多项式 差分的定义和性质 差分表示的Newton插值多项式 3.4 Newton 插值法 Newton (1624～1727) 由线性代数知识可知, 任何一个 n 次多项式都可表示成: 这n+1 个多项式的线性组合. 问：是否可以将这 n+1个多项式作为插值基函数？ 已知函数 f(x) 的插值节点 xi 及相应函数值 , 将上述线性无关的多项式取作Newton插值法的基函数,即令： Newton插值基函数 则相应的插值多项式为: 式中, 为待定参数,它们可利用插值条件 来求,即令： 可以求得： 依此类推,可求得 .为标记、推导、记忆方便,给出均差定义,可得参数 的一般表达式。 这也太复杂了吧！ 为 关于节点 的一阶均差(差商) 4-1 均差(差商) 1.均差的定义:设给定函数 在 个互异的节点 处的函数值为 ,称 为 关于节点 的二阶均差(差商) 缺倒数第二个节点 缺最后一个节点 最后一个节点－倒数第二个节点 称 可见：一个高阶差商可由两个低一阶的差商得到 缺倒数第二个节点 缺最后一个节点 称 为 关于节点 的 k 阶均差(差商) 最后一个节点－倒数第二个节点 由此定义,显然: 用归纳法可证： 将上述结果代入： 2.差商的性质 性质1：差商与函数值的关系 f(x) 关于 的 k 阶差商是 f(x) 在这些点上函数值的线性组合,即 例如: 利用对称性,可对 f(x) 关于