2011年峨眉山市事业单位公开招聘工作人员岗位表xls.ppt

数值分析 第二章 插值法 插值法 插值法是一种古老的数学方法，早在一千多年前的隋唐时期定制历法时就广泛应用了二次插值。刘焯将等距节点的二次插值应用于天文计算。 插值理论却是在17世纪微积分产生后才逐步发展起来的，Newton插值公式理论是当时的重要成果。 由于计算机的使用以及航空、造船、精密仪器的加工，插值法在理论和实践上都得到进一步发展，获得了广泛的应用。 插值法 §2.1 引言 §2.2 拉格朗日插值 §2.3 均差与牛顿插值公式 §2.4 差分与等距节点插值 §2.5 埃尔米特插值 §2.6 分段低次插值 §2.7 三次样条插值 * * 能否存在一个性能优良、便于计算的函数 一、插值问题 ------(1) 这就是插值问题, (1)式为插值条件, 其插值函数的图象如图 整体误差的大小反映了插值函数的好坏 为了使插值函数更方便在计算机上运算,一般插值函 数都使用代数多项式和有理函数。 x0 x1 x2 x3 x4 x P(x) ? f(x) 二、代数插值多项式的存在唯一性 且满足 --------(2) --------(3) --------(4) 上述方程组的系数行列式为n+1阶范德蒙行列式 由Cramer法则,线性方程组(4)有唯一解 定理1. p27 则满足插值条件 的插值多项式 存在且唯一. 反证：若不唯一，则除了Pn(x) 外还有另一 n 阶多项式 Ln(x) 满足 Ln(xi) = yi 。 考察 则 Qn 的阶数 ? n 而 Qn 有 个不同的根 n + 1 x0 … xn 注：若不将多项式次数限制为 n ，则插值多项式不唯一。 例如 也是一个插值多项式，其中 可以是任意多项式。 为了求得便于使用的简单的插值多项式P(x), 我们先讨论n=1的情形 要求线性插值多项式L1(x),使它满足： L1(x)的几何意义就是通过这两点的直线； 由两点式可以看出， L1(x)是由两个线性函数 也是线性插值多项式 称为线性插值基函数 n=2的情况，假定插值节点为 为了求出L2(x)的表达式，可采用基函数方法 同理 线性无关，作为二次插值基函数 得到二次插值多项式 考虑通过n+1个节点 …… n+1次多项式 且 (请同学们思考) 从而 n ? 1 希望找到li(x)，i = 0, …, n 使得 li(xj)=?ij ；然后令 ? = = n i i i n y x l x P 0 ) ( ) ( ，则显然有Pn(xi) = yi 。 li(x) 每个 li 有 n 个根 x0 … xi … xn ? = - = - - - = n j j ? i j i n i i i x x C x x x x x x C x l 0 0 ) ( ) )...( )...( ( ) ( ? - = = j ? i j i i i i x x C x l ) ( 1 1 ) ( Lagrange Polynomial 与 有关，而与 无关 节点 f 例1: 解: 且 在例1中,如果只给出两个节点169和225,也可以作插值 多项式,即1次Lagrange插值多项式,有两个插值基函数, 这种插值方法称为Lagrange线性插值,也可以在n+1个 节点中取相邻的两个节点作线性插值 ? 例2. 解: Lagrange插值基函数为 Lagrange线性插值多项式为 所以 Lagrange插值多项式的缺点： 插值基函数计算复杂 高次插值的精度不一定高 高次插值通常优于低次插值 ? 但绝对不是次数越高就越好，嘿嘿…… 三、插值余项Remainder 满足 不会完全成立 因此,插值多项式存在着截断误差,那么我们怎样估 计这个截断误差呢？ 令 设 其中 近似函数 误差 根据Rolle定理, 再由Rolle定理, 依此类推 由于 因此 所以 定理2. Lagrange型余项 *