第6章数值积分与数值微分-天津大学理学院数学系.ppt

(2) 将区间 [ 0, 1] 4等分, 步长 采用复化 Simpson 公式计算, 仍然利用原来 9 个分点处的函数值, 求得 这两种方法计算量基本相同, 但所得到的结果与真值π= 3.1415926 ... 比较可以看出复化 Simpson 公式求得的结果要精确得多. 6.4 变步长的梯形公式与Romberg算法 1.变步长的梯形公式 复化求积公式称为定步长的求积公式，它对提高精度是行之有效的。但对于给定的精度，要确定一个合适的步长往往难以办到。因此实际上一般常采用变步长的求积公式。即让步长逐次折半的过程中，反复使用复化求积公式进行计算，直到相邻两次计算结果之差的绝对值小于允许精度ε的要求时终止计算，这种方法称为变步长的求积方法。 例如: 对于积分 采用变步长的梯形公式进行计算. 将区间 [ a , b] n 等分 , 步长 , 按复化梯形公式 计算时, 需调用 n +1 个函数值。 现在将 h 折半, 再将上述每个区间 [ xk , xk+1 ] 对分一次, 分点增至 2n + 1 个, 设上述小子区间的中点为 在[ xk , xk+1 ] 上用复化梯形公式并求和得 上式称为变步长的梯形公式. 即在求 T2n 时, 可以利用前面已求出的结果Tn , 剩下的仅仅需要求出 n 个新分点处的函数值. 注意：h = xk+1-xk 变步长的梯形公式的算法 Step 1. 给定精度 ? 0，m为正整数，步长h =(b - a)/2m。 即将积分区间分割成2m等份。 Step 2. 计算 这里 Step 3. 计算 这里 将每一个小子区间二等分，即步长折半。 Step 4. 如果 ，则停止，输出值 , 否则，置 m = m+1，h : = h/2 ，转到Step 3。 例6.4.1 用变步长的梯形公式计算积分 解: 对于 , 定义 f (0) = 1, 首先在区间 [0 , 1] 上用梯形公式（即步长 h = 1）,求得 将 [0 , 1] 对分, 它的中点函数值 , 则有 如果 不成立，则 h := h/2 = 1/2 ，计算 （精确到10-6） 如此继续下去, 计算结果如下表 如果 不成立，则 h := h/2 = 1/4 ， 继续计算 。 k k T 2k T 2k 0 1 2 3 4 5 6 7 8 9 10 0.920 735 5 0.939 793 3 0.944 513 5 0.945 690 9 0.945 985 0 0.946 059 6 0.946 076 9 0.946 081 5 0.946 082 7 0.946 083 0 0.946 083 1 从上表可看出, 将积分区间对分了10次, 求得 I 的近似值为0.9460831 (积分精确值为0.9460831. . . ), 可见收敛速度比较缓慢。 为了加速变步长梯形公式的收敛速度，我们采用外推策略。 2. Richardson外推算法 若用一个步长为 h 的函数 I1( h ) 去逼近问题 I , 设其 截断误差可表示为 为了提高逼近的精度，选取 q 为满足 的 正数, 将上式(1)中的 h 换为 qh , 则有 其中 是与h无关的常数,并且 ， (1) 由 (1) 可知I1(h)逼近I的误差为 。 (2) (2) 式减上式 , 得 式(1)两端同乘以 得 (1) (2) 则 I2(h) 逼近I 误差降为 令 其中 是与h无关的常数,则有 ， 如此继续。 一般地, 选取 q 为满足 的正数, 由此得到序列 则 Im+1 ( h ) 逼近 I 的误差由下面的定理给出。 定理 6.4.1 设 I1 ( h ) 逼近 I 的截断误差由下式给出 则 Im+1 ( h )逼近 I 的截断误差为 其中 是与 h 无关的常数。 这种利用序列{Im+1(h)} 逐步加速去逼近 I 的方法 称为Richardson外推算法 Richardson外推公式 3. Romberg 算法 Romberg 算法是利用变步长的梯形求积序列 外推加速来逼近积分真值的算法. 考虑积分 由复化梯