[数学]Lagrange插值法1.ppt

练习: 已测得某地大气压强随高度变化的一组数据 用抛物插值计算 sin0.3367时，因为 这个结果与六位有效数字的正弦函数表完全一样，这说明查表时用二次插值精度已相当高了。插值多项式的存在性?其截断误差是多少? 所以 定理1 假设x0 ,x1,…,xn 是n+1个互异节点,函数f(x)在这组节点的值f(xj)(j=0,1,…,n)是给定的，那么存在唯一的次数不超过n 的多项式Pn (x)满足 Pn (xj)=f(xj), j=0,1,…,n 二 、 插值多项式的存在性和唯一性 证明：由插值函数的构造知n次多项式Pn (x)的存在性，由代数基本定理可证明它的唯一性。 函数f(x)的n次插值多项式Ln(x)只是在节点处有: Ln (xj)=f (xj), j=0,1,…,n 若x≠xj,一般Ln (x) ≠ f (x), 令 Rn(x)=f (x)-Ln(x) 称Rn(x)为插值多项式的余项 三 、 Lagrange插值的插值余项 (截断误差) 定理2：设Ln(x)是过点x0 ，x1 ，x2 ，…xn的 n 次插值多项式， f (x)在[a，b]上有n阶连续导数，在(a，b)内存在n+1阶导数，其中[a，b]是包含点x0 , x1 , x2 ,…, xn的任一区间，则对任意给定的x?[a，b]，总存在一点??(a，b)（依赖于x）使 其中 。 证明: Rn(x) = f(x) - Ln(x), 当x=xi时，显然有： Rn(xi) =f(xi) - Ln(xi), =0 ， ?n+1 (xi)=0 （ i=0,1,…,n）结论成立 当x≠xi时，∵Rn(xi ) =0 ，?n+1 (xi)=0 （ i=0,1,…,n） 　∴可设Rn(x)=K(x) ?n+1 (x) 需证:K(x)=f(n+1)( ?)/(n+1)! 现在[a,b]上任意固定一点x,引进辅助函数 g(t)=f(t)- Ln(t)-K(x)? n+1 (t), (*) 则g(t)在[a,b]上具有n阶连续导数，在(a,b)内存在n+1阶导数，当 t= x, x0, x1,…, xn 时，g(t)=0,，即g(t)在(a,b)内有n+2个零点,由Rolle定理知g’(t)在(a,b)内有n+1个零点,如此反复，最后可推知g(n+1)(t)在(a,b)内有1个零点, 即有g(n+1)( ?)=0, a ?b。这样，由(*)式便有 现在[a,b]上任意固定一点x,引进辅助函数 g(t)=f(t)- Ln(t)-K(x)? n+1 (t), (*) 有g(n+1)( ?)=0, a ?b。这样，由(*)式便有 因为? n+1 (t)是n+1次多项式， ? n+1 (n+1)(t)=(n+1)!,又因为Ln(t)是次数为n的多项式,因此Ln (n+1)(t) = 0 。 由此得： K(x)=f(n+1)( ?)/(n+1)! . 代入Rn(x)=K(x) ? n+1(x)，定理得证。 需证:K(x)=f(n+1)( ?)/(n+1)! 上式称为带余项的Lagrange插值公式。 当f(x)具有n+1阶导数，才有余项表达式成立，且余项为： 由此得： 那么插值多项式Ln(x)逼近f(x)的截断误差是 注：1）余项表达式只有在 f(x) 的n+1阶导数存在时才能应用。 2）? 在 （a，b）内的具体位置通常不可能给出，如果我们可以求出 例2: 已给sin0.32=0.314567,sin0.34=0.333487, sin0.36=0.352274, 用线性插值及抛物插值计算 sin0.3367的值并估计截断误差。 y2=0.352274 y1=0.333487 y0=0.314567 x2=0.36 x1=0.34 x0=0.32 解： 用线性插值计算，取 x0=0.32 及 x1=0.34 , 又由公式得 y1 - y0 sin0.3367?L1(0.3367)=y0+————(0.3367 -x0) x1 - x0 0.01892 =0.314567+ ——— (0.0167) =0.330365 . 0.02 其截断误差得 其中 其截断误差得 ?R1(0.3367)?=?sin 0.3367 –L1(0.3367)?