3.2数值积分分析.ppt

* 补充：使用用龙贝格方法计算积分 使得计算误差 * * 参考书： 数值分析及其MATLAB实现，任玉杰，高等教育出版社，2007年 数值分析（第5版）习题解答，张威，杨月婷，清华大学出版社，2010年 数值分析与计算方法，蒋勇,李建良，科学出版社，2012年 数值分析（上下册），黄明游,冯果忱，高等教育出版社，2008年 数值分析与工程应用，杨威启，李晓玲，师忠秀，国防工业出版社，2008年 梯形公式/*Trapezoidal Formula */ 1次代数精度 用梯形面积近似 - * 3次代数精度 Simpson公式 用抛物形面积近似 - * Cotes公式 5次代数精度 近似等于曲边梯形的面积 - * §3 复化求积公式/*Compound Quadrature Formula */ 思想 将积分区间 分成若干个小区间，然后在每个小区间上采用低阶的Newton—Cotes公式 一、复化梯形公式：/*Compound Trapezoidal Formula */ 将积分区间 n等分： 分点 在区间 上采用梯形公式 * 复化梯形公式 复化梯形公式的几何意义 小梯形面积之和近似 - * 复化梯形公式的余项 设 由介值定理 余项估计式 * 复化梯形公式的收敛性 其中定积分与 区间分法和 的取法无关 设 * 二、复化Simpson公式： /*Compound Simpon Formula */ 分点 在区间 上采用Simpson公式 其中 将积分区间 n等分： * 复化Simpson公式 复化Simpson公式的几何意义 小抛物面积之和近似 - * 复化Simpson公式的余项 设 由介值定理 余项估计式 * 复化Simpson公式的收敛性 类似地可以得到复化Cotes公式 * 例2：分别利用复化梯形公式、复化Simpson公式计算 积分 的近似值，要求按复化Simpson公 式计算时误差不超过 。 解： 首先来确定步长 复化Simpson公式的余项： 其中 * 本题 的求法： 由归纳法知 * 解不等式得 取区间 9个等分点，分别采用复化Simpson、梯形公式 0 1/8 1/4 3/8 1 0.997398 0.989688 0.976727 1/2 5/8 6/8 7/8 1 0.958851 0.936156 0.908858 0.877193 0.841471 * 复化梯形公式(n=8) 复化Simpson公式(n=4) * 复化梯形公式 复化Simpson公式 * （1）使用复化梯形公式、Simpson公式，首先要确定步长 ； （2）而步长要根据余项确定，这就涉及到高阶导数的估计； （3）高阶导数的估计一般比较困难，且估计值往往偏大； （4）计算机上实现起来不方便，通常采用“事后估计法”。 三、积分步长的自动选取： 注意事项： ?基本思想： 将积分区间逐次分半 ?终止法则： 前后两次近似值的误差小于已知精度 * ?具体过程（以复化梯形公式为例） 1、首先将区间 n等分： 2、再将区间 2n等分，即步长减半： * 上述条件满足，程序终止；否则，继续分半计算。 3、终止条件： 由复化梯形公式的余项知 变化不大时 由此得到近似关系式 误差控制条件 * ?对于复化Simpson公式、Cotes公式可以类似得到 * 加速 收敛 应用步长逐次减半得到的复化梯形值、复化 Simpson值、复化Cotes值与精确值的比较 * §4 Romberg积分法/*Romberg Integration Method */ ? Romberg积分思想 由上节分析知，用复化梯形公式计算积分值 的误差大约为： 令 由复化梯形公式知 * * 梯形加速公式： 利用复化梯形公式前后两次积分近似值 和 ，按 上式作出的线性组合得到了更高精度的Simpson公式。 上述公式说明： Romberg积分公式正是由此产生 * Romberg 值序列 Simpson加速公式： Cotes加速公式： 类似于梯形加速公式的处理方法，得到： * 通过上述3个积分值序列求积分近似值的方法， 称之为Romberg积分法。 4个积分值序