2012研究生数值分析课期末考试复习题及答案

一、填空

1. 设

2.3149541...x *

=，取5位有效数字，则所得的近似值x= 2.3150 .

2.设一阶差商

()()()21122114

,321f x f x f x x x x --=

=

=---，

()()()322332

615

,422f x f x f x x x x --=

=

=--

则二阶差商

()123,,______

f x x x =11/6

3. 设(2,3,1)T

X =--, 则2||||X =

14 ，

=∞||||X 3 。p49

4. 4．求方程 2

1.250x x --= 的近似根，用迭代公式 1.25x x =+，取初始

值

01

x =， 那么

1______x =。

1.5

5．解初始值问题 00

'(,)()y f x y y x y =??

=?近似解的梯形公式是

1______k y +≈。

()()[]11,,2

++++k k k k k y x f y x f h y

6、

1151A ??= ?

-??，则A 的谱半径 ＝ 6 。

7、设

2()35, , 0,1,2,... ,

k f x x x kh k =+== ，则

[]12,,n n n f x x x ++=

—————

—————3 和

[]123,,,n n n n f x x x x +++=

_______________0_____ 。

8、 若线性代数方程组AX=b 的系数矩阵A 为严格对角占优阵，则雅可比迭代和高斯-塞德尔迭代都

收敛

。

9、解常微分方程初值问题的欧拉（Euler）方法的局部截断误差为_______O(h）___。

10、为了使计算

23 123 10

1(1)(1)

y

x x x

=++-

---的乘除法运算次数尽量的少,应将

表达式改写成____________??

?

?

?

?

?

?

?

?

?

-

-

-

+

-

+

=

1

3

2

1

1

1

1

1

10

x

x

x

y_____________。

二、计算题

1、已知的满足，试问如何利用构造一个收敛的简单迭代函数，使0，1…收敛？

由()

x x

?

=，可得3()3

x x x x

?

-=-，

1

(()3)()

2

x x x x

?ψ

=--=

1

()(()3)

2

x x

ψψ

=--

’’

因，故

11

()1

22

x x

ψ?

=<<

’’

（）-3

[]

1

1

()()3 , k=0,1,....

2

k k k k

x x x x

ψ?

+

==--

故收敛。

2、试确定常数A，B，C和 a，使得数值积分公式

有尽可能高的代数精度。试问所得的数值积分公式代数精度是多少？它是否为Gauss型的？

101612

,,

995

A C

B a

====±

，该数值

求积公式具有5次代数精确度，它是Gauss型的

3、利用矩阵的LU分解法解方程组

123

123

123

2314

25218

3520

x x x

x x x

x x x

++=

?

?

++=

?

?++=

?

?

?

?

?

?

?

?

?

?

?

3

2

1

4、写出求解下列初始值问题?

?

?

=

≤

≤

-

=

2

)1(

)2

