计算方法第三章rr绪论.pptx

第三章 插值法第一节 插值多项式的基本概念;假设已经获得若干点上的函数值 即提供了一张数据表 ? 如何利用这张表求某个给定点上的函数值呢?插值方法所要研究的就是这个课题。 ;通常用多项式来作为近似函数，称为插值多项式。 ;因式定理：多项式P(x)具有r 次因式 (x-a)r 的充要条件 是P(a)= P ‘(a) = ……= P(a)(r-1) =0 最一般的插值条件： 是 重根，;设函数 y = f (x) 在闭区间 [a , b ]上有n + 1 阶导数， 满足前面的一般插值条件，且插值节点各不相同， 则插值截断误差为;证明思路：构造辅助函数，用罗尔定理。;值得注意的是在较大区间上进行插值时，误差可能会 很大！另外，一般情况下，外推不如内插好！ ;第二节 Lagrange插值公式 插值条件是 ;一次插值： 问题为求一次多项式，即一次函数，过以下两点： 容易求出，该函数为： ;二次插值： 问题为求二次多项式，即二次函数，过以下三点： 容易求出，该函数为： ;一般插值问题：求过n+1个点 的不超过n次多项式 。 称为Lagrange插值基函数，满足： ;问题:过n+1个点的Lagrange插值多项式是否唯一？;插值公式的误差为： ;计算程序 框图;第三节 逐次线性插值 函数 y = f (x)在节点 上的插值多项式记为 ，则有;Aitken（埃特肯）算法 Neville（列维尔）算法;Aitken（埃特肯）算法 ;Neville（列维尔）算法;例子：求方程 在(2 , 3)内的根 思路，用反函数;第四节 牛顿插值 ;差商 设函数 f (x) ，定义函数在两个不同点的一阶差商为 三个不同点的二阶差商为： 在点 处 K+1 阶差商为： ;给定 n +1个点的函数值，则牛顿插值公式为： ;差商的计算简表：;例子：用0、30、45、60、90五个点作出sinx 牛顿插值多项式。 做差商表;;牛顿插值的截断误差：;例子：用0、30、45、60、90五个点作出sinx 牛顿插值多项式。 做差商表 ;;差商的计算公式： 差商的对称性： 差商的线性;由于n次插值多项式是唯一的，所以牛顿插值公式与Lagrange 插值多项式一样，这意味着余项也一样，Lagrange余项为： 所以牛顿余项也一样，;差商与导数的关系 重节点差商 推论：当n个节点全为同一个点，牛顿插值变成 泰勒多项式。;差商的导数 n 次多项式的的 1 阶差商是 n-1 次多项式。 ;差分 设函数 ，定义 为该函 数在 i 点的一阶差分，记为 类似地，定义二阶差分为： K 阶差分为： 此差分称为向前差分。;类似地，向后差分定义为： 中心差分定义为： ;差商与差分的关系：等距节点时; 第五节 带导数的插值 问题的提出：如果在已知节点处不仅知道函数值， 同时还知道导数值，这样，插值多项式就要求在 已知节点处与函数值和导数值都相等。这就是所 谓埃尔米特插值，记为 H (x) ;1、牛顿插值 如果已知某个点 i 的 ，则 插值节点应视为 个相同节点 ，并注意到k+1 重节点的差商;例子：已知关于函数 y = f (x)的函数值、导数值;已知函数在n个不同的节点处的函数值和导数值： 求次数不超过2n-1次的多项式 设想该多项式具有形式： ;由条件可得： 此外，由 得：;同理： 由 ，可得： 最后，得到埃尔米特插值公式：;特别，当 n=2 时，三阶埃尔米特多项式为：;埃尔米特插值公式唯一。 误差估计，设被插值函数在插值区间上2n次连续可导，则在n个节点上的2n-1次插值多项式的余项为： 特别，对于2个节点3次插值，余项为：;例子： ;;如用距离较小的两个点插值，效果会好得多;;第六节 样条函数 由于被插值函数高阶导数未知，因此，如果高阶导数随阶数 增长出现无限增长，则由误差公式可知，高阶插值公式就不 一定无限接近被插值函数。这称为龙格现象。 所以，在进行多项式插值时，不宜进行