第5章-数值微分与数值积分.ppt

第五章 数值微分与数值积分 5.1 数值微分 5.2 数值积分 复合求积公式的误差（余项） 解：使用复化梯形公式 使用复化Simpson公式： 例 5. 3 分别用复化梯形公式、复化Simpson公式计算定积分 * 微积分中，关于导数的定义如下： 自然，而又简单的方法就是， 取极限的近似值，即差商. 5.1 数值微分 5.1.1 差商型求导公式 由Taylor展开 因此，有误差 向前差商 5.1 数值微分 5.1.1 差商型求导公式 由Taylor展开 因此，有误差 向后差商 5.1 数值微分 5.1.1 差商型求导公式 由Taylor展开 因此，有误差 中心差商 5.1 数值微分 5.1.1 差商型求导公式 由误差表达式，h越小，误差越小，但同时舍入误差增大， 所以，有个最佳步长 我们可以用事后误差估计的方法来确定 设D(h),D(h/2)分别为步长为h,h/2的差商公式。则 时的步长h/2就是合适的步长 f(x)=exp(x) -0.0032 3.1550 0.01 -0.0018 3.1600 0.06 -0.0007 3.1575 0.02 -0.0025 3.1607 0.07 -0.0001 3.1583 0.03 -0.0031 3.1613 0.08 -0.0006 3.1588 0.04 -0.0040 3.1622 0.09 -0.0008 3.1590 0.05 -0.0048 3.1630 0.10 R(x) f’(1.15) h R(x) f’(1.15) h 例： 插值是建立逼近函数的手段，用以研究原函数的性质。因此，可以用插值函数的导数近似为原函数的导数。 误差 5.1 数值微分 5.1.2 插值型求导公式 注意：为了便于估计误差，限定只能对节点上的导数值采用插值 多项式的相应导数进行近似。 1、两点公式 5.1 数值微分 5.1.2 插值型求导公式 给定两点上的函数值 这称为两点公式。 截断误差: 若给定三点上的函数值 则由 这称为三点公式，其中（5.1.9）又称为中点公式。 2、三点公式 例1：已知列表 X 2.5 2.55 2.60 2.65 2.70 Y 1.58114 1.59687 1.61245 1.62788 1.64317 解: h=0.05 二阶导数公式及误差 对其求二阶导数得 由Taylar展开可得误差估计式 5.2 数值积分 5.2.1 插值型求积公式 （1）插值型求积公式 （2）Newton-Cotes型求积公式 （3）梯形公式、Simpson公式和Cotes公式 5.2.2 复化求积公式 （1）复化梯形公式 （2）复化Simpson公式 5.2.3 Romberg积分法 （1）梯形逐步减半算法 （2）Romberg积分法 问题：如何求积分 数学分析：牛顿-莱布尼茨(Newton-Leibniz)公式： N-L公式失效的情形： （1）被积函数，诸如 等等，找不到用 初等函数表示的原函数； （2）当 是由测量或数值计算给出的一张数据表. 这时，牛顿-莱布尼茨公式也不能直接运用; 问题：点ξ的具体位置一般是不知道的，因而难以 准确算出 的值，怎么办？ 只要对平均高度 提供一种算法，相应地便可获得 一种数值求积方法. 由积分中值定理知，在积分区间 内存在一点ξ， 成立 构造数值积分公式的基本思想： （1）左矩形公式 （3）用区间中点 的“高度” 近似地取代平均 高度 ，则又可导出所谓中矩形公式 （2）右矩形公式 （4）用两端点“高度” 与 的算术平均作为平均高度 的近似值，这样导出的求积公式 是梯形公式. 一般地，可以在区间 上适当选取某些节点 ， 然后用 加权平均得到平均高度 的近似值，这样 权 仅仅与节点 的选取有关， 构造出的求积公式具有下列形式： 的具体形式. 而不依赖于被积函数 式中 称为求积节点； 称为求积系数,亦称伴随节点 的权. k A 将这种思想一般化： （1）插值型求积公式 f(xn) xn f(xi+1) xi+1 f(xi) xi f(xi-1) xi-1 --- --- f(x1) f(x0) --- --- x1 x