7.3 Romberg积分 7.3.1 Richardson外推法 外推法是用精确度较低的近似公式组合成精 确度较高的近似公式的一种方法 设h≠0是任意数,F(h是关于步长h逼近 的近似公式它们的误差估计式为 F*-F(h)=+k2 h2+k h3+ (1) 这里,k1,k2,3是一组常数 ,1,2,3
7.3 Romberg积分 设 是任意数, 是关于步长h逼近 的近似公式,它们的误差估计式为 (1) h 0 F h( ) * F * 2 3 1 2 3 F F h k h k h k h − = + + + ( ) ... 7.3.1 Richardson外推法 外推法是用精确度较低的近似公式组合成精 确度较高的近似公式的一种方法. 这里, k1 ,k2 ,k3 ,… 是一组常数
按(1)式称F(逼近F的误差为O(h)把h 的幂次称为误差的阶,例如,O(h2称为二阶误差 我们希望找到一种简便的方法,用近似公式F(h)的 组合得到误差阶较高的近似公式F(h使 F-F(h)=k,h+kh'+.(2) 此时F(h逼近F的误差为O(h 类似地,用F(h组合产生逼近F的误差为O(13) 的近似公式等下面我们给出一种具体的组合方法
我们希望找到一种简便的方法,用近似公式F(h)的 组合,得到误差阶较高的近似公式 ,使 (2) 此时, 逼近 F * 的误差为O(h2 ) 类似地,用 组合产生逼近F * 的误差 为 O(h3 ) 的近似公式等.下面我们给出一种具体的组合方法. F h( ) * F O h( ) 2 O h( ) ~ F h( ) ~ * ' 2 ' 3 2 3 F F h k h k h − = + + ( ) ... ~ F h( ) ~ F h( ) 按(1)式,称 逼近 的误差为 .把h 的幂次称为误差的阶,例如, 称为二阶误差
把(1)式改写为 F(h)+kh+k2h2+k3h+(3) 用h/2代替(3)式中的h得 h、,h,h F=F()+k1x+k2+k3+ 48(4 用2乘(4)式再减去(3)试式消去含h的项,得 h F=[F(=)+(F(=)-F(h)]+k2(--h2)+k3(-h3)+ (5) 令F{(h)=F(h)且记 F2(h)=F(=)+[F1(=)-F1(h)
把(1)式改写为 (3) 用h/2代替(3)式中的h,得 (4) 用2乘(4)式再减去(3)式,消去含h的项,得 (5) 令 ,且记 * 2 3 1 2 3 F F h k h k h k h = + + + + ( ) ... 2 3 * 1 2 3 ( ) ... 2 2 4 8 h h h h F F k k k = + + + + 2 3 * 2 3 2 3 [ ( ) ( ( ) ( ))] ( ) ( ) ... 2 2 2 2 h h h h F F F F h k h k h = + − + − + − + 1 F h F h ( ) ( ) = 2 1 1 1 ( ) ( ) [ ( ) ( )] 2 2 h h F h F F F h = + −
那么(5)式可写为 F=F2(h)-ck,h -=k,h 4(6) 这里,2(h逼近F的误差为O(h2) 再用h/2代替h,使(6)式变为 F=F(h-k,h-k,h 32°(7) 用4乘(7式减去(6式消去含h的项得 h、,F2(h/2)-F2(h)11 F=[F2()+ I+=kh'+ 2 8(8 同样记 h、F2(h/2)-F2(h) F3(h)=F2()+ 2
那么(5)式可写为 (6) 这里, 逼近 的误差为 再用 h/2 代替 h , 使(6)式变为 (7) 用4乘(7)式减去(6)式,消去含 的项,得 (8) 同样记 * 2 3 2 2 3 1 3 ( ) ... 2 4 F F h k h k h = − − − 2 F h( ) 2 O h( ) * 2 3 2 2 3 1 3 ( ) ... 8 32 F F h k h k h = − − − 2 h * 3 2 2 2 3 ( / 2) ( ) 1 [ ( ) ] ... 2 3 8 h F h F h F F k h − = + + + 2 2 3 2 ( / 2) ( ) ( ) ( ) 2 3 h F h F h F h F − = + * F
(8)式可以写为g F=F3(h)+k3h3+ 这里F(h逼近F的误差为O(h) 还是用h/2代替h代入(9)式后类似上述过 程,可以得到误差为(h的F(h) 般地,对k=2,3,…,洧有逼近的误差为 O(n递推公式 F∠1(h/2)-F(h) 22-1(10 也称为关于步长h的外推公式 表7-1列出了k=23付按(10)式产生F(h) 的计算次序,表中各列左边黑体数字表示序号
(8)式可以写为 (9) 这里 逼近 的误差为 还是用h/2代替h代入(9)式后,类似上述过 程,可以得到误差为 的 一般地,对 ,有逼近 的误差为 的递推公式 (10) 也称为关于步长h的外推公式. 表7-1列出了 时,按(10)式产生 的计算次序,表中各列左边黑体数字表示序号. * 3 3 3 1 ( ) ... 8 F F h k h = + + 3 F h( ) * F 3 O h( ) 4 O h( ) 4 F h( ) k n = 2,3,..., ( ) k O h * F 1 1 1 1 ( / 2) ( ) ( ) ( ) 2 2 1 k k k k k h F h F h F h F − − − − − = + − k = 2,3,4 ( ) F h k
表7-1 O(h) O(h) O(h) O(h) 1:F1(h)=F(h) 2:F1(=)=F(=)3:F2(h) 4:F(=F()5:F2(∞6:F3(h) h h h 7:F(n)=F()8:F2()9:F(=)10:F4(h 例1设∫(x余项的差分公式为
表7-1 例1 设 带余项的差分公式为 O h( ) 2 O h( ) 3 O h( ) 4 O h( ) 1 1: ( ) ( ) F h F h = 1 2 2 : ( ) ( )3: ( ) 2 2 h h F F F h = 1 2 3 4 : ( ) ( )5: ( )6 : ( ) 4 4 2 h h h F F F F h = 1 2 3 4 7 : ( ) ( )8: ( )9 : ( )10 : ( ) 8 8 4 2 h h h h F F F F F h = ' 0 f x( )
f(x)=[f(x+h)-f(x0-)-f(x) ch h 120 导出具有误差为O(h2的外推公式 解 F1(h)=F(h)=,f(x+h)-f(x0-h) 2h 用h2代替h得 f(x0)=F()-∫( h (5 224 1920 (12) 为消去含h的项用4乘(12)式减去(11试式,得
(11) 导出具有误差为 的外推公式. 解 令 用 h/2代替h,得 (12) 为消去含 的项,用4乘(12)式减去 (11)式,得 2 ' ''' 0 0 0 0 1 ( ) [ ( ) ( )] ( ) 2 6 h f x f x h f x h f x h = + − − − 4 (5) 0 ( ) ... 120 h − − f x 2 ( )j O h 1 0 0 1 ( ) ( ) [ ( ) ( )] 2 F h F h f x h f x h h = = + − − 2 4 ' ''' (5) 0 1 0 0 ( ) ( ) ( ) ( ) ... 2 24 1920 h h h f x F f x f x = − − − 2 h 2 h
(4 3f(x0)=4F()-F(h)+f(x) 160 从而有 f(x)=F2(h)+ (5) 480 +(3 这里 h、,F1(h/2)-F1(h) F2(h)=F(=+ 这时,F2(h近f(x的误差为O(h4 重复用M2代替h并消去含h3的项(=23…,j- (如h4,h,到逼近f(x的误差为O确的 外推公式为
从而有 (13) 这里 这时, 逼近 的误差为 . 重复用h/2代替h并消去含 的项 ,得到逼近 的误差为 的 外推公式为 (4) ' (5) 0 1 1 0 3 ( ) 4 ( ) ( ) ( ) ... 2 160 h h f x F F h f x = − + + (4) ' (5) 0 2 0 ( ) ( ) ( ) ... 480 h f x F h f x = + + 1 1 2 1 ( / 2) ( ) ( ) ( ) 2 3 h F h F h F h F − = + 2 F h( ) ' 0 f x( ) (4) O h( ) 2i h ( 2,3,..., 1) i j = − 4 6 ( , ,...) 如h h ' 0 f x( ) 2 ( )j O h
(b)=(h、F1(h12)-F2(h) )+ j=2,3,,k 4J--1 注意(14)式中第二项的分母为44-而不 是(10)式中的2-1这是由于(11式中的余项 为关于h2的幂次而不是关于h的幂次 732 Romberg求积方法 Romberg求积方法是以复化梯形公式为基 础,应用 Richardson外推法导出的数值求积方 法 回忆72.1节的复化梯形公式分别把积分区
注意(14)式中第二项的分母为 而不 是(10)式中的 .这是由于(11)式中的余项 为关于 的幂次而不是关于h的幂次. 7.3.2 Romberg求积方法 Romberg求积方法是以复化梯形公式为基 础,应用Richardson外推法导出的数值求积方 法. 回忆7.2.1节的复化梯形公式,分别把积分区 1 1 1 1 ( / 2) ( ) ( ) ( ) 2,3,..., 2 4 1 j j j j j h F h F h F h F j k − − − − − = + = − 1 4 1 j− − 1 2 1 j− − 2 h
间[a,b]分为1,2,4等分的结果列入表7-2 表72 ≈2k1 b-a h,=b-a [ f(a)+f(b) 22h2=7(b-a)2U(a)+f(b)+2/(a+b 34h=(b-a)"{f(a)+f(b)+(a+)+ f(a+2h3)+f(a+3h3)}
间[a，b]分为1,2,4等分的结果列入表7-2. 表7-2 k 1 1 2 2 3 4 1 2 k mk − = k k b a h m − = 1 h b a = − 2 1 ( ) 2 h b a = − 3 1 ( ) 4 h b a = − 1 [ ( ) ( )] 2 h f a f b + 2 2 [ ( ) ( ) 2 ( )] 2 h f a f b f a h + + + 3 3 { ( ) ( ) 2[ ( ) 2 h f a f b f a h + + + + 3 3 f a h f a h ( 2 ) ( 3 )]} + + +