第二章 插值a.pptVIP

第二章 插值 线性插值 n 1 抛物线插值 二次插值 n 2 * 已知: 求： sin 50? ？ 已知: 求： 当精确函数 y f x 非常复杂或未知时，在一系列节点 x0 … xn 处测得函数值 y0 f x0 , … yn f xn ，由此构造一个简单易算的近似函数 g x ? f x ，满足条件g xj f xj j 0, … n 。这里的 g x 称为f x 的插值函数。最常用的插值函数是 …? 多项式 x0 x1 x2 x3 x4 x g x ? f x 已知函数y f x 在区间[a, b]内一系列点xj上的函数值f xj yj , j 0,1,…,n，求一简单函数P x ，使满足条件 P xj yj , j 0,1,…,n， 称点xj为插值节点，[a, b]为插值区间，P x 为插值函数，f x 为被插函数， 上式称为插值条件，求P x 的方法为插值方法。 定义 §1 基本概念 §2 拉格朗日插值多项式 由于代数多项式的结构简单，实用中常取代数多项式 作为插值函数，称其为n次代数插值多项式，求Ln x 的问题叫做代数插值问题，求Ln x 的过程也叫做拉格朗日插值。 定理 唯一性 满足 的 n 阶插值多项式是唯一存在的。 证明： p.29-30 利用Vandermonde 行列式论证 反证：若不唯一，则除了Ln x 外还有另一 n 阶多项式 Pn x 满足 Pn xi yi 。 考察 则 Qn 的阶数 ? n 而 Qn 有 个不同的根 n + 1 x0 … xn 注：若不将多项式次数限制为 n ，则插值多项式不唯一。 例如 也是一个插值多项式，其中 可以是任意多项式。 已知： n 1 线性插值 3 2 4 x f x 用连接两点的直线近似表示f x ，则任意自变量所对应的函数值都可以求出。 n 1 已知 x0 , x1 ; y0 , y1 ，求 使得 1 1 1 0 0 1 , y x P y x P 可见 P1 x 是过 x0 , y0 和 x1, y1 两点的直线。 0 0 1 0 1 0 1 x x x x y y y x P - - - + 1 0 1 x x x x - - 0 1 0 x x x x - - y0 + y1 l0 x l1 x 称为拉格朗日插值基函数，满足条件 线性插值 n 1 抛物线插值 二次插值 n 2 而l0 x 、l1 x 、l2 x 必须满足以下条件 对于l0 x 必有形式 为满足l0 x0 1，则有 问题： 求二次式L2 x 满足条件 L2 x0 y0, L2 x1 y1, L2 x2 y2 用基函数表示 拉格朗日插值基函数 n ? 3 希望找到lk x ，k 0, …, n 使得 ，则显然有Pn xj yj。 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 然后令 条件：无重合节点，即 例： 3 1 0 2 f x 3 2 0 -1 x … … … … 解： 在插值区间[a, b]上用插值多项式Ln x 近似代替f x ，除了在插值节点上没有误差外，其他点上一般是存在误差的。若记Rn x f x -Ln x ,则Rn x 就是用Ln x 近似代替f x 时的截断误差。我们称Rn x 为插值多项式的余项，并且可根据下面定理来估计它的大小。 拉格朗日插值余项定理 注：? 通常不能确定 ?x , 而是估计 , ?x? a,b 将 作为误差估计上限。 ?当 f x 为任一个次数? n 的多项式时， , 可知 ，即插值多项式对于次数? n 的多项式是精确的。 定理 设区间[a，b]含有节点x0，x1,…,xn，而f x 在[a, b]内有直到n+1阶导数，则n次插值多项式的误差可表示为