北京科技大学计算方法试题.docxVIP

《计算方法》2008试题与答案

一、填空题（每空2分，共20分）

为了提高数值计算精度,当正数x充分大时,应将ln(x? x2?1)改写为

_?ln(x? x2?1) .

3x*的相对误差约是x*的相对误差的_1/3

3x*

?1 5 ?2?

(3).设A???2 1 0?，则A = 13 .A ???14

? ? ? 1

??3 ?8 2??

(4)已知p(x)为二次多项式，满足P(?2)?f(?2)?3,P(?1)?f(?1)?1和

P(?1)?f(?1)?1，则p(x)?f(?2)?a(x?2)?b(x?2)(x?1)，这里a= -2 ,

b= 3。

(5)设f(x)?4x3?3x2?2x?1，则差商f?0,1,2,3?= 4 f?0,1,2,3,4??_0_.

n个求积节点的求积公式的代数精确度最高为_2n?1 次.

求解初值问题y??50(y?x),y(0)?1时，若用改进欧拉方法的绝对稳定域中步长h不超过.0.04

二、（10分）用Newton法求方程x?lnx?2在区间(2,?)内的根,取x ?3，要

0

x ?xk?1xk求 ?10

x ?x

k?1

x

k

k

解此方程在区间(2,?)内只有一个根s，而且在区间（2，4）内。设

则 f(x)?1?

Newton法迭代公式为

f(x)?x?lnx?2

?1 1

?

， f(x)

x x2

x ?lnx ?2 x

(1?lnx)

kkx ?x ?

k

k

k?1 k

1?1/x

k

? k x ?1kk

， （5分）

取x ?3，得x

0 1

?3.1479184

?0.049306145

x?x1

x?x

1

x

0

0

2

?0.00054797894

x?x2

x?x

2

x

1

1

3

?010?7

x?x3

x?x

3

x

2

2

4

?10?8

x?x4x33s?

x?x

4

x

3

3

4

三.、（20分）分别用 Jacobi 迭代与高斯-赛德尔迭代法解线性方程组

?2x?x?x ??5

?x1 2 3

?

?x1

5x?x

2 3

?8 ，给出迭代格式与迭代矩阵，说明上述迭代是否收敛，并使用收敛迭

? ?x?10x ?11

1 2 3

代公式计算2步，每步结果保留4位小数，取x?0???1,1,1?T。

?x?k?1??0.5x?k??0.5x?k??2.5

?1 2 3

解：本问题的Jacobi迭代格式为?x?k?1???0.2x?k??0.2x?k??1.6

?2 1 3

??x?k?1???0.1x?k??0.1x?k??1.1

迭代矩阵为G

J

? 0

????0.2

?

3

0.5 0.5?

?0 0.2?

?

1 2

G ?0.7?1 此迭代收敛

J1

? ???0.1 ?0.1 0

? ?

取初始迭代向量为x?0???1,1,1?T

，得到x?1????1.5,1.6,0.9?T，x?2????1.25,2.08,1.09?T，

本问题的高斯-赛德尔迭代格式为

?x?k?1??0.5x?k??0.5x?k??2.5

?1 2 3

?x?k?1???0.2x?k?1??0.2x?k??1.6??0.1x?k??0.1x?k??2.1

?2 1 3 2 3

31 2 2 3??x?k?1???0.1x?k?1??0.1x?k?1??1.1??0.04x?k??0.06x?k??

3

1 2 2 3

?0 0.5 0.5?

迭代矩阵为G

??0 ?0.1 0.1?

G ?0.66?1 此迭代收敛

? ?s ? ? s1

? ?

? ?0

? ?

取初始迭代向量为x?0???1,1,1?T

，得到x?1???1.5,2.1,1.04?T，

x?2???0.93,1.994,0.9936?T，

.

四（10分）已知sin(450)?

，sin(600)? ，sin(90