1(

,

3

8

'

y

x

y

y

的欧拉迭代式，欧拉预-校迭代式及四阶龙格-库塔法迭代式。

5.设

2

1

2

S gt

=

，假定 g是准确的，而对的测量有±秒的误差，

证明当增加时的绝对误差增加，而相对误差却减少。解：

2

**2

22

11

()0.1

22

()0.10.2

()

11

22

,(),().

r

r

e S S S gt gt gt

e S gt

e S

t

gt gt

t e S e S

=-=-=

===

∴↑↑↓

6.在x

-≤≤上给出x

f x e

=的等距节点函数表，若用二次插值求e的近似值，要使截断误差不超过10-，问使用函数表的步长应取多少？解：

()4

0000

(),(),[4,4],,,, 1.

x k x

f x e f x e e x x h x x h x x th t

==≤∈--+=+≤

考察点及

(3)

2000

4

43

4

3

()

()[(()]()[()]

3!

(1)(1)

(1)(1)

3!3!

2

.(4,4).

633

f

R x x x h x x x x h

t t t

e

t h th t h e h

e

h

ξ

ξ

=----+

-+

≤+??-=

≤?∈-

则

7.已知单调连续函数y f x

=的如下数据

用插值法计算约为多少时 f x =（小数点后至少保留4位）0.2008 解：作辅助函数g x f x =-则问题转化为为多少时，g x =此时可作新

的关于

i

g x 的函数表。

由f x 单调连续知g x 也单调连续，因此可对g x 的数值进行反插。的牛顿型插值多项式为

()0.110.097345( 2.23)0.451565( 2.23)( 1.10) 0.255894( 2.23)( 1.10)(0.17)

x g y y y y y y y -==-+++++-++-

故

(0) 1.321497.x g -

== 8. 设函数f x 在区间[0,3]上具有四阶连续导数,试用埃尔米特插值法求一个次数不高于3的多项式

()

P x ， 使其满足

(0)0

P =，(1)1

P =，

'(1)3P =,(2)1

P =

并写出误差估计式。

解：由所给条件可用埃尔米特插值法确定多项式()P x ， 32357()722p x x x x =-+-

由题意可设3()()()()(1)(2)

R x f x p x k x x x x =-=--为确定待定函数

k x ，作辅助函数：3()()()()(1)(2)g t f t p t k t t t t =----

11、确定下列求积公式中的待定参数，使其代数精度尽量高，并指明所构造出的求积公式所具有的代数进度。 1）()()(0)();

h

h

f x dx A f h A f A f h --≈-++?

2）1012()()(0)();

h

h

f x dx A f h A f A f h --≈-++?

3）

1

121

(1)2()3()

();

3f f x f x f x dx --++≈

?

4）

2''0

[(0)()]

()[(0)()].

2h

h f f h f x dx ah f f h +≈

+-?

解：（1）三个参数，代入

110121102311133

3444

1324()1,,,()032

1()33()()()()33334()()(0)().333h h h h h h A h

A A A h f x x x h A A A h

h A A h A h h h h h x dx h h x dx h h h h h f x dx f h f f h -------?=???++=???

=?--=?=??????+==???

=-+≠-+∴≈-++???Q 具有三次代数精度

（2）三个参数，代入

1101211022

311123333224544

452228344()1,,,03168()33848()0()03336484816()0()53333

84()()3h h h h

h h A h

A A A h f x x x hA hA A h

h A h A h A h h h

x dx h h h h h x dx h h h h h h h f x dx f h -------?

=???++=??-?

=?-+=?=??????-+==???

=--?+==≠--?+=∴≈-+???Q 8(0)().33h h f f h +具有三次代数精度

11212121122

12221

(3)()1,()[(1)2()3()].

3

,(),2310.689900.289902310.126600.52660f x f x dx f f x f x f x x x x x x x x x x x -==-++=+===-????????+==-=????当时有两个参数令精确成立

或

1333

1211

11

1

1[123]

3

()[(1)2(0.68990)3(0.12660)]/3

()[(1)2(0.28990)3(0.52660)]/3

2.

x dx x x f x dx f f f f x dx f f f ---≠-++≈-++-≈-+-+???而 故与均具有次代数精度

200

22220

2333

20244

301

(4)()1,,1[11]0,

[0](11).2

2

(),1[0][202].212

(), [0][03]

212(),()[0][04].

212

.h

h

h

h h h

f x x dx xdx h ah f x x h x dx h ah h a h h f x x x dx h h h h f x x f x dx h h ==++=

++-==

++?-?===++-=≠++-??

???时有故令时求积公式精确成立

当时时故只有三次代数精度

12. 对线性代数方程组13423412321564833

x x x x x x x x x x x x -+=??-+=??

+-=??-+-=? 设法导出使雅可比（Jacobi ）迭代法和高斯－赛德尔（G-S ）迭代法均收敛的迭代格式，要求分别写出迭代格式，并说明收敛的理由。

13. 设线性方程组为 1122211222

2

a x a x

b a a a x a x b +=?≠?

+=?

（1） 证明用雅可比迭代法和高斯－赛德尔迭代法解此方程组要么同时收敛，

要么同时发散。 （2） 当同时收敛时，试比较其收敛速度。

14.

写出用四阶经典的龙格—库塔方法求解下列初值问题的计算公式

''3,01;,01;

(1)1)2)(0)1;

(0) 1.y y x y x y x x y y ?=<

?=??=?

