数值分析第五版全答案c h a p4第四章 数值积分与数值微分1.确定下列求积公式中的特定参数，使其代数精度尽量高，并指明所构造出的求积公式所具有的代数精度：10121012112120(1)()()(0)();(2)()()(0)();(3)()[(1)2()3()]/3;(4)()[(0)()]/2[(0)()];hh hhh f x dx A f h A f A f h f x dx A f h A f A f h f x dx f f x f x f x dx h f f h ah f f h -----≈-++≈-++≈-++''≈++-⎰⎰⎰⎰ 解：求解求积公式的代数精度时，应根据代数精度的定义，即求积公式对于次数不超过m 的多项式均能准确地成立，但对于m+1次多项式就不准确成立，进行验证性求解。

（1）若101(1)()()(0)()hh f x dx A f h A f A f h --≈-++⎰ 令()1f x =，则1012h A A A -=++令()f x x =，则110A h A h -=-+令2()f x x =，则3221123h h A h A -=+ 从而解得011431313A hA h A h -⎧=⎪⎪⎪=⎨⎪⎪=⎪⎩令3()f x x =，则3()0hhh h f x dx x dx --==⎰⎰ 101()(0)()0A f h A f A f h --++=故101()()(0)()hh f x dx A f h A f A f h --=-++⎰成立。

令4()f x x =，则4551012()52()(0)()3hh h h f x dx x dx h A f h A f A f h h ---==-++=⎰⎰故此时，101()()(0)()hh f x dx A f h A f A f h --≠-++⎰故101()()(0)()hh f x dx A f h A f A f h --≈-++⎰具有3次代数精度。

（2）若21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰令()1f x =，则1014h A A A -=++令()f x x =，则110A h A h -=-+令2()f x x =，则32211163h h A h A -=+ 从而解得011438383A h A h A h -⎧=-⎪⎪⎪=⎨⎪⎪=⎪⎩令3()f x x =，则22322()0hh h h f x dx x dx --==⎰⎰101()(0)()0A f h A f A f h --++=故21012()()(0)()hh f x dx A f h A f A f h --=-++⎰成立。

令4()f x x =，则22452264()5h h h h f x dx x dx h --==⎰⎰ 510116()(0)()3A f h A f A f h h --++= 故此时，21012()()(0)()h h f x dx A f h A f A f h --≠-++⎰因此， 21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰具有3次代数精度。

（3）若1121()[(1)2()3()]/3f x dx f f x f x -≈-++⎰ 令()1f x =，则1121()2[(1)2()3()]/3f x dx f f x f x -==-++⎰令()f x x =，则120123x x =-++令2()f x x =，则22122123x x =++从而解得120.28990.5266x x =-⎧⎨=⎩或120.68990.1266x x =⎧⎨=⎩ 令3()f x x =，则11311()0f x dx x dx --==⎰⎰ 12[(1)2()3()]/30f f x f x -++≠故1121()[(1)2()3()]/3f x dx f f x f x -=-++⎰不成立。

因此，原求积公式具有2次代数精度。

（4）若20()[(0)()]/2[(0)()]hf x dx h f f h ah f f h ''≈++-⎰ 令()1f x =，则0(),hf x dx h =⎰2[(0)()]/2[(0)()]h f f h ah f f h h ''++-=令()f x x =，则200221()21[(0)()]/2[(0)()]2hhf x dx xdx h h f f h ah f f h h ==''++-=⎰⎰令2()f x x =，则23002321()31[(0)()]/2[(0)()]22hhf x dx x dx h h f f h ah f f h h ah ==''++-=-⎰⎰故有33211232112h h ah a =-=令3()f x x =，则340024441()41111[(0)()]/2[(0)()]12244hh f x dx x dx h h f f h h f f h h h h ==''++-=-=⎰⎰ 令4()f x x =，则450025551()51111[(0)()]/2[(0)()]12236hh f x dx x dx h h f f h h f f h h h h ==''++-=-=⎰⎰ 故此时，201()[(0)()]/2[(0)()],12hf x dx h f f h h f f h ''≠++-⎰ 因此，201()[(0)()]/2[(0)()]12hf x dx h f f h h f f h ''≈++-⎰ 具有3次代数精度。

