数值分析试题(卷)与答案解析
数值分析试题
一、 填空题（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 。
x 0 1 2 y =f (x )
-2
-1
2
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 )
的二阶导数不变号，则初始点x 0
的选取依据为
f(x0)f ”(x0)>0 。
14. 使用迭代计算的步骤为建立迭代函数、 选取初值 、迭代计算。 二、判断题（10×1′）
1、 若A 是n 阶非奇异矩阵，则线性方程组AX ＝b 一定可以使用高斯消元法求解。( × )
2、 解非线性方程f (x )=0的牛顿迭代法在单根x *附近是平方收敛的。 (
)
3、 若A 为n 阶方阵，且其元素满足不等式 )
,...,2,1( 1
n i a a n
i
j j ij ii =≥∑≠=
则解线性方程组AX ＝b 的高斯——塞德尔迭代法一定收敛。 ( × ) 4、 样条插
值
一
种
分
段
插
值
。
(
)
5、 如果插值结点相同，在满足相同插值条件下所有的插值多项式是等价的。 (
)
6、 从实际问题的精确解到实际的计算结果间的误差有模型误差、观测误差、截断误差及舍
入
误
差
。
(
)
7、 解线性方程组的的平方根直接解法适用于任何线性方程组AX ＝b 。 ( × )
8、 迭代解法的舍入误差估计要从第一步迭代计算的舍入误差开始估计,直到最后一步迭
代
计
算
的
舍
入
误
差
。
( × )
9、 数值计算中的总误差如果只考虑截断误差和舍入误差，则误差的最佳分配原则是截断误差
＝
舍
入
误
差
。
(
)
10、插值计算中避免外插是为了减少舍入误差。 ( × ) 三、计算题（5×10′）
1、用列主元高斯消元法解线性方程组。
???
??=++-=+--=+-112123454 3
21321321x x x x x x x x x 解答：
（1，5，2）最大元5在第二行，交换第一与第二行：
???
??=++-=+--=+-1124 123453
21321321x x x x x x x x x L 21=1/5=,l 31=2/5= 方程化为：
???
??=--=+--=+-8.152.06.26.1 0.4 2.0123453
232321x x x x x x x （,）最大元在第三行，交换第二与第三行：
???
??-=+-=--=+-6.1 0.4 2.08.152.06.2123453
232321x x x x x x x L32==,方程化为：
???
??-==--=+-38466.00.38462 8.152.06.2123453
32321x x x x x x 回代得：
???
??-===00010.1 99999.500005.33
21x x x
2、用牛顿——埃尔米特插值法求满足下列表中插值条件的四次插值多项式P 4(x )，并写出其截断误差的表达式(设f (x )在插值区间上具有直到五阶连续导数)。
x i 0 1 2 f (x i ) 1 -1 3 f ’(x i )
1
5
解答： 做差商表 xi F(xi) F[xi,xi+1] F[++2] F[xi,xi+1,xi+2,xi+3] F[xi,xi+1,xi+2,xi+3,xi+4] 0 1 1 -1 -2 1 -1 1 3 2 3 4 3 0 2
3
5
1
-2
-1
P4(x)=1-2x-3x(x-1)-x(x-1)(x-1)(x-2) R4(x)=f(5)()/5!x(x-1)(x-1)(x-2)(x-2)
3、对下面的线性方程组变化为等价的线性方程组，使之应用雅克比迭代法和高斯——赛德尔迭代法均收敛，写出变化后的线性方程组及雅克比迭代法和高斯——赛德尔迭代法的迭代公式，并简单说明收敛的理由。
解答：
交换第二和第四个方程，使系数矩阵为严格对角占优：
雅克比迭代公式：
???????=-+-=-+=+-=+-3 38 4 65 1 2321432431421x x x x x x x x x x x x ???????=+-=-+=-+-=+-6
5 8 4 3 3 1 2431432321421x x x x x x x x x x x x ??????
?=+-=-+=-+-=+-6
5 8 4 3 3 1 2431432321421x x x x x x x x x x x x
《计算机数学基础(2)》数值分析试题
一、单项选择题(每小题3分，共15分)
1. 已知准确值x *与其有t 位有效数字的近似值x ＝…a n ×10s (a 10)的绝对误差x *－
x ( )．
(A) ×10 s －1－t (B) ×10 s －t (C) ×10s ＋1－t (D) ×10 s ＋
t 2. 以下矩阵是严格对角占优矩阵的为( )．
(A) ?
???
?????
???------21001210012100
12，
(B)?
?
???
????
???2100141101410125 (C) ?
?
???????
???--2100
14121241
0125 (D) ??
???
??
??
???-513
114120141112
4 3. 过(0，1)，(2，4)，(3，1)点的分段线性插值函数P (x )=( )
(A) ?????≤<+-≤≤+32103201
23
x x x x (B)
??
???≤<+-≤≤+32103201
2
3
2x x x x (C) ?????≤<+-≤≤-3
2103201
23
x x x x (D)
????
?≤<+-≤≤+3
24201
23
x x x x 4. 等距二点的求导公式是( )
(A) ???????
-='+-='+++)
(1)()(1)(111k k k k k k y y h x f y y h x f
(B) ???
????
-='-='+++)
(1)()(1)(111k k k k k k y y h x f y y h x f
(C) ???
????-='+-='+++)
(1)()(1)(111k k k k k k y y h x f y y h x f
(D)
5. 解常微分方程初值问题的平均形式的改进欧拉法公式是
)(2
1
1c p k y y y +=
+ 那么y p ,y c 分别为( )．
(A) ???+=+=+)
,()
,(1k k k c k k k p y x hf y y y x hf y y (B)
????
?+=+=+)
,()
,(1p k k c k k k p y x hf y y y x hf y y
(C) ?????+=+=),(),(p k k c
k k k p y x f y y y x f y y (D)
????
?+=+=+),()
,(1p k k c
k k k p y x hf y y y x hf y y 二、填空题(每小题3分，共15分)
6. 设近似值x 1,x 2满足(x 1)=，(x 2)=，那么(x 1x 2)= ．
7. 三次样条函数S (x )满足：S (x )在区间[a ,b ]内二阶连续可导，S (x k )=y k (已知)，k =0,1,2,…,n ，且满足S (x )在每个子区间[x k ,x k +1]上是 ．
8. 牛顿－科茨求积公式
∑?
=≈n k k k b
a
x f A x x f 0
)(d )(，则∑=n
k k A 0
＝ .
9. 解方程f (x )=0的简单迭代法的迭代函数(x )满足在有根区间内 ，则在有根区间内任意取一点作为初始值，迭代解都收敛．
10. 解常微分方程初值问题的改进欧拉法预报――校正公式是
预报值：),(1k k k k y x hf y y +=+，校正值：y k +1= ． 三、计算题(每小题15分，共60分)
11. 用简单迭代法求线性方程组
???
??=++=-++=+-36
12363311420238321
321321x x x x x x x x x 的X (3)．取初始值(0,0,0)T ，计算过程保留4位小数．
12. 已知函数值f (0)=6，f (1)=10，f (3)=46，f (4)=82，f (6)=212，求函数的四阶均差f (0,1,3,4,6)和二阶均差f (4，1，3)．
13.将积分区间8等分，用梯形求积公式计算定积分
?
+3
1
2d 1x x ，计算过程保留4位小数．
14. 用牛顿法求115的近似值，取x =10或11为初始值，计算过程保留4位小数． 四、证明题(本题10分)
15. 证明求常微分方程初值问题
??
?=='0
0)()
,(y x y y x f y 在等距节点a =x 0 y (x k +1) y k +1=y k + 2 h [f (x k ,y k )+f (x k +1,y k +1)] 其中h =x k +1－x k (k =0,1,2,…n －1) 《计算机数学基础(2)》数值分析试题答案 一、单项选择题(每小题3分，共15分) 1. A 2. B 3. A 4. B 5. D 二、填空题(每小题3分，共15分) 6. x 2+x 1 7. 3次多项式 8. b －a 9. (x ) r <1 10. y k + )],(),([2 11+++k k k k y x f y x f h hf (x k ＋1, 1+k y ) ． 三、计算题(每小题15分，共60分) 11. 写出迭代格式 ?????++--=+++-=+-+=+++3025.05.039090.006363.05.225.0375.00)(2)(1)1(3 ) (3)(1) 1(2 ) (3)(2)1(1k k k k k k k k k x x x x x x x x x X (0)=(0,0,0)T . ?????=++?-?-==+?++?-==+?-?+=3 30025.005.03309090.0006363.05.25.2025.00375.00) 1(3) 1(2)1(1x x x 得到X (1)＝，3，3)T ?????=++?-?-==+?++?-==+?-?+=0 000.130325.05.25.07363.2339090.005.26363.0875.25.2325.03375.00) 2(3) 2(2)2(1x x x 得到X (2)=， 7， 0)T ?????=++?-?-==+?++?-==+?-?+=6 971.0307363.225.0875.25.06045.2319090.00875.26363.04136.35.2125.07363.2375.00) 3(3) 3(2)3(1x x x 得到X (3)= 4， 6， 6)T . 12. 计算均差列给出． k x f (x k ) 一阶均差 二阶均差 三阶均差 四阶均差 0 6 1 10 4 3 46 18 14/3 4 82 36 6 1/3 6 212 65 29/3 11/15 1/15 f (0,1,3,4,6)= 15 f (4, 1, 3)=6 13. f (x )=21x +,h =25.08 2 =．分点x 0=，x 1=，x 2=，x 3=，x 4=，x 5=，x 6=，x 7=，x 8=. 函数值：f = 2，f = 8，f = 8，f = 6，f = 1，f = 2，f = 6，f = 2，f = 3． )()([2 d )(803 1 x f x f h x x f +=? ))]()()()()()()((27654321x f x f x f x f x f x f x f +++++++ (9分) = 2 25 .0×[ 2+ 3+2× 8+ 8+ 6 + 1+ 2+ 6+ 2)] =× 5+2× 3)= 1 14. 设x 为所求，即求x 2－115=0的正根．f (x )=x 2－115． 因为f (x )=2x ，f (x )=2，f (10)f (10)=(100－115)×2<0，f (11)f (11)=(121－115)×2>0 取x 0=11． 有迭代公式 x k +1=x k －)()(k k x f x f '=k k k k k x x x x x 2115 221152 + =--(k =0,1,2,…) x 1=11 2115211?+＝ 3 x 2= 3727.102115 23727.10?+＝ 8 x 3= 8 723.102115 28723.10?+＝ 8 x * 8 四、证明题(本题10分) 15. 在子区间[x k +1,x k ]上，对微分方程两边关于x 积分，得 y (x k +1)－y (x k )= ? +1 d ))(,(k k x x x x y x f 用求积梯形公式，有 y (x k +1)－y (x k )=))](,())(,([2 11+++k k k k x y x f x y x f h 将y (x k ),y (x k +1)用y k ,y k +1替代，得到 y (x k +1) y k +1=y k + 2 h [f (x k ,y k )+f (x k +1,y k +1)](k =0,1,2,…,n －1) 数值分析期末试题 一、填空题（20102=?分） （1)设???? ??????---=283012251 A ，则 =∞A ______13_______。 （2)对于方程组? ??=-=-34101522121x 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）已知?? ? ? ??=1021A ，则条件数=∞)(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 的水平直线是)(313 1 ∑==i i x f y 。 二、（10分）证明：方程组??? ??=-+=++=+-1 2112321 321321x x x x x x x x x 使用Jacobi 迭代法求解不收敛性。 证明：Jacobi 迭代法的迭代矩阵为 ?? ??? ?????---=05.05.01015.05.00J B J B 的特征多项式为 )25.1(5.05.011 5 .05.0)det(2+=---=-λλλ λλ λj B I J B 的特征值为01=λ，i 25.12=λ，i 25.13-=λ，故25.1)(=J B ρ＞1，因而迭代法不收 敛性。 三、（10分）定义内积 ? = 1 )()(),(dx x g x f g f 试在{}x Span H ,11=中寻求对于x x f =)(的最佳平方逼近元素)(x p 。 解：1)(0≡x ?，x x ≡)(1?， 1),(1 00== ? dx ??，2 1 ),(1 01= = ? xdx ??，3 1 ),(1 211= =? dx x ??，3 2),(1 0= =? dx x f ?，5 2),(1 1= = ? dx x x f ?。 法方程 ???? ??????= ??????????? ?? ???5232312 1211 10c c 解得1540= c ，15 121=c 。所求的最佳平方逼近元素为 x x p 15 12 154)(+=，10≤≤x 四、（10分）给定数据表 试用三次多项式以最小二乘法拟合所给数据。 解:3 32210)(x c x c x c c x y +++= ????? ?? ?????????----=84211111000111118421A ， ????? ? ??????= 130034003401034010001005A A T T T y A )4.14,7,2.4,9.2(= 法方程 y A Ac A T T = 的解为4086.00=c ，39167.01=c ，0857.02=c ，00833.03=c 得到三次多项式 3200833.00857.039167.04086.0)(x x x x y +++= 误差平方和为000194.03=σ 五. (10分) 依据如下函数值表 建立不超过三次的Lagrange 插值多项式，用它计算)2.2(f ，并在假设1)()4(≤x f 下，估计计算误差。 解:先计算插值基函数 147 8781)40)(20)(10()4)(2)(1()(230+-+-=------= x x x x x x x l x x x x x x x l 38 231)41)(21)(01()4)(2)(0()(231+-=------= x x x x x x x l -+-=------= 2324 5 41)42)(12)(02()4)(1)(0()( x x x x x x x l 12 1 81241)24)(14)(04()2)(1)(0()(233+-=------= 所求Lagrange 插值多项式为 12 1 445411)(3)(23)(9)()()()(2332103 03+-+- =+++==∑=x x x x l x l x l x l x l x f x L i i i 从而0683.25)2.2()2.2(3=≈L f 。 据误差公式))()()((! 4) ()(3210)4(3x x x x x x x x f x R ----= ξ及假设1)()4(≤x f 得误差估计： 0396.09504.0! 41 )42.2)(22.2)(12.2)(02.2(! 4)()()4(3=?≤ ----= ξf x R 六. (10分) 用矩阵的直接三角分解法解方程组 ? ? ??? ? ??????=????????????????????????71735 3010 3421101002014321x x x x 解 设 ????? ?????? ???????? ?????=????????????443433 24232243 42 41 3231 2102011111 3010342110100201u u u u u u l l l l l l 由矩阵乘法可求出ij u 和ij l ????? ???????=????????? ???1010121101111143 42 41 323121 l l l l l l ????? ? ??????=????????? ?? ?2121010201020 1443433 242322 u u u u u u 解下三角方程组 ????? ???????=????????????????????????7173510101211014321y y y y 有51=y ，32=y ，63=y ，44=y 。再解上三角方程组 ? ???? ???????=????????????????????????463521************x x x x 得原方程组的解为11=x ，12=x ，23=x ，24=x 。 七. (10分) 试用Simpson 公式计算积分 dx e x ? 2 1 1 的近似值, 并估计截断误差。 解： 0263.2)4(61 221 5.112 1 1 =++-≈? e e e dx e x x e x x x x f 1 5678) 4()24 36121(+++= 43.198)1()(max )4()4(2 1==≤≤f x f x 截断误差为 06890.0)(max 2880)12()4(2 15 2=-≤≤≤x f R x 八. (10分) 用Newton 法求方程2ln =-x x 在区间) ,2(∞内的根, 要求 81 10--<-k k k x x x 。 解：此方程在区间) ,2(∞内只有一个根s ，而且在区间（2，4）内。设 2ln )(--=x x x f 则 x x f 11)('- =， 21 )(''x x f = Newton 法迭代公式为 1) ln 1(112ln 1-+= - --- =+k k k k k k k k x x x x x x x x ， Λ,2,1,0=k 取30=x ，得146193221.34=≈x s 。 九. (10分) 给定数表 求次数不高于5的多项式)(5x H ，使其满足条件 ???====2 ,0 ),()(3 ,2 ,1 ,0 ),()(' 55i x f x H i x f x H i i i i 其中,1i x i +-= 3 ,2 ,1 ,0=i 。 解：先建立满足条件 )()(3i x f x p =， 3,2,1,0=i 的三次插值多项式)(3x p 。采用Newton 插值多项式 [][]))((,,)(,)()(1021001003x x x x x x x f x x x x f x f x p --+-+=+ []))()((,,,2103210x x x x x x x x x x f --- )1()1(6 1 )1()1(410-+-+-++=x x x x x x 326 1 61914x x x --+ = 再设 )2)(1()1)(()()(35--+++=x x x x b ax x p x H ，由 ???=-++==-+-+-=-1 .0)2)(()1()1(1 )6)(()1()1(' 3'5' 3'5b a p H b a p H 得 ?? ??? = +=+-6017811b a b a 解得36059- =a ，360 161 =b 。 故所求的插值多项式 )2)(1()59161(360 1 6161914)(2325---+--+ =x x x x x x x x H 数值分析试题 一、 填空题（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 卷） 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 求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. 已知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)( 二 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设，求解方程组，求雅可比迭代法与高斯-赛德尔迭代法收敛的充要条件。 解雅可比法的迭代矩阵 ， 故雅可比法收敛的充要条件是。 高斯-赛德尔法的迭代矩阵 ， 数值分析整理版试题及答案 例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 ρ=，这样，有 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、（本题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、 已知函数表 求()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 +=?=-， 数值分析模拟试卷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、讨论绝对误差（限）、相对误差（限）与有效数字之间的关系。 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 一、单项选择题（每小题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 --=?+?=-+--% 数值分析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 《数值计算方法》复习试题 一、填空题： 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.答案 数值分析期末试题 一、填空题（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 若误差限为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 时允许的相对误差限为多大（函数误差的计算） 一、（8分）用列主元素消去法解下列方程组： ??? ??=++-=+--=+-11 2123454 321321321x x x x x x x x x 二、（10分）依据下列数据构造插值多项式：y(0)=1,y(1)= —2,y '(0)=1, y '(1)=—4 三、（12分）分别用梯形公式和辛普生公式构造 复化的梯形公式、复化的辛普生公式并利用复化的梯形公式、复化的辛普生公式计算下列积分： ? 9 1dx x n=4 四、（10分）证明对任意参数t ，下列龙格－库塔方法是二阶的。 五、（14分）用牛顿法构造求c 公式，并利用牛顿法求115。保留有效数字五位。 六、（10分）方程组AX=B 其中A=????????? ?10101a a a a 试就AX=B 建立雅可比迭代法和高斯-赛德尔迭代法，并讨论a 取何值时 迭代收斂。 七、（10分）试确定常数A,B,C,a,使得数值积分公式?-++-≈2 2 ) (}0{)()(a Cf Bf a Af dx x f 有尽可能多的 代数精确度。并求该公式的代数精确度。 八、{6分} 证明： A ≤ 其中A 为矩阵，V 为向量. 第二套 一、（8分）用列主元素消去法解下列方程组： ??? ??=++=+-=+3 2221 43321 32132x x x x x x x x 二、（12分）依据下列数据构造插值多项式：y(0)=y '(0)=0, y(1)=y '(1)= 1,y(2)=1 三、（14分）分别用梯形公式和辛普生公式构造 复化的梯形公式、复化的辛普生公式,并利用复化的梯形公式、 复化的辛普生公式及其下表计算下列积分： ?2 /0 sin πxdx ????? ? ? -+-+=++==++=+1 3121231)1(,)1(() ,(),()(2 hk t y h t x f k thk y th x f k y x f k k k h y y n n n n n n n n 《计算方法》期中复习试题 一、填空题： 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 ，为了减少舍入误差，应将表达式数值分析试题及答案汇总
数值分析学期期末考试试题与答案(A)
数值分析试卷及答案
数值分析试卷及其答案
数值分析期末考试复习题及其答案.doc
数值分析试卷及答案
数值分析整理版试题及答案
数值分析试卷及其答案1
数值分析试卷及其答案2
数值分析整理版试题及答案
北航2010-2011年研究生数值分析期末模拟试卷1-3
数值分析作业思考题汇总
数值分析试题及答案
数值分析期末试卷
数值计算方法试题集及答案要点
数值分析试题及答案
数值分析期末试题
数值分析最佳习题(含答案)
数值分析期末试题
《数值计算方法》精彩试题集及问题详解1-6 2