插值法与数值微分.pptVIP

第46页，共78页，2022年，5月20日，2点39分，星期五 二、向后差分与牛顿向前后插值公式 第47页，共78页，2022年，5月20日，2点39分，星期五 【例10】已知函数表同例9，计算sin(0.58)，并估计截断误差. 因三阶向后差分接近于常数，故用三次插值进行计算,且 于是由后插公式得 第48页，共78页，2022年，5月20日，2点39分，星期五 第49页，共78页，2022年，5月20日，2点39分，星期五 定义 1 记 称 为 关于xi 的 零阶均差. 称 为 关于 xi , xi+1 的一阶均差. 称为二阶均差. 三、差商与牛顿基本插值多项式 第50页，共78页，2022年，5月20日，2点39分，星期五 一般地, k 阶均差为 均差有如下基本性质： 定理 1: (1) 均差与函数值的关系为 (2) 均差与节点的排列顺序无关, 即 第51页，共78页，2022年，5月20日，2点39分，星期五 (4) 若函数 在 上存在n 阶导数, 且节点 则 使得 第52页，共78页，2022年，5月20日，2点39分，星期五 * 三、均差的计算方法(表格法): 规定函数值为零阶均差 均差表 第53页，共78页，2022年，5月20日，2点39分，星期五 解：先构造差商表如表2-5所示。 由表可以看出牛顿基本插值多项式中各系数为 表2.5 第54页，共78页，2022年，5月20日，2点39分，星期五 故用线性插值所得的近似值为 用抛物插值所得的近似值为 第55页，共78页，2022年，5月20日，2点39分，星期五 一、拉格朗日插值公式： 问题提出： 这节就具有一般形式的代数插值问题（即已知函数 在n+1个点上的函数值 求一个n次多项式 ，并满足条件 ， ）来讨论如何构造其插值多项式 。 §2-2 拉格朗日插值多项式 第14页，共78页，2022年，5月20日，2点39分，星期五 第15页，共78页，2022年，5月20日，2点39分，星期五 第16页，共78页，2022年，5月20日，2点39分，星期五 这就是所要求的插值多项式，称为拉格朗日（Lagrange） 插值多项式。 当n=1时，就得出线形插值多项式， n=2时，就得出抛物插值多项式。 第17页，共78页，2022年，5月20日，2点39分，星期五 二、拉格朗日插值余项： 插值余项： 定理： 第18页，共78页，2022年，5月20日，2点39分，星期五 证明： 当 X为节点时，两边皆为0，显然成立。 下设 X 不为节点。作辅助函数 第19页，共78页，2022年，5月20日，2点39分，星期五 即问题得证。 这个定理所讲的余项用起来有一定的困难 ，因为实际计算时，只是给出 的一张数据表，并未给出具体的解析式子，故 并不知道，所以 也就无法得到。 第20页，共78页，2022年，5月20日，2点39分，星期五 第21页，共78页，2022年，5月20日，2点39分，星期五 【例4】在例3中分别用线性插值 和抛物插值计算了 的近似值，试估计它们的截断误差。 第22页，共78页，2022年，5月20日，2点39分，星期五 第23页，共78页，2022年，5月20日，2点39分，星期五 第24页，共78页，2022年，5月20日，2点39分，星期五 解：记 由插值多项式有 故 根据余项公式，若能估计出 的上界 ， 那么将有 第25页，共78页，2022年，5月20日，2点39分，星期五 三、插值误差的事后估计法 第26页，共78页，2022年，5月20日，2点39分，星期五 利用余项公式知： 第27页，共78页，2022年，5月20日，2点39分，星期五 稍加整理得： 这种用计算的结果来估计误差的办法，通常称为事后估计，在计算中是常用的，这种估计误差的方法，将贯穿我们计算方法这门课程的始终。 第28页，共78页，2022年，5月20日，2点39分，星期五 四、拉格朗日插值多项式的优缺点： 优点：拉格朗日插值多项式结构对称，使用方便 缺点： a.不具备递推性，当需要增加节点时需要重新计算； b.龙格（Runge）现象:高次拉格朗日插值多项式稳定性差，对于计算过程的舍入误差十分敏感，当插值节点增多时，不能保证非节点处的插值精度得到改善，有时反而误差更大。 龙格就给出了一个例子：设被插值函数 第