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

期末考试试卷（A 卷）

2007学年第二学期考试科目：数值分析考试时间：120 分钟

学号姓名年级专业

一、判断题（每小题2分，共10分）

1. 用计算机求

1000

n n

=∑时，应按照n 从小到大的顺序相加。（）

2. 为了减少误差,

进行计算。（）

3. 用数值微分公式中求导数值时，步长越小计算就越精确。（）

4. 采用龙格－库塔法求解常微分方程的初值问题时，公式阶数越高，数值解越精确。（）

5. 用迭代法解线性方程组时，迭代能否收敛与初始向量的选择、系数矩阵及其演变方式有

关，与常数项无关。（）

二、填空题（每空2分，共36分）

1. 已知数a 的有效数为0.01，则它的绝对误差限为________，相对误差限为_________.

2. 设1010021,5,1301A x -????????=-=-????????-????

则1A =_____，2x =______，Ax ∞

=_____.

3. 已

知

53()245,

f x x x x =+-则

[1,1,0]f -= ,[3,2,1,1,2,3]f ---= .

4. 为使求积公式

1231

()((0)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 产

生的向量序列{}

()k X 收敛的充分必要条件是 .

7. 使用消元法解线性方程组AX B =时，系数矩阵A 可以分解为下三角矩阵L 和上三角矩

阵U 的乘积，即.A LU = 若采用高斯消元法解AX B =，其中4221A -??

???

，则L =_______________，U =______________；若使用克劳特消元法解AX B =，则

11u =____；若使用平方根方法解AX B =，则11l 与11u 的大小关系为_____（选填：>，

<，=，不一定）。

8. 以步长为1的二阶泰勒级数法求解初值问题(0)1y x y

y '=+??=?

的数值解，其迭代公式为

___________________________.

三、计算题（第1～3、6小题每题8分，第4、5小题每题7分，共46分）

1. 以02x =为初值用牛顿迭代法求方程3

()310f x x x =--=在区间(1,2)内的根，要求

（1）证明用牛顿法解此方程是收敛的；

（2）给出用牛顿法解此方程的迭代公式，并求出这个根（只需计算12,,x x 计算结果

取到小数点后4位）。

2. 给定线性方程组

1231231230.40.410.40.820.40.83

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

++=??++=?

（1）分别写出用Jacobi 和Gauss-Seidel 迭代法求解上述方程组的迭代公式；

（2）试分析以上两种迭代方法的敛散性。

3. 已知函数()y f x =在如下节点处的函数值

（1）（2）根据后三个节点建立二阶牛顿后插公式2()P x ，并计算(1.1)y 的近似值；

（3）采用事后估计法计算（2）中近似值的截断误差（结果保留四位小数）。

4.已知如下数据表，试用最小二乘法求它的二次最小平方逼近多项式。

5. 已知函数()y f x =在以下节点处的函数值，利用差商表求(3)f '和(3)f ''的近似值。

6. 写出前进欧拉公式、后退欧拉公式，并由这两个公式构造一个预估－校正公式求解下列

常微分方程的数值解。

22(01,0.2)(0)0

y x y x h y '?=+≤≤=?

四、（8分）已知n+1个数据点(,)(0,1,2,,)i i x y i n L ，请用多种方法建立这些数据点之间的函数关系，并说明各种函数的适用条件。

期末考试答案及评分标准（A 卷）

2007学年第二学期考试科目：数值分析

一、判断题：（每小题2分，共10分）

1. ×

2. √

3. ×

4. ×

5. ×

二、填空题：（每空2分，共36分） 1. 0.005或2

0.510-? ，0.5

2. 3. 0,2 4. 1,0,1,3 5. ()A A ρ≤ 6. ()1M ρ<

7. 1042,,1,10212??

-????=????????

8. 11

()(1)2

n n n n n n y y x y x y +=+++

++或1 1.5 2.50.5,0,1,2,n n n y x y n +=++=L

三、解答题（第1～4小题每题8分，第5、6小题每题7分，共46分） 1. （1）证明：3

()31f x x x =--，由于

a) (1)30,(2)10,f f =-<=> b) 2()330

