计算方法第三章r.pptVIP

差商与导数的关系 重节点差商 推论：当n个节点全为同一个点，牛顿插值变成 泰勒多项式。 差商的导数 n 次多项式的的 1 阶差商是 n-1 次多项式。 推论：设 p(x) 是 n 次多项式，k ? n 时 k 阶差商是 n - k 次 多项式，k n 时 k 阶差商为零。 差分 设函数 ，定义 为该函 数在 i 点的一阶差分，记为 类似地，定义二阶差分为： K 阶差分为： 此差分称为向前差分。 类似地，向后差分定义为： 中心差分定义为： 差商与差分的关系： 等距节点时 第五节 带导数的插值 问题的提出：如果在已知节点处不仅知道函数值， 同时还指导导数值，这样，插值多项式就要求在 已知节点处与函数值和导数值都相等。这就是所 谓埃尔米特插值，记为 H (x) 1、牛顿插值 如果已知某个点 i 的 ，则 插值节点应视为 个相同节点 ，并注意到k+1 重节点的差商 例子：已知关于函数 y = f (x)的函数值、导数值 -1 0 0 -4 -4 0 -4 0 4 0 -4 0 3 -1 1 -2 2 2 -1 0 1 -2 5 3 1 2 1 已知函数在n个不同的节点处的函数值和导数值： 求次数不超过2n-1次的多项式 设想该多项式具有形式： 由条件可得： 此外，由 得： 同理： 由 ，可得： 最后，得到埃尔米特插值公式： 特别，当 n=2 时，三阶埃尔米特多项式为： 埃尔米特插值公式唯一。 误差估计，设被插值函数在插值区间上2n次连续可导，则在n个节点上的2n-1次插值多项式的余项为： 特别，对于2个节点3次插值，余项为： 例子： 如用距离较小的两个点插值，效果会好得多 第六节 样条函数 由于被插值函数高阶导数未知，因此，如果高阶导数随阶数 增长出现无限增长，则由误差公式可知，高阶插值公式就不 一定无限接近被插值函数。这称为龙格现象。 所以，在进行多项式插值时，不宜进行高次多项式插值。 一个解决的途径是分段低次插值。 龙格现象 样条函数插值：给定区间一个划分 如函数 S(x) 满足下面条件： （1）在每个小区间 上为m次多项式； （2） S(x) 及m－1 阶导数在整个区间上连续。 则称 S(x) 是关于该划分的 m 次样条函数，划分点称为节点，m＝3 时，就是最常用的 3 次样条函数。 3次样条函数的基本思想： 将样条函数在每一个子区间端点的二阶导数值当作参数，则用这两个二阶导数值可以将样条函数表示出来，再利用衔接条件，即每一段样条函数在相邻两个子区间端点处的二阶导数相等，建立求解二阶导数的方程组。 设S(x)在每个小区间 端点的二阶导数为： 则： 记 ， 将上式积分两次，并利用端点函数值已知，有： 我们注意到，在相邻的两个子区间 和 的共同端点处，样条函数一阶导数相等， 经过化简，最后得到： 注意到上面的方程组总共只有 N－1个方程，而未知数却共有 N 个，因此，要求解方程，还需要两个条件（即两个方程），通常有以下几种方案: 1、给出端点一阶导数值 ，这相当于增加两个方程； 2、给定端点二阶导数值，得方程： 特别，令二阶导数在端点为零，得 3、样条件函数在第一后最后一个区间上为二次多项式，即样条函数在第一和最后一个区间上的二阶导数为常数，得两个方程： 总之，可以将方程组统一写为： 4、周期性条件（这只有在给的初值满足 时才能用），此时，由周期性， ，就得到两个方程；第一个方程为 ，第二个方程为 最后一个方程为： 最后，方程组为： 小结 插值法，其目的是利用节点上的值，构造通过这些节点的多项式，从原则上说，利用n+1个节点的值，可以构造n次多项式，而且这种构造是唯一的！利用待定系数法，将节点值带入后，得到一个n+1阶线性方程，即可求出多项式。 为便于在计算机上实现，本章引入了