数值分析第五版全答案chap4.docxVIP

. . 精选文档 精选文档 第四章 数值积分与数值微分 确定下列求积公式中的特定参数，使其代数精度尽量高，并指明所构造出的求积公式所具有的代数精度： h (1) f ( x)dx A f ( h) 1h 1 2h A0 f (0) A1 f ( h); (2) f (x)dx A f ( h) 12h 1 1 A0 f (0) A1 f (h); (3) f ( x)dx 1 [ f ( 1) 2 f (x1) 3 f ( x2 )]/ 3; h (4) 0 解： f ( x)dx h[ f (0) f (h)]/ 2 ah2[ f (0) f (h)]; 求解求积公式的代数精度时，应根据代数精度的定义，即求积公式对于次数不超过 m 的 多项式均能准确地成立，但对于 m+1 次多项式就不准确成立，进行验证性求解。 h 若 (1) f ( x)dx A 1 f ( h) h A0 f (0) A1 f ( h) 令 f (x) 1 ，则 2 h A 1 A0 A1 令 f (x) x ，则 0 A 1h A1h 令 f (x) 2 x ，则 2 h3 3 h2 A h2 A 1 1 1 1 0A 4 h 0 3 1A 1 h 1 3 1A 1 h 1 3 令 f (x) x3 ，则 h f (x)dx h x3dx 0 h h A 1 f ( h) A0 f (0) A1 f ( h) 0 h 故 f ( x)dx A h 1 f ( h) A0 f (0) A1 f (h) 成立。 令 f (x) 4 x ，则 h f ( x)dx x4dx 2 h5 hh h 5 h A f ( h) A f (0) A f (h) 2 h5 1 0 1 3 故此时， h f (x)dx A 1 f ( h) h A0 f (0) A1 f (h) h 故 f ( x)dx A h 1 f ( h) A0 f (0) A1 f (h) 具有 3 次代数精度。 若 2 h 1f ( x)dx A f ( h) 1 2 h A0 f (0) A1 f ( h) 令 f (x) 1 ，则 4 h A 1 A0 A1 令 f (x) x ，则 0 A 1h A1h 令 f (x) x2 ，则 16 h3 3 从而解得 A0 h2 A h2 A 1 14 1 1 3 1A 8 h 1 3 1A 8 h 1 3 令 f (x)  3 x ，则 2 h f ( x)dx 2 h 2 h x3dx 0 2 h A 1 f ( h) A0 f (0) A1 f ( h) 0 2 h 1故 f ( x)dx A f ( h) 1 2h A0 f (0) A1 f ( h) 成立。 令 f (x) x4 ，则 2 h f ( x)dx 2 h 2 h 64 h5 64 h5 5 (h) 16 4 2 h A f ( h) A f (0) A f h5 31 0 1 3 故此时， 2 h 1f ( x)dx A f ( h) 1 2 h A0 f (0) A1 f (h) 因此， 2 h 1f ( x)dx A f ( h) 1 2 h A0 f (0) A1 f (h) 具有 3 次代数精度。 若 1 f (x)dx 1  [ f ( 1) 2 f  ( x1 ) 3 f  (x2 )]/ 3 令 f (x) 1 ，则 1 f ( x)dx 1 2 [ f ( 1) 2 f ( x1 ) 3 f (x2 )]/ 3 令 f (x) x ，则 0 1 2 x1 3x2 令 f (x) x2 ，则 2 1 2 x2 3 x2 1 2 从而解得 x1 0.2899 x1 或 x2 0.5266 x2 0.6899 0.1266 令 f (x) x3 ，则 1 f ( x)dx x3dx 0 11 1 1 [ f ( 1) 2 f (x1) 3 f ( x2 )]/ 3 0 1 故 f (x)dx 1 [ f ( 1) 2 f ( x1 ) 3 f (x2 )]/ 3 不成立。 因此，原求积公式具有 2 次代数精度。 h 2 若 f (x)dx h[ 0 f (0) f (h)]/ 2 ah [ f (0) f (h)] 令 f (x) 1 ，则 h f (x)dx h, 0 h[ f (0) f (h)]/ 2 ah2[ f (0) f (h)] h 令 f (x) h f (x)dx x ，则 h xdx  1 h2 0 0 2 h[ f  (0)  f (