34Newton插值.docVIP

§4　Newton插值 在前面我们知道不同的插值基函数，会有不同形式的插值多项式，本质是一样的，计算方法不同，可举两种基函数为例，还有其它基函数吗？(或从时的点斜式出发)从Taylor展开的形式出发，引导出下列基函数： 对于给定的插值节点，以及节点处的函数值，设法找出个系数，使次多项式 成为函数的插值多项式，即 具体求出 ， ， ， ， 从中引出均差的概念 4－1均差 定义4.1　设在互异的节点上的函数值，称 为关于的一阶均差（差商），称 为关于的二阶均差（差商）。一般地，称 为关于的阶均差（差商）。 两个方面：1　前面插值多项式的系数可用均差的概念来表达；2　对于均差本身来说，有什么性质？定义式有什么问题？如何计算？下面先解决第二个问题。 定理　（1）:关于的阶均差是在这些点上的函数值的线性组合，即 （2）阶均差关于是对称的，即与节点的排列次序无关。 且 可知在定义阶均差时，两个均差的选取可以是任意的：（第1项缺，第2项缺） 或有 均差表 一阶差商 二阶差商 三阶差商 四阶差商 例：已知节点及其函数值如表（注意分母），计算各阶差商。 一阶差商 二阶差商 三阶差商 四阶差商 五阶差商 0.40 0.41075 1.116 0.55 0.57815 0.28 1.186 0.1973333 0.65 0.69675 0.358933 0.031238095 1.275733 0 2.9304029×10-4 0.80 0.88811 0.433466 0.031428571 1.3841 0 0.90 1.02652 0.524933 1.515333 1.05 1.25382 4－2　Newton插值公式及其余项 现在回答第1个问题：对于给定的插值节点，以及节点处的函数值，设法找出个系数，使次多项式 成为函数的插值多项式，由插值条件，知 ，得 ， 得 …… 得 可知应取 ， （记）于是得到满足插值条件的次多项式： 称为的次Newton插值多项式。所需的差商（）在表中上对角线上。 用Newton插值法也应当避免使用高次插值多项式，类似Lagrange插值多项式，可用分段插值方法来减少误差。 例1：写出满足上例中给出的插值条件的5次Newton插值多项式，用分段三次Newton插值多项式计算的近似值。 解：5次插值多项式为 计算的近似值时，选取节点，得三次插值多项式为 选取节点，得三次插值多项式为 求Newton插值多项式时，节点的顺序无关，若把节点视为按的次序排列，则牛顿插值多项式可写为 其中系数（均差）在表中次对角线上。 例：利用上例数据，用三次Newton插值多项式计算的近似值。 或 故 讨论Newton插值多项式与Lagrange插值多项式的异同：两者本质是一样的，但当增加一个节点时，Lagrange插值多项式需全部重新计算，而从上面计算差商的过程来看，Newton插值多项式只需增加一项，原来的计算结果还可利用。 下面讨论Newton插值多项式的余项，由于插值多项式的唯一性，Lagrange插值余项与Newton插值余项应相等，故可利用导数界来估计误差。能否利用差商来估计误差呢？ 注意到直接演算得到 称为Newton插值多项式的余项，应有 （导数存在时） 从而有 ，在之间 把视为一个节点，上述式子可推广到一般情形 上述余项公式用于估计误差不可行，因为用到待求值，但若阶差商接近于一个常数时，高于阶均差就会接近于零，从而可用下式作为误差估计式 4－3　差分 在实际应用Newton插值多项式时，经常遇到插值节点是等距的，即个插值节点： ，其中为步长。于是在差商中，分母部分将变得简单，计算量主要集中在分子（两节点处函数值的差）。分析差商的形式，引入差分概念 定义　设在等距节点上的函数值为，称 ， 分别为在处的一阶向前差分和一阶向后差分。符号分别称为向前差分算子和向后差分算子。 类似地，可定义在处的阶向前差分和阶向后差分 ， 可具体展开2阶差分增加感性认识（与哪些有关）。 向后差分算子和向前差分算子的关系如下： 易知，反复应用该式，不难得到 以为例，下面给出计算向前差分的差分表