龙贝格求积公式.ppt

用龙贝格方法计算积分的步骤为： （1）：准备初值，先用梯形公式计算积分近 似值： （2）：按变步长梯形公式计算积分近似值： 令 计算： （3）：按加速公式求积分（为便于编程，写 为下列形式） 梯形加速： 辛普生加速： 柯特斯加速： （4）：精度控制 当 时终止计算 并取 为近似值，否则，将步长折半，转 （2）执行。实际计算时的加工流程图7.4.1示： 例7.4.2 用龙贝格算法加工例7.4.1得到的近似值 解：用算法图7.4.2计算结果见表7.4.1,表中 k表示二分次数 ： * * * （1）使用复化梯形公式、Simpson公式，首先要确定步长 ； （2）而步长要根据余项确定，这就涉及到高阶导数的估计； （3）高阶导数的估计一般比较困难，且估计值往往偏大； （4）计算机上实现起来不方便，通常采用“事后估计法”。 三、积分步长的自动选取： 注意事项： ?基本思想： 将积分区间逐次分半 ?终止法则： 前后两次近似值的误差小于已知精度 ?具体过程（以复化梯形公式为例） 1、首先将区间 n等分： 2、再将区间 2n等分，即步长减半： 上述条件满足，程序终止；否则，继续分半计算。 3、终止条件： 由复化梯形公式的余项知 变化不大时 由此得到近似关系式 误差控制条件 ?对于复化Simpson公式、Cotes公式可以类似得到 ?对于复化梯形公式 加速 收敛 应用步长逐次减半得到的复化梯形值、复化 Simpson值、复化Cotes值与精确值的比较 Romberg积分法/*Romberg Integration Method */ ? Romberg积分思想 由上节分析知，用复化梯形公式计算积分值 的误差大约为： 令 由复化梯形公式知 梯形加速公式： 利用复化梯形公式前后两次积分近似值 和 ，按 照上式作出的线性组合得到了具有更高精度的积分值。 上述公式说明： Romberg积分公式正是由此思想产生 Romberg 值序列 Simpson加速公式： Cotes加速公式： 类似于梯形加速公式的处理方法，得到： 通过上述3个积分值序列求积分近似值的方法， 称之为Romberg积分法。 4个积分值序列： 梯形值序列 Simpson值序列 Romberg值序列 Cotes值序列 Romberg积分法的一般公式 其中 Romberg积分表 0.9460831 0.9460830 0.9460831 0.9461459 0.9460869 0.9460834 0.9207355 0.9397933 0.9445135 0.9456909 0 1 2 3