第五章-插值法.pptVIP

三次样条插值函数 在每一个小区间上是3次的多项式(4个系数), 在整个插值区间上有4n个系数. 且有4n-2个约束: 内节点 要确定4n个系数，还需附加2个约束条件. 常用的约束条件有以下三类： ②.弯矩边界条件 特别的称 为自然边界条件. ①.转角边界条件 2. 三次样条函数造法 记 , 基本步骤如下： ①.取 为待定参数，并用S(x)的插值条件写出 的表达式。 ④.代入S(x)的表达式，得各个区间上的表达式。 ②.用 在内节点 的连续条件及边界条件导出关于 的方程组。 ③.求解后得到 。 加以整理后可得 由条件 由于以上两式相等,得 基本方程组 其中 如果问题要求满足第一类(一阶)边界条件: 基本方程组化为n-1阶方程组 将上式化为矩阵形式 这是一个三对角方程组 矩阵形式 如果问题要求满足第二类(二阶自然)边界条件: 都是三对角方程组,并且都严格对角占优 可以使用追赶法求解,并且解是唯一的 3. 差商的计算 根据插值多项式的存在唯一性知，如果f(x)充分光滑，则有估计 不足： 对函数的光滑性要求高； 需估计导函数的最值； 偏保守。 导数型误差估计 三、Newton插值余项 差商型误差估计 导数和差商的关系 差商型误差估计特点：对被插值函数光滑性要求不高；但不适用于实际计算。 四、例题 解 1）建立差商表 1.0 1.5 2.0 0.8415 0.9975 0.9093 0.312 -0.1764 -0.4884 2）插值 Newton插值多项式适用于节点任意分布的情形。但当节点等距分布时，可以简化Newton插值公式。 §5.3 等距节点插值 设a=x0x1… xn=b，yi=f(xi)为等距节点xi=x0+h(i=0,1,…,n)上的函数值，其中h=(b-a)/n称为步长。 在此基础上我们先定义差分，用差分表示Newton插值多项式，从而得到等距节点的插值公式。 一、差分的定义与性质 定义：称 ?yi=yi+1-yi (i=0,1,…,n-1) 为f(x)在xi处以h为步长的一阶向前差分。 ?2yi＝?yi＋1-?yi =yi+2-2yi+1+yi (i=0,1,…,n-2) 称为f(x)在xi处以h为步长的二阶向前差分。 一般地，?myi＝?m-1yi＋1-?m-1yi (i=0,1,…,n-m) 称为f(x)在xi处以h为步长的 m 阶向前差分。 差分的性质 性质1：各阶差分可用函数值线性表示，其计算公式为： 其中 性质2：差分与差商满足下述关系： 证明：利用数学归纳法 当k＝1时，有 即结论成立。 设k＝m-1时结论成立，即 则当k＝m时，有 由数学归纳法知，结论成立。 证毕 Remark：类似地可以定义向后差分与中心差分： 性质3：差分与导数满足关系： 证明：利用差商与导数、差分的关系，有： 证毕 二、Newton向前插值公式 令x=x0+th，由xi=x0+ih(i=0,1,…,n)得： x-xi=(t-i)h，则有: 将差商与差分的关系式 带入Newton插值多项式，得: 从而可得Newton向前插值多项式及其余项为： 三、差分表 Newton向前插值公式，又称表初公式，它利用差分表的最上面一个斜行的数值进行计算。 四、例题 解 # 五、Newton向后插值公式(不讲) 类似于向前差分，也可以得到差商与向后差分的关系： 将插值节点从大到小排列，即 类似于向前插值公式，可得到Newton向后插值公式，又称表末公式，它利用差分表的最下面一个斜行的数值进行计算。 同样，还可以利用中心差分，构造插值公式，称为贝塞尔（Bessel）插值公式。 这一类插值问题为埃尔米特(Hermite)插值问题。其几何意义是在插值点上插值曲线与被插值曲线有公共切线。由这2n+2个条件可以唯一确定一个2n+1次的插值多项式。具体我们采用基函数的方法来确定。 §5.4 埃尔米特插值 一、问题 1.辅助问题及Hermit插值 二、一般情形