15. 证明对任意参数，下列龙格－库塔公式是二阶的

1

231

2131();2

(,);

(,);((1),(1)).n n n n n n n n h y y K K K f x y K f x th y thK K f x t h y t hK +?

=++???=??=++?

=+-+-?

?

16. 证明0sin 1=--x x 在]1,0[内仅有一个根，若用二分法求误差不大于

4

?的根，求需要迭代的次数。

5.0-

10

答案：https://www.wendangku.net/doc/b413980225.html,/upload/shuzhifenxi/jxxt.htm

数值分析试题及答案汇总

数值分析试题 一、 填空题（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 )

数值分析试卷及答案

二 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设，求解方程组，求雅可比迭代法与高斯-赛德尔迭代法收敛的充要条件。 解雅可比法的迭代矩阵 ， 故雅可比法收敛的充要条件是。 高斯-赛德尔法的迭代矩阵 ，

2014级硕士研究生数值分析上机实习报告

2014级硕士研究生数值分析上机实习（第一次） 姓名：学号：学院： 实习题目：分别用二分法和Newton迭代法求方程x3■ 2x210x-20=0的根.实习目的：掌握两种解法，体会两种解法的收敛速度. 实习要求：用C程序语言编程上机进行计算，精确到8位有效数字. 报告内容： 1.确定实根的个数以及所在区间 2.将最后两次计算结果填入下表（保留8位数字）： 3.实习过程中遇到哪些问题？如何解决？有何心得体会?

4.两种解法的计算程序（此页写不下时可以加页）：

2014级硕士研究生数值分析上机实习（第二次）姓名：学号：学院： 实习题目：计算8阶三对角矩阵A=tridiag（0.235, 1.274, 0.235）的行列式.实习目的：掌握计算行列式的方法. 实习要求：首先选择一种算法，然后用C程序语言编程上机进行计算.报告内容： 1.简单描述所采用的算法: 2?计算结果： A 3.实习过程中遇到哪些问题？如何解决？有何心得体会?

4.写出C语言计算程序（此页写不下时可以加页）：

2014级硕士研究生数值分析上机实习（第三次） 姓名：学号：学院： 分别用Jacobi迭代法和Gauss-Seidel迭代法求解线性方程组实习题目： 2lx + 9.8y+ 3.4z= 6.7 <2.7x + 1.8y+ 7.2z= 2.4 8.6x + 1.5y + 3.4z = 1.9 实习目的：感受两种迭代法的收敛速度. 首先构造收敛的Jacobi迭代法和Gauss-Seidel迭代法，然后用实习要求： C程序语言编程上机进行求解，初始值均取为0,精确到4位小 数. 报告内容： 1.写出收敛的Jacobi迭代法和Gauss-Seidel迭代法:

数值分析试卷及答案

二 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对称正定，，故也对称正定。 设迭代矩阵的特征值为，为相应的特征向量，则与做内积，有 因正定，故，从而，格式收敛。

数值分析试卷及其答案

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、 已知函数表 x -1 1 2 ()f x -3 0 4 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解： （1）k x -1 1 2 k y -3 0 4 插值基函数分别为 ()()()()()()()()()() 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 ---===-----= = =----=== --- k x ()k f x 一阶 二阶 -1 -3 1 0 3/ 2 2 4 4 5/6 故所求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 ρ=，这样，有

数值分析试卷及其答案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 迭代格式为

研究生数值分析试卷