2.分别用梯形公式和辛普森公式计算下列积分：12012101(1),8;4(1)(2),10;(3),4;(4),6;x x dx n x e dx n x n n ϕ-=+-===⎰⎰⎰ 解：21(1)8,0,1,,()84x n a b h f x x=====+ 复化梯形公式为781[()2()()]0.111402k k h T f a f x f b ==++=∑ 复化辛普森公式为7781012[()4()2()()]0.111576k k k k h S f a f x f x f b +===+++=∑∑ 121(1)(2)10,0,1,,()10x e n a b h f x x --===== 复化梯形公式为9101[()2()()] 1.391482k k h T f a f x f b ==++=∑ 复化辛普森公式为99101012[()4()2()()] 1.454716k k k k h S f a f x f x f b +===+++=∑∑(3)4,1,9,2,()n a b h f x =====复化梯形公式为341[()2()()]17.227742k k h T f a f x f b ==++=∑ 复化辛普森公式为3341012[()4()2()()]17.322226(4)6,0,,,()636k k k k h S f a f x f x f b n a b h f x ππ+===+++======∑∑ 复化梯形公式为 561[()2()()] 1.035622k k h T f a f x f b ==++=∑ 复化辛普森公式为5561012[()4()2()()] 1.035776k k k k h S f a f x f x f b +===+++=∑∑ 3。

直接验证柯特斯教材公式（2。

4）具有5交代数精度。

证明：柯特斯公式为01234()[7()32()12()32()7()]90ba b a f x dx f x f x f x f x f x -=++++⎰ 令()1f x =，则01234()90[7()32()12()32()7()]90ba b a f x dx b a f x f x f x f x f x b a -=-++++=-⎰令()f x x =，则2222012341()()21[7()32()12()32()7()]()902b ba a f x dx xdxb a b a f x f x f x f x f x b a ==--++++=-⎰⎰ 令2()f x x =，则23333012341()()31[7()32()12()32()7()]()903bb a a f x dx x dx b a b a f x f x f x f x f x b a ==--++++=-⎰⎰ 令3()f x x =，则34444012341()()41[7()32()12()32()7()]()904bb a a f x dx x dx b a b a f x f x f x f x f x b a ==--++++=-⎰⎰ 令4()f x x =，则45555012341()()51[7()32()12()32()7()]()905bb a a f x dx x dx b a b a f x f x f x f x f x b a ==--++++=-⎰⎰ 令5()f x x =，则56666012341()()61[7()32()12()32()7()]()906bb a a f x dx x dx b a b a f x f x f x f x f x b a ==--++++=-⎰⎰ 令6()f x x =，则012340()[7()32()12()32()7()]90hb a f x dx f x f x f x f x f x -≠++++⎰ 因此，该柯特斯公式具有5次代数精度。

4。

用辛普森公式求积分10x e dx -⎰并估计误差。

解：辛普森公式为[()4()()]62b a a b S f a f f b -+=++ 此时，0,1,(),x a b f x e -===从而有1121(14)0.632336S e e --=++= 误差为4(4)04()()()1802110.00035,(0,1)1802b a b a R f f e ηη--=-≤⨯⨯=∈5。

推导下列三种矩形求积公式：223()()()()();2()()()()();2()()()()();224b a b a ba f f x dxb a f a b a f f x dx b a f b b a a b f f x dx b a f b a ηηη'=-+-'=---''+=-+-⎰⎰⎰证明：(1)()()()(),(,)f x f a f x a a b ηη'=+-∈两边同时在[,]a b 上积分，得()()()()()bba a f x dxb a f a f x a dx η'=-+-⎰⎰ 即2()()()()()2(2)()()()(),(,)ba f f x dxb a f a b a f x f b f b x a b ηηη'=-+-'=--∈⎰ 两边同时在[,]a b 上积分，得()()()()()bba a f x dxb a f a f b x dx η'=---⎰⎰ 即22()()()()()2()(3)()()()()(),(,)22222ba f f x dxb a f b b a a b a b a b f a b f x f f x x a b ηηη'=---''++++'=+-+-∈⎰ 两连边同时在[,]a b 上积分，得2()()()()()()()22222b b b a a a a b a b a b f a b f x dx b a f f x dx x dx η''++++'=-+-+-⎰⎰⎰ 即 3()()()()();224ba ab f f x dx b a f b a η''+=-+-⎰ 6。