《数值分析》习题1

习题1

1． 以下各表示的近似数，问具有几位有效数字？并将它舍入成有效数。

（1）*

1x ＝451.023， 1x ＝451.01；

（2）*

2x ＝－0.045 113， 2x ＝－0.045 18；

（3）*

3x ＝23.421 3， 3x ＝23.460 4；

（4）*4x

＝

3

1

， 4x ＝0.333 3； （5）*

5x ＝23.496， 5x ＝23.494； （6）*

6x ＝96×510， 6x ＝96.1×510； （7）*

7x ＝0.000 96， 7x ＝0.96×310-； （8）*

8x ＝－8 700， 8x ＝－8 700.3。

解：(1) =*

1x 451.023 =1x 451.01

=-1*1

x x 0.0131102

1

-?≤，1x 具有4位有效数字。→1x 451.0

(2) -=*

2x 0.045 113 -=2x 0.045 18

=-

1x x 0.045 18045113.0-=0.000 06731021-?<

2x 具有2位有效数字，045.02-→x

(3)=*

3x 23.4213 =3x 23.4604

=-3*3

x x =-4604.234213.23=-4213.234604.231102

1

0391.0-?≤

3x 具有3位有效数字，4.233→x (不能写为23.5)

(4) =

*

4x 3

1

，=4x 0.3333

=-4*4x x 4102

1

000033.0-?< ,4x 具有4位有效数字，=4x 0.3333

(5) =*

5

x 23.496，=5x 23.494 =-5*

5x x =-494.23496.232102

1002.0-?<

5x 具有4位有效数字， →5x 23.50 (不能写为23.49)

(6) =*

6

x 51096?71096.0?= =6x 5101.96?710961.0?= =-6*

6x x 710001.0-?7210102

1--??≤

6x 具有2位有效数字，57610961096.0?=?=x

(7) =*

7

x 0.00096 371096.0-?=x 3*71096.0-?=x =-7*

7x x 0 7x 精确 (8) 8700*

8

-=x 8x 3.8700-= 8*

8

x x -0102

1

3.0?≤= 8x 具有4位有效数字，8x 8700-=精确

2.以下各数均为有效数字： (1) 0.1062 + 0.947; (3)2.747?6.83;

(2)23.46―12.753; (4)1.473 / 0.064 。

问经过上述运算后，准确结果所在的最小区间分别是什么？ 解：(1) 1x =0.1062，2x =0.947，1x +2x =1.0532

)(1x e 4

1021-?≤，)(2x e 3102

1-?≤

)()()(2121x e x e x x e +≈+≤

+≤)()(21x e x e 34102

1

1021--?+? =0.00055

*2*1x x +∈[1.053200055.0-，1.0532+0.00055]=[1.05265，1.05375]

(2) 1x =23.46， -=2x 12.753 =-21x x 10.707

)(1x e 2

1021-?≤，)

(2x e 3102

1-?≤ )()()(2121x e x e x x e -≈-≤)()(21x e x e +

32

102

11021--?+?≤=0.0055 ∈-*

2*1x x [10.7070055.0-, 10.707+0.0055]=[10.7015,10.7125]

(3) =1x 2.747 =2x 6.83 =21x x 18.76201,

≤)(1x e 3

1021-?, )(2x e 2102

1-?≤

)()()()()(2112211221x e x x e x x e x x e x x x e +≤+≈

??=??+??≤---223102

1

1021747.2102183.6(0.683+2.747)=0.01715

]77916.18,74486.18[]01715.076201.18,01715.076201.18[*

2*1=+-∈x x

(4) =1x 1.473 , =2x 0.064 , =21x x 23.015625

)(1x e 31021-?≤, )(2x e 31021-?≤ )()(1)(22211221x e x x x e x x x

e -≈ =+≤

)()(1)(22211221x e x x x e x x x e 3

231021064

.0473.11021064.01--??+?? =0.187622

∈*2

*

1

x x [23.015625187622.0-, 23.015625+0.187622]

=[22.828003 , 23.203247]

3.对一元2次方程01402=+-x x ，如果≈39919.975具有5位有效数字，求其具有5位有效数字的根。

解：01402=+-x x

399400402=+-x x

=*1x 39920+ , =*2x 399

20139920+=

-

记 =*x 399 ,=x 19.975 )(x e 3102

1

-?≤

=1x 20x +=20+19.975=39.975 )()(21x e x e =3102

1

-?≤

∴

1x 具有5位有效数字。

=2x 025*******.0975.391

975.19201201==+=+x

-

≈)(2x e 2

)

20()(x x e + ，

≤+≈22)20()()(x x e x e 6623102110313.0975.391021---?

121x x 025*******.0975

.391

= 21

12)

()(x x e x e -≈ 2

62112975.391021)()(-?≤≈x x e x e 4.若937.01≈x 具有3位有效数字，问1x 的相对误差限是多？设

x x f -=1)(，求)(1x f 的绝对误差限和相对误差限。

解：=1x 0.937 )(1x e 3102

1

-?≤

≤=111)()(x x e x e r 33

10534.0937

.0102

1

--?=?

x x f -=1)( ,x

x f --=

'121

)(

='≈)()()(x e x f f e 21-

)(11x e x

-? , ≈

))((1x f e 21)(1111x e x -?3310996.0102

1937.01121--?=??-?≤ -≈=f

f e f e r )()()(11

21x e x -?

， ≈

))((1x f e r 21≤-?)(1111x e x 3102

1937.01121-??-? =31097.300397.0-?= 5.取

42.101.2≈，41.100.2≈试按00.201.2-=A 和

)00.201.2(01.0+=A 两种算法求A 的值，并分别求出两种算法所得A 的近似值的绝对误差限和相对误差限，问两种结果各至少具有几位有效

数字？

解：1) 记 =*1x 01.2 ，=1x 1.42 ，*2x =00.2 ，=2x 1.41

则 ≤

)(1x e 21021-? ，)(2x e 2102

1

-?≤ *A 01.041.142.100.201.2=-≈-= 1A 01.041.142.1=-=

)()()()(21211x e x e x x e A e -≈-=

)()()()()(21211x e x e x e x e A e +≤-≈22210102

1

1021---=?+?=

≤

=1

11)()(A A e A e r 101.0102

=- 不能肯定所得结果具有一位有效数字。 2 ) *A =)00.201.2(01.0+， =2A 00353356.083.201.0)41.142.1(01.0==+ )()

