【精选】5.3 Romberg算法.pdf

第4章 数值积分 § 4.3 Romberg算法 k 1,2,3,  § 4.3 Romberg算法 综合前几节的内容,我们知道 梯形公式,Simpson公式,Cotes公式的代数精 度分别为1次,3次和5次 复化梯形、复化Simpson、复化Cotes公式的收 敛速度依次加快 无论从代数精度还是收敛速度,复合梯形公式 都是较差的 有没有办法改善梯形公式呢? 一、复合梯形公式的递推化 b 将定积分I  f (x)dx的积分区间[a,b]分割为n等份 a 各节点为xj a jh , j 0,1, ,n h b a n 复合梯形(Trapz)公式为 b a n1 T [f (a) 2f (x ) f (b)] n 2n j (1) j 1 如果将[a,b]分割为2n等份，而h (b a)/n不变, 则 b a n 1 n1 T [f (a) 2f (x ) 2f (x ) f (b)] 2 n j 1 4n j  j 1 j 0 2 (2) 1 1 其中x 1 xj  h a  (j  )h j  2 2 2 b a n 1 n1 T [f (a) 2f (x ) 2f (x ) f (b)] 2 n j 1 4n j  j 1 j 0 2 b a n1 b a n1 [f (a) 2f (x ) f (b)]  2f (x ) j 1 4n 4n j  j 1