2005~2006学年第一学期硕士研究生期末考试试题（A 卷） 科目名称：数值分析 学生所在院： 学号： 姓名： 注意：所有的答题内容必须答在答题纸上，凡答在试题或草稿纸上的一律无效。 一、（15分）设求方程 0cos 2312=+-x x 根的迭代法 k k x x cos 3 2 41+=+ (1) 证明对R x ∈?0,均有*lim x x k k =∞ →,其中*x 为方程的根. (2) 此迭代法收敛阶是多少? 证明你的结论. 二、（12分）讨论分别用Jacobi 迭代法和Gauss-Seidel 迭代法求解下列方程组的收敛性。 ??? ??=++-=++=-+. 022,1, 122321 321321x x x x x x x x x 三、（8分）若矩阵??? ? ? ??=a a a a A 000002，说明对任意实数0≠a ，方程组b AX =都是非病态的。（范数用∞?） 四、（ 求)(x f 的Hermite 插值多项式)(3x H ，并给出截断误差)()()(3x H x f x R -=。 五、（10分）在某个低温过程中，函数 y 依赖于温度x （℃）的试验数据

为 已知经验公式的形式为 2bx ax y += ，试用最小二乘法求出 a ，b 。 六、（12分）确定常数 a ，b 的值，使积分 [ ] dx x b ax b a I 2 1 1 2 ),(?--+= 取得最小值。 七、（14分）已知Legendre(勒让德)正交多项式)(x L n 有递推关系式： ?? ? ? ???=+-++===-+),2,1()(1)(112)()(, 1)(1110 n x L n n x xL n n x L x x L x L n n n 试确定两点的高斯—勒让德（G —L ）求积公式 ? -+≈1 1 2211)()()(x f A x f A dx x f 的求积系数和节点，并用此公式近似计算积分 ?=2 11 dx e I x 八、（14分）对于下面求解常微分方程初值问题 ?????==0 0)() ,(y x y y x f dx dy 的单步法： ??? ? ??? ++==++=+) ,() ,()2 121(1 21211 hk y h x f k y x f k k k h y y n n n n n n

数值分析试卷及其答案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分

数值分析整理版试题及答案

例1、 已知函数表 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解： （1） 故所求二次拉格朗日插值多项式为 （2）一阶均差、二阶均差分别为 例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 ρ=，这样，有 所以，法方程为

011231261192 34a a ??????????=?????????? ?????????? ，经过消元得012311 62110123a a ??? ???????=???????????????????? 再回代解该方程，得到14a =，011 6 a = 故，所求最佳平方逼近多项式为* 111()46S 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 ?=，这样，有 所以，法方程为 解法方程，得到00.8732a =，1 1.6902a =， 故，所求最佳平方逼近多项式为 例4、 用4n = 的复合梯形和复合辛普森公式计算积分1 ? 。 解： （1）用4n =的复合梯形公式 由于2h =，( )f x =()121,2,3k x k k =+=，所以，有 （2）用4n =的复合辛普森公式 由于2h =，( )f x =()121,2,3k x k k =+=，()12 220,1,2,3k x k k + =+=，所以，有 例5、 用列主元消去法求解下列线性方程组的解。 解：先消元 再回代，得到33x =，22x =，11x = 所以，线性方程组的解为11x =，22x =，33x = 例6、 用直接三角分解法求下列线性方程组的解。 解： 设 则由A LU =的对应元素相等，有 1114u = ，1215u =，1316u =， 2111211433l u l =?=，3111311 22 l u l =?=， 2112222211460l u u u +=?=-，2113232311 545l u u u +=?=-，

2008级研究生数值分析试题

太原科技大学 2008级硕士研究生08/09学年第一学期 《数值分析》考试试卷 说明：1、Legendre 正交多项式)(x L n 有三项递推关系式： ?? ?? ???=+-++===-+ ,2,1)(1)(112)()(,1)(1110n x L n n x xL n n x L x x L x L n n n 2、Chebyshev 多项式)(x T n 有三项递推关系式： ?? ? ??=-===-+ ,2,1)()(2)()(,1)(1110n x T x xT x T x x T x T n n n 一、填空题：（每题4分，共20分） 1、设??? ? ??-=1511A ，则=∞)(A Cond 2、为提高数值计算精度，当x 充分小时，应将 x x sin cos 1-改写为 3、设)5()(2 -+=x a x x ?，要使)(1k k x x ?=+局部收敛到5* = x ，则a 的取值范围为 4、近似数235.0* =x 关于真值229.0=x 有 位有效数字。 5、设,1)(3 -+=x x x f 则差商=]3,2,1,0[f 二、（本题满分10分）用数值积分的方法建立求解初值问题b x a y a y y x f y a ≤≤==',)(),,(的Simpson 公式： )4(3 1111-+-++++=n n n n n f f f h y y 其中1,,1),,(+-==n n n i y x f f i i i ，11-+-=-=n n n n x x x x h . 三、（本题满分15分）设要用Gauss-Seidel 迭代法求解下列线性方程组

数值分析试题及答案

一、单项选择题（每小题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 --=?+?=-+--%

数值分析习题集及答案Word版

数值分析习题集 （适合课程《数值方法A 》和《数值方法B 》） 长沙理工大学 第一章 绪 论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2％,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指 出它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式(3.3)求下列各近似值的误差限: ********12412324(),(),()/,i x x x ii x x x iii x x ++其中**** 1234 ,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1％,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 1n n Y Y -=…) 计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字27.982). 8. 当N 充分大时,怎样求2 1 1N dx x +∞+?? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 212S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对 误差增加,而相对误差却减小. 11. 序列 {}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字), 计算到 10y 时误差有多大?这个计算过程稳定吗? 12. 计算6 1)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?

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

