第二节 复化求积式和龙贝格求积公式.ppt

* 第三节 复化求积公式 将积分区间 分成若干个小区间，然后在每个小区间上采用低阶的牛顿-柯特斯公式。然后将所有小区间的计算结果加起来。 一、复化梯形公式： 将积分区间 n等分： 分点 在区间 上采用梯形公式 复化梯形公式 复化梯形公式的几何意义 小梯形面积之和近似 复化梯形公式的余项 设 ，则余项估计式为： 的误差是二阶的。 由上式可知，误差与 同阶，此时称复化梯形公式 误差的阶数越高，精度越好！ 二、复化辛蒲生公式： 分点 在区间 上采用辛蒲生公式 其中 将积分区间 n等分： 复化辛蒲生公式 复化辛蒲生公式的几何意义 小抛物面积之和近似 复化辛蒲生公式的余项 设 ，则有余项估计式 复化辛蒲生公式中“半点”的处理 可将整个区间等分成偶数个小区间，每两个小区间 合并起来视为复化辛蒲生公式中的一个小区间。 类似地，可以得到复化柯特斯公式 它的余项为 例2：将[0,1]区间八等分，根据如下函数值表，利用复 化梯形公式、复化辛蒲生公式计算积分 的近似值。 0.976727 0.989688 0.997398 1 3/8 1/4 1/8 0 0.841471 0.877193 0.908858 0.936156 0.958851 1 7/8 6/8 5/8 1/2 解： 分别采用复化辛蒲生公式、复化梯形公式 复化梯形公式 复化辛蒲生公式 复化梯形公式(n = 8)， 复化辛蒲生公式(n = 4)， （1）使用复化梯形公式、辛蒲生公式，首先要确定步长 ； （2）而步长要根据余项确定，这就涉及到高阶导数的估计； （3）高阶导数的估计一般比较困难，且估计值往往偏大； （4）计算机上实现起来不方便，通常采用 “事后估计法”。 三、积分步长的自动选取： 注意事项： ?基本思想： 将积分区间逐次分半，比较前后两次的近似值 ?终止法则： 前后两次近似值的误差小于已知精度 ?具体过程（以复化梯形公式为例） 1、首先将区间 n等分： 2、再将区间 2n等分，即步长减半： 上述条件满足，程序终止；否则，继续分半计算。 3、终止条件： 由复化梯形公式的余项知 变化不大时 由此得到近似关系式 误差控制条件 例3：根据如下函数值表，利用复化梯形公式计算积分 的近似值，要求误差不超过 。 0.976727 0.989688 0.997398 1 3/8 1/4 1/8 0 0.841471 0.877193 0.908858 0.936156 0.958851 1 7/8 6/8 5/8 1/2 解： 先在整个区间上用梯形公式 然后将区间二等分，利用递推公式求出 递推公式 进一步二分积分区间，类似可求出 如此不断二分并利用递推公式，可得下表中的结果 k 表示二分次数，区间数 1 2 3 4 5 6 7 8 9 0.9397933 0.9445135 0.9456909 0.9459850 0.9460596 0.9460769 0.9460815 0.9460827 0.9460830 由表中可以看出，对分8次和对分7次之间的差 因而 是满足精度要求的解。 误差控制条件 收敛速度慢 ?对于复化辛蒲生公式、柯特斯公式可以类似得到 不足 加速 收敛 应用步长逐次减半得到的复化梯形值、复化 辛蒲生值、复化柯特斯值与精确值的比较 第四节 龙贝格求积公式 ? 龙贝格积分思想 由上节分析知，用复化梯形公式计算积分值 的误差大约为： 令 由复化梯形公式知