数值分析整理版试题及答案
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
例1、 已知函数表
求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解:
(1)
插值基函数分别为
()()()()()()()()()()
1200102121()1211126
x x x x x x l x x x x x x x ----=
==--------
()()()()()()()()
()()021*******
()1211122x x x x x x l x x x x x x x --+-=
==-+---+-
()()()()()()()()()()0122021111
()1121213
x x x x x x l x x x x x x x --+-=
==-+--+-
故所求二次拉格朗日插值多项式为
()
()()()()()()()()()()2
20
2()11131201241162314
121123537623k k k L x y l x x x x x x x x x x x x x ==⎡⎤
=-⨯
--+⨯-+-+⨯+-⎢⎥⎣⎦=---++-=+-∑
(2)一阶均差、二阶均差分别为
[]()()[]()()[][][]010*********
011201202303
,11204
,412
3
4,,5
2,,126
f x f x f x x x x f x f x f x x x x f x x f x x f x x x x x ---===
-----=
==----===
---
故所求Newton 二次插值多项式为
()()[]()[]()()
()()()20010012012,,,35
311126537623P x f x f x x x x f x x x x x x x x x x x x =+-+--=-+
+++-=+-
例2、 设2
()32f x x
x =++,[0,1]x ∈,试求()f x 在[0, 1]上关于()1x ρ=,{}
span 1,x Φ=的最佳平方逼近多项式。
解:
若{}span 1,x Φ=,则0()1x ϕ=,1()x x ϕ=,且()1x ρ=,这样,有
()()()()()()()()1
1
200110
1
1
2011000
1
210
1
,11,
,3
1
23
,,,
,3226
9,324
dx x dx xdx f x x dx f x x x dx ϕϕϕϕϕϕϕϕϕϕ====
====++=
=++=
⎰⎰⎰⎰⎰ 所以,法方程为
01123126119234a a ⎡⎤⎡⎤⎢⎥⎢⎥⎡⎤=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥
⎣⎦⎢⎥⎢⎥⎣
⎦⎣
⎦,经过消元得012311
62110123a a ⎡⎤⎡
⎤⎢⎥⎢
⎥⎡⎤=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎣⎦⎣⎦ 再回代解该方程,得到14a =,011
6
a =
故,所求最佳平方逼近多项式为*
111
()46
S x x =
+ 例3、 设()x
f x e =,[0,1]x ∈,试求()f x 在[0, 1]上关于()1x ρ=,{}span 1,x Φ=的最佳
平方逼近多项式。 解:
若{}span 1,x Φ=,则0()1x ϕ=,1()x x ϕ=,这样,有
()()()()()()1
0001
21101
01100
1
001
10,111,3
1
,,2
, 1.7183,1
x x dx x dx xdx f e dx f xe dx ϕϕϕϕϕϕϕϕϕϕ====
=======⎰⎰⎰⎰⎰
所以,法方程为
0111 1.718321112
3a a ⎡⎤
⎢⎥⎡⎤⎡⎤=⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎢⎥⎣⎦
解法方程,得到00.8732a =,1 1.6902a =, 故,所求最佳平方逼近多项式为
*1()0.8732 1.6902S x x =+
例4、 用4n =
的复合梯形和复合辛普森公式计算积分1
⎰
。
解:
(1)用4n =的复合梯形公式
由于2h =,(
)f x =()121,2,3k x k k =+=,所以,有
()()(
)
4
1
3
1
[129]22
22
17.2277
k k T h
f f x f =≈=++=+⨯=⎰∑
(2)用4n =的复合辛普森公式
由于2h =,(
)f x =()121,2,3k x k k =+=,()12
220,1,2,3k x
k k +
=+=,所以,有