《数值计算方法》复习试题 一、填空题： 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的改进的欧拉公

数值分析试题及答案

数值分析试题及答案 一、单项选择题（每小题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.答案

2012数值分析试卷答案

昆明理工大学2012级硕士研究生试卷 科目： 数值分析 考试时间： 出题教师： 集体 考生姓名： 专业： 学号： 考试要求：考试时间150分钟；填空题答案依顺序依次写在答题纸上，填在试卷卷面上的不予计分；可带计算器。 一、 填空题（每空2分，共40分） 1．设*0.231x =是真值0.228x =的近似值，则*x 有 位有效数字，*x 的相对误差限 为 。 2．设 133)(47+++=x x x x f ，则=]2,,2,2[710 f ，=]2,,2,2[810 f 。 3. 过点)0,2(),0,1(-和)3,1(的二次拉格朗日插值函数为 )(2x L = ， 并计 算=)0(2L 。 4．设 32()3245f x x x x =+-+在[]1,1-上的最佳二次逼近多项式为 ， 最佳二次平方逼近多项式为 。 5．高斯求积公式 )()()(1101 0x f A x f A dx x f x +≈? 的系数0A = ， 1A = ，节点0x = ， 1x = 。 6．方程组 b Ax =，,U L D A --=建立迭代公式f Bx x k k +=+)()1(，写出雅可比迭代法和 高斯-赛德尔迭代法的迭代矩阵， =Jacobi B ，=-Seidel Gauss B 。 7．0 0100A ??? =? ???，其条件数2()Cond A = 。 8．设?? ? ???=2113A ，计算矩阵A 的范数，1||||A = ， 2||||A = 。

9．求方程 ()x f x =根的牛顿迭代格式是 。 10．对矩阵??? ? ? ??=513252321A 作LU 分解，其L=________________， U= __________________。 二、计算题（每题10分，共50分） 1. 求一个次数不高于4次的多项式P (x ), 使它满足:1)1(,0)0(,0)0('===p p p ,1)1(,'=p ,1)2(=p 并写出其余项表达式（要求有推导过程）。 2. 若用复合梯形公式计算积分 dx e x ? 1 ，问区间[0, 1]应分成多少等分才能使截断误差不超过 5102 1 -?? 若改用复合辛普森公式，要达到同样的精度区间[0, 1]应该分成多少等份? 由下表数据，用复合辛普森公式计算该积分的近似值。 3. 线性方程组b Ax =，其中???? ??????=18.04.08.014.04.04.01A ，T b ]3,2,1[=，（1）建立雅可比迭代法和 高斯-赛德尔迭代法的分量形式。（2）问雅可比迭代法和高斯-赛德尔迭代法都收敛吗 ? 4. 已知如下实验数据4,,1,0),,( =i y x i i , 用最小二乘法求形如x a a y 10+=的经验公式，并 计算最小二乘法的误差。 5. 用改进的欧拉公式(预估-校正方法)，解初值问题0)0(,10022=+=y y x dx ，取步长,1.0=h 计算到2.0=x （保留到小数点后四位） 。 三、证明题（共10分） 1． 如果 A 是对称正定矩阵，则A 可唯一地写成T LL A =，其中L 是具有正对角元的下三角 阵。

研究生《数值分析》教学大纲

研究生《数值分析》教学大纲 课程名称：数值分析 课程编号：S061005 课程学时：64 学时 课程学分： 4 适用专业：工科硕士生 课程性质：学位课 先修课程：高等数学，线性代数，计算方法，Matlab语言及程序设计 一、课程目的与要求 “数值分析”课是理工科各专业硕士研究生的学位课程。主要介绍用计算机解决数学问题的数值计算方法及其理论。内容新颖，起点较高，并加强了数值试验和程序设计环节。通过本课程的学习，使学生熟练掌握各种常用的数值算法的构造原理和过程分析，提高算法设计和理论分析能力，并且能够根据数学模型，提出相应的数值计算方法编制程序在计算机上算出结果。力求使学生掌握应用数值计算方法解决实际问题的常用技巧。 二、教学内容、重点和难点及学时安排： 第一章? 数值计算与误差分析( 4学时) 介绍数值分析的研究对象与特点，算法分析与误差分析的主要内容。 第一节数值问题与数值方法 第二节数值计算的误差分析 第三节数学软件工具----MATLAB 语言简介 重点：误差分析 第二章? 矩阵分析基础( 10学时) 建立线性空间、赋范线性空间、内积空间的概念，为学习以后各章打好基础。矩阵分解是解决数值代数问题的常用方法，掌握矩阵的三角分解、正交分解、奇异值分解，并能够编写算法程序。 第一节? 矩阵代数基础

第二节? 线性空间 第三节? 赋范线性空间 第四节? 内积空间和内积空间中的正交系 第五节矩阵的三角分解 第六节矩阵的正交分解 第七节矩阵的奇异值分解 难点：内积空间中的正交系。矩阵的正交分解。 重点：范数，施密特(Schmidt) 正交化过程，正交多项式，矩阵的三角分解, 矩阵的正交分解。 第三章? 线性代数方程组的数值方法( 12学时) 了解研究求解线性代数方程组的数值方法分类及直接法的应用范围。高斯消元法是解线性代数方程组的最常用的直接法，也是其它类型直接法的基础。在此方法基础上加以改进，可得选主元的高斯消元法、按比例增减的高斯消元法，其数值稳定性更高。掌握用列主元高斯消元法解线性方程组及计算矩阵的行列式及逆，并且能编写算法程序。掌握矩阵的直接三角分解法：列主元LU 分解，Cholesky分解。了解三对角方程组的追赶法的分解形式及数值稳定性的充分条件。掌握矩阵条件数的定义，并能利用条件数判别方程组是否病态以及对方程组的直接方法的误差进行估计。 迭代解法是求解大型稀疏方程组的常用解法。熟练掌握雅可比迭代法、高斯- 塞德尔迭代法及SOR 方法的计算分量形式、矩阵形式，并能在计算机上编出三种方法的程序用于解决实际问题。了解极小化方法：最速下降法、共轭斜量法。迭代法的收敛性分析是研究解线性代数方程组的迭代法时必须考虑的问题。对于上述常用的迭代法，须掌握其收敛的条件。而对一般的迭代法，掌握其收敛性分析的基本方法和主要结果有助于进一步探究新的迭代法。 第一节求解线性代数方程组的基本定理 第二节高斯消元法及其计算机实现 第三节矩阵分解法求解线性代数方程组 第三节? 误差分析和解的精度改进 第四节? 大型稀疏方程组的迭代法 第五节? 极小化方法 难点：列主元高斯消元法，直接矩阵三角分解。迭代法的收敛性，雅可比迭代法，高斯-塞德尔迭代法，SOR 迭代法。

数值分析最佳习题(含答案)

第一章 绪论 姓名 学号 班级 习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。 1 若误差限为5105.0-?,那么近似数有几位有效数字（有效数字的计算） 解：2*103400.0-?=x ，325*102 1102 1---?=?≤-x x 故具有3位有效数字。 2 14159.3=π具有4位有效数字的近似值是多少（有效数字的计算） 解：10314159.0?= π，欲使其近似值*π具有4位有效数字，必需 41*102 1 -?≤-ππ，3*3102 1102 1--?+≤≤?-πππ，即14209.314109.3*≤≤π 3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +， b a ?有几位有效数字（有效数字的计算） 解：3*1021 -?≤-a a ，2*102 1-?≤-b b ，而1811.2=+b a ，1766.1=?b a 2123****102 1 10211021)()(---?≤?+?≤ -+-≤+-+b b a a b a b a 故b a +至少具有2位有效数字。 2123*****102 1 0065.01022031.1102978.0)()(---?≤=?+?≤ -+-≤-b b a a a b b a ab

故b a ?至少具有2位有效数字。 4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差（误差的计算） 解：已知δ=-* *x x x ，则误差为 δ=-= -* **ln ln x x x x x 则相对误差为 * * ** * * ln ln 1ln ln ln x x x x x x x x δ = -= - 5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5*=， 已知cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v 2π=的绝对误差 限与相对误差限。（误差限的计算） 解：*2******2),(),(h h r r r h r r h v r h v -+-≤-ππ 绝对误差限为 πππ252.051.02052)5,20(),(2=??+????≤-v r h v 相对误差限为 %420 1 20525) 5,20() 5,20(),(2 ==??≤ -ππv v r h v 6 设x 的相对误差为%a ,求n x y =的相对误差。（函数误差的计算） 解：%* *a x x x =-， )%(* **** *na x x x n x x x y y y n n n =-≤-= - 7计算球的体积，为了使体积的相对误差限为%1，问度量半径r 时允许的相对误差限为多大（函数误差的计算）

硕士研究生数值分析试卷

数值分析（研究生，2008-12-15） 1.（10分）求函数???≤≤++<≤-+=1 0,101,1sin )(2x x x x x x f 在区间[-1，1]上的最佳平方逼近式 x e a x a a x 210)(++=φ。 2．（15分）利用乘幂法计算下列矩阵的主特征值和相应的特征向量 ???? ??????----110141012，初始向量为T x ]0,0,1[0=（要求结果有三位有效数字）。同时计算该矩阵的1-条件数和谱条件数。

3.（15分）已知函数x x f sin )(=在36.0,3 4.0,32.0210===x x x 处的值分别为352274.0,333487.0,314567.0210===y y y 。用Lagrange 插值多项式对3167.0=x 的函数值进行近似计算，并估计近似计算的误差界。

4.（15分）用Newton 迭代法求方程0ln 2=+x x 在区间（0，2 π）内的解，选择你认为合适的初始点，计算方程的根，使得近似解具有四位有效数字。请从理论上估计达到所需精度所需的迭代次数。

5.（15分）用Gauss-Seidel 迭代法解方程组 ?????? ????-=????????????????????---542834*********x x x 取初始近似向量0[0,0,0]T x =，估计达到4位有效数字需要的迭代次数，并实际计算之。就该具体问题分析计算过程中总的乘除法计算量。

6. （10分）应用拟牛顿法解非线性方程组 ?????=-+=-+. 12,2322112221x x x x x x 取T x ]1,0[)0(= ，终止容限210-=ε。 7.（10分） 求解矛盾方程组 ???????=++=++=++=++2 32328.12221321321 321321x x x x x x x x x x x x

数值分析版试题及答案

例1、已知函数表 求() f x的Lagrange二次插值多项式和Newton二次插值多项式。 解： （1）由题可知 插值基函数分别为 故所求二次拉格朗日插值多项式为 （2）一阶均差、二阶均差分别为 均差表为

故所求Newton 二次插值多项式为 例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 ρ=，这样，有 所以，法方程为 011231261192 34a 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 ?=，这样，有 所以，法方程为 解法方程，得到00.8732a =，1 1.6902a =， 故，所求最佳平方逼近多项式为 例4、 用4n =的复合梯形和复合辛普森公式计算积分1?。 解： （1）用4n =的复合梯形公式 由于 2h =，()f x =，()121,2,3k x k k =+=，所以，有 （2）用4n =的复合辛普森公式 由于2h =，()f x =，()121,2,3k x k k =+=，()1 2 220,1,2,3k x k k +=+=，所以，有 例5、 用列主元消去法求解下列线性方程组的解。 解：先消元 再回代，得到33x =，22x =，11x =