计算方法第3章数据近似.docVIP

第三章 数据近似 已知函数关系的数据点：，求一个适当的函数，它是已知函数族——称为基函数——的线性组合，即：寻求系数，使 ， 与按某种标准相近似；通常按下述标准之一，使之取得极小： 若取向量形式： ,则上述三数分别为： ； 3.1 多项式插值 一般，若有个数据点，上述近似标准要求 即：，则称近似函数为满足插值条件(3.2)的插值函数，而点称为插值节点。若满足插值条件(3.2)的函数的基函数取 ，则是多项式 称满足插值条件(3.2) 的多项式为插值多项式。 代数学基本定理：n次多项式有且仅有n个根（包括重根） 插值多项式 根据插值条件(3.2)，可形成有个未知量的个方程： ， 这个方程组的系数矩阵称为Verdermonde矩阵，由其性质可知，当 时，该方程组的解存在且唯一。因此有 定理3.1：对数据点，存在满足条件(3.3)的唯一的插值多项式； 强调：插值多项式唯一 —— 形式可以不同 3.1.2 Lagrange(形式)插值多项式 由两点直线方程谈起，； 若有不超过次的多项式： ，可形成如下表： … … 1 0 0 … 0 0 … 0 0 0 0 1 … 0 0 … 0 0 … … … … … … … … … … … … 0 0 0 0 … 1 … 0 0 … … … … … … … … … … 0 0 0 0 … 0 0 … 1 … … 称为Lagrange基本插值多项式，显然它有个根：；由于次多项式有且仅有个根，因此 注意到 ，可知： ； 若记 便有 不难证明 为所求的插值多项式，它被称为Lagrange（形式）插值多项式。 缺点：增加插值节点后，Lagrange插值多项式发生大变化。 举例， 3.1.3 Newton(形式)插值多项式 定义：差商 一阶差商： ； 二阶差商： ； n阶差商： ； 考虑逐次增加插值节点后，插值多项式的变化：记满足插值条件 的（不超过）次插值多项式为，并记，即，注意：是次多项式；由 故有 ；由于次多项式有且仅有个根，因此 ； 确定系数： 同样的方法可得： 由此形成的插值多项式称作Newton插值多项式，记作： 差商-性质1、对称性：对的任意排列，有 证明：前方法的Newton插值多项式的导出过程是：从仅满足一个插值条件的插值多项式开始，逐步增加插值条件，直到最后的，形成的，它的次项的系数就是；同样，若从仅满足一个插值条件的插值多项式开始，逐步增加插值条件，直到最后的，形成的的最高次（次）项的系数就必定是。由插值多项式的唯一性可知，这两个插值多项式是同一的，从而它们的次项的系数也是同一的，由此便得结论。 差商-性质2、设有直到n阶导数，则 证明：设为插值节点，则对应的Newton插值多项式 ， 令，显然，它有个零点：；由定理可知有n个零点，同理可知有n-1个零点，依此不难推出在 中至少有一个零点，记为，即：，而，由此得结论。 例：多项式，且与点的选择无关； 例：设 3.1.4 带导数的插值多项式 若在点可导，利用差商和导数的定义，有 ， 类似，可得到 ； 上述结论也可以从差商性质2推出；例如，当时可获得上述第一个结论，当 时 可得到后一个结论。 说明：由此可知，Newton插值多项式是Taylor展开式的推广 例：求满足下述插值条件的插值多项式， 解：建立差商表： 由此可写出插值多项式 ； 注：利用Newton插值多项式，重节点差商公式，以及差商性质2，容易看到：Taylor展开式实质上是当时的Newton插值多项式。 3.1.5 插值公式的余项 若将以为插值节点插值多项式 ， 再增加一个插值节点，得插值多项式 ，由插值性质，有 ， 据此，并由所取 的任意性，可将改作，便得插值公式的余项公式： 利用差商性质2，及符号，又可以有余项的另一种表达形式 Runge现象： 问题： 答1：若有界，有界， 则可。 例：，显然 ，取区间，插值点等距分布：等分 区间，步长 ，节点 逐步向外扩充，比较：使 最小：是区间中点，等分区间 等分，则 而当或，则 所以，， 从而 若插值点非等距分布，例如 可证明：，因此当 时， 另一方面，即使 有界，当 无界时，仍可能无界。例：， 显然，当． 经典的例子：函数，区间等分份，节点；（参考后图） 注意：通常插值公式用到. 附表：的函数插值多项式值 0.84 0.86 0.88 0.90 0.92 0.95 0.96 2.6744 4.0691 3.9451 0.0471 -10.3345 -39.9524 -50