第5章 数值积分与数值微分11.ppt

数值积分/Numerical Integration / 5.5 高斯(Gauss)求积公式 所以有 由此得 一方面，若|C2n-Cn|63? ,则有近似误差|I-C2n|?. 另一方面,(64C2n-Cn)/63应比Cn和C2n的近似程度更好. 记(64C2n-Cn)/63=Rn,称为Romberg求积公式. 用Tm(k)(m=1,2,3,4)分别表示把区间2k等分的复化梯形公式,复化Simpson公式, 复化Cotes公式和Romberg求积公式. 而且，要使|I*-Tm(k)|? ,只要|Tm(k)-Tm(k-1)|(4m-1)? (m=1,2,3,4). 则有 若对Romberg求积公式作组合也有 一般有： Romberg 序列 ? Romberg 算法： ? ? ? ? ? ? … … … … … … ? T1 = ) 0 ( 1 T ? T8 = ) 3 ( 1 T ? T4 = ) 2 ( 1 T ? T2 = ) 1 ( 1 T ? S1 = ) 0 ( 2 T ? R1 = ) 0 ( 4 T ? S2 = ) 1 ( 2 T ? C1 = ) 0 ( 3 T ? C2 = ) 1 ( 3 T ? S4 = ) 2 ( 2 T 实际计算可按下表顺序进行 T4(0) T4(1) ? T3(0) T3(1) T3(2) ? T2(0) T2(1) T2(2) T2(3) ? T1(0) T1(1) T1(2) T1(3) T1(4) ? 1 2 4 8 16 ? 0 1 2 3 4 ? Romberg公式 T4(k) Cotes公式 T3(k) Simpson公式 T2(k) 梯形公式 T1(k) 区间等分数 n=2k k 例 利用Romberg积分公式计算积分 解 按递推公式计算,结果如下 可见,由于|T1(4)-T1(3)|=0.0019531,应有|I-T1(4)|0.000651033. 3.1415858 3.1415926 3.1421177 3.1415941 3.1415926 3.1333333 3.1415687 3.1415925 3.1415926 3.0000000 3.1000000 3.1311765 3.1389885 3.1409416 1 2 4 8 16 0 1 2 3 4 T4(k) T3(k) T2(k) T1(k) n=2k k 由于|T2(3)-T2(2)|=0.0000001,应有|I-T2(3)|0.00000000666. 由于|T3(2)-T3(1)|=0.0000015,应有|I-T3(2)|0.0000000238. 由于|T4(1)-T4(0)|=0.0000068,应有|I-T4(1)|0.0000002666. ? 理查德森外推法 /* Richardson’s extrapolation */ 利用低阶公式产生高精度的结果。 设对于某一 h ? 0，有公式 T0(h) 近似计算某一未知值 I。由Taylor展开得到： T0(h) ? I = ?1 h + ?2 h2 + ?3 h3 + … ?i 与 h 无关 现将 h 对分，得： ( ) ( ) ( ) ... ) ( 3 2 3 2 2 2 2 1 2 0 + + + = - h h h h I T a a a Q：如何将公式精度由 O(h) 提高到 O(h2) ？ ... 4 3 2 1 1 2 ) ( ) ( 2 3 3 2 2 0 2 0 - - - = - - - h h I h T T h a a 即： 前面介绍的 n+1个节点的 Newton -Cotes求积公式， 其特征是节点是等距的。这种特点使得求积公式便于 构造，复化求积公式易于形成。但同时也限制了公式 的精度。 n是偶数时，代数精度为n+1， n是奇数时， 代数精度为n 。 我们知道 n+1个节点的插值型求积公式的代数精 确度不低于n 。设想：能不能在区间[a,b]上适当选择 n+1个节点 x0x1,x2,……,xn ,使插值求积公式的代数精 度高于n？ 答案是肯定的，适当选择节点，可使公式的精度 最高达到2n+1，这就是本节所要介绍的高斯求积公式。 构造具有2n+1次代数精度的求积公式 将节点 x0 … xn 以及系数 A0 … An 都作为待定系数。令 f (x) = 1, x, x2, …, x2n+1 代入可求解，得到的公式具有2n+1 次代数精度。这样的节点称为Gauss 点，公式称为Gauss 型求积公式。 例