计算方法课件 第五章 插值法.ppt

三、几何意义、内插法、外插法 四、多项式插值问题 对于不同的函数族Φ的选择，得到不同的插值问题 当Φ为一些三角函数的多项式集合时:三角插值; 当Φ为一些有理分式集合时：有理插值； 当Φ为一些多项式集合时：多项式插值（代数插值） 五、插值多项式的存在唯一性 分析 对于多项式插值问题，插值条件（1）等价于确定多项式的系数，使得满足如下的线性方程组 六、插值余项 引理 已知函数f x 在[a,b]上具有m-1阶连续导函数，且在（a,b）上存在m阶导数。 若它在该区间上有m+1个零点，则它的m阶导函数在 a,b 内至少存在一个零点。 分析： 二、Lagrange型插值公式 例题 四、反插值法 问题求解 Lagrange 插值公式的特点： 形式对称 通常用于理论分析 当增加插值节点时，在计算实践中不方便 1.Lagrange插值多项式间的关系 2.Newton型插值公式 3. 差商的计算 差商型误差估计特点：对被插值函数光滑性要求不高；但不适用于实际计算。 解 1）建立差商表 一、差分的定义与性质 二、Newton向前插值公式 令x x0+th，由xi x0+ih i 0,1,…,n 得： x-xi t-i h，则有: 三、差分表 四、例题 1.辅助问题及Hermit插值 2.辅助问题的求解 3.Hermite插值问题解的存在唯一性 4.插值余项 定理证明 三、特殊情形－带不完全导数的插值问题举例 差商型误差估计 导数和差商的关系 四、例题 1.0 1.5 2.0 0.8415 0.9975 0.9093 0.312 -0.1764 -0.4884 2）插值 Newton插值多项式适用于节点任意分布的情形。但当节点等距分布时，可以简化Newton插值公式。 §5.3 等距节点插值 设a x0 x1 … 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插值多项式，得: 从而可得Newton向前插值多项式及其余项为： Newton向前插值公式，又称表初公式，它利用差分表的最上面一个斜行的数值进行计算。 解 # 五、Newton向后插值公式 类似于向前差分，也可以得到差商与向后差分的关系： 将插值节点从大到小排列，即 类似于向前插值公式，可得到Newton向后插值公式，又称表末公式，它利用差分表的最下面一个斜行的数值进行计算。 同样，还可以利用中心差分，构造插值公式，称为贝塞尔（Bessel）插值公式。 这一类插值问题为埃尔米特 Hermite 插值问题。其几何意义是在插值点上插值曲线与被插值曲线有公共切线。由这2n+2个条件可以唯一确定一个2n+1次的插值多项式。具体我们采用基函数的方法来确定。 §5.4 埃尔米特插值 一、问题 二、一般情形 存在性： 唯一性： 0 分析： 函数 零点（从小到大） …… …… 至少2n+1个零点 …… 至少1个零点ξ 证毕 分析（方法1）： 误差： # 方法2：（用带有重节点的差商表） # # 第五章 插值法 §5.0 插值问题 §5.1 拉格朗日插值 §5.2 牛顿插值 §5.3 等距节点插值 §5.4 埃尔米特插值 §5.5 三次样条插值 1 函数表达式过于复杂不便于计算, 而又需要计算许多点处的函数值 2 仅有几个采样点处的函数值, 而又需要知道非采样点处的函数值 …… 上述问题的一种解决思路：建立复杂函数或者未知函数的一个便于计算的近似表达式. 解决方法－插值法 §5.0 插值问题 一、问题提出 二、插值问题定义 求插值函数? x 的问题称为插值问题。 内插 外插 定理1