原创力文档- 1、本文档共2页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
（五）课后习题
4.1 对于积分
a
f ( x)dx ，以 x0
a, x1
0, x2
a 为节点，构造形如
a
a
f ( x)dx
A0 f ( x0 ) A1 f ( x1 ) A2 f (x2 )
a
的插值型求积公式，并讨论所得公式的代数精度。
解答： A0
a ( x
x1 )( x x2 )
a
( x 0)( x a)
1
a (x0
x1 )( x0 x2 )
dx
dxa
a ( a 0)( a a)
3
A1
a
(x
x0 )( x x2 )
a (x
a)( x
a)
4
a ( x1
x0 )( x1
dx
a (0
a)(0
dx
a
x2 )
a)
3
A2
a
( x
x0 )( x
x1 ) dx
a ( x
a)( x
0) dx
1 a
a (x2
x0 )( x2
x1 )
a (a
a)( a
0)
3
易知为 Simpson 公式 ,因此代数精度为 3
4.2 确定 下列求积公式中的待定参数，使其代数精度尽量高，并指出所得公式的代数精度。
2
（ 1）
f ( x)dx
A0
f (0)
A1 f (1)
A2 f ( 2)
0
h
h[ f (0)
f (h)]
h 2 [ f (0)
f ( h)]
（ 2）
f ( x) dx
0
2
解答 ：（ 1）令 f (x)
1, x, x2
，假定求积公式均准确成立，从而有：
2
dx
2
A0
A1
A2
0
2
2
A0
0
A1
1
A2 2
xdx
0
2
8
A0020 A112
A2 22
x 2 dx
0
3
1
4
A0
解以上三元线性方程组从得：
A2,A1
，显然仍为 Simpson 公式 ,因此代数精度为 3
3
3
（ 2）求积公式中只含一个待定参数
，当 f ( x)
1, x 时，有
h
h
[1 1] 0，
h
h [ 0 h]h
2(1 1)
dx
xdx
0
2
0
2
故令 f ( x)
x 2 时求积公式准确成立，即
h
h [ 0 h2 ]
h 2 [2
0 2h] ，解得
1
x 2 dx
0
2
12
将 f ( x)
x3 代入上述确定的求积公式，有：
h
2
3 dx
h [ 0 h3 ]
h [0
3h 2 ] ，这说明求积公式至少有
f ( x) x4
x
3 次代数精度，再令
，代入
0
2
12
求积公式时有：
h
2
4 dx
h [ 0 h4 ]
h [ 0 4h3 ]
x
0
2
12
h
h [ f (0)
h
2
[ f (0)
f (h)]
f ( x)dx
f (h)]
故所建求积公式为
0
2
2
4.3 对于
f ( x)
sin x ，利用下表数据，计算
x
n 4,8 时的复合梯形公式
T4 , T8 ，以及
n
4 复合
Simpson
公式
S4 的值。
k
xk
f ( xk )
k
xk
f (xk )
0
0
1.000 000 0
5
0.625
0.936 155 6
1
0.125
0.997 397 8
6
0.75
0.908 851 6
2
0.25
0.989 615 8
7
0.875
0.877 192 5
3
0.375
0.976 726 7
8
1
0.841 470 9
4
0.5
0.958 851 0
1
解： T4
4 [ f (0)
2
( f (0.25)
f (0.5) f (0.75))
f (1)] 0.9445135
2
1 [ f (0)
7
T8
2
f ( xk ) f (1) ] 0.9456909
16
k
1
S4
1 [ f (0)
4
( f (0.125)
f (0.375)
f (0.625) f (0.875))
24
2 ( f (0.25)
f (0.5)
f (0.75))
f (1)]
0.9460832 （或利用 S4
4 T8
1T4）
3
3
4.4 如果用复化梯形公式计算定积分
e x dx ，要求截断误差不超过
－
4，试问 n 至少取多少？
0.5× 10
解答 ；复化的梯形公式的截断误差为
M 2
max f ( x)
max(e x )
1
0 x 1
0 x 1
RN ( f )
b a h M
.
， n=40.8
，取 n 41。
n
文档评论（0）