计算方法第五讲精要.pptVIP

* §5.3 复化求积公式 (1) 如果使用低阶的Newton-Cotes公式，误差将会较大； 所谓复化求积法，就是 (3) 最后将每个小区间上的积分的近似值相加。 为了提高公式的精度,又使算法简单易行,往往使用复化求积法. 一、问题 (1) 首先把整个积分区间分成若干个子区间（通常采用等分) ； (2) 然后在每个小区间上使用低阶Newton-Cotes公式; 如果使用高阶的Newton-Cotes公式，则系数复杂， 且公式的稳定性和收敛性也不能保证。 二、复化求积公式 各节点为 复化梯形公式 复化辛普森公式 复化柯特斯公式 三、复化求积公式的截断误差 我们上一节知道，三个求积公式的截断误差分别为 单纯的求积公式 复化求积公式在每个小区间 复化梯形公式的截断误差 由于 即有 复化simpson公式的截断误差 复化cotes公式的截断误差 比较三种复化公式的的余项 例1 解: =3 例2 解: 各节点的值如下面表格 各节点的值如下表 0 1 0.125 0.9973978 0.25 0.9896158 0.375 0.9767267 0.5 0.9588510 0.625 0.9361556 0.75 0.9088516 0.875 0.8771925 1 0.8414709 分别由复化梯形、复化辛普森。复化柯特斯公式有 原积分的精确值为 精度最高 精度次高 精度最低 三个公式的结果比较 三种方法的工作量基本相同。收敛速度一个比一个快。实际应用时，较多地应用辛普森公式。 四、自动选取积分步长 (1) 高阶导数的估计往往是很困难的； “事后估计法”的基本思想是 (2) 利用前后两次的计算结果来判断误差是否满足精度要求， 从而确定n. 为了改正上述缺点，实际常采用“事后估计法” 事前确定步长的问题 (1) 求数值积分时，将区间逐次分半； 下面以复化梯形公式为例来介绍这种步长逐次减半求积法 (2) 这种估计往往是很保守的，得到的n往往偏大。 将区间[a,b] n等分，步长 ，复化梯形公式为 再将步长减半，增加新分点 ， 公式变为 如何根据Tn和T2n来确定误差是否满足要求？ 则有 即 如果二阶导数在区间[a,b]上变化不大 类似地可以推得 例3 解: 0.0000003 0.9460830 512 0.0002941 0.9459850 16 0.0000012 0.9460827 256 0.0011774 0.9456909 8 0.0000046 0.9460815 128 0.0047202 0.9445135 4 0.0000183 0.9460769 64 0.0190578 0.9397933 2 0.0000736 0.9460586 32 0.9207355 1 |T2n-Tn| T2n 2n |T2n-Tn| T2n 2n §4.3 龙贝格 （Romberg）算法 一、问题 复化梯形公式简单，但计算精度较差 有没有办法改善复化梯形公式的精度呢? 由上节关于复化梯形公式的误差公式 可得 此式可视作用误差对 作一种补偿，可以提高精度。