Chap03数值积分演示课件文.pptVIP

例如: 对于积分 采用变步长的梯形公式进行计算. 将区间 [ 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。 例2.3 用变步长的梯形公式计算积分 解: 对于 , 定义 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. . . ), 可见收敛速度比较缓慢。 为了加速变步长梯形公式的收敛速度，我们采用外推策略。 3. Richardson外推算法 若用一个步长为 h 的函数 I1( h ) 去逼近问题 I , 设其 截断误差可表示为 为了提高逼近的精度，选取 q 为满足 的 正数, 将上式(1)中的 h 换为 qh , 则有 其中 是与h无关的常数,并且 ， (2.1) 由 (1) 可知I1(h)逼近I的误差为 。 (2.2) (2.2) 式减上式 , 得 式(1)两端同乘以 得 (2.1) (2.2) 则 I2(h) 逼近I 误差降为 令 其中 是与h无关的常数,则有 ， 如此继续。 一般地, 选取 q 为满足 的正数, 由此得到序列 则 Im+1 ( h ) 逼近 I 的误差由下面的定理给出。 定理 2.1 设 I1 ( h ) 逼近 I 的截断误差由下式给出 则 Im+1 ( h )逼近 I 的截断误差为 其中 是与 h 无关的常数。 这种利用序列{Im+1(h)} 逐步加速去逼近 I 的方法 称为Richardson外推算法 Richardson外推公式 4. Romberg 算法 Romberg 算法是利用变步长的梯形求积序列 外推加速来逼近积分真值的算法. 考虑积分 由复化梯形公式有 现在将 Tn 记为 T1( h ), 则T2n 为 T1(h/2) 且 设 f (x) 在区间 [a , b] 上任意次可微, 根据 Euler-Maclaurin公式有 其中 是与 h 无关的常数。 因为Pm = 2m , 带入上式整理后得 易知 Tm+1( h ) 逼近 I 的误差为 O ( h2(m+1) ) ，这种算法称为 Romberg算法。 ， 。 则有 选取 利用 Richardson 外推公式, 知T2 ( h ) = Sn， 当m =1