习 题 三 解 答
1、用高斯消元法解下列方程组。
(1)1231231
22314254
27x x x x x x x x -+=??
++=??+=?①②③
解:?4②+(-)①2,1
2
?③+(-)①消去第二、三个方程的1x ,得:
1232323231425313222
x x x x x x x ?
?-+=?
-=???-=?④⑤⑥ 再由5
2)4
?⑥+(-⑤消去此方程组的第三个方程的2x ,得到三角方程组:
1232332314272184x x x x x x ?
?-+=?
-=???-=
?
回代,得:
36x =-,21x =-,19x = 所以方程组的解为
(9,1,6)T x =--
注意:
①算法要求,不能化简。化简则不是严格意义上的消元法,在算法设计上就多出了步骤。实际上,由于数值计算时用小数进行的,化简既是不必要的也是不能实现的。无论是顺序消元法还是选主元素消元法都是这样。
②消元法要求采用一般形式,或者说是分量形式,不能用矩阵,以展示消元过程。
要通过练习熟悉消元的过程而不是矩阵变换的技术。 矩阵形式错一点就是全错,也不利于检查。 一般形式或分量形式:
1231231
2231425427x x x x x x x x -+=??
++=??+=?①②③
矩阵形式
123213142541207x x x -?????? ??? ?= ??? ? ??? ???????
向量形式
123213142541207x x x -???????? ? ? ? ?++= ? ? ? ? ? ? ? ?????????
③必须是方程组到方程组的变形。三元方程组的消元过程要有三个方程组,不能变形出单一的方程。
④消元顺序12x x →→L ,不能颠倒。按为支援在方程组中的排列顺序消元也是存储算法的要求。实际上,不按顺序消元是不规范的选主元素。 ⑤不能化简方程,否则系数矩阵会变化,也不利于算法设计。
(2)12312312311323
23110
221x x x x x x x x x --=??
-++=??++=-?
①②③
解:?23②+(
)①11,1
11
?③+(-)①消去第二、三个方程的1x ,得: 123232311323523569111111252414111111x x x x x x x ?
--=??
?
-=??
?
+=-??④⑤⑥ 再由25
11)5211
?⑥+(-⑤消去此方程组的第三个方程的2x ,得到三角方程组:
123233113235235691111111932235252x x x x x x ?
?--=?
?
-=?
?
?
=-??
回代,得:
32122310641
,,193193193
x x x =-
==, 所以方程组的解为 41106223(,,)193193193T
x =-
2、将矩阵
102001112011001
1A ??
?
?= ?-
???
作LU 分解。 解:设
111213
1421222324313231
324142434410001020100001111000201110
000011u u u u l u u u A LU l l u u l l l u ?????? ??? ? ???
?=== ??? ?- ??? ?
????
??
根据矩阵乘法,先求U 的第一行,由11j j a u =,得
111213141,0,2,0u u u u ====。
再求L 的第一列,由矩阵乘法,因为1111i i a l u =,所以1
111
i i a l u =,而111u =,所以11i i l a =,所以2131410,2,0l l l ===。 再求U 的第二行,得 21122211l u u ?+?=,则
22211211001u l u =-?=-?=, 21132333101l u u u ?+?+?=,则 23211311021u l u =-?=-?=, 21142434441001l u u u u ?+?+?+?=,则 24211411001u l u =-?=-?=,
再求L 的第二列,得
3112322210000l u l u ?+?+?+?=,则 32311200200l l u =-?=-?=
41124222430000l u l u l ?+?+?+?=,则 42411200000l l u =-?=-?=
再求U 的第三行,得
311332233311l u l u u ?+?+?=-,则
33311332231122015u l u l u =--?-?=--?+?=- 311432243444101l u l u u u ?+?+?+?=,则 34311432241120011u l u l u =-?-?=-?-?= 再求L 的第三列,得
411342234333101l u l u l u ?+?+?+?=,则
4311(10201)55l =-?-?+?=-
再求U 的第四行,得
4114422443344411l u l u l u u ?+?+?+?=,则
4441144224433416
110001(1)55
u l u l u l u =-?+?+?=-?-?--?=
所以,矩阵A 的LU 分解为:
1000102001000111,201000511600100055L U ???? ? ? ? ?== ? ?- ? ?
? ?- ? ????
?
指出:
用分数而表示元素,不能化成近似小数也不化成小数表示。 3、用LU 分解紧凑格式分解法解方程组。 123457910168109171087157651x x x x ?????? ? ? ? ? ? ?= ? ? ? ? ? ??????? 解一,用一般格式求解: 将系数矩阵作LU 分解得:
100057
91062410003555,71171000552213
000101105L U ???? ? ? ? ?--- ? ?
? ?==- ?-- ? ? ? ? ? ? ?
????
Ly=b 方程组为 12341000611005171
1105213
1015y y y y ?? ????? ? ? ? ? ? ? ?= ? ?-
? ? ? ????? ? ???
解之得
123411512310y y y y ?? ??? ?- ? ? ? ?= ?- ? ? ??? ?
???
同样地,解方程组Ux=y 得 1234201253x x x x ???? ? ?- ? ?= ? ?- ? ?????。
解二,用LU 紧凑格式分解法求解: 对增广矩阵三角分解:
57910157
910
1579101662
418109136810915555
5710871771
108718715
5257651176511065157
9101579101624163555557117155222310
515???? ? ???
? ?----
? ? ? ?→→ ? ?
?- ? ?
?
? ?
?? ? ????
?
?? ? ?----- ? ?→→----
? ? ? ???24135557117
15522231
31051010?? ? ?--- ? ?---- ? ? ? ??
?
原方程组化成同解的上三角方程组为: 1234234344579101241
3555171522131010x x x x x x x x x x +++=??
?---=-???--=-??
?=
??
回代得
(20,12,5,3)T x =--。 指出:
紧凑格式是直接应用公式进行计算,计算结果保存在A 的相应元素位置。从算法的角度,紧凑格式实际体现在数据的存储方法上。
由于紧凑格式计算时不再需要A 的前面的元素,因此可以进行。
4、 用列主元的三角分解法解线性方程组。
1231231
232213472320
x x x x x x x x x -+-=-??
-+=??--=? 解一,列选主元素消元法:
先选第一列主元为213a =,将第一个方程与第二个方程交换,消去1x 得:
123232334752433371414
333x x x x x x x ?
?-+=?
?
-=?
?
?--=-?
?
再选第二列主元为3273
a =-,交换第二、三两个方程,消去2x 得三角形方程组:
12323334771414333126
33x x x x x x ?
?-+=?
?
--=-??
?-=-?
?
回代求得方程组的解312x =,21x =,12x =
所以方程组的解为
1
(2,1,)2
T x =。
解二,列主元素三角分解法:
21233247122132471(,)324712212213232023202
32033247324722714
143203333
31152212133
7A b r r r r ?? ?
-----???? ?
? ? ?=-?---→--- ? ? ?
? ?---- ????? ?--?????
? ? -- ? ? ?--→--- ? ? ? -----
--????u u u u u u u r u u u u u u u r 32472714
14333
3154
23
7?
?? ?-? ?? ?→---? ?? ?? ?--
--??
同解的三角形方程组为
1232333477141433342x x x x x x -+=??
?
--=-??
?-=-?
回代求得方程组的解312x =,21x =,12x =
所以方程组的解为
1
(2,1,)2
T x =。
说明:
用矩阵讨论中,矩阵元素进行了化简。 5.用追赶法解方程组
2111210,121
01210120A b -???? ? ?-- ? ?
? ?==-- ? ?
-- ? ?
? ?-????。 分析: 三对角矩阵
112
21n n n A αβγαβγα-??
? ?= ?
? ???
O O O 可以分解如下形式的两个矩阵:
11
22
2
331111
,1n n
n u l u L l U u l u βββ-????
?
? ? ? ? ?==
? ?
?
? ? ????
?
O
O O
O
。 即
111122
2
2233
111
111n n n n n
n u l u l u l u βαββγαββγα--????
??
??
?
? ??? ? ???= ?
??? ?
??? ??
?
????
??
?
O O
O O O O O
由矩阵乘法规则,有
111
1(2,3,)(2,3,,)
i i
i i i i i u l i n u u l i n αγαβ--?=?
?
==???=-=?L L , 这样可以求出矩阵L 和U 的所有元素。 设有系数矩阵为A 的方程组:
12,(,,)T n Ax b b b b b ==L ,
这样的方程组称为三对角方程组。 三对角方程组经LU 分解分解为
,Ly b Ux y ==,
求解之
11
1,2,3,,i
i i i y b y b l y i n -=??
=-=?L , 1(),1,2,,1n n n
i
i i i i x y u x y x i n n β+=??
=-=--?L 这就是所谓追赶法。 解:由公式 11221
22213323332443444355455542,
11,22
13
2()(1)22
12332
24
2()(1)33
13443
35
2()(1)44
1455
4
46
2()(1)55
u l u u l l u u l l u u l l u u l αγαβγαβγαβγαβ==-=
=
=-=-=--?-=
-===-
=-=--?-=
-===-
=-=--?-=
-===-
=-=--?-=
由此得下三角方程组 123451111202
1033
0140415y y y y y ?? ? ?-???? ? ? ? ? ? ?- ? ? ?= ? ? ? ? ? ?- ? ? ???
?? ? ?- ???
和上三角方程组
1122334454213124
135
1465x y x y x y x y x y -?? ?
?-???? ? ? ? ? ? ?- ? ? ?= ? ? ? ? ? ?-
? ? ???
?? ?
? ?
?
?
解上三角方程组
11223344551
111112202
1103
330114
4041155y y y y y y y y y y ???? ? ? ? ?-?????? ? ? ? ? ? ? ? ? ? ?- ? ? ? ? ?=?= ? ? ? ? ? ? ? ? ? ?- ? ?
? ? ??????? ? ? ? ?- ? ??
???
代入并解上三角方程组
112233445521131122234
114335
1514466155x x x x x x x x x x -???? ? ?
? ?-?????? ? ? ? ? ? ? ? ? ? ?- ? ? ? ? ?=?= ? ? ? ? ? ? ? ? ? ?- ? ? ? ? ???
???? ? ? ? ?
? ?
????
6、用改进的Cholesky 分解法解方程组
12312312342410
21710341097
x x x x x x x x x --=??
-++=??-++=-?
解:设此方程组的系数矩阵为A ,右端向量为b ,则
4241021710,341097A b --???? ? ?=-= ? ? ? ?-????
矩阵A 是对称正定矩阵,可以进行乔累斯基分解。 设
111112131222222313233333424217104109u u u u u u u u u u u u --??????
? ???
-= ? ??? ? ???-??????
由矩阵乘法得
1112132223332,1,2,4,2,1u u u u u u ==-=-===
由
11
1122221323
3331037u y u u y u
u u y ??????
??? ?= ??? ? ??? ???????
得
1235,2,1y y y ===-,
再由
1112
13122
232333521u u u x u u x u x ??????
??? ?= ??? ? ??? ?-?
?????
得
1232,1,1x x x ===-。
7、用改进的Cholesky 分解法解方程组
123441107131081152400246x x x x -?????? ? ? ?- ? ? ?= ? ? ?--- ? ? ??????? 解:
11121314222324333444u u u u u u u P u u u ?? ?
?= ? ?
??
1120220222
1155P ??- ? ?
?
- ?
?= ? ? ? ? ?
?
?
第一章 误差 1. 试举例,说明什么是模型误差,什么是方法误差. 解: 例如,把地球近似看为一个标准球体,利用公式2 4A r π=计算其表面积,这个近似看为球体的过程产生 的误差即为模型误差. 在计算过程中,要用到π,我们利用无穷乘积公式计算π的值: 12 222...q q π=? ?? 其中 11 2,3,... n q q n +?=?? ==?? 我们取前9项的乘积作为π的近似值,得 3.141587725...π≈ 这个去掉π的无穷乘积公式中第9项后的部分产生的误差就是方法误差,也成为截断误差. 2. 按照四舍五入的原则,将下列各数舍成五位有效数字: 816.956 7 6.000 015 17.322 50 1.235 651 93.182 13 0.015 236 23 解: 816.96 6.000 0 17.323 1.235 7 93.182 0.015 236 3. 下列各数是按照四舍五入原则得到的近似数,它们各有几位有效数字? 81.897 0.008 13 6.320 05 0.180 0 解: 五位 三位 六位 四位 4. 若1/4用0.25表示,问有多少位有效数字? 解: 两位 5. 若 1.1062,0.947a b ==,是经过舍入后得到的近似值,问:,a b a b +?各有几位有效数字? 解: 已知4311 d 10,d 1022 a b --, 又0.2053210a b +=?, ()4332111 10100.551010222 d a b da db da db ----+=+≤+=?+?=?, 所以a b +有三位有效数字; 因为0.1047571410a b ?=?, ()4332111 0.94710 1.1062100.600451010222 d a b b da a db ----?=+=??+??=? 所以a b ?有三位有效数字.
2.1 用二分法求方程013=--x x 在[1, 2]的近似根,要求误差不超过3102 1-?至少要二分多少? 解:给定误差限ε=0.5×10-3,使用二分法时,误差限为 )(211*a b x x k k -≤-+ 只要取k 满足ε<-+)(2 11 a b k 即可,亦即 96678.912lg 10lg 35.0lg 12lg lg )lg(=-+-=---≥εa b k 只要取n =10. 2.3 证明方程1 -x –sin x =0 在区间[0, 1]内有一个根,使用二分法求误差不超过 0.5×10-4的根要二分多少次? 证明 令f (x )=1-x -sin x , ∵ f (0)=1>0,f (1)=-sin1<0 ∴ f (x )=1-x -sin x =0在[0,1]有根.又 f '(x )=-1-c os x<0 (x ∈[0.1]),故f (x ) 在[0,1]单调减少,所以f (x ) 在区间 [0,1]内有唯一实根. 给定误差限ε=0.5×10-4,使用二分法时,误差限为 )(211*a b x x k k -≤-+ 只要取k 满足ε<-+)(211 a b k 即可,亦即 7287.1312 lg 10lg 45.0lg 12lg lg )lg(=-+-=---≥εa b k 只要取n =14. 2.4 方程0123=--x x 在x =1.5附近有根,把方程写成四种不同的等价形式,并建立相应的迭代公式: (1)211x x +=,迭代公式2111k k x x +=+ (2)231x x +=,迭代公式3211k k x x +=+ (3)112-=x x ,迭代公式111-=+k k x x (4)13-=x x ,迭代公式131-=+k k x x 试分析每种迭代公式的收敛性,并选取一种收敛迭代公式求出具有四位有效数字的近似根。 解:(1)令211)(x x f + =,则3 2)(x x f -=',由于 159.05.112)(33<≈≤='x x f ,因而迭代收敛。 (2)令321)(x x f +=,则322)1(3 2)(-+='x x x f ,由于
【 数值计算方法试题一 一、 填空题(每空1分,共17分) 1、如果用二分法求方程043=-+x x 在区间]2,1[内的根精确到三位小数,需对分( )次。 2、迭代格式)2(2 1-+=+k k k x x x α局部收敛的充分条件是α取值在( )。 3、已知?????≤≤+-+-+-≤≤=31)1()1()1(211 0)(2 33x c x b x a x x x x S 是三次样条函数, 则 a =( ), b =( ), c =( )。 4、)(,),(),(10x l x l x l n 是以整数点n x x x ,,,10 为节点的Lagrange 插值基函数,则 ∑== n k k x l 0)(( ), ∑== n k k j k x l x 0 )(( ),当2≥n 时 = ++∑=)()3(20 4x l x x k k n k k ( )。 ; 5、设1326)(2 47+++=x x x x f 和节点,,2,1,0,2/ ==k k x k 则=],,,[10n x x x f 和=?07 f 。 6、5个节点的牛顿-柯特斯求积公式的代数精度为 ,5个节点的求积公式最高代数精度为 。 7、{}∞ =0)(k k x ?是区间]1,0[上权函数x x =)(ρ的最高项系数为1的正交多项式族,其中1)(0=x ?,则?= 1 4)(dx x x ? 。 8、给定方程组?? ?=+-=-2211 21b x ax b ax x ,a 为实数,当a 满足 ,且20<<ω时,SOR 迭代法收敛。 9、解初值问题 00 (,)()y f x y y x y '=?? =?的改进欧拉法 ??? ??++=+=++++)],(),([2),(] 0[111] 0[1n n n n n n n n n n y x f y x f h y y y x hf y y 是 阶方法。
习 题 四 解 答 1、设010,1x x ==,写出()x f x e -=的一次插值多项式1()L x ,并估计插值误差。 设插值函数为1()L x ax b =+,由插值条件,建立线性方程组为 1 01 1a b a b e -?+=???+=? 解之得11 1a e b -?=-?=? 则11()(1)1L x e x -=-+ 因为(),()x x y x e y x e --'''=-= 所以,插值余项为 (1)(2) (2)011 ()()()()() (1)! 1()()2!1 ()()()2!1 (0)(1)((0,1))2n r x f x p x f x n f x f x x x x e x x ξξπξπξξ+-=-=+= =--=--∈ 所以 01 0101 ()max max (1) 2111248x r x e x x e ξξ-≤≤≤≤-≤-=??=。 2选用合适的三次插值多项式来近似计算f 和f 。 解:设三次插值多项式为230123()f x a a x a x a x =+++,由插值条件,建立方程组为 23012323 012323 01232301 23(0.1)(0.1)(0.1)0.9950.30.30.30.995 0.70.70.70.7651.1 1.1 1.10.454 a a a a a a a a a a a a a a a a ?+?-+?-+?-=?+?+?+?=??+?+?+?=??+?+?+?=?
即 012301230123 123012312301230.10.010.0010.9950.10.010.0010.9950.30.090.0270.9950.40.080.02800.70.490.3430.7650.80.480.344 1.761.1 1.21 1.3310.454a a a a a a a a a a a a a a a a a a a a a a a a a a -+-=-+-=??+++=++=??? +++=++=??+++=?12301231232330.40.720.9880.3110.10.010.0010.9950.40.080.02800.320.288 1.760.384 3.831a a a a a a a a a a a a a ??????++=-? -+-=??++=??? +=? ?-=-? 解之得 01 230.416.293.489.98 a a a a =??=-?? =-??=? 则所求的三次多项式为23()0.41 6.29 3.489.98f x x x x =--+。 所以 2323 (0.2)0.41 6.290.2 3.480.29.980.20.91 (0.8)0.41 6.290.8 3.480.89.980.8 1.74f f =-?-?+?=-=-?-?+?=- 3、设(0,1,2,,)i x i n =L 是 n+1个互异节点,证明: (1)0()(0,1,2,,)n k k i i i x l x x k n ===∑L ; (2)0 ()()0(0,1,2,,)n k i i i x x l x k n =-==∑L 。 证明: (1)由拉格朗日插值定理,以x 0,x 1,x 2,…x n 为插值节点,对y=f(x)=x k 作n 次插值,插值多项式为 0()()n n i i i p x l x y ==∑, 而y i =x i k , 所以0 ()()()n n k n i i i i i i p x l x y l x x ====∑∑ 同时,插值余项 (1)(1)11 ()()()()()()0(1)!(1)! n k n k n r x x p x f x x x n n ξξππ++=-= ==++ 所以0 ()n k k i i i l x x x ==∑ 结论得证。 (2)取函数()(),0,1,2,,k f x x t k n =-=L 对此函数取节点(0,1,2,,)i x i n =L ,则对应的插值多项式为
《计算方法》习题答案 第一章 数值计算中的误差 1.什么是计算方法?(狭义解释) 答:计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算,以便在计算机上编程上机,求出问题的数值解,并对算法的收敛性、稳定性和误差进行分析、计算。 2.一个实际问题利用计算机解决所采取的五个步骤是什么? 答:一个实际问题当利用计算机来解决时,应采取以下五个步骤: 实际问题→建立数学模型→构造数值算法→编程上机→获得近似结果 4.利用秦九韶算法计算多项式4)(5 3 -+-=x x x x P 在3-=x 处的值,并编程获得解。 解:400)(2 3 4 5 -+?+-?+=x x x x x x P ,从而 所以,多项式4)(5 3 -+-=x x x x P 在3-=x 处的值223)3(-=-P 。 5.叙述误差的种类及来源。 答:误差的种类及来源有如下四个方面: (1)模型误差:数学模型是对实际问题进行抽象,忽略一些次要因素简化得到的,它是原始问题的近似,即使数学模型能求出准确解,也与实际问题的真解不同,我们把数学模型与实际问题之间存在的误差称为模型误差。 (2)观测误差:在建模和具体运算过程中所用的一些原始数据往往都是通过观测、实验得来的,由于仪器的精密性,实验手段的局限性,周围环境的变化以及人们的工作态度和能力等因素,而使数据必然带有误差,这种误差称为观测误差。 (3)截断误差:理论上的精确值往往要求用无限次的运算才能得到,而实际运算时只能用有限次运算的结果来近似,这样引起的误差称为截断误差(或方法误差)。 (4)舍入误差:在数值计算过程中还会用到一些无穷小数,而计算机受机器字长的限制,它所能表示的数据只能是一定的有限数位,需要把数据按四舍五入成一定位数的近似的有理数来代替。这样引起的误差称为舍入误差。 6.掌握绝对误差(限)和相对误差(限)的定义公式。 答:设* x 是某个量的精确值,x 是其近似值,则称差x x e -=* 为近似值x 的绝对误差(简称误差)。若存在一个正数ε使ε≤-=x x e * ,称这个数ε为近似值x 的绝对误差限(简称误差限或精度)。 把绝对误差e 与精确值* x 之比* **x x x x e e r -==称为近似值x 的相对误差,称
数值计算方法试题一 一、填空题(每空1分,共17分) 1、如果用二分法求方程在区间内的根精确到三位小数,需对分()次。 2、迭代格式局部收敛的充分条件是取值在()。 3、已知是三次样条函数,则 =( ),=(),=()。 4、是以整数点为节点的Lagrange插值基函数,则 ( ),( ),当时( )。 5、设和节点则 和。 6、5个节点的牛顿-柯特斯求积公式的代数精度为,5个节点的求积公式最高代数精度为。 7、是区间上权函数的最高项系数为1的正交多项式族,其中,则。 8、给定方程组,为实数,当满足,且时,SOR迭代法收敛。 9、解初值问题的改进欧拉法是 阶方法。 10、设,当()时,必有分解式,其中为下三角阵,当其对角线元素满足()条件时,这种分解是唯一的。 二、二、选择题(每题2分) 1、解方程组的简单迭代格式收敛的充要条件是()。(1), (2) , (3) , (4) 2、在牛顿-柯特斯求积公式:中,当系数是负值时,公式的稳定性不能保证,所以实际应用中,当()时的牛顿-柯特斯求积公式不使用。 (1),(2),(3),(4), (1)二次;(2)三次;(3)四次;(4)五次 4、若用二阶中点公式求解初值问题,试问为保证该公式绝对稳定,步长的取值范围为()。 (1), (2), (3), (4)
三、1、 2、(15 (1)(1) 试用余项估计其误差。 (2)用的复化梯形公式(或复化 Simpson公式)计算出该积分的近似值。 四、1、(15分)方程在附近有根,把方程写成三种不同的等价形式(1)对应迭代格式;(2)对应迭代格式;(3)对应迭代格式。判断迭代格式在的收敛性,选一种收敛格式计算附近的根,精确到小数点后第三位。选一种迭代格式建立Steffensen迭代法,并进行计算与前一种结果比较,说明是否有加速效果。 2、(8分)已知方程组,其中 , (1)(1)列出Jacobi迭代法和Gauss-Seidel迭代法的分量形式。 (2)(2)求出Jacobi迭代矩阵的谱半径,写出SOR 迭代法。 五、1、(15分)取步长,求解初值问题用改进的欧拉法求的值;用经典的四阶龙格—库塔法求的值。 2、(8分)求一次数不高于4次的多项式使它满足 ,,,, 六、(下列2题任选一题,4分) 1、1、数值积分公式形如 (1)(1)试确定参数使公式代数精度尽量高;(2)设,推导余项公式,并估计误差。 2、2、用二步法 求解常微分方程的初值问题时,如何选择参数使方法阶数尽可能高,并求局部截断误差主项,此时该方法是几阶的。 数值计算方法试题二 一、判断题:(共16分,每小题2分) 1、若是阶非奇异阵,则必存在单位下三角阵和上三角阵,使唯一成立。()
1、解:将)(x V n 按最后一行展开,即知)(x V n 是n 次多项式。 由于 n i i i n n n n n i n x x x x x x x x x x V ...1...1... ......... ...... 1 )(21110 20 0---= ,.1,...,1,0-=n i 故知0)(=i n x V ,即110,...,,-n x x x 是)(x V n 的根。又)(x V n 的最高 次幂 n x 的系数为 )(...1...1... ...... .........1),...,,(101 1 21 11 2 2221 02001101j n i j i n n n n n n n n n n n x x x x x x x x x x x x x x V -== ∏-≤<≤-----------。 故知).)...()()(,...,,()(1101101------=n n n n x x x x x x x x x V x V 6、解:(1)设 .)(k x x f =当n k ,...,1,0=时,有.0)()1(=+x f n 对 )(x f 构造Lagrange 插值多项式, ),()(0 x l x x L j n j k j n ∑== 其 0)()! 1() ()()()(1)1(=+=-=++x w n f x L x F x R n n n n ξ, ξ介于j x 之间,.,...,1,0n j = 故 ),()(x L x f n =即 .,...,1,0,)(0 n k x x l x k j n j k j ==∑= 特别地,当0=k 时, 10) (=∑=n j x j l 。 (2) 0)()1(1) ()1()()(0000=-=??? ? ??-??? ? ??-=--=-===∑∑∑∑k j j i j i k j k i i j i i k j n j k i i j k n j j x x x x i k x l x x i k x l x x )利用(。 7、证明:以b a ,为节点进行线性插值,得 )()()(1 b f a b a x a f b a b x x P --+--= 因 0)()(==b f a f ,故0)(1=x P 。而 ))()(("2 1 )()(1b x a x f x P x f --= -ξ,b a <<ξ。 故)("max )(8 122)("max )(max 2 2 x f a b a b x f x f b x a b x a b x a ≤≤≤≤≤≤-=??? ??-≤。 14、解:设 ))...()(()(21n n x x x x x x a x f ---=, k x x g =)(,记)() (1 ∏=-=n j j n x x x w ,则 ),()(x w a x f n n =).()(' j n n j x w a x f = 由差商的性质知 [])! 1()(1,..,,1) (' 1 )(')('1 211 11 -== ==-===∑∑∑ n g a x x x g a x w x a x w a x x f x n n n n n j j n k j n n j j n n k j n j j k j ξ, ξ介于n x x ,...,1之间。 当20-≤≤ n k 时,0)()1(=-ξn g , 当 1-=n k 时,)!1()(1-=-n g n ξ, 故 ???-=-≤≤=-= --=∑1,,20,0)!1()(1) ('1 11 n k a n k n g a x f x n n n n j j k j ξ 16、解:根据差商与微商的关系,有 [] 1! 7! 7!7)(2,...,2,2)7(7 10===ξf f , [ ] 0! 80 !8)(2,...,2,2)8(8 1 ===ξf f 。 ( 13)(47+++=x x x x f 是7次多项式, 故 ,!7)()7(=x f 0)()8(=x f )。 25、解:(1) 右边= [][]dx x S x f x S dx x S x f b a b a ??-+-)(")(")("2)(")("2 = [] d x x S x f x S x S x S x f x f b a ?-++-)("2)(")("2)(")(")("2)(" 222 = [] d x x S x f b a ?-)(")(" 22 = [][]dx x S dx x f b a b a 2 2 )(")("??- =左边。 (2)左边= ? -b a dx x S x f x S ))(")(")(("
《计算方法》练习题一 练习题第1套参考答案 一、填空题 1. 14159.3=π的近似值3.1428,准确数位是( 2 10- )。 2.满足d b f c a f ==)(,)(的插值余项=)(x R ( ))((!2) (b x a x f --''ξ ) 。 3.设)}({x P k 为勒让德多项式,则=))(),((22x P x P (5 2 )。 4.乘幂法是求实方阵(按模最大 )特征值与特征向量的迭代法。 5.欧拉法的绝对稳定实区间是( ]0,2[-)。 二、单选题 1.已知近似数,,b a 的误差限)(),(b a εε,则=)(ab ε(C )。 A .)()(b a εε B.)()(b a εε+ C.)()(b b a a εε+ D.)()(a b b a εε+ 2.设x x x f +=2 )(,则=]3,2,1[f ( A )。 A.1 B.2 C.3 D.4 3.设A=?? ? ? ??3113,则化A为对角阵的平面旋转=θ( C ) . A. 2π B.3π C.4π D.6 π 4.若双点弦法收敛,则双点弦法具有(B )敛速. A.线性 B.超线性 C.平方 D.三次 5.改进欧拉法的局部截断误差阶是( C ). A .)(h o B.)(2 h o C.)(3 h o D.)(4 h o 三、计算题 1.求矛盾方程组:??? ??=-=+=+2 42321 2121x x x x x x 的最小二乘解。 2 212 212 2121)2()42()3(),(--+-++-+=x x x x x x x x ?, 由 0,021=??=??x x ? ?得:???=+=+9 629232121x x x x , 解得14 9 ,71821== x x 。
数值计算方法习题一(2) 习题二(6) 习题三(15) 习题四(29) 习题五(37) 习题六(62) 习题七(70) 2009.9,9
习题一 1.设x >0相对误差为2%4x 的相对误差。 解:由自变量的误差对函数值引起误差的公式: (())(())'()()()() f x x f x f x x f x f x δδ?= ≈得 (1)()f x = 11 ()()*2%1% 22x x δδδ≈ ===; (2)4 ()f x x =时 44 4 ()()'()4()4*2%8%x x x x x x δδδ≈ === 2.设下面各数都是经过四舍五入得到的近似数,即误差不超过最后一位的半个单位,试指出他们各有几位有效数字。 (1)12.1x =;(2)12.10x =;(3)12.100x =。 解:由教材9P 关于1212.m n x a a a bb b =±型数的有效数字的结论,易得上面三个数的有效 数字位数分别为:3,4,5 3.用十进制四位浮点数计算 (1)31.97+2.456+0.1352; (2)31.97+(2.456+0.1352) 哪个较精确? 解:(1)31.97+2.456+0.1352 ≈2 1 ((0.3197100.245610)0.1352)fl fl ?+?+ =2 (0.3443100.1352)fl ?+ =0.3457210? (2)31.97+(2.456+0.1352) 2 1 (0.319710(0.245610))fl fl ≈?+? = 21 (0.3197100.259110)fl ?+? =0.34562 10? 易见31.97+2.456+0.1352=0.3456122 10?,故(2)的计算结果较精确。 4.计算正方形面积时,若要求面积的允许相对误差为1%,测量边长所允许的相对误差限为多少?
第一章绪论 习题一 1.设x>0,x*的相对误差为δ,求f(x)=ln x的误差限。解:求lnx的误差极限就是求f(x)=lnx的误差限,由公式(1. 2.4)有 已知x*的相对误差满足,而 ,故 即 2.下列各数都是经过四舍五入得到的近似值,试指出它们有几位有效数字,并给出其误差限与相对误差限。 解:直接根据定义和式(1.2.2)(1.2.3)则得 有5位有效数字,其误差限,相对误差限 有2位有效数字, 有5位有效数字, 3.下列公式如何才比较准确? (1) (2)
解:要使计算较准确,主要是避免两相近数相减,故应变换所给公式。 (1) (2) 4.近似数x*=0.0310,是 3 位有数数字。 5.计算取,利用:式计算误差最小。 四个选项: 第二、三章插值与函数逼近 习题二、三 1. 给定的数值表 用线性插值与二次插值计算ln0.54的近似值并估计误差限. 解:仍可使用n=1及n=2的Lagrange插值或Newton插值,并应用误差估计(5.8)。线性插值时,用0.5及0.6两点,用Newton插值 误差限,因
,故 二次插值时,用0.5,0.6,0.7三点,作二次Newton插值 误差限 ,故 2. 在-4≤x≤4上给出的等距节点函数表,若用二次插值法求的近似值,要使误差不超过,函数表的步长h 应取多少? 解:用误差估计式(5.8), 令 因 得 3. 若,求和.
解:由均差与导数关系 于是 4. 若互异,求 的值,这里p≤n+1. 解:,由均差对称性 可知当有 而当P=n+1时 于是得 5. 求证. 解:解:只要按差分定义直接展开得 6. 已知的函数表
习 题 三 解 答 1、用高斯消元法解下列方程组。 (1)1231231 22314254 27x x x x x x x x -+=?? ++=??+=?①②③ 解:?4②+(-)①2,1 2 ?③+(-)①消去第二、三个方程的1x ,得: 1232323231425313222 x x x x x x x ? ?-+=? -=???-=?④⑤⑥ 再由5 2)4 ?⑥+(-⑤消去此方程组的第三个方程的2x ,得到三角方程组: 1232332314272184x x x x x x ? ?-+=? -=???-= ? 回代,得: 36x =-,21x =-,19x = 所以方程组的解为 (9,1,6)T x =-- 注意: ①算法要求,不能化简。化简则不是严格意义上的消元法,在算法设计上就多出了步骤。实际上,由于数值计算时用小数进行的,化简既是不必要的也是不能实现的。无论是顺序消元法还是选主元素消元法都是这样。 ②消元法要求采用一般形式,或者说是分量形式,不能用矩阵,以展示消元过程。 要通过练习熟悉消元的过程而不是矩阵变换的技术。 矩阵形式错一点就是全错,也不利于检查。 一般形式或分量形式: 1231231 22314254 27x x x x x x x x -+=?? ++=??+=?①②③ 矩阵形式 123213142541207x x x -?????? ??? ?= ??? ? ??? ???????
向量形式 123213142541207x x x -???????? ? ? ? ?++= ? ? ? ? ? ? ? ????????? ③必须是方程组到方程组的变形。三元方程组的消元过程要有三个方程组,不能变形出单一的方程。 ④消元顺序12x x →→L ,不能颠倒。按为支援在方程组中的排列顺序消元也是存储算法的要求。实际上,不按顺序消元是不规范的选主元素。 ⑤不能化简方程,否则系数矩阵会变化,也不利于算法设计。 (2)1231231231132323110 221x x x x x x x x x --=?? -++=??++=-? ①②③ 解:?23②+( )①11,1 11 ?③+(-)①消去第二、三个方程的1x ,得: 123232311323523569111111252414111111x x x x x x x ? --=?? ? -=? ? ? +=-??④⑤⑥ 再由25 11)5211 ?⑥+(-⑤消去此方程组的第三个方程的2x ,得到三角方程组: 123233113235235691111111932235252x x x x x x ? ?--=? ? -=?? ? =-?? 回代,得: 32122310641 ,,193193193 x x x =- ==, 所以方程组的解为 41106223(,,)193193193T x =- 2、将矩阵 1020011120110011A ?? ? ?= ?- ???
练习题与答案 练习题一 练习题二 练习题三 练习题四 练习题五 练习题六 练习题七 练习题八 练习题答案 练习题一 一、是非题 1.*x=–1 2.0326作为x的近似值一定具有6位有效数字,且其误差限 ≤ 4 10 2 1 - ? 。( ) 2.对两个不同数的近似数,误差越小,有效数位越多。( ) 3.一个近似数的有效数位愈多,其相对误差限愈小。( )
4. 用 2 12x -近似表示cos x 产生舍入误差。 ( ) 5. 3.14和3.142作为π的近似值有效数字位数相同。 ( ) 二、填空题 1. 为了使计算()()2334912111y x x x =+-+---的乘除法次数尽量少,应将该表达式改写为 ; 2. *x =–0.003457是x 舍入得到的近似值,它有 位有效数字,误差限 为 ,相对误差限为 ; 3. 误差的来源是 ; 4. 截断误差为 ; 5. 设计算法应遵循的原则 是 。 三、选择题 1.*x =–0.026900作为x 的近似值,它的有效数字位数为( ) 。 (A) 7; (B) 3; (C) 不能确定 (D) 5. 2.舍入误差是( )产生的误差。 (A) 只取有限位数 (B) 模型准确值与用数值方法求得的准确值 (C) 观察与测量 (D) 数学模型准确值与实际值
3.用 1+x 近似表示e x 所产生的误差是( )误差。 (A). 模型 (B). 观测 (C). 截断 (D). 舍入 4.用s *=21 g t 2表示自由落体运动距离与时间的关系式 (g 为重力加速度),s t 是 在时间t 内的实际距离,则s t s *是( )误差。 (A). 舍入 (B). 观测 (C). 模型 (D). 截断 5.1.41300作为2的近似值,有( )位有效数字。 (A) 3; (B) 4; (C) 5; (D) 6。 四、计算题 1. 3.142,3.141,22 7分别作为π的近似值,各有几位有效数字? 2. 设计算球体积允许的相对误差限为1%,问测量球直径的相对误差限最大为多少? 3. 利用等价变换使下列表达式的计算结果比较精确: (1)1||,11211<<+-++x x x x , (2) 1||1112<<+?+x dt t x x (3) 1||,1<<-x e x , (4) 1)1ln(2>>-+x x x 4.真空中自由落体运动距离s 与时间t 的关系式是s =21 g t 2,g 为重力加速度。现设g 是精确的,而对t 有0.1±秒的测量误差,证明:当t 增加时,距离的绝对误差增加,而相对误差却减少。
《计算方法》期中复习试题 一、填空题: 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 d )(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 ,为了减少舍入误差,应将表达式1999 2001-
数值分析作业答案 插值法 1、当x=1,-1,2时,f(x)=0,-3,4,求f(x)的二次插值多项式。 (1)用单项式基底。 (2)用Lagrange插值基底。 (3)用Newton基底。 证明三种方法得到的多项式是相同的。 解:(1)用单项式基底 设多项式为: , 所以: 所以f(x)的二次插值多项式为: (2)用Lagrange插值基底 Lagrange插值多项式为: 所以f(x)的二次插值多项式为: (3) 用Newton基底: 均差表如下: xk f(xk) 一阶均差二阶均差 1 0 -1 -3 3/2 2 4 7/ 3 5/6 Newton插值多项式为: 所以f(x)的二次插值多项式为: 由以上计算可知,三种方法得到的多项式是相同的。 6、在上给出的等距节点函数表,若用二次插值求ex的近似值,要使截断误差不超过10-6,问使用函数表的步长h应取多少? 解:以xi-1,xi,xi+1为插值节点多项式的截断误差,则有 式中 令得 插值点个数
是奇数,故实际可采用的函数值表步长 8、,求及。 解:由均差的性质可知,均差与导数有如下关系: 所以有: 15、证明两点三次Hermite插值余项是 并由此求出分段三次Hermite插值的误差限。 证明:利用[xk,xk+1]上两点三次Hermite插值条件 知有二重零点xk和k+1。设 确定函数k(x): 当或xk+1时k(x)取任何有限值均可; 当时,,构造关于变量t的函数 显然有 在[xk,x][x,xk+1]上对g(x)使用Rolle定理,存在及使得 在,,上对使用Rolle定理,存在,和使得 再依次对和使用Rolle定理,知至少存在使得 而,将代入,得到 推导过程表明依赖于及x 综合以上过程有: 确定误差限: 记为f(x)在[a,b]上基于等距节点的分段三次Hermite插值函数。在区间[xk,xk+1]上有 而最值 进而得误差估计: 16、求一个次数不高于4次的多项式,使它满足,,。
比较详细的数值分析课后习题答案
0.1算法 1、 (p.11,题1)用二分法求方程013 =--x x 在[1,2]的近似根,要求误差不超过 10-3. 【解】 由二分法的误差估计式31 1*102 1 2||-++=≤=-≤ -εk k k a b x x ,得到100021≥+k .两端取自然对数得96.812ln 10 ln 3≈-≥ k ,因此取9=k ,即至少需 2、(p.11,题2) 证明方程210)(-+=x e x f x 在区间[0,1]有唯一个实根;使用二 分法求这一实根,要求误差不超过2102 1 -?。 【解】 由于210)(-+=x e x f x ,则)(x f 在区间[0,1]上连续,且 012010)0(0<-=-?+=e f ,082110)1(1>+=-?+=e e f ,即0)1()0(+=x e x f ,即)(x f 在区间[0,1]上是单调的,故)(x f 在区间[0,1]有唯一实根.
由二分法的误差估计式21 1*1021 2 12||-++?=≤=-≤-εk k k a b x x ,得到1002≥k .两端取自然对数得6438.63219.322 ln 10 ln 2=?≈≥ k ,因此取7=k ,即至少需二分 0.2误差 1.(p.12,题8)已知e=2.71828…,试问其近似值7.21=x ,71.22=x ,x 2=2.71, 718.23=x 各有几位有效数字?并给出它们的相对误差限。 【解】有效数字: 因为111021 05.001828.0||-?= <=- x e ,所以7.21=x 有两位有效数字; 因为1 2102105.000828.0||-?=<=- x e ,所以71.22=x 亦有两位有效数字; 因为3 3102 10005.000028.0||-?=<=- x e ,所以718.23=x 有四位有效数字; %85.17.205 .0||111=<-= x x e r ε; %85.171 .205 .0||222=<-= x x e r ε;
《数值计算方法》复习试题 一、填空题: 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的改进的欧拉公
roots([1 -1 -1]) x=linspace(0,2*pi,10); y=sin(x); xi=linspace(0,2*pi,100); y1=interp1(x,y,xi); y2=interp1(x,y,xi,'spline'); y3=interp1(x,y,xi,'cublic'); plot(x,y,'o',xi,y1,xi,y2,xi,y3) x=[0 300 600 1000 1500 2000]; y=[0.9689 0.9322 0.8969 0.8519 0.7989 0.7491]; xi=linspace(0,2000,20); yi=1.0332*exp(-(xi+500)/7756); y1=interp1(x,y,xi,'spline'); subplot(2,1,1);plot(x,y,'o',xi,yi,xi,y1,'*') p=polyfit(x,y,2); y2=polyval(p,xi); subplot(2,1,2);plot(x,y,'o',xi,yi,xi,y2,'*') x=[0 300 600 1000 1500 2000]; y=[0.9689 0.9322 0.8969 0.8519 0.7989 0.7491]; xi=linspace(0,2000,20); y1=interp1(x,y,xi,'spline'); subplot(2,1,1);plot(x,y,'-o', xi,y1,'-*') p=polyfit(x,y,2); y2=polyval(p,xi); subplot(2,1,2);plot(x,y,'-o',xi,y2,'-*')
数值分析 (p11页) 4 试证:对任给初值x 0, 0)a >的牛顿迭代公式 112(),0,1 ,2,......k a k k x x x k +=+= 恒成立下列关系式: 2112(1)(,0,1,2,.... (2)1,2,...... k k k x k x x k x k +-=≥= 证明: (1 )(2 1122k k k k k k x a x x x x +-??=+= =? ?? (2) 取初值00>x ,显然有0>k x ,对任意0≥k , a a x a x x a x x k k k k k ≥+??? ? ??-=???? ??+=+2 12121 6 证明: 若k x 有n 位有效数字,则n k x -?≤ -1102 1 8, 而() k k k k k x x x x x 28882182 1-=-???? ? ?+=-+ n n k k x x 21221102 1 5.22104185 .28--+?=??<-∴>≥ 1k x +∴必有2n 位有效数字。 8 解: 此题的相对误差限通常有两种解法. ①根据本章中所给出的定理: (设x 的近似数* x 可表示为m n a a a x 10......021*?±=,如果* x 具有l 位有效数字,则其相对误差限为 ()11 * *1021 --?≤ -l a x x x ,其中1a 为*x 中第一个非零数)
则7.21=x ,有两位有效数字,相对误差限为 025.0102 21 111=??≤--x x e 71.22=x ,有两位有效数字,相对误差限为 025.0102 21 122=??≤--x x e 3 2.718x =,有两位有效数字,其相对误差限为: 00025.0102 21 333=??≤--x e x ②第二种方法直接根据相对误差限的定义式求解 对于7.21=x ,0183.01<-e x ∴其相对误差限为00678.07 .20183.011≈<-x e x 同理对于71.22=x ,有 003063 .071 .20083 .022≈<-x e x 对于718.23=x ,有 00012.0718 .20003 .033≈<-x e x 备注:(1)两种方法均可得出相对误差限,但第一种是对于所有具有n 位有效数字的近似数都成立的正确结论,故他对误差限的估计偏大,但计算略简单些;而第二种方法给出较好的误差限估计,但计算稍复杂。 (2)采用第二种方法时,分子为绝对误差限,不是单纯的对真实值与近似值差值的四舍五入,绝对误差限大于或等于真实值与近似值的差。 11. 解: ......142857.3722≈,.......1415929.3113 255≈ 2102 1 722-?≤-∴ π,具有3位有效数字