((1,2)),f x x x '=-≠∈

()60((1,2)),f x x x ''=>∈ 即()f x ''在(1,2)上不变号，

d) 对于初值02x =，满足(2)(2)0,f f ''>

所以用牛顿迭代法求解此方程是收敛的。

………………………………………4分

（2）解：牛顿迭代法的迭代公式为

()31

()33

n n n n n n n n f x x x x x x f x x +--=-=-'- ………………………………………2分

取初值02x =进行迭代，得

1 1.8889,x =

………………………………………1分

2 1.8795.x =

………………………………………1分

2. 解：（1）Jacobi 迭代公式为

(1)()()

123(1)()()

13(1)()()3

120.40.410.40.820.40.83

k k k k k k k k k x x x x x x x x x +++?=--+?=--+??=--+? ……………………………2分 Gauss-Seidel 迭代公式为

(1)()()

123(1)(1)()

13(1)(1)(1)3

120.40.410.40.820.40.83

k k k k k k k k k x x x x x x x x x ++++++?=--+?=--+??=--+?……………………………2分（2）Jacobi 迭代矩阵的特征方程为0.40.4

0.4

0.800.40.8λ

λλ

=，展开得

30.960.2560λλ-+=，即

(0.8)(0.40.40λλλ-+++-=，

从而得 123-1.0928,0.8000,0.2928λλλ===，（或由单调性易判断必有一个大于1

的特征根，）因此迭代矩阵的谱半径等于必大于1，所以Jacobi 迭代法发散。

……………………………2分

Gauss-Seidel 迭代矩阵的特征方程为0.40.4

0.40.800.40.8λ

λλ

=，展开得

2(0.8320.128)0λλλ-+=，解得1230,0.628,0.204,λλλ=≈≈迭代矩阵的谱半径

小于1，所以Gauss-Seidel 迭代法收敛。

……………………………2分

（2）建立牛顿后插公式为

2232

022********

()()()()!!

()()()P x x x x x x x x =-

----=-----=-+ 则所求近似值为

211279(.).P =

………………………………………3分

（3）根据前三个节点建立牛顿后插公式为

12214

31112312124

()()()()!!

()()P x x x x x x x x x =-

---=----=-++ 则 1211268()

(.).P =

根据事后误差估计法

12222

09091

()()(.)(.)x R x P P x -??≈

-??+ 故截断误差

209

112792680047121

.(.)(..)..R -≈

?-≈- ………………………………………3分

4. 解：设所求二次最小平方逼近多项式为2

2012().P x a a x a x =++ 根据已知数据，得

01211111002,,11151

240a M A a Y a -????

?????????

???===????

??????????????

……………………………2分

则

4268268,468186M M M Y ????????''==????

????????

……………………………1分

建立法方程组为

0124268268468186a a a ????????????=???????????

??????? ……………………………2分

解得

0123.5, 1.5, 1.5.a a a ===-

……………………………1分

从而得所求一次最小平方逼近多项式为2

1() 3.5 1.5 1.5.P x x x =+-

……………………………1分

5. 解：设2()P x 为已知节点数据的插值二次多项式。构造如下差商表：

2分

因为二次多项式的二阶差商为常数，又2()P x 是

()f x 的插值函数，故有

225[4,3,3][3,3,3]2

P P ==

……………………………2分

而

22[3,3]75

[4,3,3]342

P P -=

=-，

因此得

[3,3]2

P =

， ……………………………1分

由于

L 1442443

()()![,,,,]k n k f x k P x x x x +≈，

从而得

3332

()[,],f P '==

2323335()![,,].f P ''==

……………………………2分

6. 解：前进欧拉公式：

1(,)0.20.2n n n n n n n y y h f x y y x y +=+?=++…………1分

后退欧拉公式：

2211111(,)0.20.2n n n n n n n y y h f x y y x y +++++=+?=++ ……1分

预估时采用欧拉公式

*22

10.20.2n n n n y y x y +=++

……………………………1分

校正时采用后退欧拉公式

()

10.20.2n n n n y y x

+++=++

……………………………1分

由初值000002,,.x y h ===知，节点分别为0.2,(1,2,3,4,5)i x i i ==

当1

0.2,x =

*2210000.20.20,y y x y =++=

()2

2101102020008*...y y x y =++=，

……………………………1分

当2

0.4,x =

*2221110.20.20.0160,y y x y =++≈

()

212

020200401*...y y x y

=++≈.

……………………………1分

当3

0.6,x =

*2232220.20.20.0724,y y x y =++≈

()

323

020201131*...y y x y

=++≈.

……………………………1分

当4

0.8,x =

*2243330.20.20.1877,y y x y =++≈

()

434

020202481*...y y x y

=++≈.

……………………………1分

当5

1.0,x =

*2254440.20.20.3884,y y x y =++≈

()2

25455020204783*...y y x y =++≈.

四、（8分）

答：1、可以建立插值函数：（1）Newton 基本差商公式

00100121001110()()()[,]()()[,,]

()()()[,,,]

L L L n n n P x f x x x f x x x x x x f x x x x x x x x x f x x x -=+-+--++---

……………………………1分

（2）Lagrange 插值多项式

0011()()()()()L L n i i n n L x a f x a f x a f x a f x =+++++

其中01101101()()()()

,(,,,)()()()()

L L L L L i i n i i i i i i i n x x x x x x x x a i n x x x x x x x x -+-+----=

=----.

……………………………1分

这两类插值函数的适用条件是：n 不太大；而且要求函数严格通过已知数据点。

……………………………2分

2、可以建立拟合函数：

2012()m m m P x a a x a x a x =++++L

……………………………1分

其中系数012,,,,n a a a a L 满足法方程组M MA M Y ''=,

00000021111112()1()1,,()1

m m m n n n

n n a f x y x x x a f x y x x x M A Y a f x y x x x ????????

?????????

???????====??

????????????????????????

K K

L L L K K K K

K K

……………………………1分

拟合函数的适用条件是：n 比较大，而且并不要求函数严格通过已知数据点，或者已知数据点本身的误差较大。

……………………………2分

数值分析试题及答案汇总

数值分析试题一、填空题（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 》和《数值方法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 =按递推公式 11783 100 n n Y Y -=- ( n=1,2,…) 计算到100Y .若取783≈27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字(783≈27.982). 8. 当N 充分大时,怎样求 2 11N dx x +∞+?? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设2 12S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加, 而相对误差却减小. 11. 序列{}n y 满足递推关系1101 n n y y -=-(n=1,2,…),若02 1.41y =≈(三位有效数字),计算到10 y 时误差有多大?这个计算过程稳定吗? 12. 计算6 (21)f =-,取 2 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 6 3 11,(322), ,9970 2. (21) (322) --++ 13. 2 ()ln(1)f x x x =- -,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式 2 2 ln(1)ln(1)x x x x - -=-+ + 计算,求对数时误差有多大? 14. 试用消元法解方程组{ 10 10 12121010; 2. x x x x +=+=假定只用三位数计算,问结果是否可靠? 15. 已知三角形面积 1sin , 2 s ab c = 其中c 为弧度, 02c π << ,且测量a ,b ,c 的误差分别为,,.a b c ???证明面积的误差s ?满足 . s a b c s a b c ????≤ ++ 第二章插值法 1. 根据( 2.2)定义的范德蒙行列式,令

数值分析试卷及答案

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

数值分析学期期末考试试题与答案(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、已知函数表 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 ρ=，这样，有

数值分析试卷及答案

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

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

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

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

数值分析试题

《计算机数学基础(下)》数值分析试题 2000、8 之六(2002、7已用) 一、单项选择题(每小题3分，共15分) 1.数值x *的近似值x =0.1215×10－ 2，若满足≤-*x x ( )，则称x 有4位有效数字. (A) 21×10－3 (B) 21×10－4 (C) 21×10－5 (D) 2 1×10－6 2. 设矩阵A ＝?? ?? ? ?????------52111021210，那么以A 为系数矩阵的线性方程组A X ＝b 的雅可比迭代矩阵为( ) (A)??????????04.02.01.002.01.02.00 (B) ???? ? ?? ???14.02 .01.012.01.02.01 (C) ??????????------04.02.01.002.01.02.00 (D) ???? ??????021 102120 3. 已知y =f (x )的均差f (x 0,x 1,x 2)=314，f (x 1,x 2,x 3)=315，f (x 2,x 3,x 4)=15 91，f (x 0,x 2,x 3)=318 , 那么均差f (x 4,x 2,x 3)=( ) (A) 315 (B) 318 (C) 1591 (D) 3 14 4. 已知n =4时牛顿－科茨求积公式的科茨系数,15 2,4516,907)4(2)4(1) 4(0===C C C 那么 )4(3C ＝( ) 90 39 152********)D (152)C (4516)B (907)A (=--- 5.用简单迭代法求方程的近似根，下列迭代格式不收敛的是( ) (A) e x －x －1＝0，[1,1.5]，令x k +1=1e -k x (B) x 3－x 2－1=0，[1.4,1.5], 令211 1k k x x +=+ (C) x 3－x 2－1=0，[1.4,1.5], 令32 11k k x x +=+ (D) 4－2x =x ，[1,2], 令)4(log 21x x k -=+ 二、填空题(每小题3分，共15分) 6.sin1有2位有效数字的近似值0.84的相对误差限是 . 7.设矩阵A 是对称正定矩阵，则用迭代法解线性方程组A X =b ，其迭代解数列一定收敛. 8. 已知f (1)=1,f (2)=3,那么y =f (x )以x =1，2为节点的拉格朗日线性插值多项式为 . 9. 用二次多项式2210)(x a x a a x ++=?,其中a 0, a 1, a 2是待定参数，拟合点(x 1,y 1),(x 2,y 2),…,(x n ,y n ). 那么参数a 0, a 1, a 2是使误差平方和取最小值的解. 10. 设求积公式 ∑?=≈n k k k b a x f A x x f 0 )(d )(，若对的多项式积分公式

数值分析习题集及答案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、（本题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 迭代格式为

数值分析试题集

2 A J ：；［则 || A 「一— 仙二 ------------- 'a+1 2 3 设「_1 J ，当a 满足条件时，A 可作LU 分解。 (试卷一) 一 (10 分)已知％ =1.3409, x 2 =1.0125都是由四舍五入产生的近似值，判断x-i x 2及x 1 - x 2 有几位有效数字。二 ( 1 多项式三(15分)设f(x)? C 4［a,b ］，H (x )是满足下列条件的三次多项式 H (a)二 f (a) , H (b)二 f (b) , H (c) = f (c) , H (c)二 f (c) ( a ：：： c ：： b ) 求f (x) -H(x)，并证明之。 1 四(15分)计算， : =10』。 o 1 +X 五(15分)在［0，2］上取X 。= 0, X 1 = 1, X 2 = 2，用二种方法构造求积公式，并给出其公式的代数精度。六(10分)证明改进的尢拉法的精度是 2阶的。七(10分)对模型y ■ = ■?y , ■：■ 0，讨论改进的尢拉法的稳定性。八(15分)求方程x 3 4x 2 - 7x - 1 = 0在-1.2附近的近似值，；=10 "。 (试卷二) 一填空(4*2分) 1 { k (x) }k￡是区间［0，1］上的权函数为'(x)=x 2的最高项系数为1的正交多项式族，其中 1 （x ） =1，贝y . X 0（ x ）dx = ------------ ， 1（X ）工 ------- 数值分析试题集

3 2 * * * 4设非线性方程f (x)二(x -3x - 3x -1)(x ? 3) = 0，其根& = -3 ,他 =-1，则求为的近似值时，二阶局部收敛的牛顿迭代公式是 -------------------------------------- 。广1 —0.5 a ' 二(8 分)方程组AX=b，其中A= — 0.5 2 -0.5，X, R3 l -a -0.5 1 』 1试利用迭代收敛的充要条件求出使雅可比迭代法收敛的a的取值范围，a取何值时雅可比迭代收敛最快？ 2选择一种便于计算的迭代收敛的充要条件，求出使高斯-塞德尔迭代法收敛的a的取值范围。 "V " = f(X y) 三(9分)常微分方程初值问题丿'的单步法公式为y n* = y n」+2hf (x n, y n)，求该、、y°= y(x°) 公式的精度。四(14分)设A X =b为对称正定方程组 1求使迭代过程X k 1二X k ?〉(b-A?X k)收敛的数〉的变化范围；『2 -1 -1、、 1、『0 、 2用此法解方程组-12 0-X2=1 L1 0

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

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____.

数值分析作业思考题汇总

￥数值分析思考题1 1、讨论绝对误差（限）、相对误差（限）与有效数字之间的关系。 2、相对误差在什么情况下可以用下式代替 3、查阅何谓问题的“病态性”，并区分与“数值稳定性”的不同点。 4、取，计算，下列方法中哪种最好为什么（1）(3 3-，（2）(2 7-，（3） ()3 1 3+ ，（4） ()6 1 1 ，（5）99- , 数值实验数值实验综述：线性代数方程组的解法是一切科学计算的基础与核心问题。求解方法大致可分为直接法和迭代法两大类。直接法——指在没有舍入误差的情况下经过有限次运算可求得方程组的精确解的方法，因此也称为精确法。当系数矩阵是方的、稠密的、无任何特殊结构的中小规模线性方程组时，Gauss消去法是目前最基本和常用的方法。如若系数矩阵具有某种特殊形式，则为了尽可能地减少计算量与存储量，需采用其他专门的方法来求解。 Gauss消去等同于矩阵的三角分解，但它存在潜在的不稳定性，故需要选主元素。对正定对称矩阵，采用平方根方法无需选主元。方程组的性态与方程组的条件数有关，对于病态的方程组必须采用特殊的方法进行求解。数值计算方法上机题目1 1、实验1. 病态问题实验目的：算法有“优”与“劣”之分，问题也有“好”和“坏”之别。所谓坏问题就是问题本身的解对数据变化的比较敏感，反之属于好问题。希望读者通过本实验对此有一个初步的体会。数值分析的大部分研究课题中，如线性代数方程组、矩阵特征值问题、非线性方程及方程组等都存在病态的问题。病态问题要通过研究和构造特殊的算法来解决，当然一般要付出一些代价（如耗用更多的机器时间、占用更多的存储空间等）。 $ r e x x e x x ** * ** - == 141 . ≈)61

数值分析习题集及答案

数值分析习题集（适合课程《数值方法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 -=( n=1,2,…) 计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字27.982). 8. 当N 充分大时,怎样求 2 11N dx x +∞ +? ? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 2 12S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加,而相对误差却减小. 11. 序列{}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y ≈(三位有效数字),计算到10y 时误差有多大?这个计算过程稳定吗?

数值分析试题及答案

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

数值分析试题及答案

数值分析试题及答案一、单项选择题（每小题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的改进的欧拉公

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

数值分析历年考题

数值分析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 =___________，此时该方法是几阶的：___________

《数值计算方法》精彩试题集及问题详解1-6 2

《计算方法》期中复习试题一、填空题： 1、已知3.1)3(,2.1)2(,0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得 ?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。答案：2.367，0.25 2、1)3(,2)2(,1)1(==-=f f f ，则过这三点的二次插值多项式中2 x 的系数为，拉格朗日插值多项式为。答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 3、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 4、设)(x f 可微,求方程)(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 ),=]4,3,2,1,0[f ( 0 )； 6、计算方法主要研究( 截断 )误差和( 舍入 )误差； 7、用二分法求非线性方程f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为 ( 1 2+-n a b )； 8、已知f (1)＝2，f (2)＝3，f (4)＝5.9，则二次Newton 插值多项式中x 2系数为 ( 0.15 )； 11、两点式高斯型求积公式?1 0d )(x x f ≈(?++-≈1 )] 321 3()3213([21d )(f f x x f )，代数精度为( 5 )； 12、为了使计算 32)1(6 )1(41310-- -+-+ =x x x y 的乘除法次数尽量地少，应将该表达式改写为 11 ,))64(3(10-= -++=x t t t t y ，为了减少舍入误差，应将表达式