Chapter6_2_RombergGauss积分例析.ppt

? 理查德森外推法 /* 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 即： Romberg求积算法 以复合梯形公式算法为例介绍： 将[a,b]n等分，h为步长，复合梯形公式为 若将[a,b]2n等分，即将求积区间[xk,xk+1]再二分一次，只增加一个分点xk+1/2=(xk+xk+1)/2，用复合梯形公式求得该区间的积分值为： 分析误差： 假定 则 则有 即 若T2n与Tn接近，则T2n误差很小。 这种以计算结果估计误差的方法称事后估计法。 若用T2n的误差作为T2n的一种补偿，得 可能是更好的结果。 由 复合Simpson公式 即T2n与Tn作线性组合，可得Simpson公式的值Sn。 考察Simpson方法，类似推导可得 即 复合Cotes公式 重复上述过程，可得Romberg公式： 上述讨论说明，由梯形公式出发，将[a,b]逐次二分可提高精度。 设 若记Tn=T(h)，当[a,b]2n等分时，T2n=T(h/2)，且有 将T(h)展开成h2的幂级数形式： （其中αk与h无关） 则 记 则T1(h)与I的近似阶为O(h4)，且序列T1(h),T1(h/2),……即Simpson序列Sn,S2n,……。 又 若令 则又可进一步消去h4项。 记为 序列T2(h),T2(h/2),……，即Cotes序列Cn,C2n,……，近似阶为O(h6)。 继续下去，每加速一次，误差量级提高2阶。 一般地，若记T0(h)=T(h)，则有 经过m次加速后， 上述处理方法称理查森（Richardson）外推加速方法。 设以T0(k)表示二分k次后求得的梯形值，且以Tm(k)表示序列{T0(k)}的m次加速值，则以外推公式 (m=1,2,…,k) 称Romberg求积算法。 得 T表： （1）取k=0，h=b-a，求T0(0)=h/2[f(a)+f(b)]； 由1→k（k为[a,b]二分次数）计算： （2）求梯形值T0(k)， （3）求加速值。用Romberg求积公式逐个求出T表的第k行其余元素Tj(k-j)，(j=1,2,…,k)； （4）若| Tk(0)－Tk-1(0) |ε，则中止计算，并取Tk(0)≈I；否则，令k+1→k，转(2)继续。 计算过程： 例1：用Romberg求积算法计算定积分 解： k T0(k) T1(k) T2(k) T3(k) T4(k) T5(k) 0 1 2 3 4 5 0.9207 0.9398 0.9445 0.9457 0.9460 0.9461 0.9461 0.9461 0.9461 0.9461 0.9461 0.9461 0.9461 0.9461 0.9461 0.9461 0.9461 0.9461 0.9461 0.9461 0.9462 例2：用Romberg求积算法计算定积分 解： k T0(k) T1(k) T2(k) T3(k) T4(k) T5(k) 0 1 2 3 4 5 0.5 0.426777 0.407018 0.401812 0.400463 0.400118 0.400002 0.400002 0.400002 0.400002 0.400002 0.400014 0.400009 0.400009 0.400009 0.400077 0.400054 0.400050 0.400432 0.400302 0.402369 (I=0.4) (取| Tk(0)－Tk-1(0) |10-5) § 高斯型积分 /* Gaussian Quadrature */ 构造具有2n+1次代数精度的求积公式 将节点 x0 … xn 以及系数 A0 … An 都作为待定系数。令 f (x) = 1, x,