1、第 4 章 复习与思考题 习题 1、给出计算积分的梯形公式及中矩形公式，说明它们的几何意义。 答： 用两端点的算术平均值作为 ()f 的近似值，这样导出的求积公式 ( ) d ( ) ( ) 2aabaf x x f a f b ，就是梯形求积公式 。 而如果改用区间中点 2bac 近似取代 ()f ，则导出中矩形公式 ( ) d ( ) ( )2aaabf x x b a f 几何意义的图形，略 。 2、什么是求积公式的代数精确度？梯形公式及中矩形公式的代数精确度是多少？ 答：如果某个求积公式对次数不超过 m 的多项式均能准确成立，但对于 m+1 次多项式就不准确成立，则称该求积公式具有 m
2、 次代数精度 梯形公式和中矩形公式的代数精度为 1. 3、对给定求积公式的节点，给出两种计算求积系数的方法。 由于是给定求积公式的节点，因此，不能使用高斯型求积公式 由于未说明是等距节点，因此不能用牛顿 -科特斯求积公式。 未找到明确的资料 . 答： 插值型求积公式和 4、 什么是牛顿 -柯特斯求积？它的求积节点如何分布？它的代数精确度是多少？ 答：设积分区间 a,b划分为 n 等份，步长 h=（ b-a） /n,选取等距节点 kx a kh 构造出的插值型求积公式 0( ) ( )n kn n kkI b a C f x 成为牛顿 -柯特斯求积 公式 ，式中 knC 称为柯特斯系数。 其节点
3、是等距分布的，代数精度为 节点数 n-1 次 。 5、 什么是辛普森求积公式？它的余项是什么？它的代数精确度是多少？ 当 n=2 时，牛顿 -柯特斯求积公式即为辛普森求积公式，其余项为 4 ( 4 )( ) ( ) ( )180 2b a b aR f f 其，代数精度为 3. 6、 什么是复合求积法？给出复合梯形公式及其余项表达式。 答：为了提高计算精度，通常把积分区间分成若干子区间（通常是等分），再在每个子区间上使用低阶求积公式。这种方法称为复合求积法。 复合梯形公式为 111012( ) ( ) ( a ) ( ) ( ) 22( ) ( ) , a , b32nnk k kkknnhh
4、T n f x f x f f x f bb a hR f I T f 7、 给出复合辛普森公式及其余项表达式。如何估计它的截断误差？ 复合辛普森公式为 111 / 2014( 4 ) ( a ) 4 ( ) 2 ( ) ( ) 6( ) ( ) , , b180 2nnkkkknnhSn f f x f x f bb a hR f I S f a 8、 什么是龙贝格求积？它有什么优点？ 龙贝格求积公式也称为逐次分半加速法。它是在梯形公式、 辛普森 公式和柯特斯公式之间的关系的基础上，构造出一种加速计算积分的方法。 使用理查森 外推算法 , 它在不增加计算量的前提下提高了误差的精度 . 在等距
5、基点的情况下，用计算机计算积分值通常都采用把 区间 逐次分半的方法进行。这样，前一次分割得到的函数值在分半以后仍可被利用 ,且易于 编程 。 龙贝格算法公式 ( ) ( 1 ) ( )1141 , 1 , 2 , 3 . . .4 1 4 1mk k km m mmmT T T k 9、 什么是高斯型求积公式？它的求积节点是如何确定的？它的代数精确度是多少？为何称它是具有最高 代数精确度的求积公式？ 答： 如果求积公式 0( ) ( ) d ( )an nkkka f x x x A f x 具有 2n-1（ n 为秋季节点数）次代数精度，则称其节点 kx 为高斯点，求积公式为高斯型求积公式。
6、 可以使用证明的方法求证，插值公式的代数精度不超过 2n-1。 即回到了最后一个问题。 根据老师的讲课，给出证明的方法 。 10、 牛顿 -柯特斯求积和高斯求积的节点分布有什么不同？对同样数目的节点，两种求积方法哪个更精确？为什么？ 牛顿 -柯特斯求积节点等距分布 高斯求积的节点分布是插值型多项式的零件。 对同样数目的节点，高斯求积更精确。 11、 描述自动求积的一般步骤。怎样得到所需的误差估计？ 答：如果求积 区间中被积函数变化很大，有的部分函数值变化剧烈，需要使用小不长，另一部分函数值变化平缓，可以使用大步长，针对被积函数在区间上的不同情形采用不同的步长，使得在满足精度前提下积分计算工作量
7、尽可能小，针对这类问题的算法技巧是在不同区间上预测被积函数变化的剧烈程度确定响应步长。就是自动求积的一般步骤。 12、 怎样利用标准的一维求积公式计算矩形域上的二重积分 基本原则：累次积分。 多重积分的辛普森公式： 110 1 / 201( , ) ( , ) 4 ( , ) 2 ( , ) ( , ) 6bdacMMb b b bi i Ma a a aiif x y dy dxk f x y dx f x y dx f x y dx f x y dx 对每一个积分再次利用辛普森公式 111 / 201 ( ) 4 ( ) 2 ( ) ( ) ( 6)nnkkka kb h f a f xf
8、 f x f bx d x 13、对给定函数，给出两种近似求导的方法。若给定函数值有扰动，在你的方法中怎样处理这个问题？ 14、 判断如下命题是否正确： （ 1）如果被积函数在区间 a , b 上连续，则它的黎曼（ Riemann）积分一定存在。 （ 2）数值求积公式计算总是稳定的。 （ 3）代数精确度是衡量算法稳定性的一个重要指标。 （ 4） n + 1 个点的插值型求积公式的代数精确度至少是 n 次，最多可达到 2n + 1 次。 （ 5）高斯求积公式只能计算区间 -1, 1上的积分。 （ 6）求积公式的阶数与所依据的插值多项式的次数一样。 （ 7）梯形公式与两点高斯公式精度一样。 （ 8
9、）高斯求积公式系数都是正数，故计算总是稳定的。 （ 9）由于龙贝格求积节点与牛顿 -柯特斯求积节点相同，因此它们的精度相同。 （ 10）阶数不同的高斯求积公式没有公共节点。 1） 正确 2） 错误 3） 错误，是衡量计算准确度的一个指标 4） 正确 5） 错误，可以通过变化使得计算时区间在 -1,1上。 6） 错误，典型的例子是，当 n 为偶数时，牛顿 -柯斯特公式至少为 n+1 阶代数精度。 7） 错误。梯形公式，代数精度为 1，两点高斯公式代数精度为 3 8） 正确 9） 错误。龙贝格精度为 2n，牛顿 -柯特斯精度最大为 n+1 10） 错误。 习题 1、 确定下列求积公式中的待定参数，
10、使其代数精度尽量高，并指明所构造出的求积公式所具有的代数精度。 1） )()0()()(101 hfAfAhfAdxxfhh ； 2） )()0()()(1012 2 hfAfAhfAdxxfhh ； 3） 3/)(3)(2)1()(211 1 xfxffdxxf ； 4） )()0(2/)()0()( 20 hffahhffhdxxfh 。 1） )()0()()(101 hfAfAhfAdxxfhh 由于需要确定 3 个 未知量，因此，需要给定 3 个方程。 设 2( ) 1, ,f x x x ，有 1 0 1113 2 2112023h A A AhA hAh h A h A 解出1
11、0 11012134313h A A AAhAhAh 令 3()f x x ， 有 3 3 3 4 4 34 1 1( ) 0 0 03 3 3 3 3 hhh h hh h h h x dx 令 4()f x x ， 有 4 4 4 5 4 5 54 2 1 2( ) 0 03 3 3 3 5 5h hhhh h hh h h x dx x h 因此，具有 3次代数精度。 2） )()0()()(1012 2 hfAfAhfAdxxfhh 由于需要确定 3 个 未知量，因此，需要给定 3 个方程。 设 2( ) 1, ,f x x x ，有 1 0 1113 2 21140163h A A
12、AhA hAh h A h A 解出1 0 11014834383h A A AAhAhAh 令 3()f x x ， 有 3 3 3 4 4 38 4 8 8 8( ) 0 0 03 3 3 3 3 hhh h hh h h h x dx 令 4()f x x ， 有 24 4 4 5 4 5 528 4 8 16 1 64( ) 0 03 3 3 3 5 5h hhhh h hh h h x dx x h 因此，具有 3次代数精度。 3） 1 121 ( ) ( 1 ) 2 ( ) 3 ( ) / 3f x d x f f x f x 需要确定 2个待定参数，因此，令 设 2( ) 1,
13、,f x x x ，有 1222122 1 2 3 / 30 1 2 3 / 32 1 2 3 / 33xxxx 解出0.68990.28990.68990.51212 266xxxx 令 3()f x x ， 有1143113310 ( )4 ( 1 ) 2 ( - 0. 68 99 ) 3 ( 0. 28 99 ) / 3 1 2 - 0. 68 99 3 0. 28 99 / 3- 0. 52 78 8x x dx f x dxf f f 因此，具有 2次代数精度。 4) )()0(2/)()0()( 20 hffahhffhdxxfh 需要确定 2个待定参数，因此，令 设 2( ) 1
14、, ,f x x x ，有 223330022/ 2 23hhhhh h ah 解出 112a ， h 为任意常数 令 3()f x x ，有 4 3 4 4 40011( ) / 24 4 4hhh x d x f x d x h h h 令 4()f x x ，有 5 4 5 5 50011( ) / 25 3 6hhh x d x f x d x h h h 所以代数精度为 3. 2、 分别用梯形公式和辛普森公式计算下列积分： (1)1 20 ,84 x dx nx (2)91 ,4xdx n (3) /6 20 4 sin , 6dn (1)1 20 ,84 x dx nx 梯形公式
15、11 ( ) 2 ( ) ( ) 2nnkkhT f a f x f b 8n ,所以 , 0 , 1 , 2 , 3 , 4 , 5 , 6 , 7, 88k kxk 0123456780 0.03110.06150.09060.11760.1423()()()()()()() .16440( 0.18360)() .200fxfxfxfxfxfxfxfxfx所以有 18 1 ( 0 ) 2 ( ) (1 ) 20 . 1 1 1 4nkkhT f f x f 辛普森公式 1 1 11 / 20 0 1 ( ) 4 ( ) 2 ( ) ( ) 6n n nn k kk k khS f a f
16、 x f x f b 4n ,所以 , 0 ,1, 2 , 3, 44k kxk 1 / 2 1 , 0 ,1, 2 , 384k kxk 所以 1 1 14 1 / 20 0 1 ( 0 ) 4 ( ) 2 ( ) ( 1 ) 60 . 1 1 1 5 7n n nkkk k khS f f x f x f (2)91 ,4xdx n 梯形公式 11 ( ) 2 ( ) ( ) 2nnkkhT f a f x f b 4n ,所以 1 2 , 0 ,1, 2 , 3 , 4kx k k 1 . 0 0 0 0 1 . 7 3 2 1 2 . 2 3 6 1 2 . 6 4 5 8 3 . 0
17、 0 0 , 0 , 1 . .0 .( ) 4k kfx 14 1 ( ) 2 ( ) ( ) 21 7 .2 2 7 8nkkhT f a f x f b 辛普森公式 1 1 11 / 20 0 1 ( ) 4 ( ) 2 ( ) ( ) 6n n nn k kk k khS f a f x f x f b 2n ,所以 1 4 , 0 ,1, 2 ,kx k k 1 / 2 3 4 , 0 ,1kx k k 所以 1 1 12 1 / 20 0 1 ( ) 4 ( ) 2 ( ) ( ) 61 7 . 3 3 2 1n n nkkk k khS f a f x f x f b (3) /
18、6 20 4 sin , 6dn 梯形公式 11 ( ) 2 ( ) ( ) 2nnkkhT f a f x f b 6n ,所以 , 0 , 1 , 2 , 3 , 4 , 5 , 636kx k k 2 . 0 0 0 0 1 . 9 9 8 1 1 . 9 9 2 4 1 . 9 8 3 2 1 . 9 7 0 5 1 . 9 5 4 8 1 . 9 ,( ) 0 , 1 . . .35 66k kfx 16 1 ( ) 2 ( ) ( ) 21 .0 3 5 6nkkhT f a f x f b 辛普森公式 1 1 11 / 20 0 1 ( ) 4 ( ) 2 ( ) ( ) 6n
19、n nn k kk k khS f a f x f x f b 3n ,所以 , 0 ,1, 2 , 318kx k k 1 / 2 , 0 , 1 , 236 18kx k k 所以 1 1 13 1 / 20 0 1 ( ) 4 ( ) 2 ( ) ( ) 61 . 3 5 7 7n n nkkk k khS f a f x f x f b 3、 直接验证柯特斯公式（ 2.4）具有 5 次代数精度。 验证以下公式在 5()f x x 时，等式成立 ，在 6()f x x 时，等式不成立 0 1 2 3 4( ) 7 ( ) 3 2 ( ) 1 2 ( ) 3 2 ( ) 7 ( ) 90b
20、a baf x d x f x f x f x f x f x 验证过程同题 1。 2345( ) 1, , , , ,f x x x x x x 时， （以 5()f x 为例） badxxabbbaabaabbabbababaaabbbabbababaababbababaababbababaaaabbbababaaab56642245432234555432234554322345543223455555556)2)(6)151515151515(907)2434052709015(321)510105(83)1590270405243(3217907)43(32)2(12)43(3279
21、06()f x x 时， badxxabbabbabababaaabbbabbabababaababbabababaababbabababaaaabbbababaaab677654233245666542332456654233245665423324566666666)648 2 5324 0 51 2 81 7 1 0161 9 51 2 81 7 1 0324 0 5648 2 5(907)7 2 91 4 5 81 2 1 55 4 01 3 518(1 2 81)61520156(163)181 3 55 4 01 2 1 51 4 5 87 2 9(1 2 817907)43(32
22、)2(12)43(32790从而此求积公式最高具有 5次代数精度。 4、 用辛普森公式求积分 10 dxex并估计误差。 解： 辛普森公式 ( ) 4 ( ) ( ) 62b a a bS f a f f b 0.51(0) 1(0.5)(1)ffefe所以 0. 5 11 ( 0) 4 ( 0.5 ) (1 ) 61 1 4 60.63233S f f fee 余项为 44( ) ( ) ( ) , ( a , b)180 2b a b aR f f 40( ) 1f e e 0 .0 61|( 2 5 1 |= 3 .4 7 e - 4) | | 180Rf 5、 推导下列三种矩形求积公式
23、： 2)(2 )()()( abfabafdxxfba ； 2)(2 )()()( abfabbfdxxfba ；3)(24 )()(2()( abfabbafdxxfba 解：本题有 2 中证明方式。 解法 1， 利用插值型求积公式的余项公式证明。 1 11() ( ) ( ) ( )( 1 ) !mbb mmmaa fR f R x d x w x d x Kfm 其中 2 2 1011 ()( 1 ) ! 2 nm m mkkkK b a A xmm 由于只用到 1个插值点 x=a，此时 m=0， n=0，此时插值基函数为 1. 00 ( ) 1 ( )bbaaA l x d x d x
24、 b a 所以 有2 2 1022211 ()( 1 ) ! 21 ( ) ( )21 ()2nm m mkkkK b a A xmmb a b a aba 同理 当用到 1个插值点 x=b，此时 m=0， n=0，此时插值基函数为 1. 00 ( ) 1 ( )bbaaA l x d x d x b a 所以 有2 2 1022211 ()( 1 ) ! 21 ( ) ( )21 ()2nm m mkkkK b a A xmmb a b a bba （ 1）和（ 2）即可证得。 但对于（ 3）， 3)(24 )()(2()( abfabbafdxxfba 使用此种方法 需要调整 。 用 到
25、1个插值点 2bax ，此时 m=0， n=0，此时插值基函数为 1. 00 ( ) 1 ( )bbaaA l x d x d x b a 所以 有2 2 102211 ()( 1 ) ! 21 ( ) ( ) ( )220nm m mkkkK b a A xmmbab a b a 说明对 m=0时，求积公式没有误差。 令 m=1， n=0，此时插值基函数为 1. 00 ( ) 1 ( )bbaaA l x d x d x b a 所以 有2 2 103 3 23 3 2 2311()( 1 ) ! 211( ) ( ) ( )2 3 21 1 1 10. 5 ( )12 12 4 41()2
26、4nm m mkkkK b a A xmmbab a b ab a a b abba 即可证得。 解法 2， 利用微分中值定理证明。 由微分中值定理有： )()()( axfafxf ，从而 22)(2 )()()(2 )()()()()(abfabafaxfxafdxaxfafdxxf bababa 再由微分中值定理有： )()()( bxfbfxf ，从而 22)(2 )()()(2 )()()()()(abfabbfbxfxbfdxbxfbfdxxf bababa 。 由微分中值定理有： 2)2(2 )()2)(2()2()( baxfbaxbafbafxf ，从而33322)(24)(
27、)(2(4)(6)()(2()2(6)()2)(2(21)2()2(2)()2)(2()2()(abfabbafabfabbafbaxfbaxbafxbafdxbaxfbaxbafbafdxxfbababa 6、若用复合梯形公式计算积分 10 xI e dx，问区间 0,1应分多少等份才能使截断误差不超过 51 102 ？若改用辛普森公式，要达到同样精度区间 0,1应分多少等份 ？ 本题 为书上的 例题 。 考核的是 复合 梯形公式和 复合 辛普森公式的误差与节点的关系。 解： 复合梯形公式 余项 2( ) ( )12n baR f h f ( ) e , 0,1f 2 2 51 1 1| (
28、 ) | | ( ) | | ( ) | 1012 12 2n baR f h f en 解得： 213 212.849n 复合辛普森公式余项 4 ( 4 )( ) ( )2880n baR f h f ( 4 ) ( ) e , 0,1f 4 ( 4 ) 4 51 1 1| ( ) | | ( ) | | ( ) | 10288 0 288 0 2n baR f h f en 解得： 4 3.71n ，需要 8 个节点。 7、 如果 ( ) 0fx ，证明用梯形公式计算积分 10 xI e dx所得结果比准确值 I 大，并说明其几何意义 。 解： ( ) 0fx ，说明 ()fx在 0,1区
29、间上是一个内凹的函数。 本题，选择梯形公式余项证明： 11 20 0( ) ( ) ( ) ( ) 012n k kkbaR f f x d x A f x h f 所以 110 0( ) ( )kkkI f x d x A f x，即 梯形公式计算积分 10 xI e dx所得结果比准确值 I大 。 几何解释 8、 用龙贝格方法计算下列积分，要求误差 不超过 510 。 （ 1） 102 xe dx （ 2） 20 sinx xdx（ 3） 3 20 1x x dx解： 龙贝格 求积算法 （本题还不会，特别是求 T0,还需要多看书） ( ) ( 1 ) ( )11414 1 4 1mk k
30、km m mmmT T T 计算 T0 时的递推公式 1210 21 ()22nnn kkhT T f x （ 1） 102 xe dx T 表如下 K h T0 T1 T2 T3 T4 0 1 0.6839397 1 0.5 2 0.25 3 0.125 4 0.0625 9、用辛普森公式的自适应积分计算 1.5 21 lnx xdx，允许误差 310 。 本题考核 自适应积分 , 计算量较大。 关于步长的确定还需要再复习一下 。 10、 是构造高斯型求积公式 1 0 0 1 101 ( ) ( ) ( )f x d x A f x A f xx 本题同例 9，考核高斯节点的求法 和高斯公式
31、的构造方法 。 解：具有最高代数精度的求积公式是高斯型求积公式，其节点为关于权函数 1()xx 的正交多项式的零点 01,xx，设 201( ) ( ) ( )w x x x x x x b x c ，由于正交性 ()wx 与 1,x 带权正交，即得 ： 1100( ) 0 , ( ) 0w x d x x w x d xxx于是有 1 201 201 2 2( ) 2 0532 2 2( ) 07 5 3x b x c d x b cxx x b x c d x b c 解得 -6/73/35bc所以， 2 63() 7 3 5w x x x 解得： 010.11560.7416xx 由于两
32、个节点的高斯型求积公式具有 3 次代数精度，故公式对 23( ) 1, , ,f x x x x 成立， 为方便计算 A0， A1 取前两个。 当 ( ) 1fx 时， 101 0110 0 1 1 0 1 001 2120.11 56 0.71 463A A dxxA x A x A A x dx x dxx 解得： 011.273.7270AA 所以高斯型求积公式为 101 ( ) 1 . 2 7 3 (0 . 1 1 5 6 ) 0 . 7 2 7 (0 . 7 4 1 6 )f x d x f fx 11、 用 2,3n 的高斯 -勒让德公式计算积分 31 sinxe xdx 本题考查
33、积分区间的变化，高斯 -勒让德公式。 解： 高斯 -勒让德积分区间为 -1,1,因此需要进行区间变化， 则 31 2111 21si n si n( 2) ( 2)si n( 2)xxxe x dx e x d xe x dx 即，变化区间后 2( ) sin( 2 )xf x e x 当 n=2 时， 为 3 点的高斯型求积公式， 由表 4-7，高斯节点为 0120.774596700.7745967xxx系数为 0210 .5 5 5 5 5 5 60AAA 所以高斯 -勒让德即公式为 1 21 s in ( 2 ) 0 . 5 5 5 5 5 5 6 ( 0 . 7 7 4 5 9 6
34、7 ) 0 . 8 8 8 8 8 8 9 ( 0 ) 0 . 5 5 5 5 5 5 6 ( 0 . 71 0 . 9 4 7 4 5 9 6 7 )84xe x d x f f f 同理可得 3 点的高斯 -勒让德计算 。 不再计算。 12、 卫星轨道是一个椭圆，椭圆周长的计算公式是 2/2 204 1 s i ncS a da ，这里 a是椭圆的半长轴， c 是地球中心与轨道中心（椭圆中心）的距离，记 h 为近地点距离， H 为远地点距离， 6371R 公里为地球半径，则 (2 ) / 2a R H h ， ( ) / 2c H h 。我国第一颗人造卫星近地点距离 439h 公里，远地点
35、距离为 2384H 公里，试求卫星轨道的周长。 本题属于沿线积分求周长的题型 解： 5.77822 43923846371222 hHRa ， 5.9722 4392 3 8 42 hHc 22/ 2 / 222009 7 2 . 54 1 s i n 3 1 1 3 0 1 s i n7 7 8 2 . 5cS a d da 则有 2 29 7 2 .51 s i n7782) .5(f 使用复合梯形公式或复合辛普森公式求解。 设梯形公式节点为 n=6，辛普森公式 n=3，所以 50 , , , , , ,12 6 4 3 12 2k 1 . 0 0 0 0 0 . 9 9 9 5 0 .
36、9 9 8 0 0 . 9 9 6 1 0 . 9 9 4 1 0 . 9 9 2 7 0 .() 9922kf 所以 1143 1 1 0 ( 0 ) ( ) ( ) 224 . 8 7 0 7 1 0nkhS f f f 使用辛普森公式 1/ 20, , ,6 3 25,12 4 12kk 1 . 0 0 0 0 0 . 9 9 8 0 0 . 9 9 4 1 0 . 9 2) 2( 9kf 1 / 2 0 .9 9 9 5 0 .9 9 6 1 0() .9 9 2 7 kf 111 / 20143 1 1 0 ( 0 ) 4 ( ) 2 ( ) ( ) 624 . 8 7 0 8 1
37、0nnkkkkhS f f f f 两者计算精度基本相当 。 13、 证明等式 3524s in 3 ! 5 !n n n n ，试依据 si n ( 3, 6 ,12)nnn 的值，用外推算法求 的近似值。 14、 用下列方法计算积分 311dyy ，并比较结果。 1）龙贝格方法；（ 2）三点及五点高斯公式； 3）将积分区间分为四等分，用复化两点高斯公式。 三点及五点高斯公式； 需要将区间化为 -1,1 用复化两点高斯公式 , 需要将区间化为 -1,1 15、用 n=2 的高斯 -拉盖尔求积公式计算积分 20 1xxe dxe 本题需要先构造高斯 -拉盖尔求积公式的形式，随后查表即可求得 2
38、200 111x xxxe d x e d xee ，所以21() 1 xfx e 。 查表 4-8， n=2 时有 0120.4157745572.2942803606.289945083xxx0120.7110930100.2785177340.010389257AAA则： 0120.696680.989931xxx有： 0 0 1 1 2 2200 . 7 8 1 51 (1) ( ) ( )1x xe d x A f x A f x A f xe 16、用辛普森公式（取 N=M=2）计算二重积分 0.5 0.500 yxe dydx 本题考核多重积分的求法，考核多重积分的辛普森公式，对
39、其中的单个积分再次应用辛普森公式求积。 解： 多重积分的辛普森公式： 110 1 / 201( , ) ( , ) 4 ( , ) 2 ( , ) ( , ) 6bdacMMb b b bi i Ma a a aiif x y dy dxk f x y dx f x y dx f x y dx f x y dx 对每一个积分再次利用辛普森公式 111 / 201 ( ) 4 ( ) 2 ( ) ( ) ( 6)nnkkka kb h f a f xf f x f bx d x 解：取 N=M=2，即 h=k=0.25 得 0 . 5 0 . 5000.5 0 . 1 2 5 0 . 3 7 5
40、 0 . 2 5 0 . 500 . 5 0 . 5 0 . 5 0 . 5 0 . 50 . 1 2 5 0 . 3 7 5 0 . 2 5 0 . 50 0 0 0 01 4 4 2 241 4 4 2 24yxb b b bx x x x xa a a ax x x x xe d y d xe d x e d x e d x e d x e d xe d x e e d x e e d x e e d x e e d x 对 0.50 xe dx再次使用辛普森求积公式有 0 .500 .1 2 5 0 .3 7 5 0 .2 5 0111 / 201.5 ( ) 4 ( ) 2 ( )
41、( ) 61 1 4( 2 240.)3935nnkkkxkh f a f x f x fe dxe e e eb 所以 0 . 5 0 . 5000 . 50 . 1 2 5 0 . 3 7 5 0 . 2 5 0 . 500 . 50 . 1 2 5 0 . 3 7 5 0 . 2 5 0 . 500 . 1 2 5 0 . 3 7 5 0 . 2 5 0 . 51 4 4 2 241( 1 4 4 2 )241( 1 4 4 2 ) 0 .30 .293245535yxb b b bx x x x xa a a axe d y d xe d x e d x e d x e d x e d
42、 xe e e e e d xe e e e 17、 确定数值微分公式的截断误差表达式 0 0 0 01( ) 4 ( ) 3 ( ) ( 2 ) 2f x f x h f x f x hh 此为三点求导 前插公式 ，根据式 8.6，有 截断误差 2 ()3h f 18、用三点公式求21() 1fx x 在 1.0,1.1,1.2x 处的导数值，并估计误差。 ()fx的值由下表给出： X 1 1.1 1.2 f（ x） 0.25 0.2268 0.2066 本题给出了数表，应属于考核插值型求导公式。 题目明确需要用三点公式求解。 三点求导公式如下： 0 0 1 21 0 22 0 1 21( ) 3 ( ) 4 ( ) ( ) 21( ) ( ) ( ) 21( ) ( ) 4 ( ) 3 ( ) 2f x f x f x f xhf x f x f xhf x f x f x f xh 解： 1(1 ) 3 (1 .0 ) 4 (1 .1 )
