龙贝格积分实验报告总结计划.docxVIP

二、 Romberg 积分法 1.变步长 Romberg 积分法的原理 复化求积方法对于提高精度是行之有效的方法， 但复化公式的一个主要缺点在于要事先估计出部长。若步长过大，则精度难于保证；若步长过小，则计算量 又不会太大。 而用复化公式的截断误差来估计步长， 其结果是步长往往过小， 而 且 f ( x) 和 f (4) ( x) 在区间 [ a, b] 上的上界 M 的估计是较为困难的。 在实际计算中 通常采用变步长的方法，即把步长逐次分半（也就是把步长二等分） ，直到达到某种精度为止，这种方法就是 Romberg 积分法的思想。 在步长的逐步分半过程中，要解决两个问题： 在计算出 TN 后，如何计算 T2N , 即导出 T2 N 和 TN 之间的递推公式； 在计算出 TN 后，如何估计其误差，即算法的终止的准则是什么。 首先推导梯形值的递推公式，在计算 TN 时，需要计算 N 1个点处的函数值 在计算出 TN 后，在计算 T2N 时，需将每个子区间再做二等分，共新增 N 个节点。 为了避免重复计算，计算 T2 N 时，将已计算的 N 1个点的数值保留下来，只计算 新增 N 个节点处的值。为此，把 T2N 表示成两部分之和，即 1 2 N 1 T2 N h2N [ f ( a) f (b) 2 f (a kh2N )] 2 k 1 1 h2 N [ f (a) N 1 N f (b) 2 f (a kh2 N ) 2 f (a (2k 1)h2 N )] 2 k 1 k 1 1 hN [ f (a) N 1 N f (b) 2 f (a kh2N )] h2N f (a (2k 1)h2 N )] 2 2 k 1 k 1 由此得到梯形值递推公式 1 TN N f (a (2 k 1)h2 N )L L L L L L L L T2N h2 N 1 2 k 1 因此 h1 b a,T1 h1 [ f ( a) f (b)], 2 h2 1 h1, T2 1 T1 h2 f ( a h2 ) 2 2 由复化梯形公式的截断误差有 I TN b a hN2 f ( 1), a 1b 12 I T b a h2 f ( 2 ), a 2 b 2 N 12 2 N 若 f ( x) 变化不大时，即 f ( 1) f ( 2 ) ，则有 I 4T2N TN T2 N 1(T2N TN)LLLLLLLL2 4 1 3 式（ 2）表明，用 T2N 作为定积分 I 的近似值，其误差大致为 1 (T2N TN) ， 3 因此其终止条件为 T2N TN 其中 是预先给定的精度。 2.Romberg 积分公式 将上述方法不断推广下去， 可以得到一个求积分的序列， 而且这个序列很快收敛到所求的定积分。记 TN(0)  TN  ，将区间  N 等分的梯形值。  TN(1)  SN ， 将区间  N 等分的  Simpson TN(2)  C N  ，将区间  N 等分的  Cotes。 TN(3)  RN  ，将区间  N 等分的  Romberg。 由其可构造一个序列 { TN(k) } ，次序列称为 Romberg 序列，并满足如下递推关系： (0) b a (0) 1 (0)b a N b a T1 [ f ( a) f (b)], T2N TN 2N f ( a (2 k 1) ), 2 2 k 1 2N k 4k T ( k 1) T ( k 1) TN 2 N N , k 1,2,L 4k 1 以上递推公式就是 Romberg 积分递推公式。 3.Romberg 积分程序 1. 置 N 1 ，精度要求 ， h1 b a ； 2. 计算 T1(0) b a [ f (a) f (b)] ； 2 3. 置 h2N hN (0) 1 (0) b a N f ( a b a 2 ，并计算 T2N TN 2N (2 k 1) ) ； 2 k 1 2N 4. 置 M N,N 2N,K 1; k ( k 1) ( k 1) 5. 计算 TMk 4 T2 M k TM ； 4 1 6. 若 M 1 ， 则转（ 7）；否则置 M k k 1转（ ）； ， M 2 5 7. 若 T1 ( k ) T1(k 1) ， 则停止计算（输出 T1(k ) ），否则转（ 3）。 4.Romberg 积分法的应用 function [T,n] = romb(f,a,b,eps) double R; if nargin4,eps=1e-8; end h=b-a;R(1,1)=(h/2)*(feval(f,a)+feva