计算方法复合求积公式Romberg算法ch04br讲述.ppt

第四章数值积分与数值微分 本讲内容 复合求积公式 复合梯形公式 复合 Simpson 公式 举例 举例 举例 举例 Romberg 算法 梯形法递推公式 梯形法递推公式 梯形法递推公式 举例 梯形法的加速 梯形法的加速 举例 Romberg 算法 举例 作业 * * 计算方法 —— 复合求积公式 Romberg 算法 复合求积公式 复合梯形公式 复合 Simpson 公式 梯形法的递推化计算 Romberg 算法基本思想: 外推技巧 Romberg 算法: 计算过程 Romberg (龙贝格) 算法 提高积分计算精度的常用两种方法 用 复合公式 用 非等距节点 将积分区间分割成多个小区间 在每个小区间上使用低次牛顿－科特斯求积公式 复合求积公式 将 [a, b] 分成 n 等分 [xi , xi+1] ，其中 (i = 0, 1, …, n) 复合梯形公式 余项 复合 Simpson 公式 余项 性质：复合梯形公式和复合Simpson 公式都是收敛的，也都是稳定的。 例：设 ，利用下表中的数据分别用复合梯形公式和复合simpson公式计算定积分 ，并估计误差。 0.936 5/8 0.909 6/8 0.990 2/8 0.977 3/8 0.841 0.877 0.954 0.997 1 f (xi ) 1.0 7/8 1/8 0 4/8 xi 解： 误差估计 例：计算定积分 用复合梯形公式和复合simpson公式时，n 分别取多大时才能使得误差不超过 0.5 ? 10-5 解： 要使误差不超过 0.5 ? 10-5 ，需要 取 n=213 213 等分 复合梯形公式 复合simpson公式 要使误差不超过 0.5 ? 10-5 ，需要 故取 n=4 8 等分 太 大 利用复合梯形公式、复合simpson公式、复合Cotes公式等计算定积分时，如何选取步长 h ？ 计算精度难以保证 太 小 增加额外的计算量 解决办法：采用 变步长算法 通常采取将区间不断对分的方法，即取 n = 2k ，反复使用复合求积公式，直到所得到的计算结果满足指定的精度为止。 步长折半：[xi , xi+1/2] ， [xi +1/2 , xi+1] 将 [a, b] 分成 n 等分 [xi , xi+1] ， xi xi +1 xi +1/2 记 解： 例：用梯形法的递推公式计算定积分 , 要求计算精度满足 ex42.m 0.946082975 9 0.946082687 8 0.946081539 7 0.946076943 6 0.946058561 5 0.945985030 4 0.945690864 3 0.944513522 2 0.939793285 1 0.920735492 0 0.946083046 10 T (k) k 梯形法递推公式算法简单，编程方便 梯形法的加速－－龙贝格 (Romberg) 算法 但收敛速度较 慢 定理：设 f(x)?C ?[a, b], 记 Tn = T (h), 则有 Richardson 外推算法 例：计算定积分 ex43.m 0.946082975 9 0.946082687 8 0.946081539 7 0.946076943 6 0.946058561 5 0.945985030 4 0.945690864 3 0.944513522 2 0.939793285 1 0.920735492 0 0.946083046 10 T0(k) k ① T1 =T0 (0) ② T2 =T0 (1) ③ S1 =T1 (0) ④ T4 =T0 (2) ⑤ S2 =T1 (1) ⑥ C1 =T2 (0) ⑦ T8 =T0 (3) ⑧ S4 =T1 (2) ⑨ C2 =T2 (1) ⑩ R1 =T3 (0) 记: : k 次等分后梯形公式计算所得的近似值 : m 次加速后所得的近似值 Romberg 算法是收敛的 例：用 Romberg 算法计算定积分 , 要求计算精度满足 ex44.m 000000003 0004 000000.4004