数值分析整理版试题及答案
例1、 已知函数表
求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。 解:
(1)
插值基函数分别为
()()()()()()()()()()
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 ---===
-----=
==----===
---
故所求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
1
200110
1
1
2011000
1
210
1
,11,
,3
1
23
,,,
,3226
9,324
dx x dx xdx f x x dx f x x x dx ??????????====
====++=
=++=
????? 所以,法方程为
01123126119234a 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 ?=,这样,有
()()()()()()1
0001
21101
01100
1
001
10,111,3
1
,,2
, 1.7183,1
x x dx x dx xdx f e dx f xe dx ??????????====
=======?????
所以,法方程为
0111 1.718321112
3a a ??
??????=????????????????
解法方程,得到00.8732a =,1 1.6902a =, 故,所求最佳平方逼近多项式为
*1()0.8732 1.6902S x x =+
例4、 用4n =
的复合梯形和复合辛普森公式计算积分1
?
。
解:
(1)用4n =的复合梯形公式
由于2h =,(
)f x =()121,2,3k x k k =+=,所以,有
()()(
)
4
1
3
1
[129]22
22
17.2277
k k T h
f f x f =≈=++=+?=?∑
(2)用4n =的复合辛普森公式
由于2h =,(
)f x =()121,2,3k x k k =+=,()12
220,1,2,3k x
k k +
=+=,所以,有
()()(
)
4
1
3
310
1
2[1429]61
[1423]
3
17.3321
k k k k S h
f f x f x f +==≈??
?=+++ ??
?=+?+?
+=?∑∑
例5、 用列主元消去法求解下列线性方程组的解。
123123123123315183156x x x x x x x x x -+=??
-+-=-??++=?
解:先消元
()1212331518311511161831151233151116r r A b ?-??
??=---??
????---???????→-??
????
212131312332322
,31,186
,71831
1501750
7171831183
115076171831601735183107617100m m m m r r m m =-?-+→=-?-+→?=-?-+→---??
???????????????→-?????
?
---???????→??
??-??
--?????????????→第1行()第2行第2行
第1行()第3行第3行第2行()第3行第3行
158********-??????
????
再回代,得到33x =,22x =,11x =
所以,线性方程组的解为11x =,22x =,33x =
例6、 用直接三角分解法求下列线性方程组的解。
1231231231
11945611183451
282x x x x x x x x x ?++=??
?++=???++=??
解: 设
1112
1321222331323311
145610011110034510
011
22
u u u A l u u LU l l u ??
??????
???
????
?===????????????????
??????
则由A LU =的对应元素相等,有
1114u =
,1215u =,131
6u =, 2111211433l u l =?=,3111311
22
l u l =?=,
2112222211460l u u u +=?=-,2113232311
545
l u u u +=?=-,
3112322232136l u l u l +=?=-,31133223333313
215l u l u u u ++=?=
因此,
1
1
11004
564
11100
360452361130
15A LU ????????????
?
?==--?????
???-?
??????
?
解Ly b =,即1231
094
108382361y y y ??
????
???????
?=????????????????-??
,得19y =,24y =-,3154y =- 解Ux y =,即1
231
1
14569110
460451541300
15x x x ??
???????
???????--=-???
???????-???????????
?
,得3177.69x =-,2476.92x =,1227.08x =- 所以,线性方程组的解为1227.08x =-,2476.92x =,3177.69x =-
1、若A 是n n ?阶非奇异阵,则必存在单位下三角阵L 和上三角阵U ,使
LU A =唯一成立。 ( )
2、当8≥n 时,Newton -cotes 型求积公式会产生数值不稳定性。( )
3、形如)
()(1i n
i i b
a x f A dx x f ∑?=≈的高斯(Gauss )型求积公式具有最高代数精
确度的次数为12+n 。 ( )
4、矩阵??
??? ??=210111012A 的2-范数2A =9。( )
5、设??
??? ??=a a a a A 000002,则对任意实数0≠a ,方程组b Ax =都是病态的。(用
∞
?) ( )
6、设n n R A ?∈,n
n R
Q ?∈,且有I Q Q T
=(单位阵),则有22QA A =。
( )
7、区间[]b a ,上关于权函数)(x W 的直交多项式是存在的,且唯一。( )1、( Ⅹ ) 2、( ∨ ) 3、( Ⅹ ) 4、( ∨ ) 5、( Ⅹ ) 6、( ∨ )7、( Ⅹ ) 8、( Ⅹ )
一、判断题(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、插值计算中避免外插是为了减少舍入误差。 ( × )
1. 用计算机求
1000
1000
1
1
n n
=∑时,应按照n 从小到大的顺序相加。 ( )
2. 为了减少误差,
进行计算。 ( 对 )
3. 用数值微分公式中求导数值时,步长越小计算就越精确。 ( )
4. 用迭代法解线性方程组时,迭代能否收敛与初始向量的选择、系数矩阵及其演变方式有关,与
常数项无关。 ( )
复习试题
一、填空题:
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 ,则过这三点的二次插值多项式中2
x 的系数为 ,
拉格朗日插值多项式为 。
答案:-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 );
10、已知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、 解线性方程组A x =b 的高斯顺序消元法满足的充要条件为(A 的各阶顺序主子式均
不为零)。
13、 为了使计算
32)1(6
)1(41310--
-+-+
=x x x y 的乘除法次数尽量地少,应将该表
达式改写为
11
,))64(3(10-=
-++=x t t t t y ,为了减少舍入误差,应将表达式
19992001-改写为
199920012
+ 。
14、 用二分法求方程01)(3
=-+=x x x f 在区间[0,1]内的根,进行一步后根的所在区间
为 0.5,1 ,进行两步后根的所在区间为 0.5,0.75 。 15、 计算积分?1
5
.0d x
x ,取4位有效数字。用梯形公式计算求得的近似值为 0.4268 ,
用辛卜生公式计算求得的近似值为 0.4309 ,梯形公式的代数精度为 1 ,辛卜生公式的代数精度为 3 。
16、 求解方程组??
?=+=+042.01532121x x x x 的高斯—塞德尔迭代格式为
?????-=-=+++20/3/)51()1(1)1(2)(2)1(1
k k k k x x x x ,该迭
代格式的迭代矩阵的谱半径)(M ρ= 121
。
17、 设46)2(,16)1(,0)0(===f f f ,则=)(1x l )2()(1--=x x x l ,)(x f 的二次牛顿
插值多项式为 )1(716)(2-+=x x x x N 。
18、 求积公式
?∑=≈b
a k n
k k x f A x x f )(d )(0
的代数精度以( 高斯型 )求积公式为最高,具
有( 12+n )次代数精度。
19、 已知f (1)=1,f (3)=5,f (5)=-3,用辛普生求积公式求?5
1
d )(x
x f ≈( 12 )。
20、 设f (1)=1, f (2)=2,f (3)=0,用三点式求≈')1(f ( 2.5 )。
21、如果用二分法求方程043
=-+x x 在区间]2,1[内的根精确到三位小数,需对分( 10 )
次。
23、)(,),(),(10x l x l x l n 是以整数点n x x x ,,,10 为节点的Lagrange 插值基函数,则
∑==
n
k k
x l
)(( 1 ),∑==
n
k k j
k x l
x 0
)((
j
x ),当2≥n 时=
++∑=)()3(20
4
x l x x
k k n
k k ( 32
4++x x )。
26、改变函数f x x x ()=+-1 (x >>1)的形式,使计算结果较精确
()x x x f ++=
11
。
27、若用二分法求方程()0=x f 在区间[1,2]内的根,要求精确到第3位小数,则需要对分 10
次。
29、若用复化梯形公式计算?
10
dx
e x ,要求误差不超过6
10-,利用余项公式估计,至少用 477
个求积节点。
30、写出求解方程组
??
?=+-=+2
4.016.12121x x x x 的Gauss-Seidel 迭代公式
()()
()() ,1,0,4.026.111112211=???+=-=+++k x x x x k k k k ,迭代矩阵为
????
??--64.006.10,此迭代法是否收敛 收敛 。
31、设
A =?? ?
??
5443,则=∞A 9 。
32、设矩阵
482257136A ????=??
????的A LU =,则U = 4820161002U ??
????=??
??-???? 。 33、若
4
321()f x x x =++,则差商2481632[,,,,]f = 3 。 34、数值积分公式1
12
18019()[()()()]f x dx f f f -'≈-++?的代数精度为 2 。
35、
线性方程组121015112103x ????
????????=?????
???????的最小二乘解为
11??
?
?? 。
36、设矩阵
321204135A ??
??=??
????分解为A LU =,则U = 32141003321002??
????
??-??????
?
? 。 二、单项选择题:
1、 Jacobi 迭代法解方程组b x =A 的必要条件是( C )。 A .A 的各阶顺序主子式不为零 B . 1)( 2、设 ?? ??? ?????--=700150322A ,则)(A ρ为( C ). A . 2 B . 5 C . 7 D . 3 3、三点的高斯求积公式的代数精度为( B )。 A . 2 B .5 C . 3 D . 4 4、求解线性方程组A x =b 的LU 分解法中,A 须满足的条件是( B )。 A . 对称阵 B . 正定矩阵 C . 任意阵 D . 各阶顺序主子式均不为零 5、舍入误差是( A )产生的误差。 A. 只取有限位数 B .模型准确值与用数值方法求得的准确值 C . 观察与测量 D .数学模型准确值与实际值 6、3.141580是π的有( B )位有效数字的近似值。 A . 6 B . 5 C . 4 D . 7 7、用 1+x 近似表示e x 所产生的误差是( C )误差。 A . 模型 B . 观测 C . 截断 D . 舍入 8、解线性方程组的主元素消去法中选择主元的目的是( A )。 A .控制舍入误差 B . 减小方法误差 C .防止计算时溢出 D . 简化计算 9、用1+3x 近似表示3 1x +所产生的误差是( D )误差。 A . 舍入 B . 观测 C . 模型 D . 截断 10、-324.7500是舍入得到的近似值,它有( C )位有效数字。 A . 5 B . 6 C . 7 D . 8 11、设f (-1)=1,f (0)=3,f (2)=4,则抛物插值多项式中x 2的系数为( A )。 A . –0.5 B . 0.5 C . 2 D . -2 12、三点的高斯型求积公式的代数精度为( C )。 A . 3 B . 4 C . 5 D . 2 13、( D )的3位有效数字是0.236×102。 (A) 0.0023549×103 (B) 2354.82×10-2 (C) 235.418 (D) 235.54×10-1 14、用简单迭代法求方程f(x)=0的实根,把方程f(x)=0表示成x=?(x),则f(x)=0的根是 ( B )。 (A) y=?(x)与x 轴交点的横坐标 (B) y=x 与y=?(x)交点的横坐标 (C) y=x 与x 轴的交点的横坐标 (D) y=x 与y=?(x)的交点 15、用列主元消去法解线性方程组??? ??-=+--=-+-=+-1 340921433 21321321x x x x x x x x x ,第1次消元,选择主元为 ( A ) 。 (A) -4 (B) 3 (C) 4 (D)-9 16、拉格朗日插值多项式的余项是( B ),牛顿插值多项式的余项是( C ) 。 (A) f(x,x0,x1,x2,…,xn)(x -x1)(x -x2)…(x -xn -1)(x -xn), (B) )!1() ()()()()1(+= -=+n f x P x f x R n n n ξ (C) f(x,x0,x1,x2,…,xn)(x -x0)(x -x1)(x -x2)…(x -xn -1)(x -xn), (D) ) ()!1() ()()()(1)1(x n f x P x f x R n n n n +++=-=ωξ 17、等距二点求导公式f '(x1) ≈( A )。 1011 0101 0010 101)()() D ()()() C ()()() B () ()() A (x x x f x f x x x f x f x x x f x f x x x f x f +--+---- 18、用牛顿切线法解方程f(x)=0,选初始值x0满足( A ),则它的解数列{xn}n=0,1,2,… 一定收敛到方程f(x)=0的根。 )()()D (0 )()()C (0 )()()B (0 )()()A (0000<'<''>'>''x f x f x f x f x f x f x f x f 19、为求方程x3―x2―1=0在区间[1.3,1.6]内的一个根,把方程改写成下列形式,并建 立相应的迭代公式,迭代公式不收敛的是(A )。 (A) 1 1:,1 1 12-=-=+k k x x x x 迭代公式 (B)21211:,11k k x x x x +=+ =+迭代公式 (C) 3 /12123) 1(:,1k k x x x x +=+=+迭代公式 (D) 11:,12 2 1 2 3+++==-+k k k k x x x x x x 迭代公式 21、解方程组b Ax =的简单迭代格式g Bx x k k +=+)()1(收敛的充要条件是( )。 (1)1)(A ρ, (4) 1)(>B ρ 22、在牛顿-柯特斯求积公式: ? ∑=-≈b a n i i n i x f C a b dx x f 0 )() ()()(中,当系数)(n i C 是负值时,公式的 稳定性不能保证,所以实际应用中,当( )时的牛顿-柯特斯求积公式不使用。 (1)8≥n , (2)7≥n , (3)10≥n , (4)6≥n , (1)二次; (2)三次; (3)四次; (4)五次 251732.≈计算4 1)x =,下列方法中哪种最好?( ) (A)28-; (B)24(-; (C ) ; (D) (A); (B)4; (C) ; (D ) 2。 28、形如112233()()()() b a f x dx A f x A f x A f x ≈++? 的高斯(Gauss )型求积公式的代数精度为 ( ) (A)9; (B)7; (C ) 5; (D) 3。 29的Newton 迭代格式为( ) (A) 132k k k x x x += +;(B )1322k k k x x x +=+;(C) 122k k k x x x +=+;(D) 133k k k x x x +=+。 30、用二分法求方程32 4100x x +-=在区间12[,]内的实根,要求误差限为3 1102ε-=?,则对分 次数至少为( ) (A )10; (B)12; (C)8; (D)9。 32、设()i l x 是以019(,,,)k x k k ==为节点的Lagrange 插值基函数,则9 ()i k kl k == ∑( ) (A)x ; (B )k ; (C )i ; (D )1。 33、5个节点的牛顿-柯特斯求积公式,至少具有( )次代数精度 (A )5; (B)4; (C)6; (D)3。 35、已知方程3250x x --=在2x =附近有根,下列迭代格式中在0 2x =不收敛的是( ) (A)1k x +=; (B)1k x += (C )315k k k x x x +=--; (D) 3 1225 32k k k x x x ++=-。 (A ) 4; (B)2; (C)1; (D)3。 37、5个节点的Gauss 型求积公式的最高代数精度为( ) (A)8; (B )9; (C)10; (D)11。 三、是非题(认为正确的在后面的括弧中打√,否则打?) 1、已知观察值)210()(m i y x i i ,,,, , =,用最小二乘法求n 次拟合多项式)(x P n 时,)(x P n 的次数n 可以任意取。 ( ) 2、用1-22 x 近似表示cos x 产生舍入误差。 ( ) 3、))(() )((210120x x x x x x x x ----表示在节点x 1的二次(拉格朗日)插值基函数。 ( √ ) 4、牛顿插值多项式的优点是在计算时,高一级的插值多项式可利用前一次插值的结果。 ( √ ) 5、矩阵A =? ???? ? ?-521352113具有严格对角占优。 ( ) 四、计算题: 1、用高斯-塞德尔方法解方程组 ??? ??=++=++=++225218241124321321321x x x x x x x x x ,取T )0,0,0()0(=x ,迭代四次(要求按五位有效数字计算)。 答案:迭代格式 ??? ??? ???--=--=--=++++++)222(51) 218(41)211(41)1(2)1(1)1(3)(3)1(1)1(2) (3)(2)1(1k k k k k k k k k x x x x x x x x x 2、求A 、B 使求积公式?-+-++-≈1 1)]21()21([)]1()1([)(f f B f f A dx x f 的代数精度尽量 高,并求其代数精度;利用此公式求 ? =2 1 1dx x I (保留四位小数)。 答案:2 ,,1)(x x x f =是精确成立,即 ??? ??=+=+32212222B A B A 得98,91==B A 求积公式为)]21 ()21([98)]1()1([91)(1 1f f f f dx x f +-++-=?- 当3 )(x x f =时,公式显然精确成立;当4 )(x x f =时,左=52,右=31 。所以代 数精度为3。 69286.0140 97 ] 3 21132/11[98]311311[9131111322 1 ≈= +++-++++-≈+=??--=dt t dx x x t 3、已知 分别用拉格朗日插值法和牛顿插值法求)(x f 的三次插值多项式)(3x P ,并求)2(f 的近似值(保留四位小数)。 答案: )53)(43)(13() 5)(4)(1(6 )51)(41)(31()5)(4)(3(2 )(3------+------=x x x x x x x L )45)(35)(15() 4)(3)(1(4 )54)(34)(14()5)(3)(1(5 ------+------+x x x x x x 差商表为 ) 4)(3)(1(41 )3)(1()1(22)()(33---+----+==x x x x x x x N x P 5.5)2()2(3=≈P f 6、已知x sin 区间[0.4,0.8]的函数表 如用二次插值求63891.0sin 的近似值,如何选择节点才能使误差最小?并求该近似值。 答案:解: 应选三个节点,使误差 |)(|!3|)(|33 2x M x R ω≤ 尽量小,即应使|)(|3x ω尽量小,最靠近插值点的三个节点满足上述要求。即取节点 }7.0,6.0,5.0{最好,实际计算结果 596274.063891.0sin ≈, 且 4 1055032.0)7.063891.0)(6.0963891.0)(5.063891.0(! 31 596274 .063891.0sin -?≤----≤ - 7、构造求解方程0210=-+x e x 的根的迭代格式 ,2,1,0),(1==+n x x n n ?,讨论其收敛 性,并将根求出来,4 110||-+<-n n x x 。 答案:解:令 010)1(, 02)0(,210e )(>+=<-=-+=e f f x x f x . 且010e )(>+='x x f )(∞+-∞∈?, 对x ,故0)(=x f 在(0,1)内有唯一实根.将方程0)(=x f 变形为 )e 2(101 x x -= 则当)1,0(∈x 时 )e 2(101 )(x x -= ?, 1 10 e 10e |)(|<≤-='x x ? 故迭代格式 )e 2(101 1n x n x -= + 收敛。取5.00=x ,计算结果列表如下: 且满足 6671095000000.0||-<≤-x x .所以008525090.0*≈x . 8﹑利用矩阵的LU 分解法解方程组 ? ?? ??=++=++=++2053182521432321321321x x x x x x x x x 。 答案:解: ?? ????????--??????????-==244132 11531 21LU A 令b y =L 得T )72,10,14(--=y ,y x =U 得T )3,2,1(=x . 9﹑对方程组 ??? ??=-+=--=++8 41025410151023321321321x x x x x x x x x (1) 试建立一种收敛的Seidel 迭代公式,说明理由; (2) 取初值T )0,0,0()0(=x ,利用(1)中建立的迭代公式求解,要求 3)()1(10||||-∞+<-k k x x 。 解:调整方程组的位置,使系数矩阵严格对角占优 ??? ??=++=-+=--15 1023841025410321321321x x x x x x x x x 故对应的高斯—塞德尔迭代法收敛.迭代格式为 ??? ??? ???+--=++-=++=++++++)1523(101)842(101)54(101)1(2)1(1)1(3) (3)1(1)1(2 ) (3)(2)1(1k k k k k k k k k x x x x x x x x x 取T )0,0,0() 0(=x ,经7步迭代可得: T )010000.1,326950999.0,459991999.0()7(*=≈x x . 10、已知下列实验数据 试按最小二乘原理求一次多项式拟合以上数据。 解:当0 ,则 e )(≤''x f ,且x x d e 1 0?有一位整数. 要求近似值有5位有效数字,只须误差 4) (11021 )(-?≤ f R n . 由 )(12)()( 2 3 ) (1ξf n a b f R n ''-≤,只要 4 22) (1102112e 12e ) e (-?≤≤≤n n R x n ξ 即可,解得 ???=?≥ 30877.67106e 2n 所以 68=n ,因此至少需将 [0,1] 68等份。 11、用列主元素消元法求解方程组 ??????????--=??????????????????? ?--11124112345111321x x x 。 解: ??? ?? ?????----???→??????????? ?----111124111123451111212345411121r r ?? ????? ???????? ?-----???→?????????? ??? ?????------???→?-585 25 10 57951513 012345579515 130585251 0123455 2 51 321312r r r r r r ?? ???? ? ? ??????? ?----?? ?→?+135 1350579515 13 0123 45131 23r r 回代得 3,6,1123==-=x x x 。 12、取节点1,5.0,0210===x x x ,求函数x x f -=e )(在区间[0,1]上的二次插值多项式 )(2x P ,并估计误差。 解: )15.0)(05.0() 1)(0()10)(5.00()1)(5.0()(5.002----? +----? =--x x e x x e x P )5.0(2)1(4)1)(5.0(2) 5.01)(01() 5.0)(0(15.01-+----=----? +---x x e x x e x x x x e 又 1 |)(|max ,)(,)(] 1,0[3='''=-='''=∈--x f M e x f e x f x x x 故截断误差 |)1)(5.0(|!31 |)(||)(|22--≤ -=-x x x x P e x R x 。 14、给定方程 01e )1()(=--=x x x f 1) 分析该方程存在几个根; 2) 用迭代法求出这些根,精确到5位有效数字; 3) 说明所用的迭代格式是收敛的。 解:1)将方程 01e )1(=--x x (1) 改写为 x x -=-e 1 (2) 作函数1)(1-=x x f ,x x f -=e )(2的图形(略)知(2)有唯一根)2,1(*∈x 。 2) 将方程(2)改写为 x x -+=e 1 构造迭代格式 ?? ?=+=-+5.1e 101x x k x k ),2,1,0( =k 计算结果列表如下: 3) x x -+=e 1)(?,x x --='e )(? 当]2,1[∈x 时,]2,1[)]1(),2([)(?∈???x ,且 1e |)(|1<≤'-x ? 所以迭代格式 ),2,1,0()(1 ==+k x x k k ?对任意]2,1[0∈x 均收敛。 15、用牛顿(切线)法求3的近似值。取x 0=1.7, 计算三次,保留五位小数。 解:3是03)(2 =-=x x f 的正根,x x f 2)(=',牛顿迭代公式为 n n n n x x x x 23 2 1 -- =+, 即 ) ,2,1,0(2321 =+=+n x x x n n n 取x 0=1.7, 列表如下: 16、已知f (-1)=2,f (1)=3,f (2)=-4,求拉格朗日插值多项式)(2x L 及f (1,5)的近似值,取五位小数。 解: )12)(12() 1)(1(4)21)(11()2)(1(3)21)(11()2)(1(2)(2-+-+? --+-+?+------? =x x x x x x x L )1)(1(34 )2)(1(23)2)(1(32-+--+---= x x x x x x 04167 .0241 )5.1()5.1(2≈=≈L f 17、n =3,用复合梯形公式求x x d e 10?的近似值(取四位小数),并求误差估计。 解: 7342.1]e )e e (2e [3201d e 1210 310 ≈+++?-= ≈?T x x x x x f x f e )(,e )(=''=,10≤≤x 时,e |)(|≤''x f 05.0025.0108e 312e |e |||2 3≤==?≤ -= T R x 至少有两位有效数字。 18、用Gauss-Seidel 迭代法求解线性方程组 ????? ? ?--411131103????? ??321x x x =? ???? ??--815, 取x (0)=(0,0,0)T ,列表计算三次,保留三位小数。 解:Gauss-Seidel 迭代格式为: 数值分析试题 一、 填空题(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、(本题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设,求解方程组,求雅可比迭代法与高斯-赛德尔迭代法收敛的充要条件。 解雅可比法的迭代矩阵 , 故雅可比法收敛的充要条件是。 高斯-赛德尔法的迭代矩阵 , 故高斯-赛德尔法收敛的充要条件是。 9 设求解方程组的雅可比迭代格式为,其中,求证:若,则相应的高斯-赛德尔法收敛。证明由于是雅可比法的迭代矩阵,故 又,故, 即,故故系数矩阵A按行严格对角占优,从而高斯-赛德尔法收敛。 10设A为对称正定矩阵,考虑迭代格式 求证:(1)对任意初始向量,收敛; (2)收敛到的解。 证明(1)所给格式可化为 这里存在是因为,由A对称正定,,故也对称正定。 设迭代矩阵的特征值为,为相应的特征向量,则与做内积,有 因正定,故,从而,格式收敛。 第四章 数值积分与数值微分 1.确定下列求积公式中的特定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度: 101210121 12120 (1)()()(0)(); (2)()()(0)(); (3)()[(1)2()3()]/3; (4)()[(0)()]/2[(0)()]; h h h h h f x dx A f h A f A f h f x dx A f h A f A f h f x dx f f x f x f x dx h f f h ah f f h -----≈-++≈-++≈-++''≈++-?? ?? 解: 求解求积公式的代数精度时,应根据代数精度的定义,即求积公式对于次数不超过m 的多项式均能准确地成立,但对于m+1次多项式就不准确成立,进行验证性求解。 (1)若101(1) ()()(0)()h h f x dx A f h A f A f h --≈-++? 令()1f x =,则 1012h A A A -=++ 令()f x x =,则 110A h Ah -=-+ 令2 ()f x x =,则 3 221123 h h A h A -=+ 从而解得 01 1431313A h A h A h -?=?? ?=?? ?=?? 令3 ()f x x =,则 3()0h h h h f x dx x dx --==? ? 101()(0)()0A f h A f A f h --++= 故 101()()(0)()h h f x dx A f h A f A f h --=-++? 成立。 令4 ()f x x =,则 数值分析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求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.确定下列求积公式中的特定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度: 101210121 12120 (1)()()(0)(); (2)()()(0)(); (3)()[(1)2()3()]/3; (4)()[(0)()]/2[(0)()]; h h h h h f x dx A f h A f A f h f x dx A f h A f A f h f x dx f f x f x f x dx h f f h ah f f h -----≈-++≈-++≈-++''≈++-?? ?? 解: 求解求积公式的代数精度时,应根据代数精度的定义,即求积公式对于次数不超过m 的多项式均能准确地成立,但对于m+1次多项式就不准确成立,进行验证性求解。 (1)若101(1) ()()(0)()h h f x dx A f h A f A f h --≈-++? 令()1f x =,则 1012h A A A -=++ 令()f x x =,则 110A h Ah -=-+ 令2 ()f x x =,则 3 221123 h h A h A -=+ 从而解得 011431313A h A h A h -?=?? ? =?? ?=?? 令3 ()f x x =,则 3()0h h h h f x dx x dx --==? ? 101()(0)()0A f h A f A f h --++= 令4()f x x =,则 455 1012()5 2 ()(0)()3 h h h h f x dx x dx h A f h A f A f h h ---== -++=? ? 故此时, 101()()(0)()h h f x dx A f h A f A f h --≠-++? 故 101()()(0)()h h f x dx A f h A f A f h --≈-++? 具有3次代数精度。 (2)若 21012()()(0)()h h f x dx A f h A f A f h --≈-++? 令()1f x =,则 1014h A A A -=++ 令()f x x =,则 110A h Ah -=-+ 令2 ()f x x =,则 3 2211163 h h A h A -=+ 从而解得 1143 8383A h A h A h -?=-?? ? =?? ?=?? 令3 ()f x x =,则 22322()0h h h h f x dx x dx --==? ? 101()(0)()0A f h A f A f h --++= 一、单项选择题(每小题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 --=?+?=-+--% 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.设0x >,x 的相对误差为δ,求ln x 的误差。 解:近似值*x 的相对误差为* **** r e x x e x x δ-=== 而ln x 的误差为()1ln *ln *ln **e x x x e x =-≈ 进而有(ln *)x εδ≈ 2.设x 的相对误差为2%,求n x 的相对误差。 解:设()n f x x =,则函数的条件数为'()||() p xf x C f x = 又1'()n f x nx -= , 1 ||n p x nx C n n -?∴== 又((*))(*)r p r x n C x εε≈? 且(*)r e x 为2 ((*))0.02n r x n ε∴≈ 3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指 出它们是几位有效数字:*1 1.1021x =,*20.031x =, *3385.6x =, *456.430x =,*5 7 1.0.x =? 解:*1 1.1021x =是五位有效数字; *20.031x =是二位有效数字; *3385.6x =是四位有效数字; *456.430x =是五位有效数字; *57 1.0.x =?是二位有效数字。 4.利用公式(2.3)求下列各近似值的误差限:(1) *** 124x x x ++,(2) ***123x x x ,(3) **24 /x x . 其中****1234,,,x x x x 均为第3题所给的数。 解: *4 1*3 2*13*3 4*1 51 ()102 1()102 1()102 1()102 1()102x x x x x εεεεε-----=?=?=?=?=? ***124***1244333 (1)() ()()() 111101010222 1.0510x x x x x x εεεε----++=++=?+?+?=? ***123*********123231132143 (2)() ()()() 1111.10210.031100.031385.610 1.1021385.610222 0.215 x x x x x x x x x x x x εεεε---=++=???+???+???≈ **24****24422 *4 33 5 (3)(/)()() 110.0311056.430102256.43056.430 10x x x x x x x εεε---+≈??+??=?= 5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少? 解:球体体积为343V R π= 则何种函数的条件数为 2 3 '4343 p R V R R C V R ππ=== (*)(*)3(*)r p r r V C R R εεε∴≈= 又(*)1r V ε= 数值分析整理版试题及答案 例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. 设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. 利用公式求下列各近似值的误差限: ********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 .(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字. 8. 当N 充分大时,怎样求 2 11N dx x +∞ +? ? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 212S gt = 假定g 是准确的,而对t 的测量有±秒的误差,证明当t 增加时S 的绝对误 差增加,而相对误差却减小. 11. 序列 {}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字), 计算到 10y 时误差有多大?这个计算过程稳定吗? 12. 计算61)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式 ln(ln(x x =- 计算,求对数时误差有多大? 第一章 绪论(12) 1、设0>x ,x 的相对误差为δ,求x ln 的误差。 [解]设0*>x 为x 的近似值,则有相对误差为δε=)(*x r ,绝对误差为**)(x x δε=,从而x ln 的误差为δδεε=='=* ****1)()(ln )(ln x x x x x , 相对误差为* * ** ln ln ) (ln )(ln x x x x r δ εε= = 。 2、设x 的相对误差为2%,求n x 的相对误差。 [解]设*x 为x 的近似值,则有相对误差为%2)(*=x r ε,绝对误差为**%2)(x x =ε,从而n x 的误差为n n x x n x n x x n x x x ** 1 *** %2%2) ()()()(ln * ?=='=-=εε, 相对误差为%2) () (ln )(ln *** n x x x n r == εε。 3、下列各数都是经过四舍五入得到的近似数,即误差不超过最后一位的半个单位,试指出它们是几位有效数字: 1021.1*1=x ,031.0*2=x ,6.385*3=x ,430.56*4=x ,0.17*5 ?=x 。 [解]1021.1*1 =x 有5位有效数字;0031.0* 2=x 有2位有效数字;6.385*3=x 有4位有效数字;430.56* 4 =x 有5位有效数字;0.17*5?=x 有2位有效数字。 4、利用公式(3.3)求下列各近似值的误差限,其中* 4*3*2*1,,,x x x x 均为第3题所给 的数。 (1)* 4*2*1x x x ++; [解]3 334* 4*2*11** *4*2*1*1005.1102 1 10211021)()()()()(----=?=?+?+?=++=? ??? ????=++∑x x x x x f x x x e n k k k εεεε; (2)* 3*2 *1x x x ; ()()()()()()()()()收敛较慢 代入上式得:将解: 收敛速度次并分析该迭代公式的迭代的根求方程 取试用迭代公式∴≠<<*'*+++-='∴+*+*=*∴=+?+?? ? ??===++= =∴++= ==-++=++=++014.01022220||10 2202613381013202132020 132010212010220. 2.0 20102110220 4.1222 222212012123021x x x x x x x x x x x x x x x x x x x x x x x x k k k k k k k ?????? )))()()()[]()()[])49998.0cos 215.0cos 2 1,022,00cos 2 102 12,0210,2,0.cos 2 10sin 2 11,cos 2 113cos 2 12; 1.0cos 2 12.4120101==== ==->-=<-=-=>+='-===-+x x x x x x x f f x x x f x x f x x x f x x x x k k 则 取上有一个根在所以上在为单调递增函数故则令解: 位有效数字求出这些根,精确到用迭代公式分析该方程有几个根给定方程ππππ 500 .0105.0102.0||3412≈*?=---x x x 所以方程的根 41444444466666.6663.4k k S S S S s +=+=++++++=+故迭代公式为可知: 由解: 动点迭代公式:导出下列连根公式的不Λ ΛΛΛ()()()()()()()()()()()()))()))() )()?得到的是什么迭代公式步迭代时选取第?得到的是什么迭代公式选取使收敛速度快; 选取的单根附近收敛; ,使迭代在选取值写出迭代公式是参数其中的迭代公式 给定方程不收敛 解: 都不收敛于迭代则对任何初值都有数证明:如果对于任何实为一实数设k k k k k k k k k k k k k k k k k k k x f k x f x f x f x x x x x x x G x G x x x x x G x G x x x G x G x G x Gx x x x x x G x x x G x '='==-=∴*-≥≥-≥-=*-∴-≥-∴≥--='*=*≠≥'**=*+++++++++1514302. 1. , 1.45.41 ,,1.,4.40101111111100λλλλλλΛΛ 《数值计算方法》复习试题 一、填空题: 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.答案 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) 插值基函数分别为 ()()()()()()()()()() 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 ---=== -----= ==----=== --- 故所求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 1 200110 1 1 2011000 1 210 1 ,11, ,3 1 23 ,,, ,3226 9,324 dx x dx xdx f x x dx f x x x dx ??????????==== ====++= =++= ????? 所以,法方程为 01123126119234a 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 ?=,这样,有 数值分析复习试题 第一章 绪论 一. 填空题 1.* x 为精确值 x 的近似值;() **x f y =为一元函数 ()x f y =1的近似值; ()**,*y x f y =为二元函数()y x f y ,2=的近似值,请写出下面的公式:**e x x =-: *** r x x e x -= ()()()*'1**y f x x εε≈? ()() () ()'***1**r r x f x y x f x εε≈ ? ()()()() ()* *,**,*2**f x y f x y y x y x y εεε??≈?+??? ()()()()() ** * *,***,**222r f x y e x f x y e y y x y y y ε??≈ ?+??? 2、 计算方法实际计算时,对数据只能取有限位表示,这时所产生的误差叫 舍入误 差 。 3、 分别用2.718281,2.718282作数e 的近似值,则其有效数字分别有 6 位和 7 位;又取 1.73≈-21 1.73 10 2 ≤?。 4、 设121.216, 3.654x x ==均具有3位有效数字,则12x x 的相对误差限为 0.0055 。 5、 设121.216, 3.654x x ==均具有3位有效数字,则12x x +的误差限为 0.01 。 6、 已知近似值 2.4560A x =是由真值T x 经四舍五入得 到,则相对误差限为 0.0000204 . 7、 递推公式,??? ? ?0n n-1y =y =10y -1,n =1,2, 如果取0 1.41y ≈作计算,则计算到10y 时,误 差为 81 10 2 ?;这个计算公式数值稳定不稳定 不稳定 . 8、 精确值 14159265.3* =π,则近似值141.3*1=π和1415.3*2=π分别有 3 第一章 绪论 p19 2.设x 的相对误差为2%,求n x 的相对误差。 解:设()n f x x =,则函数的条件数为'() | |() p xf x C f x = 又 1 '()n f x nx -=, 1 ||n p x nx C n n -?∴== 又 ((*))(*)r p r x n C x εε≈? 且(*)r e x 为2% ((*))0.02n r x n ε∴≈ 5计算球体积要使相对误差限为1,问度量半径R 时允许的相对误差限是多少? 解:球体体积为343 V R π= 则何种函数的条件数为 2 3'4343 p R V R R C V R ππ=== (*)(*)3(*)r p r r V C R R εεε∴≈= 又 (*)1r V ε= 故度量半径R 时允许的相对误差限为1 (*)10.333 r R ε= ?≈ 7.求方程2 5610x x -+=的两个根,使它至少具有427.982 =)。 解:2 5610x x -+= , 故方程的根应为1,228x =故 128 2827.98255.982x = ≈+= 1x ∴具有5位有效数字 211 280.0178632827.98255.982 x =-= ≈ =≈+ 2x 具有5位有效数字 9.正方形的边长大约为了100cm ,应怎样测量才能使其面积误差不超过2 1cm ? 解:正方形的面积函数为2 ()A x x = p7 当*100x =时,若(*)1A ε≤, 则21 (*)102 x ε-≤ ? 故测量中边长误差限不超过0.005cm 时,才能使其面积误差不超过2 1cm 第二章 插值法p48 1.当1,1,2 x =-时,()0,3,4f x =-, 分别用单项式基底、拉格朗日基底、牛顿基底求() f x 的二次插值多项式。 解: 0120121200102021101201220211,1,2, ()0,()3,()4;()()1 ()(1)(2)()()2()()1 ()(1)(2) ()()6 ()()1 ()(1)(1) ()()3 x x x f x f x f x x x x x l x x x x x x x x x x x l x x x x x x x x x x x l x x x x x x x ==-===-=--==-+-----==------= =-+-- 则二次拉格朗日插值多项式为 2 20 ()()k k k L x y l x ==∑ 0223()4() 14 (1)(2)(1)(1)23537623 l x l x x x x x x x =-+=---+-+=+- 2.给出()ln f x x =的数值表 用线性插值及二次插值计算的近似值。 解:由表格知,数值分析试题及答案汇总
数值分析试卷及其答案2
数值分析试卷及答案
数值分析第4章答案
最新数值分析历年考题
数值分析试卷及答案
数值分析第四章数值积分与数值微分习题答案
数值分析试题及答案
数值分析试卷及其答案
数值分析第五版全答案chap1
数值分析整理版试题及答案
数值分析第四版习题及答案
最新数值分析课程第五版课后习题答案(李庆扬等)1
数值计算第四章课后习题答案
数值计算方法试题集及答案要点
数值分析试题及答案
数值分析试卷及其答案1
数值分析整理版试题及答案
数值分析计算方法试题集及答案
数值分析第五版答案