(101.0))

(01

.0()(212

21212x x e x x x x e A e ++?

-=+=

)10211021()

41.142.1(1

01.0)(2

222--?+??+?

≤A e 44

102

1

10

12486.0--?

∴ 具有2位有效数字。

≤

=2

22)()(A A e A e r 24

103533547.000353356.01012486.0--?=? 3) 2*121*A A A A A A -+-=- 2*121*A A A A A A ---≥-

24102

1

006.0102101.000353356.0--?>=?--=

∴ 1A 无有效位数。

6.计算球的体积所产生的相对误差限为1%。若根据所得体积的值推算球的半径，问相对误差为多少？

解： π3

4

=V 3R ,π4=dV dR R 2

323

44R dR R V dV ππ==3R dR

)(3

1

)(V e R e r r ≈

由 )(V e r =210- 知 2103

1

)(-?≤R e r

7.有一圆柱，高为25.00 cm ，半径为20.00±0.05 cm 。试求按所给数据计算这个圆柱的体积和圆柱的侧面积所产生的相对误差限。 解：1) h R R V 2)(π=

)(2)(2)()()(2R e R e h

R R

hR R e V R R V V e r r r r =?=?

'≈ππ ≤≈)(2)(R e V e r r 005.02001

2005.02==?

2) π2)(=R S Rh

)()(22)()()(R e R e Rh R h R e S R R S S e r r r r =?=?

'≈ππ )()(R e S e r r ≈0025.020

05

.0=≤

答 计算体积的相对误差限为0.005，计算侧面积的相对误差限为0.0025 9.试改变下列表达式，使计算结果比较精确：

(1) 2

1

cos 1cos 1??

? ??+-x x , 当1<

(2)

x x -+1, 当1>>x 时；

(3) x

x

x +--+11211, 当1<

x

x

sin cos 1-, 当1<

解 : (1) 22cos 22sin 2cos 1cos 12

1222

1

x tg x x x x =?????

? ??=?

?

? ??+- (2)

x

x x x ++=

-+11

1 (3) )1)(21()

21()1()1)(21()21)(1()1(112112x x x x x x x x x x x x x ++-+-+=+++--+=+--+ =)

1)(21(22x x x ++

(4)

22

cos 2sin 22

sin 2sin cos 12

x

tg x x x x x ==-

10.若1个计算机的字长3=n ，基数10=β，阶码22≤≤-p ，问这台计算机能精确表示几个实数。 解：3=n , 0=β, 2-=L , 2=U

所能精确表示的实数个数为

9001)122(10921)1()1(2121=++???+=+--+-L U n ββ

11.给定规格化的浮点数系F:2=β，3=n ,1-=L ,1=U ,求F 中规格化的浮点数的个数，并把所有的浮点数在数轴上表示出来。 解：2=β, 4=n , 1-=L , 1=U

所有规格化浮点数个数为

49)111(2121)1()1(2131=++???+=+--+-L U n ββ

机器零 0

p=1 121000.0?±, 121001.0?±, 121010.0?±, 121011.0?±

121100.0?±, 121101.0?±, 121110.0?±, 121111.0?±

p=0 021000.0?±, 021001.0?±, 021010.0?±, 021011.0?±

021100.0?±, 021101.0?±, 021110.0?±, 021111.0?±

p=1- 121000.0-?±, 121001.0-?±, 121010.0-?±, 121011.0-?± 121100.0-?±, 121101.0-?±, 121110.0-?±, 121111.0-?± 12.设有1计算机：3=n ，2==-U L ,10=β,试求下列各数的机器近似值（计算机舍入装置）： (1) 41.92; (2) 328.7 (3) 0.0483 (4) 0.918; (5)0.007 845; (6)98 740; (7) 1.82310?; (8) 4.71610-?; (9)6.644 52110?; (10) 3.8791010-?;

(11) 3.19610010-?;

(12) 13.6549910? 。

解：3=n , 2-=L , 2=U , 10=β (1) 41.92 (2) 328.7 (3) 0.0483 (4) 0.918 (5) 0.007845 (6) 98740 (7) 1.82310? (8) 4.71610-? (9) 6.64452110? (10) 3.8791010-?

(11) 3.19610010-? (12) 13.6549910?

210419.0)92.41(?=fl

)7.328(fl 溢出

110483.0)0483.0(-?=fl 010918.0)918.0(?=fl 210785.0)007845.0(-?=fl

)98740(fl 溢出

)1082.1(3?fl 溢出 )1071.4(6-?fl 溢出

)106445.6(21?fl 溢出 )10879.3(10-?fl 溢出

)10196.3(100-?fl 溢出 )10654.13(99?fl 溢出

16.考虑数列1，31，91，271

，81

1， 。设10=p ，则用递推公式

13

1

-=n n p p (n =2，3， )

可以生成上述序列。试考察计算n p 的算法的稳定性。 解：13

1

-=

n n p p , ,3,2,1=n 。若0p 有误差，则实际按如下递推 1~

~

3

1-=n p p n

1~1~

3131---=-n n n p p p n p =)(3

1

1~1---n n p p

记 ~

n n n p p e -=, 则有

013

1

31e e e n n n ===-

03

1e e n

n =

13

1

-=n n e e ,误差逐步缩小，数值稳定

17.考虑数列1，31，91，271，811， 。设10=p ，3

1

1=p ，则用递推公式

213

10

---=n n n p p p (n =2，3， ) 可以生成上述序列。试问计算的上述公式是稳定的吗？ 解：213

10

---=

n n n p p p , ,3,2=n 。若0p 和1p 有误差，则实际按如下

递推：

,3

102~1~

~

---=n n n p p p ,3,2=n 。

记 ~

n n n p p e -= ,则有 213

10

---=

n n n e e e ， ,3,2=n )3

1

(3)31(333101121211e e e e e e e e n n n n n n n -=-=-=------- (A)

)3(31

)3(3131301121211e e e e e e e e n n n n n n n -=-=-=------- (B)

9(A)-(B) 得

??

?

???---=-+)3(31)31(381011011e e e e e n n n

只需 03

1

01≠-e e , 则

∞=∞

→n n e lim 因而递推过程不稳定

18.已知47311230125)(2

3

5

-+-+=x x x x x p ，用秦九韶法求)5(p 。 解： 125 0 230 11- 3 47- 5 625 3125 16775 83820 419115

125 625 3355 16764 83823 419068

419068)5(=p

数值分析试题及答案汇总

数值分析试题 一、 填空题（2 0×2′） 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值，则x 有 2 位 有效数字。 2. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= 1 ， f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设，‖A ‖∞＝___5 ____，‖X ‖∞＝__ 3_____， ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ，则使用该迭代 函数的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商 公式的 前插公式 ，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的 后插公式 ；如果要估计结果的舍入误差，应该选用插值公式中的 拉格朗日插值公式 。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=n i i x a 0)( 1 ；所以当 系数a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%，至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收 敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。 11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。 12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的，其中的残差 r i ＝ (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ，(i =0,1,…,n )。 13. 在非线性方程f (x )=0使用各种切线法迭代求解时，若在迭代区间存在唯一解，且f (x )

数值分析学期期末考试试题与答案(A)

期末考试试卷（A 卷） 2007学年第二学期 考试科目： 数值分析 考试时间：120 分钟 学号 姓名 年级专业 一、判断题（每小题2分，共10分） 1. 用计算机求 1000 1000 1 1 n n =∑时，应按照n 从小到大的顺序相加。 （ ） 2. 为了减少误差,进行计算。 （ ） 3. 用数值微分公式中求导数值时，步长越小计算就越精确。 （ ） 4. 采用龙格－库塔法求解常微分方程的初值问题时，公式阶数越高，数值解越精确。（ ） 5. 用迭代法解线性方程组时，迭代能否收敛与初始向量的选择、系数矩阵及其演变方式有 关，与常数项无关。 （ ） 二、填空题（每空2分，共36分） 1. 已知数a 的有效数为0.01，则它的绝对误差限为________，相对误差限为_________. 2. 设1010021,5,1301A x -????????=-=-????????-???? 则1A =_____，2x =______，Ax ∞ =_____. 3. 已知5 3 ()245,f x x x x =+-则[1,1,0]f -= ,[3,2,1,1,2,3]f ---= . 4. 为使求积公式 1 1231 ()()(0)33 f x dx A f A f A f -≈- ++? 的代数精度尽量高，应使1A = ，2A = ，3A = ，此时公式具有 次的代数精度。 5. n 阶方阵A 的谱半径()A ρ与它的任意一种范数A 的关系是 . 6. 用迭代法解线性方程组AX B =时，使迭代公式(1) ()(0,1,2,)k k X MX N k +=+=产 生的向量序列{ }() k X 收敛的充分必要条件是 . 7. 使用消元法解线性方程组AX B =时，系数矩阵A 可以分解为下三角矩阵L 和上三角矩

数值分析试卷及其答案

1、（本题5分）试确定7 22 作为π的近似值具有几位有效数字，并确定其相对误差限。 解 因为 7 22 =3.142857…=1103142857 .0-? π=3.141592… 所以 312102 11021005.0001264.0722--?=?=<=- π （2分） 这里，3,21,0=-=+-=n n m m 由有效数字的定义可知7 22 作为π的近似值具有3位有效数字。 （1分） 而相对误差限 3102 1 0005.00004138.0001264.07 22-?= <≈= -= π π πε r （2分） 2、（本题6分）用改进平方根法解方程组：??? ?? ??=????? ??????? ??--654131*********x x x ； 解 设???? ? ??????? ? ?????? ??===????? ??--11111 1 131321112323121 32 132 31 21 l l l d d d l l l LDL A T 由矩阵乘法得： 5 7,21,215 27 ,25,2323121321- ==-== -==l l l d d d （3分） 由y D x L b Ly T 1 ,-==解得 T T x y )9 23 ,97,910(,)563, 7,4(== （3分） 3、（本题6分）给定线性方程组???????=++-=+-+=-+-=-+17 7222382311387 510432143213 21431x x x x x x x x x x x x x x 1）写出Jacoib 迭代格式和Gauss-Seidel 迭代格式； 2）考查Jacoib 迭代格式和Gauss-Seidel 迭代格式的敛散性； 解 1）Jacoib 迭代格式为

数值分析试卷及答案

二 1 求A的LU分解，并利用分解结果求 解由紧凑格式 故 从而 故 2求证：非奇异矩阵不一定有LU分解 证明设非奇异，要说明A不一定能做LU分解，只需举出一个反例即可。现考虑矩阵，显然A为非奇异矩阵。若A有LU分解，则 故，而，显然不能同时成立。这矛盾说明A不能做LU分解，故只假定A非奇异并不能保证A能做LU分解，只有在A的前阶顺序主子式时才能保证A一定有LU分解。 3用追赶法求解如下的三对角方程组 解设有分解 由公式 其中分别是系数矩阵的主对角线元素及其下边和上边的次对角线元素，故有 从而有 故，，， 故，，， 4设A是任一阶对称正定矩阵，证明是一种向量范数 证明（1）因A正定对称，故当时，，而当时， （2）对任何实数，有 （3）因A正定，故有分解，则 故对任意向量和，总有 综上可知，是一种向量范数。 5 设，，已知方程组的精确解为 （1）计算条件数； （2）若近似解，计算剩余； （3）利用事后误差估计式计算不等式右端，并与不等式左边比较，此结果说明了什么？解（1） （2）

（3）由事后误差估计式，右端为 而左端 这表明当A为病态矩阵时，尽管剩余很小，误差估计仍然较大。因此，当A病态时，用大小作为检验解的准确度是不可靠的。 6矩阵第一行乘以一数成为，证明当时，有最小值 证明设，则 又 故 从而当时，即时，有最小值，且 7 讨论用雅可比法和高斯-赛德尔法解方程组时的收敛性。如果收敛，比较哪一种方法收敛较快，其中 解对雅可比方法，迭代矩阵 ， 故雅可比法收敛。 对高斯-赛德尔法，迭代矩阵 ，故高斯-赛德尔法收敛。 因=故高斯-赛德尔法较雅可比法收敛快。 8设，求解方程组，求雅可比迭代法与高斯-赛德尔迭代法收敛的充要条件。 解雅可比法的迭代矩阵 ， 故雅可比法收敛的充要条件是。 高斯-赛德尔法的迭代矩阵 ， 故高斯-赛德尔法收敛的充要条件是。 9 设求解方程组的雅可比迭代格式为，其中，求证：若，则相应的高斯-赛德尔法收敛。证明由于是雅可比法的迭代矩阵，故 又，故， 即，故故系数矩阵A按行严格对角占优，从而高斯-赛德尔法收敛。 10设A为对称正定矩阵，考虑迭代格式 求证：（1）对任意初始向量，收敛； （2）收敛到的解。 证明（1）所给格式可化为 这里存在是因为，由A对称正定，，故也对称正定。 设迭代矩阵的特征值为，为相应的特征向量，则与做内积，有 因正定，故，从而，格式收敛。

数值分析期末考试复习题及其答案.doc

数值分析期末考试复习题及其答案 1. 已知325413.0,325413* 2* 1==X X 都有6位有效数字，求绝对误差限。（4分） 解： 由已知可知,n=6 5.01021 ,0,6,10325413.0016*1=?= =-=?=ε绝对误差限n k k X 2分 620* 21021,6,0,10325413.0-?=-=-=?=ε绝对误差限n k k X 2分 2. 已知?????=001A 220 - ???? ?440求21,,A A A ∞ （6分） 解： {},88,4,1max 1==A 1分 {},66,6,1max ==∞A 1分 () A A A T max 2λ= 1分 ?????=001A A T 420 ?? ?? ? -420?????001 220 - ?????440=?????001 080 ???? ?3200 2分 {}3232,8,1max )(max ==A A T λ 1分 24322==A 3. 设3 2 )()(a x x f -= （6分） ① 写出f(x)=0解的Newton 迭代格式 ② 当a 为何值时，)(1k k x x ?=+ （k=0,1……）产生的序列{}k x 收敛于2 解： ①Newton 迭代格式为： x a x x x a x a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(2 2 32 1 += +=---=-=+? 3分

②时迭代收敛即当222,112 10)2(',665)('2<<-<-=-=a a x a x ?? 3分 4. 给定线性方程组Ax=b ，其中：? ??=1 3A ??? 22，??????-=13b 用迭代公式)()()()1(k k k Ax b x x -+=+α（k=0,1……）求解Ax=b ，问取什么实数α，可使迭代收 敛 （8分） 解： 所给迭代公式的迭代矩阵为?? ? --? ??--=-=ααααα21231A I B 2分 其特征方程为 0) 21(2)31(=----= -αλα ααλλB I 2分 即，解得αλαλ41,121-=-= 2分 要使其满足题意，须使1)(

数值分析试卷及其答案2

1、（本题5分）试确定7 22作为π的近似值具有几位有效数字，并确定其相对误差限。 解 因为 7 22=3.142857…=1103142857.0-? π=3.141592… 所以 3 12 10 2 110 21005.0001264.07 22--?= ?= <=- π （2分） 这里，3,21,0=-=+-=n n m m 由有效数字的定义可知7 22作为π的近似值具有3位有效数字。 （1分） 而相对误差限 3 10 2 10005.00004138.0001264.07 22-?= <≈= -= π π πε r （2分） 2、（本题6分）用改进平方根法解方程组：???? ? ??=????? ??????? ??--654131321 112321x x x ； 解 设???? ? ? ?????? ? ?????? ??===????? ? ?--11 1 11113 1321 11232312132 1 32 31 21 l l l d d d l l l LDL A T 由矩阵乘法得： 5 7,21,21527,25,2323121321- == - == -==l l l d d d （3分） 由y D x L b Ly T 1 ,-==解得 T T x y )9 23,97,910( ,)5 63, 7,4(== （3分） 3、（本题6分）给定线性方程组??? ? ? ??=++-=+-+=-+-=-+17722238231138751043214321 321431x x x x x x x x x x x x x x 1）写出Jacoib 迭代格式和Gauss-Seidel 迭代格式； 2）考查Jacoib 迭代格式和Gauss-Seidel 迭代格式的敛散性； 解 1）Jacoib 迭代格式为

数值分析试卷及答案

二 1求A的LU分解，并利用分解结果求 解由紧凑格式 故 从而 故 2求证：非奇异矩阵不一定有LU分解 证明设非奇异，要说明A不一定能做LU分解，只需举出一个反例即可。现考虑矩阵，显然A为非奇异矩阵。若A有LU分解，则 故，而，显然不能同时成立。这矛盾说明A不能做LU分解，故只假定A非奇异并不能保证A能做LU分解，只有在A的前阶顺序主子式 时才能保证A一定有LU分解。

3用追赶法求解如下的三对角方程组 解设有分解 由公式 其中分别是系数矩阵的主对角线元素及其下边和上边的次对角线元素，故有 从而有 故，，， 故，，，

4设A是任一阶对称正定矩阵，证明是一种向量范数 证明（1）因A正定对称，故当时，，而当时， （2）对任何实数，有 （3）因A正定，故有分解，则 故对任意向量和，总有 综上可知，是一种向量范数。 5 设，，已知方程组的精确解为 （1）计算条件数； （2）若近似解，计算剩余； （3）利用事后误差估计式计算不等式右端，并与不等式左边比较，此结果说明了什么？解（1） （2） （3）由事后误差估计式，右端为 而左端

这表明当A为病态矩阵时，尽管剩余很小，误差估计仍然较大。因此，当A病态时，用大小作为检验解的准确度是不可靠的。 6矩阵第一行乘以一数成为，证明当时，有最小值 证明设，则 又 故 从而当时，即时，有最小值，且 7讨论用雅可比法和高斯-赛德尔法解方程组时的收敛性。如果收敛，比较哪一种方 法收敛较快，其中 解对雅可比方法，迭代矩阵 ， 故雅可比法收敛。 对高斯-赛德尔法，迭代矩阵

，故高斯-赛德尔法收敛。 因=故高斯-赛德尔法较雅可比法收敛快。 8设，求解方程组，求雅可比迭代法与高斯-赛德尔迭代法收敛的充要条件。 解雅可比法的迭代矩阵 ， 故雅可比法收敛的充要条件是。 高斯-赛德尔法的迭代矩阵 ，

昆明理工大学—数值分析各年考试题及答案

昆明理工大学数值分析考试题 （07） 一．填空（每空3分，共30分） 1． 设A 0.231x =是真值0.229T x =的近似值，则A x 有 位有效数字。 2． 若74()631f x x x x =+++，则017[2,2,...2]f = ，018[2,2,...2]f = 。 3． A=1031?? ? ?-?? ,则1 A = ； A ∞ = ； 2 A = 2()cond A = 。 4． 求方程()x f x =根的牛顿迭代格式是 。 5．设105%x =± ，则求函数()f x =的相对误差限为 。 6．A=2101202a a ?? ? ? ??? ,为使其可分解为T L L （L 为下三角阵，主对角线元素>0），a 的取值范 围应为 。 7．用最小二乘法拟合三点A(0,1),B(1,3),C(2,2)的直线是 。 （注意：以上填空题答案标明题号答在答题纸上，答在试卷上的不给予评分。） 二．推导与计算 （一）对下表构造f(x)的不超过3次的插值多项式，并建立插值误差公式。（12分） （二）已知()x x =Φ和()x 'Φ满足∣()x 'Φ-3∣<1。请利用()x Φ构造一个收敛的简单迭代函数()x ψ，使1(),0,1,......k k x x k +=ψ=收敛。（8分）

（三）利用复化梯形公式计算2 1 x I e dx -=?，使其误差限为60.510-?，应将区间[0，1] 等份。（8分） （四）设A= 1001005a b b a ?????????? ，detA ≠0，推导用a ，b 表示解方程组AX=f 的Seidel(G-S) 迭代法收敛的充分必要条件。（10分） （五）确定节点及系数，建立如下 GAUSS 型求积公式 1 11220 ()()dx A f x A f x ≈+? 。（10分） （六）对微分方程初值问题'00(,) ()y f x y y x y ?=?=? （１） 用数值积分法推导如下数值算法： 1111(4)3 n n n n n h y y f f f +-+-=+ ++，其中(,)i i i f f x y =，(1,,1)i n n n =-+。（8分） （２） 试构造形如 1011011(),n n n n n y a y a y h b f b f +--=+++ 的线形二步显格式差分格式，其中111(,),(,)n n n n n n f f x y f f x y ---==。试确定系数0101,,,a a b b ，使 差分格式的阶尽可能高，写出其局部截断误差主项，并指明方法是多少阶。（14分） （考试时间2小时30分钟）

数值分析试题及答案

一、单项选择题（每小题3分，共15分） 1. 3.142和3.141分别作为π的近似数具有（ ）和（ ）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．4和4 2. 已知求积公式 ()()2 1 121 1()(2)636f x dx f Af f ≈ ++? ，则A ＝（ ） A ． 16 B ．13 C ．12 D ．2 3 3. 通过点 ()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足（ ） A ． ()00l x ＝0， ()110l x = B ． ()00l x ＝0， ()111l x = C ．() 00l x ＝1，()111 l x = D ． () 00l x ＝1，()111 l x = 4. 设求方程 ()0 f x =的根的牛顿法收敛，则它具有（ ）敛速。 A ．超线性 B ．平方 C ．线性 D ．三次 5. 用列主元消元法解线性方程组 1231231 220223332 x x x x x x x x ++=?? ++=??--=? 作第一次消元后得到的第3个方程（ ）. A ． 232 x x -+= B ．232 1.5 3.5 x x -+= C ． 2323 x x -+= D ． 230.5 1.5 x x -=- 单项选择题答案 1.A 2.D 3.D 4.C 5.B 得 分 评卷人 二、填空题（每小题3分，共15分）

1. 设T X )4,3,2(-=, 则=1||||X ，2||||X = . 2. 一阶均差 ()01,f x x = 3. 已知3n =时，科茨系数()()() 33301213,88C C C ===，那么 () 33C = 4. 因为方程()420 x f x x =-+=在区间 []1,2上满足 ，所以()0f x =在区间 内有根。 5. 取步长0.1h =，用欧拉法解初值问题 ()211y y y x y ?'=+?? ?=? 的计算公式 . 填空题答案 1. 9和29 2. ()() 0101 f x f x x x -- 3. 1 8 4. ()()120 f f < 5. ()12 00.1 1.1,0,1,210.11k k y y k k y +???? ?=+? ?=+???? =??L 得 分 评卷人 三、计算题（每题15分，共60分） 1. 已知函数 21 1y x = +的一组数据： 求分 段线性插值函数，并计算 () 1.5f 的近似值. 计算题1.答案 1. 解 []0,1x ∈， ()1010.510.50110x x L x x --=?+?=---% []1,2x ∈，()210.50.20.30.81221x x L x x --=?+?=-+--%

2014-2015数值分析考试试题卷

太原科技大学硕士研究生 2014/2015学年第1学期《数值分析》课程试卷 一、填空题（每空4分，共32分） 1、设?????≤≤-++<≤+=2 1,1321 0,)(2 323x x bx x x x x x s 是以0,1,2为节点三次样条函数，则b=__-2___ 2、解线性方程组12312312388 92688 x x x x x x x x x -++=-?? -+=??-+-=? 的Jacobi 迭代格式（分量形式）为 ?? ???+--=++-=++=+++)(2)(1)1(3) (3)(1)1(2) (3)(2)1(1882/)96(88k k k k k k k k k x x x x x x x x x ，其相应的迭代矩阵为??????????-0812/102/9810。 3、方程03 =-a x 的牛顿法的迭代格式为__3 12 3k k k k x a x x x +-=-__________，其收敛的阶为 2 。 4、已知数x 的近似值0.937具有三位有效数字，则x 的相对误差限是310534.0-? 解：x 1≈0.937， 31102 1 )(-?≤ x ε 3 31111 10(x )2 (x )0.53410x 0.937 r εε--?=≤=? 5、用列主元高斯消去法解线性方程组 ??? ??=--=++=++2333220221 321321x x x x x x x x 作第1次消元后的第2，3个方程分别为? ? ?=+--=-5.35.125 .15.03232x x x x 6、设???? ??-=3211A ，则=∞)(A Cond __4____.

华南理工大学数值分析试题-14年下-C

华南理工大学研究生课程考试 《数值分析》试卷C （2015年1月9日） 1. 考前请将密封线内各项信息填写清楚； 所有答案请按要求填写在本试卷上； 课程代码：S0003004； 4. 考试形式：闭卷； 5. 考生类别：硕士研究生； 本试卷共八大题，满分100分，考试时间为150分钟。 一、（12分）解答下列问题： 1）设近似值0x >，x 的相对误差为δ，试证明ln x 的绝对误差近似为δ。 2）利用秦九韶算法求多项式 542()681p x x x x x =-+-+ 在3x =时的值（须写出计算形式），并统计乘法次数。 （12分）解答下列问题： 1）设()235f x x =+，求[]0,1,2f 和[]0,1,2,3f 。 2）利用插值方法推导出恒等式： 33220,0[]j j i i x j i x i j =≠=-=-∑∏ 。

（1）设{}∞ =0)(k k x q 是区间[]1,0上带权1=ρ而最高次项系数为1的正交多项式族，其中1)(0=x q ，求1()q x 和2()q x 。 （2）求形如2y a bx =+的经验公式，使它与下列数据拟合： 四、（14分）对积分()10I f x dx = ?，试 （1）构造一个以012113,,424 x x x ===为节点的插值型求积公式； （2）指出所构造公式的代数精度； （3）用所得数值求积公式计算积分1 203x dx ?的精确值； （4）指出所得公式与一般的Newton-Cotes 型公式在形式上的重要区别。

（1）设?? ????=4321A ,计算1A 、()Cond A ∞和()A ρ。 （2）用列主元Gauss 消去法解方程组： 12312315410030.112x x x ????????????=????????????-?????? 六、（13分）对2阶线性方程组 11112212112222 a x a x b a x a x b +=??+=? （11220a a ≠ ） （1）证明求解此方程组的Jacobi 迭代与Gauss-Seidel 迭代同时收敛或同时发散； （2）当同时收敛时，试比较它们的收敛速度。

北航2010-2011年研究生数值分析期末模拟试卷1-3

数值分析模拟试卷1 一、填空（共30分，每空3分） 1 设??? ? ??-=1511A ，则A 的谱半径=)(a ρ______，A 的条件数=________. 2 设 ,2,1,0,,53)(2==+=k kh x x x f k ，则],,[21++n n n x x x f ＝________, ],,[321+++n n n n x x x x f ，＝________. 3 设?????≤≤-++≤≤+=2 1,121 0,)(2 323x cx bx x x x x x S ，是以0，1，2为节点的三次样条函数，则b=________,c=________. 4 设∞=0)]([k k x q 是区间［0，1］上权函数为x x =)(ρ的最高项系数为1的正交多项式族，其中1)(0=x q ，则 ?=1 )(dx x xq k ________，=)(2 x q ________. 5 设???? ??????=11001a a a a A ，当∈a ________时，必有分解式，其中L 为下三角阵，当 其对角线元素)3,2,1(=i L ii 满足条件________时，这种分解是唯一的. 二、（14分）设4 9,1,41,)(2102 3 === =x x x x x f , （1）试求)(x f 在]4 9,41[上的三次Hermite 插值多项式)(x H 使满足 2,1,0),()(==i x f x H i i ，)()(11x f x H '='. （2）写出余项)()()(x H x f x R -=的表达式. 三、（14分）设有解方程0cos 2312=+-x x 的迭代公式为n n x x cos 3 2 41+ =+， （1） 证明R x ∈?0均有? ∞ →=x x n x lim （? x 为方程的根）； （2） 取40=x ，用此迭代法求方程根的近似值，误差不超过，列出各次迭代值； （3）此迭代的收敛阶是多少？证明你的结论. 四、(16分) 试确定常数A ，B ，C 和，使得数值积分公式 有尽可能高的代数精度. 试问所得的数值积分公式代数精度是多少？它是否为Gauss 型的？

数值分析试卷及其答案1

1. 已知325413.0,325413*2*1==X X 都有6位有效数字，求绝对误差限。（4分） 解： 由已知可知6 5.0102 1 ,0,6,10325413.0016*1=?= =-=?=ε绝对误差限n k k X 2分 620*2102 1 ,6,0,10325413.0-?= -=-=?=ε绝对误差限n k k X 2分 2. 已知?? ???=0 01 A 220- ?????440求21,,A A A ∞ （6分） 解： {}, 88,4,1max 1==A 1分 {}, 66,6,1max ==∞A 1分 () A A A T max 2λ= 1分 ?????=0 1 A A T 4 2 ???? ? -420?????0 01 2 20 - ???? ?440= ?????0 01 80 ???? ?3200 2分 {}32 32,8,1max )(max ==A A T λ

1分 24322==A 3. 设32)()(a x x f -= （6分） ① 写出f(x)=0解的迭代格式 ② 当a 为何值时，)(1k k x x ?=+ （0,1……）产生的序列{}k x 收敛于 2 解： ①迭代格式为： x a x x x a x a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(2 2 32 1 += +=---=-=+? 3 分 ②时迭代收敛即当222,112 10)2(',665)('2<<-<-=-= a a x a x ?? 3分 4. 给定线性方程组，其中：?? ?=13A ?? ?2 2，?? ? ???-=13b 用迭代公式 )()()()1(k k k Ax b x x -+=+α（0,1……）求解，问取什么实数α ，可使 迭代收敛 （8分） 解： 所给迭代公式的迭代矩阵为?? ? --???--=-=ααααα21231A I B 2分

数值分析期末试卷

数值分析2006 — 2007学年第学期考试 课程名称：计算方法 A 卷 考试方式：开卷[] 闭卷[V ] 半开卷[] IV 类 充要条件是a 满足 二、（18分）已知函数表如下 1?设 f(0) = 0, f (1) =16 , f( 2) =46，则 f [0,1]= ,f[0,1,2]二 2 ?设 AJ <2 -3 -1 ，则X ，A := A 1 1 j — 3 ?计算积分 xdx ，取4位有效数字。用梯形公式求得的近似值为 "0.5 （辛普森）公式求得的近似值为 ，用 Spsn 4?设f （x ）二xe x -3，求方程f （x ） =0近似根的牛顿迭代公式是 ，它的收 敛阶是 5 ?要使求积公式 1 1 [f (x)dx 拓一(0) + A , f (x 1)具有2次代数精度，则 捲= _________________ , 0 4 6 ?求解线性方程组 x 1 ax 2 = 4 , 12_3 （其中a 为实数）的高斯一赛德尔迭代格式收敛的 10 11 12 13 In x 2.3026 2.3979 2.4849 2.5649

三、(20分)构造如下插值型求积公式，确定其中的待定系数，使其代数精度尽可能高, 并指出所得公式的代数精度。 2 f (x)dx ： A o f (0) A f (1) A2f(2) o

X 2 4 6 8 y 2 11 28 40 五、（14分）为求方程X ’ -X 2 -1 =0在X o =1.5附近的一个根，将方程改写为下列等价 形式，并建立相应的迭代公式: 试问上述两种迭代公式在 x 0 =1.5附近都收敛吗？为什么？说明理由。 （1）X =1 ?丄，迭代公式 X 1 X k 1 = 1 - X k （2） X 2二1 ，迭代公式 X —1 2 （X k ）； X k 1

数值分析试题及答案

数值分析试题及答案 一、单项选择题（每小题3分，共15分） 1. 3.142和3.141分别作为的近似数具有（）和（）位有效数字. A．4和3 B．3和2 C．3和4 D．4和4 2. 已知求积公式，则＝（） A． B．C．D． 3. 通过点的拉格朗日插值基函数满足（） A．＝0，B．＝0， C．＝1，D．＝1， 4. 设求方程的根的牛顿法收敛，则它具有（）敛速。 A．超线性B．平方C．线性D．三次 5. 用列主元消元法解线性方程组作第一次消元后得到的第3个方程（）. A．B． C．D． 单项选择题答案 1.A 2.D 3.D 4.C 5.B 得分评卷 人 二、填空题（每小题3分，共15分） 1. 设, 则， . 2. 一阶均差 3. 已知时，科茨系数，那么 4. 因为方程在区间上满足，所以在区间内有根。 5. 取步长，用欧拉法解初值问题的计算公式.填空题答案

1. 9和 2. 3. 4. 5. 得分评卷 人 三、计算题（每题15分，共60分） 1. 已知函数的一组数据：求分段线性插值函数，并计算的近似值. 计算题1.答案 1. 解， ， 所以分段线性插值函数为 2. 已知线性方程组 （1）写出雅可比迭代公式、高斯－塞德尔迭代公式； （2）对于初始值，应用雅可比迭代公式、高斯－塞德尔迭代公式分别计算（保留小数点后五位数字）. 计算题2.答案 1.解原方程组同解变形为 雅可比迭代公式为 高斯－塞德尔迭代法公式 用雅可比迭代公式得 用高斯－塞德尔迭代公式得 3. 用牛顿法求方程在之间的近似根 （1）请指出为什么初值应取2？ （2）请用牛顿法求出近似根，精确到0.0001. 计算题3.答案

数值计算方法试题集及答案要点

《数值计算方法》复习试题 一、填空题： 1、 ?? ??? ?????----=410141014A ，则A 的LU 分解为 A ? ???????? ???=????????? ?? ?。 答案： ?? ????????--??????????--=1556141501 4115401411A 2、已知3.1)3(,2.1)2(, 0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求 得?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。 答案：2.367，0.25 3、1)3(,2)2(, 1)1(==-=f f f ，则过这三点的二次插值多项式中2x 的系数 为 ，拉格朗日插值多项式为 。 答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 4、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 5、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )； 答案 )(1)(1n n n n n x f x f x x x '--- =+ 6、对 1)(3++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )； 7、计算方法主要研究( 截断 )误差和( 舍入 )误差； 8、用二分法求非线性方程f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为( 1 2+-n a b )； 9、求解一阶常微分方程初值问题y '= f (x ,y )，y (x 0)=y 0的改进的欧拉公

数值分析历年考题

数值分析A 试题 2007.1 第一部分：填空题10?5 1.设3112A ?? = ??? ,则A ∞=___________ 2()cond A =___________ 2.将4111A ??= ??? 分解成T A LL =，则对角元为正的下三角阵L =___________ ,请用线性最小二乘拟合方法确定拟合函数()bx f x ae =中的参数：a = ___________ b =___________ 4.方程13 cos 2044x x π--=在[0,1]上有 个根，若初值取00.95x =，迭代方法 113 cos 244 k k x x π+=-的收敛阶是 5.解方程2 210x x -+=的Newton 迭代方法为___________，其收敛阶为___________ 6.设()s x = 323 2 323,[0,1]31,[1,2] ax x x x x x bx x +-+∈--+∈为三次样条函数，则a = ___________ b =___________ 7.要想求积公式： 1 121 ()(()f x dx A f f x -≈+? 的代数精度尽可能高，参数1A = ___________ 2x =___________此时其代数精度为：___________ 8.用线性多步法2121(0.50.5)n n n n n y y h f f f ++++-=-+来求解初值问题 00'(,),(),y f x y y x y ==其中(,)n n n f f x y =，该方法的局部截断误差为___________，设 ,0,f y μμ=?其绝对稳定性空间是___________ 9.用线性多步法 2121()n n n n n y ay by h f f ++++-+=-来求解初值问题 00'(,),(),y f x y y x y ==其中(,)n n n f f x y =，希望该方法的阶尽可能高，那么a = ___________ b =___________，此时该方法是几阶的：___________

数值分析期末试题

数值分析期末试题 一、填空题（20102=?分） （1)设??? ? ? ??? ??---=28 3 012 251A ，则=∞ A ______13_______。 （2)对于方程组?? ?=-=-3 4101522121x x x x ，Jacobi 迭代法的迭代矩阵是=J B ?? ? ? ??05.25.20。 （3)3*x 的相对误差约是*x 的相对误差的 3 1倍。 （4）求方程)(x f x =根的牛顿迭代公式是) ('1)(1n n n n n x f x f x x x +-- =+。 （5）设1)(3 -+=x x x f ，则差商=]3,2,1,0[f 1 。 （6）设n n ?矩阵G 的特征值是n λλλ,,,21 ，则矩阵G 的谱半径=)(G ρi n i λ≤≤1max 。 （7）已知?? ? ? ??=1021 A ，则条件数=∞ )(A Cond 9 （8）为了提高数值计算精度，当正数x 充分大时,应将)1ln(2 -- x x 改写为 )1ln(2 ++ -x x 。 （9）n 个求积节点的插值型求积公式的代数精确度至少为1-n 次。 （10）拟合三点))(,(11x f x ，))(,(22x f x ，))(,(33x f x 的水平直线是)(3 1 3 1 ∑== i i x f y 。 二、（10分）证明：方程组? ?? ??=-+=++=+-1 211 2321321321x x x x x x x x x 使用Jacobi 迭代法求解不收敛性。 证明：Jacobi 迭代法的迭代矩阵为 ???? ? ?????---=05 .05 .01015.05.00J B J B 的特征多项式为

数值分析试题1

数值分析试卷1 一、填空题（每空2分，共30分） 1. 近似数231.0=*x 关于真值229.0=x 有____________位有效数字； 2. 设)(x f 可微，求方程)(x f x =根的牛顿迭代格式是_______________________________________________； 3. 对1)(3++=x x x f ，差商=]3,2,1,0[f _________________； =]4,3,2,1,0[f ________； 4. 已知??? ? ??-='-=1223,)3,2(A x ，则=∞||||Ax ________________，=)(1A Cond ______________________ ； 5. 求解线性方程组?????=+=+045 11532121x x x x 的高斯—赛德尔迭代格式为_______________________________________；该迭代格式迭代矩阵的谱半径=)(G ρ_______________； 二、（12分）(1)设LU A =，其中L 为下三角阵，U 为单位上三角阵。已知 ?????? ? ??------=2100121001210012A ，求L ，U 。 (2)设A 为66?矩阵，将A 进行三角分解：LU A =，L 为单位下三角阵，U 为上三角阵，试写出L 中的元素65l 和U 中的元素56u 的计算公式。 三、给定数据表如下 x 0.20.40.60.81 1.2f(x)212523202124 （1） 用三次插值多项式计算f ( 0.7 ) 的近似值； （2） 用二次插值多项式计算f ( 0.95 ) 的近似值： （3） 用分段二次插值计算 f ( x ) )2.12.0(≤≤x 的近似值能保证有几位有

数值分析试题A卷10.1

中国石油大学（北京）2009--2010学年第一学期 研究生期末考试试题A (闭卷考试) 课程名称：数值分析 注：计算题取小数点后四位 一、填空题（共30分，每空3分） 1、 已知x =是由准确数a 经四舍五入得到的近似值，则x 的绝对误差 界为_______________。 2、数值微分公式()() '()i i i f x h f x f x h +-≈ 的截断误差为 。 3、已知向量T x =，求Householder 变换阵H ，使(2,0)T Hx =-。 H = 。 4、利用三点高斯求积公式 1 1 ()0.5556(0.7746)0.8889(0)0.5556(0.7746)f x dx f f f -≈-++? 导出求积分 4 0()f x dx ?的三点高斯求积公式 。 5、4 2 ()523,[0.1,0.2,0.3,0.4,0.5]_____.f x x x f =+-= 若则 6、以n + 1个互异节点x k ( k =0,1,…,n ),（n >1）为插值节点的 Lagrange 插值基函数为l k (x)( k =0,1,…,n )，则 (0)(1)__________.n k k k l x =+=∑ 7、已知3()P x 是用极小化插值法得到的cos x 在[0,4]上的三次插值多项式，则3()P x 的 截断误差上界为3()cos ()R x x P x =-≤_________.

8、已知向量(3,2,5)T x =-，求Gauss 变换阵L ，使(3,0,0)T Lx =。L =_________. 9、设3 2 ()(7)f x x =-, 给出求方程()0f x =根的二阶收敛的迭代格式_________。 10、下面M 文件是用来求解什么数学问题的________________________. function [x,k]=dd （x0） for k=1:1000 x=cos (x0); if abs(x-x0)<, break end x0=x; end 二、（15分）已知矛盾方程组Ax=b ，其中11120,1211A b ???? ????==???????????? ， （1）用施密特正交化方法求矩阵A 的正交分解，即A=QR 。 （2）用此正交分解求矛盾方程组Ax=b 的最小二乘解。 三、（10分）已知求解线性方程组Ax=b 的分量迭代格式 1 (1) (1) ()1 +1 /, 121,,i n k k k i i ij j ij j ii j j i x b a x a x a i n n -++===-- =-∑∑（），， （1）试导出其矩阵迭代格式及迭代矩阵； （2）若11a A a ?? = ??? ，推导上述迭代格式收敛的充分必要条件。 四、（15分）（1）证明对任何初值0x R ∈，由迭代公式11 1sin ,0,1,2, (2) k k x x k +=+ = 所产生的序列{}0k k x ∞ =都收敛于方程1 1sin 2 x x =+ 的根。 （2）迭代公式11 21sin ,0,1,2, (2) k k k x x x k +=-- =是否收敛。 五、（15分）用最小二乘法确定一条经过原点(0,0)的二次曲线，使之拟合下列数据