《计算方法》习题答案
第一章 数值计算中的误差
1.什么是计算方法?(狭义解释)
答:计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算,以便在计算机上编程上机,求出问题的数值解,并对算法的收敛性、稳定性和误差进行分析、计算。
2.一个实际问题利用计算机解决所采取的五个步骤是什么?
答:一个实际问题当利用计算机来解决时,应采取以下五个步骤: 实际问题→建立数学模型→构造数值算法→编程上机→获得近似结果 4.利用秦九韶算法计算多项式4)(53-+-=x x x x P 在3-=x 处的值,并编程获得解。 解:400)(2345-+?+-?+=x x x x x x P ,从而 1 0 -1 0 1 -4 -3 -3 9 -24 72 -219
1
-3
8
-24
73
-223
所以,多项式4)(53-+-=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 的相对误差,称
*
x ε
η=
为近似值x 的相对误差限η≤r e ,由于真值*
x 是未知的,所以常常用
x
e x x x e r =-=*来表示相对误差,于是相对误差可以从绝对误差求出。
7.近似值的规格化表示形式如何?
答:一般地,对于一个精确值*
x ,其近似值x 的规格化形式为m p x x x x 10.021?±= ,
其中{}),2,1(9,2,1,0,01p i x x i =∈≠,p 为正整数,m 为整数。
8.有效数字的概念是什么?掌握有效数字与误差的关系。
答:若近似值x 的(绝对)误差限是它的某一位的半个单位,也就是说该近似值准确到这一位,且从该位起直到前面第一个非零数字为止的所有数字都称为有效数字。 若近似值x 的(绝对)误差限为n m x x e -?≤-=102
1
*
,则称x 为具有n 位有效数字的有效数,或称它精确到n
m -10
位,其中的每一位数字n x x x ,,21都是x 的有效数字。
设精确值*
x 的近似值x 的规格化形式为m
p x x x x 10.021?±= ,若x 具有n 位有效数字,则其相对误差限为n r x e -?≤
11
1021
;反之,若x 的相对误差限为n r x e -?+≤
1110)
1(21
,则x 至少有n 位有效数字。
9.下列各数都是对真值进行四舍五入后获得的近似值,试分别写出它们的绝对误差限,相对误差限和有效数字的位数。
(1)024.01=x (2)4135.02=x (3)50.573=x (4)600004=x (5)55108?=x ;
解:(1)0005.01*
1
≤-=x x e ;0021.0*≤=-=x
e x x x e r ;有三位有效数字。
(2)00005.02*
2≤-=x x e ;000121
.0*≤=-=x
e
x x x e r ;有四位有效数字。 (3)005.03*
3≤-=x x e ;000087
.0*≤=-=x
e
x x x e r ;有四位有效数字。 (4)5.04*
4≤-=x x e ;0000084
.0*≤=-=x
e
x x x e r ;有五位有效数字。
(5)5.05*
5≤-=x x e ;000000625
.0*≤=-=x
e
x x x e r ;有六位有效数字。 10.为了使19的相对误差≤0.1%,问至少应取几位有效数字?
解:由19的首位数是4.设近似数*
x 有n 位有效数字,由定理4.1可知,相对误差
001.0104
21
)(1*≤??≤
-n r x e ,解得097.3≥n ,即取4位有效数字,近似数的相对误差不超过0.1%。
11.已知33,3100,1150)(*
2
==-+==x x x x x P y ,计算)3
100
(*p y =及)33(P y =,并求x 和y 的相对误差。 解: 55555.51150)3
100
()3100()3100(
2*
-≈-+==p y 281150)33()33()33(2-=-+==P y
333
.0)(*
≈-=x x x e 0101.0)
()(≈=
x
x e x e r
44444.22)(*
≈-=y y y e 801587
.0)
()(≈=
y
y e y e r 12.写出误差估计的一般公式(以二元函数),(y x f z =为例)。 解:二元函数),(y x f z =的绝对误差: )(|)(|)(),(),(y e y
f
x e x f z e y x y x ???+???≈
二元函数的相对误差: z y e y f z x e x f z z e z e y x y x r )
(|)(|)()(),(),(???+???≈=
)(|)(|),(),(y e y
f
z y x e x f z x r y x r y x ????+????=
13.用电表测得一个电阻两端的电压和流过的电流范围分别为V V 2220±=,A I 1.010±=,求这个电阻的阻值R ,并估算其绝对误差和相对误差。
解:2)(≤V e ,1.0)(≤I e ,又2,1,I
V I R I V R I V R -=??=??=
。所以: 42.01.0100
2202101)(|)(|)(|)(|)(),(),(),(),(=?+?=???+???≤???+???≈
I e I R V e V R I e I
R
V e V R R e I V I V I V I V
21099.1)
()(-?≈=
R
R e R e r 。 14.若01.045.0,01.003.1*
2*1±=±=x x ,计算2
2
12
1x e x y +
=的近似值,并估计)(y e 及其上界。
解:45.02
2
1)03.1(e y +
≈ )(2
1))(()21()21()(2*22*21*
11*11*1*x x x x e e x x x x e x e x y y y e -++-=+-+=-=
),(,01.02
11006.2)(21))((*
2221*1
1*
12*2x x e e e x x x x x x ∈??+?=-++-≤-ξξ
15.已测得某场地长为m l 110=,宽d 的值为m d 80=,已知m l l l e 2.0)(*
≤-=,
m d d d e 1.0)(*≤-=,试求面积ld s =的绝对误差限和相对误差限。
解:由ld s =,
l d
s d l s =??=??,,m l l l e 2.0)(*≤-=,m d d d e 1.0)(*≤-=。可得:30
1.080
2.0110)
(|)(|)(|)(|)(),(),(),(),(=?+?=???+???≤???+???≈
d e d
s l e l s d e d s l e l s s e d l d l d l d l 3104.3)
()(-?≈=
s
s e s e r 。
16.掌握二元函数的加、减、乘、除和开方运算的绝对误差和相对误差估计公式。 解:(1)加、减运算: 由于()1/=?+?x y x ()()(),1/,1/,1/-=?-?=?-?=?+?y y x x y x y y x ,所以
()()()()()()()()()()()()()()()()()()()()()|
||/|||/|||,//,,//,y e y x y x e y x x y x e y e y x y x e y x x y x e y e x e y x e y e y x y x e y x x y x e y e x e y x e r r r r r r r r r ?-+?-≤-?--?-≈--≈-?++?+≈++≈+从而有
(2)乘法运算: 由于
()(),x y
xy y x xy =??=??,所以()()()()()()y e x e xy e y xe x ye r r r +≈+≈,xy e ,从而()()()|||||||||y e x x e y xy e ?+?≤
(3)除法运算: 由于2)(,1)
(y
x y y x y x y
x
-=??=
??,所以)()(1)(2y e y
x
x e y y x
e -≈
,)()()(y e x e y
x
e r r r -≈
(4)乘方及开方运算:
由于()
1-=??n n
nx x
x ,所以()()()()x ne x e x e nx x e r n r n n ≈≈-,1 17.求方程01562
=+-x x 的两个根,使它至少具有4位有效数字(982.27783≈)。
解:782.55982.27281
21
14)56(5621=+≈???--+=x
017863
.0782
.551
12≈==
x c x 19.求方程01162
=+-x x 的较小正根,要求有3位有效数字。
解:937.15937.781
21
14)16(1621=+≈???--+=x
062747.0937
.15112≈==
x c x 所以较小正根为062747.02≈x 。 20.设4
1
10
,,2,1,0, ==
?n dx e x I x n n 。
(1)证明:4
110,,2,1,0, =-=-n nI e I n n ;
(2)给出一个数值稳定的算法,并证明算法的稳定性。 (1)证明:11
11
1
---=-===???
n x n x
n x
n n nI e x d e nx e e d x dx e x I
(2))(1
1n n I e n
I -=
-
设n n n I I e -=*
,则
n n
n
n n n n n n n e n I I e e n
I I e e n I I e 1
11
0*
0022
*221*
11=
-==-==-=------
当n 无限大时,n e 越小,所以该算法稳定。
21.用递推算法计算积分?=+=1
010,2,1,0,41 n dx x
x I n
n ,并验证算法的数值稳定性。 解:1101
101101141
41)41(4141441------=+-=+-+=???n n n n n n n I n dx x x dx x dx x x x x I 设0*
00I I e -=,则
1010*
1010022*
2201*11414141e I I e e I I e e I I e =-==-==
-=
所以该算法是稳定的。 22.设计一个计算362412
163)(x x x
x f ++=的最小计算量的算法。
解:24121212442362412
163163)(x x x x x x x x x x x x
x f ??+??+????=++=
23.什么是数值稳定的算法?数值计算应遵循的六条规则是什么?
答:一个算法如果原始数据有误差(扰动),而计算过程中舍入误差不增长或增长可以控制,则称此算法是数值稳定的。否则,称此算法是数值不稳定的。
数值计算应遵循的六条规则是:
(1) 选用数值稳定的算法(计算公式); (2) 尽量避免两个相近数相减;
(3) 尽量避免用绝对值很大的数作乘数; (4) 尽量避免用绝对值很小的数作除数; (5) 防止大数“吃掉”(或“淹没”)小数(即合理安排运算顺序); (6) 简化计算步骤,减少运算次数。
第二章 非线性方程的数值解法
1.叙述零点定理的内容。
答:设函数)(x f 在闭区间],[b a 上连续且0)()(
0)(*=x f ,即)(x f 在区间),(b a 内存在实的零点,称区间],[b a 为方程的有根区间。
2.方程求根的两个步骤是什么?确定方程有根区间的方法有哪些? 答:第一步 确定方程0)(=x f 的有根区间。
第二步 近似根的精确化。
确定方程有根区间的方法有两种:作图法和逐步搜索法。 3.利用作图法确定方程01)(3=--=x x x f 的有根区间。 解
:
由于,05128)2(,01)0(>=--=<-=f f 于是,在区间)2,0(内至少有一个根,取步长5.0=h 向右进行根的搜索,即计算)5.1(),0.1(),5.0(f f f 的值得到
0)5.1(,0)0.1(,0)5.0(>< 4.利用逐步搜索法确定方程0343)(2 3 =++-=x x x x f 的有根区间。 解:由于,05)1(,03)0(<-=->=f f 于是,方程在)0,1(-内至少有一个实根,所以,从1-=x ,取步长5.0=h 向右进行根的搜索,即计算)5.0(-f 得到08 1 )5.0(>= -f ,从y 1-1 1 )(3--=x x x f 而,原方程的有根区间缩小为)2 1,1(--。 5.确定方程01042 3 =-+x x 的有根区间。 解:由于函数104)(23-+=x x x f 的定义域为()+∞∞-,,用逐步搜索法:由于 014)2(,010)0(>=<-=f f ,于是,方程在)2,0(内至少有一个实根,所以,从0=x , 取步长5.0=h 向右进行根的搜索,即计算)5.1(),0.1(),5.0(f f f 的值得到 0)5.1(,0)1(,0)5.0(>< 6.二分发的基本思想是什么? 解:二分发的基本思想是将方程0)(=x f 的有根区间逐步分半,通过判别)(x f 在端 点的符号以及零点定理来缩小有根区间,使在足够小的区间内使方程0)(=x f 有且仅有一个根,并满足给定的精度要求为止。 7.以方程0)(=x f 的有根区间为[]b a ,为例)0)(,0)((> 解:第一步:将有根区间[]b a ,分半,用区间[]b a ,的中点2 b a +将[] b a ,分为两个相等区间,计算中点的函数值)2(b a f +。若0)2(=+b a f ,则2 *b a x +=就是方程0)(=x f 的根;否则,若0)2( <+b a f ,由于)(x f 在左半区间?? ????+2,b a a 内不变号,所以方程的有根区间变为???? ??+b b a ,2。同理,若0)2(>+b a f ,则方程的有根区间变为??? ?? ?+2,b a a ,从而将新的有根区间记为[]11,b a ,且区间[]11,b a 的长度仅为区间[]b a ,的一半,即2 11a b a b -=-。 第二步:对压缩了的有根区间[]11,b a 又可施行同样的方法,即用中点 2 1 1b a +将区间[]11,b a 再分为两半,然后通过根的搜索判定所求的根位于哪半个区间,从而又确定一个新的有根区间[]22,b a ,该区间的长度是区间[]11,b a 的一半。 如此反复可得出一系列有根区间且具有关系[][][] ????k k b a b a b a ,,,11,其中后一个区间长是前一个区间长的一半,因此区间[]k k b a ,的长度k k k a b a b 2-= -,当∞→k 时,区间[]k k b a ,的长度必趋于零,即这些区间最终收缩于一点*x ,显然*x 就是方 程0)(=x f 的根。 8.以方程0)(=x f 的有根区间为[]b a ,,精度要求为ε,试写出利用二分法求该方程的近似根所需二分次数k 的计算公式。 解:若事先给定的精度要求为0>ε,则只需ε<-≤ -+1 * 2k k a b x x ,即2 ln 2ln εa b k ->,此时k x 就是满足给定精度要求的近似值,k 为二分法的次数。 9.用二分法求下列方程在给定的有限区间及精度要求下的近似值及二分次数k (编程) (1)2)(-=x xe x f []1,5.0 0001.0=JD 解:852600.0=k x 12=k (2)343)(23-+-=x x x x f []5.1,1 00001.0=JD 解:499992.1=k x 15=k (3)104)(23-+=x x x f []2,1 0005.0=JD 解:364746.1=k x 10=k (4)1)(3 --=x x x f []5.1,1 00005.0=JD 解:324707.1=k x 13=k 10.若应用二分法求方程02 sin =--x e x π在区间[]1,0上误差不超过 5 21 的近似值,应二分多少次? 解:其近似根为437500.0,应分5=k 次。 11.迭代法的基本思想是什么? 解:迭代法是一种逐次逼近法,首先给定方程0)(=x f 的一个粗糙的初始近似根0x , 然后用一个固定公式反复校正这个根的近似值使之逐步精确化,直到满足预先给定的精度要求为止。 12.迭代法的具体做法如何? 解:(1)将方程0)(=x f 改写成等价形式)(x x ?=,在根*x 的附近任取一个初始近似根 0x 。 (2)构造近似根序列:将0x 代入)(x ?计算得到)(01x x ?=,一般01x x ≠,再把1x 作 为新的近似根代入)(x ?得到)(12x x ?=,重复上述步骤即可。 13.迭代法的几何意义是什么? 答:方程()x x ?=的求根问题在几何上就是确定曲线()x y ?=与直线x y =交点*p 的横坐标* x 。设迭代初值为0x ,曲线()x y ?=上以0x 为横坐标的点为0p ,()0x ?为0p 点的 纵坐标,过0p 点引平行于x 轴的直线,并与直线x y =相交于0P ',其横坐标为()01x x ?=,然后过点0P '引平行线于y 轴的直线,并与曲线()x y ?=的交点记作1p ,重复上述过程可得点列,,,,,21 k p p p 他们横坐标依次由迭代公式()k k x x ?=+1, 1,0=k 所确定。如果点列,,,,,21 k p p p 逐步逼近* p ,则迭代过程收敛,否则迭代过程发散。 14.叙述迭代过程收敛定理的内容。 解:假设迭代函数满足下列两个条件 (1)对任意的[]b a x ,∈有b x a ≤≤)(?; (2)存在正数1 则(1)对任意初值[]b a x ,0∈迭代过程)(1k k x x ?=+均收敛于方程)(x x ?=的根* x , 即)(lim *∞→=k x x k 。 (2)误差事后估计公式为k k k x x L x x --≤ -+1* 11 。 15.试构造收敛的迭代公式求解下列方程: (1)4 sin cos x x x += ; (2)x x 24-=。 解:(1)将方程4sin cos x x x +=改写为4) 4sin (2π += x x ,从而得到迭代公式 0,1,2,k 4 ) 4sin(21=+ = +π k k x x 。 (2)将方程x x 24-=改写为)4l n (x x -=,从而得到迭代公式 ,2,1,0 )4ln(1=-=+k x x k k 。 16.判断迭代法解方程0)2ln()(=+-=x x x f 在[]2,0内的根时所用的迭代过程的收敛性。 解:将方程0)2ln(=+-x x 改写为)2l n (+=x x ,从而得到迭代公式 ,2,1,0),2ln(1=+=+k x x k k 。则)2ln()(+=x x ?为迭代函数。由12 1 )(<+= 'x x ?,由定理3.2可得该迭代法是收敛的。 17.用迭代法计算4 4446666 ++++= s 的近似值。 19.牛顿法的基本思想是什么?具体做法如何? 解:基本思想:牛顿迭代法实质上是一种线性化的方法,其基本思想是将非线性方程 0)(=x f 逐步归结为某种线性方程来求解的方法。 具体做法:设已知方程0)(=x f 有近似根k x ,将)(x f 在k x 作一阶泰勒展开,于是方程0)(=x f 可近似地表示为0))(()(=-'+k k k x x x f x f 是一个线性方程,设0)(≠'k x f ,则)()(k k k x f x f x x '- =,于是就有牛顿迭代公式 ,2,1,0,) () (1='-=+k x f x f x x k k k k 。 20.牛顿法的几何意义是什么? 解:牛顿迭代法实质上是用过点))(,(k k x f x 的切线与x 轴交点的横坐标1+k x 来逐步逼 近曲线)(x f y =与x 轴交点的横坐标* x ,所以牛顿法又叫切线法。 22.试证:用牛顿法求方程0)3()2(2 =+-x x 在[]3,1内的根2* =x 是线性收敛的。 证明:由牛顿迭代公式 ,2,1,0,) () (1='- =+k x f x f x x k k k k ,可得,4 36 32)()()(2+++= '-=x x x x f x f x x ?,显然,0)2(≠'?,所以该迭代过程是线性收敛的。 23.用牛顿法求方程03 =-a x ,导出求立方根3a 的迭代公式,并讨论其收敛性。 解:设0)(3 =-=a x x f ,得牛顿迭代公式为 ,1,0,32 3 1 ==--=+k x a x x x k k k k ,牛顿迭代函数2332)(x a x x +=?,3 3322)(x a x x -='?,10)(3 <='a ?,所以该迭代公式收敛。 26.正割迭代法的基本思想是什么?具体做法如何?几何意义是什么? 解:基本思想:用过两点))(,(k k x f x ,))(,(11--k k x f x 的直线的斜率这个差商来代替牛堆迭代公式中的倒数)(k x f '。 具体做法:对方程0)(=x f 经过k 次迭代后得到近似根1-k x ,k x ,从而取 ) ())()(()(11----≈ 'k k k k k x x x f x f x f ,于是牛顿迭代公式变为 )() ()() (111--+--- =k k k k k k k x x x f x f x f x x ,此公式为正割法迭代公式。 几何意义:正割迭代法是用过两点))(,(k k x f x A ,))(,(11--k k x f x B 的直线与x 轴交点的横坐标1+k x 来逐步逼近曲线)(x f 与x 轴交点的横坐标* x ,因此正割迭代法又叫割线法。 27.简述正割迭代法与牛顿迭代法的区别。 解:牛顿迭代法在计算时只需要一个初值0x ,在计算1+k x 只用到前一步的值Xk ,但要 计算)(k x f ';而正割法在计算时需要两个初值0x ,1x ,在计算1+k x 时要用到前两次的迭代值1-k x ,k x ,但不用计算导数。 30.使迭代法加速的方法有哪些?并分别写出它们的迭代公式。 答:使迭代法加速的方法有艾特肯加速公式和斯蒂芬森方法: 艾特肯加速公式: 校正:()K +K X =X ?1 再校正:() 12+K +K X =X ? () ,2,1,0=K 改进:()2 1 2 1 2 221 +K +K K +K +K X +X -X X -X +K +K -X =X 斯蒂芬森方法: 迭代:()()();,2,1,0, =K Y =Z X =Y K K K K ?? 加速:()().,2,1,0212 =K -X =X K K K K K X +Y -Z X -Y K +K 第三章 线性方程组的数值解法 1.线性方程组的数值解法有哪两大类?并简述他们的概念。 答:线性方程组的数值解法有两大类: (1)直接法 :直接法就是在没有舍入误差的情况下,经过有限步算术运算可求得 方程组精确解的算法。 (2)迭代法:迭代法就是用某种极限过程去逐步逼近线性方程组精确解的方法,即先给定一个初始解向量,然后按新的迭代公式逐步求出解的更准确值的方法。 2.高斯消去法的基本思想是什么? 答:高斯消去法的基本思想是用逐次消去未知量的方法把原来方程组b AX =化为与其同解的三角形方程组,而求解三角形方程组就容易了。 3.高斯主元素消去法是在何种情况下提出来的? 答:用高斯消去法解线性方程组b AX =的消元过程中,可能会出现以下两种情况:第一是主元素全是0的情形,致使消元过程无法进行下去;第二即使主元素不为0,但其绝对值很小时作除数可能会导致其他元素数量级的严重增长和舍入误差的传播,使计算结果不可靠。所以对于一般矩阵来说,最好每一步选取系数矩阵中绝对值大的元素作为主元素。 4.用高斯顺序消去法,完全主元素消去法和列主元素消去法解下列方程组,并写出高斯顺序消去法的程序。 (1)?????-=--=+-=++44385522321321321x x x x x x x x x ; (2)??? ??=---=-+-=+-0 232122743321 321321x x x x x x x x x 。 解:(1)将方程组的增广矩阵进行初等变化,并利用高斯顺序消去法得: ?????=-==→????? ??→????? ??→????? ??211x 4 2 0 09 8 7 05 2 1 213- 10- 7- 09- 8- 7- 05 2 1 24- 4- 3- 18 1 1- 55 2 1 23 21x x ; 利用完全主元素消去法得: ??? ?? ??→????? ??→????? ??→????? ??4 2 3 09 7 8 08 1 - 1 54 3 2 09 8 7 08 1 1- 54- 4 - 3- 15 2 1 28 1 1- 54- 4- 3- 18 1 1- 55 2 1 2 ??? ??-===→????? ??→121x 5- 5 0 09 7 8 08 1 - 1 52 31x x ; 利用列主元素消去法得: ?????=-==→????? ??→????? ??→????? ??→????? ??2 11x 10 5 0 09 8 7 08 1 1- 54 3 2 09 8 7 08 1 1- 54- 4 - 3- 15 2 1 28 1 1- 54- 4- 3- 18 1 1- 55 2 1 2321x x (2)将方程组的增广矩阵进行初等变化,并利用高斯顺序消去法得: ??? ? ??? ===→????? ??→????? ??→????? ??21 12x 1 2 0 04 2- 5 07 4 1- 32 2 1 04 2- 5 07 4 1- 30 2- 3- 21- 2- 2 1-7 4 1- 33 21x x ; 利用完全主元素消去法得: ???? ? ??→????? ??→????? ??→????? ??1 1 1- 05 1 3 07 3 1- 41 1 1- 05 1 3 07 3 1- 40 2 3- 2-1- 1- 2 2-7 3 1- 40 2- 3- 21- 2- 2 1-7 4 1- 3 ??? ???? === →????? ??→2 121x 8 4 0 05 1 3 07 3 1- 4123x x ; 利用列主元素消去法得: ??? ???? ===→????? ??→????? ??→????? ??21 12x 1 2 0 04 2- 5 07 4 1- 32 2 1 04 2- 5 07 4 1- 30 2- 3- 21- 2- 2 1-7 4 1- 33 21x x 。 5.用矩阵的三角分解法解下列方程组,并掌握三角分解法的编程思路。 (1)??????????=????????????????????785 20- 2 6- 16- 18 4-8 4 2-321x x x ; (2)?????? ????=????????????????????201814 5 1 3 2 5 2 3 2 1321x x x 。 解:(1)对系数矩阵A 作如下的三角分解: ???? ????????????????=??????????332322131211323121 u 0 0u u 0u u 1 l l 0 1 l 0 0 1 20- 2 6- 16- 18 4-8 4 2-u 。 根据矩阵的乘法可得: 2211111-=?-=?u u ,4411212=?=?u u ,8811313=?=?u u ; 24211121=?-=?l u l ,1018122221221=?=?+?u u u l , 3216123231321-=?-=?+?u u u l ;36311131=?-=?l u l , 123222321231-=?=?+?l u l u l , 76201333323321331-=?-=?+?+?u u u l u l 。 于是有LU A =?? ?? ??????-??????????=76- 0 032- 10 08 4 21 1- 30 1 20 0 1,则原方程组可表示为 ?? ???? ????=????????????????????-??????????78576- 0 032- 10 08 4 21 1- 30 1 20 0 1321x x x 。解方程组b Ly =,即????? ?????=????????????????????7851 1- 30 1 20 0 1321y y y ,得 ??????????--=1025y 。解方程组y Ux =,即??????????--=????????????????????-102576- 0 032- 10 08 4 2321x x x ,得?? ?? ?? ? ? ? ?????????-=385 9521 190291x 。 (2)对系数矩阵A 作如下的三角分解: ???? ? ???????????? ???=??????????332322131211323121u 0 0u u 0u u 1 l l 0 1 l 0 0 1 5 1 3 2 5 23 2 1u 。 根据矩阵的乘法可得: 1111111=?=?u u ,2211212=?=?u u ,3311313=?=?u u ; 22211121=?=?l u l ,151********=?=?+?u u u l , 42123231321-=?=?+?u u u l ;33311131=?=?l u l , 513222321231-=?=?+?l u l u l , 2451333323321331-=?=?+?+?u u u l u l 。 于是有LU A =?? ?? ? ???????????????=24- 0 04- 1 03 2 11 5- 30 1 20 0 1,则原方程组可表示为 ?? ????????=??????????????????????????????20181424- 0 04- 1 03 2 11 5- 30 1 20 0 1321x x x 。解方程组b Ly =,即????? ?????=????????????????????2018141 5- 30 1 20 0 1321y y y ,得??????????--=721014y 。解方程组y Ux =,即?????? ????--=????????????????????721014 24- 0 04- 1 03 2 1321x x x ,得?? ??? ?????=321x 。 6.用追赶法解下列方程组。 (1)???? ????????=????????????????????? ???0001 2 1- 0 0 1- 2 1- 00 1- 2 1-0 0 1- 24321x x x x ; (2) ????? ? ??????-=????????????????????????1 21 6 5 3- 0 0 3- 4 2- 00 1- 3 1-0 0 1- 24321x x x x 。 解:(1)由LU A =得: ?? ???? ??????????? ?? ?????=????????? ??? 1 0 0 0 1 0 00 1 00 0 1 0 00 00 0 0 0 0 2 1- 0 0 1- 2 1- 00 1- 2 1-0 0 1- 23 214433221βββαγαγαγα。 于是有,2 3 2,1,211,2221221111=?=+?-=-=?-=?=ααβγγββαα 1,1,3 4 2,1,32133433233222-=?-==?=+?-=-=?-=?βαγααβγγββα 4 5 2,4344343=?=+?-=?ααβγβ。 从而f Ly =为 ????????????=??????????????????????????????? ?0 00145 1- 0 00 34 1- 00 0 23 1-0 0 0 243 21y y y y ,解得??????? ???? ???????????=51413121y 。y Ux =为 ??????????? ???????????=??????????????????????????????? ?51413121 1 0 0 043- 1 0 00 3 2- 1 00 0 21- 14321x x x x ,解得??? ??????? ????????????=51525354x 。 (2)由LU A =得: ?? ????????? ???????? ?????=????????? ??? 1 0 0 0 1 0 00 1 00 0 1 0 00 00 0 0 0 0 5 3- 0 0 3- 4 2- 00 1- 3 1-0 0 1- 23 214433221βββαγαγαγα。 于是有,2 5 3,1,211,2221221111=?=+?-=-=?-=?=ααβγγββαα 3,3,5 16 4,2,52 133433233222-=?-== ?=+?-=-=?-=?βαγααβγγββα 16 355,161544343=?=+?- =?ααβγβ。 从而f Ly =为 ? ? ?? ? ???????-=????????????? ????????? ??????????1 21 6 1635 3- 0 00 516 2- 00 0 25 1-0 0 0 2432 1 y y y y ,解得????????? ???????????=353483583y 。y Ux =为 ????????? ???????????=????????????????????????????????353483583 1 0 0 01615- 1 0 00 52- 1 00 0 21- 14321x x x x ,解得?? ?????? ? ? ? ???????????=353479357435142x 。 7.设()n T n R x x x x ∈=,,,21 ,掌握常用向量范数的定义式P x x x x ,,,21∞。 解:∑==+++=n i i n x x x x x 1 211 ; {}n i x x x x x ,,,max max 21 ==∞; (又叫最大范数) 2 11 2222 2 1 2 )(∑==+++=n i i n x x x x x ; P n i P i P x x 11 )(∑==。 8.已知()T x 0,3,1,1-=,计算P x x x x ,,,21∞。 解:∑===+++=n i i n x x x x x 1 211 5 ; {}3,,,max max 21===∞ n i x x x x x ; 11)(21 122 22212 ==+++=∑=n i i n x x x x x ; P P P n i P i P x x 111 )32()(+==∑=。 9.设n n n n ij R a A ??∈=)(,掌握常用矩阵范数的定义式F A A A A ,,,21∞。 解:∑==n i ij a A 1 1max ; ∑=∞ =n j ij a x 1 max ; )(max 2A A A T λ=; 2 11 ,2) (∑==n j i ij F a A 。 10.已知? ? ? ? ??=0 31- 2A ,计算F A A A A ,,,21∞。 解:5max 1 1==∑=n i ij a A ; 3max 1 ==∑=∞ n j ij a x ; 1027)(max 2+==A A A T λ; 14)(211 ,2==∑=n j i ij F a A 。 12.解线性方程组的迭代法有哪三种方法? 答:(1)雅可比迭代法(Jacobi ) (2)高斯--赛德尔迭代法(G--S ) (3)超松弛迭代法(SOR ) 13.设有方程组??? ??=+-=++--=++3 103220241225321 321321x x x x x x x x x (1)写出用Jacobi 迭代法解此方程组迭代公式的分量形式和矩阵形式。 (2)Jacobi 迭代法是否收敛?为什么? 解:该方程组可化为:??? ??++-=+-=---=3 .03.02.055.025.04.22.04.0213 312321x x x x x x x x x ,从而得到Jacobi 迭代法的公式: ?? ???++-=+-=---=+++3.03.02.055.025.04.22.04.0)(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 ,其矩阵形式为:b D X L U D X k k 1)(1) 1()(--+++-=,其中:????? ??=10 0 00 4 00 0 5D ,????? ??=0 0 02 0 01 2 0U ,????? ??=0 3- 20 0 1-0 0 0L ,??? ? ? ??-=32012b (2)用Jacobi 迭代法解此方程组是收敛的。因为系数矩阵???? ? ??=10 3- 22 4 1-1 2 5A 是严格对角占 优阵,所以Jacobi 迭代法收敛。 14.设有方程组??? ??=+-=++=++30 1532128243220321 321321x x x x x x x x x (1)写出用G-S 迭代法解此方程组迭代公式的分量形式。 (2)G-S 迭代法是否收敛?为什么? 解:该方程组可化为:??? ??++-=+-=+--=2 2.013 3.05.1125.0125.02.115.01.0213 312321x x x x x x x x x ,从而得到G-S 迭代法的公式: ?????++-=+-=---=++++++22.0133.05.1125.0125.02.115.01.0) 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)用G-S 迭代法解此方程组是收敛的。因为系数矩阵???? ? ??=15 3- 21 8 1 3 2 20A 是严格对角占优 阵,所以G-S 迭代法收敛。 15.设有方程组?? ? ??=--=+?+-=?++-7 9 890708321321321x x x x x x x x x 怎样改变方程的顺序使Jacobi 迭代法和G-S 迭代法均收敛。 解:将方程组变化成?? ? ??=+?+-=?++-=--890 70879321321321x x x x x x x x x ,此时系数矩阵????? ??=9 0 1-0 8 1-1- 1- 9A 为严格对角占优矩阵,所以Jacobi 迭代法和G-S 迭代法均收敛。 16.设方程组???=+=+22221121212111b x a x a b x a x a )0(2211≠?a a ,迭代公式为??? ? ???-=-=--) (1)(1)1(121222)(2) 1(212111)(1k k k k x a b a x x a b a x ),2,1( =k 。求证:由上述迭代公式产生的向量序列{} )(k x 收敛的充要条件是 122 1112 21 证明:由题设知:???? ??+???? ??+???? ? ?=???? ??=0 0 00 0 0 00 1221221122211211a a a a a a a a A ,所以: 迭代矩阵? ????? ? ?? -=???? ?????? ??-=-0 - 00 0 00 2221111221121 2211a a a a a a a a B ,22 112112a a a a B ±=λ,所以由迭代法收敛的充要条件1)( (k x 收敛的充要条件 是 122 1112