《计算方法》习题答案

《计算方法》习题答案

第一章数值计算中的误差

1.什么是计算方法？（狭义解释）

答：计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算，以便在计

算机上编程上机，求出问题的数值解，并对算法的收敛性、稳定性和误差进行分析、计算。

2.一个实际问题利用计算机解决所采取的五个步骤是什么？

答：一个实际问题当利用计算机来解决时，应采取以下五个步骤：

实际问题一建立数学模型一构造数值算法一编程上机一获得近似结果

4.利用秦九韶算法计算多项式尸（x）=x-/+》5-4在x=—3处的值，并编程获得解。

解：P（x）=x5+0-x4-x3+0-x2+x-4,从而

10-101-4

-3-39-2472-219

1-38-2473-223

所以，多项式尸（x）=x—/+/-4在x=-3处的值尸（-3）=-223。

5.叙述误差的种类及来源。

答：误差的种类及来源有如下四个方面：

（1）模型误差：数学模型是对实际问题进行抽象，忽略一些次要因素简化得到的，

它是原始问题的近似，即使数学模型能求出准确解，也与实际问题的真解不同，我们把数学

模型与实际问题之间存在的误差称为模型误差。

（2）观测误差：在建模和具体运算过程中所用的一些原始数据往往都是通过观测、

实验得来的，由于仪器的精密性，实验手段的局限性，周围环境的变化以及人们的工作态度

和能力等因素，而使数据必然带有误差，这种误差称为观测误差。

（3）截断误差：理论上的精确值往往要求用无限次的运算才能得到，而实际运算时

只能用有限次运算的结果来近似，这样引起的误差称为截断误差（或方法误差）。

（4）舍入误差：在数值计算过程中还会用到一些无穷小数，而计算机受机器字长的

限制，它所能表示的数据只能是一定的有限数位，需要把数据按四舍五入成一定位数的近似

的有理数来代替。这样引起的误差称为舍入误差。

6.掌握绝对误差（限）和相对误差（限）的定义公式。

答：设X*是某个量的精确值，X是其近似值，则称差2=%*-X为近似值X的绝对误差

（简称误差）。若存在一个正数£使|e|=|X*—X称这个数£为近似值X的绝对误

差限（简称误差限或精度）。

*

.pX—Y

把绝对误差e与精确值x之比二二=J之称为近似值x的相对误差，称

xx

〃为近似值x的相对误差限，由于真值x*是未知的，所以常常用

*

e,=±二£=刍来表示相对误差，于是相对误差可以从绝对误差求出。

xx

7.近似值的规格化表示形式如何？

答：一般地，对于一个精确值X*,其近似值x的规格化形式为'=±0.七X23》0*10"',

其中占#0,x,e{0,l,2,…9}(i=l,2,…p),p为正整数，机为整数。

8.有效数字的概念是什么？掌握有效数字与误差的关系。

答：若近似值x的(绝对)误差限是它的某一位的半个单位，也就是说该近似值准确到

这一位，且从该位起直到前面第一个非零数字为止的所有数字都称为有效数字。

若近似值x的(绝对)误差限为忖=k*-Xw；xlO'"",则称x为具有〃位有效

数字的有效数，或称它精确到10"'一"位，其中的每一位数字项,它，…猫都是》的有效数字。

w

设精确值x*的近似值x的规格化形式为》=±0.x,x2■■■xpx10,若x具有〃位有

效数字，则其相对误差限为WKjxlOi;反之，若X的相对误差限为

led<---xIO1'",则x至少有〃位有效数字。

9.下列各数都是对真值进行四舍五入后获得的近似值，试分别写出它们的绝对误差限，相

对误差限和有效数字的位数。

5

(1)项=0.024(2)x2=0.4135(3)x3=57.50(4)x4=60000(5)x5=8x10；

♦

解：(1)|e1=10.0005；|e,|=匚'=240.0021；有三位有效数字。

x-x=|-|<0.000121；有四位有效数字。

(2)|e|=|x2*-x2|<0.00005；|er|=

x-x

(3)[e|=|x3*-x3|<0.005；\et\<0.000087；有四位有效数字。

x-x=|-|<0.0000084；有五位有效数字。

XX

(4)|e|=|4*-4|<0.5；\er\

x-x

(5)|e|=|x5*-X5|<0.5；\er\==-<0.000000625；有六位有效数字。

X

10.为了使M的相对误差40.1%,问至少应取几位有效数字？

解：由J历的首位数是4.设近似数x*有〃位有效数字，由定理4.1可知，相对误差

|e,,（x*）|<^xl0'-"<0.001,解得〃23.097,即取4位有效数字，近似数的相对误差

2x4

不超过0.1%。

11.已知y=P(x)=/+x—1150,x*=等,x=33,计算y*=p(与)及y=P(33),

并求x和）;的相对误差。

2

解：y*=^(122)=(122)+(122)-1150®-5.55555•••

y=尸(33)=(33y+»3)_1150=_28

|e(x)|--x|«0.333…

\er(x)|=«0.0101---

|e(y)|=\y*-y卜22.44444•••

e(y)

|/(刈=。0.801587…

y

12.写出误差估计的一般公式（以二元函数2=/（XJ）为例）。

解:二元函数z=/（x,y）的绝对误差:

二元函数的相对误差:

"**•吸加号

洛小）匕3+谭心他

13.用电表测得一个电阻两端的电压和流过的电流范围分别为修=220±2/，

/=10±0.14,求这个电阻的阻值R,并估算其绝对误差和相对误差。

解:|e(K)|<2,|e(/)|<0.1,又火=匕，空=,您=一与。所以：

1111IdvIdiI2

匕⑻核防g)，e(P)+行|(々)・e(/)<

dR./〃、dR./八1.

方k“)・e(/)+防l(“)-e(Z)=—x2-

e,.(R)=^^a1.99x10-2。

计算丁=；；的近似值，并估计及

14.若x：=1.03±0.01,x；=0.45±0.01,x+*e(y)

其上界。

解：ya(L03)2+；e°45

|e(M=|y*-y\=(x；+；*)—(项+；*)=(x：—X|)(x：+X|)+；(/2_*)

2

<|(x：-X|)(x：+x1)|+_e*2)=2.06x10-+gx/xO.Ol,^e(/,x；)

15.已测得某场地长为/=110加，宽d的值为d=80加，已知上(/)|=/一/卜0.2”，

|e(d)|=|4*-440.1加，试求面积s=/d的绝对误差限和相对误差限。

解：由5=〃，5=d,：=/,|e(/)|=『—/]<0.2加，|e(d)|=|d*—4〈0.1加。

合(a,

可得：|e(s)|®—1(/.，/)+—l(/,(y)e(d)<—|(M)e(/)+—|(/>rf)-e(d)

=110x0.2+80x0.1=30

er(.$,)=e")«3.4x10'。

s

16.掌握二元函数的加、减、乘、除和开方运算的绝对误差和相对误差估计公式。

解：(1)力口、减运算：

由于d(x+y)/dx=l5(x+y)/dy=\,5(x-y)ldx=1,3(x-y)/dy=-\,,所以

e(x+力e(x)+e(y),e,(x+力x/(x+y)xe,.(x)+y/(x+y)xer(y),e(x-力e(x)-e(y\

e,(x-y)«x/(x-y)xer(x)-y/(x-y)xer(y\从而有|e,(x-y)|<|x/(x-y)\-er(x)|+

|^/(x-y)|-|er(y)|

(2)乘法运算:

由于。(叶)=乂=x，所以e(xy)«ye(x)+xe(y),er(xy)«er(x)+er(y),从而

drdy

|e®)|w3，|e(x)Hx|“e(y)|

(3)除法运算：

由于5(—)/dx=—,5(—)/dy=---,所以e(—)«—e(x)

y/yy/yyyy

e,(与Re,(x)-e,O)

y

(4)乘方及开方运算：

n

由于="x"T,所以e(x")xnx~'e(x),er(x"/ner(x)

17.求方程/-56x+l=0的两个根，使它至少具有4位有效数字(a27.982)。

56+J(-56)2-4x1x1

解：再»28+27.982=55.782

2x1

1

»0.017863

55.782

19.求方程x2-16x+l=0的较小正根，要求有3位有效数字。

16+7(-16)2-4x1x1

解:X]«8+7.937=15.937

2x1

1

»0.062747

15.937

所以较小正根为/«0.062747o

20.设/〃=fx"e?x,〃=0,l,2,…JO'

J0

(1)证明：I,,=e-=0,1,2,­••,104；

(2)给出一个数值稳定的算法，并证明算法的稳定性。

xnxn

(1)证明：/“=£x"edx=£xde=e-(nx-'e'dx=e-nln_y

(2)/,-=L(e—/“)

n

设同=，:一小则

le»-i|=|Ci射|="l

|^-2|=\C-2

匕。|=|/；-/。卜/⑸I

当〃无限大时，e,越小，所以该算法稳定。

21.用递推算法计算积分/“=pfl―X」右,〃=0,1,2…,10,并验证算法的数值稳定性。

J。1+4%

=-(['x,-'dx-f=—

解：4〃一

4Jol+4x4J。J(M+4x4”41

设⑷=，；—//则

团=，；-小；跖|

|e2|T，j|=*o|

|修。|=/：。-。|=/高|

所以该算法是稳定的。

22.设计一个计算/（x）=x12+3x24+16/6的最小计算量的算法。

解:/(x)=x12+3x24+16x36=x-x-x2-x4-x4+3-x12-x]2+16-x12-X24

23.什么是数值稳定的算法？数值计算应遵循的六条规则是什么？

答：一个算法如果原始数据有误差（扰动），而计算过程中舍入误差不增长或增长可以

控制，则称此算法是数值稳定的。否则，称此算法是数值不稳定的。

数值计算应遵循的六条规则是：

(1)选用数值稳定的算法（计算公式）；

(2)尽量避免两个相近数相减；

(3)尽量避免用绝对值很大的数作乘数；

(4)尽量避免用绝对值很小的数作除数；

(5)防止大数“吃掉”（或“淹没”）小数（即合理安排运算顺序）;

(6)简化计算步骤，减少运算次数。

第二章非线性方程的数值解法

1.叙述零点定理的内容。

答：设函数/(x)在闭区间口,切上连续且/(。>/3)<0,则存在X*w(q,b)使

/(x*)=0,即/(x)在区间(a,b)内存在实的零点，称区间［凡切为方程的有根区间。

2.方程求根的两个步骤是什么？确定方程有根区间的方法有哪些？

答：第一步确定方程/(x)=0的有根区间。

第二步近似根的精确化。

确定方程有根区间的方法有两种：作图法和逐步搜索法。

3.利用作图法确定方程/(x)=/一x-1=0的有根区间。

由于/(0)=—1<0,/(2)=8-2—1=5>0,于是，在区间(0,2)内至少有一个根，取步

长/?=0.5向右进行根的搜索，即计算/(0.5),/(1.0),/(1.5)的值得到

/(0.5)<0,/(1.0)<0,/(1.5)>0,从而,原方程的有根区间缩小为。,1.5)

4.利用逐步搜索法确定方程/(X)=X3-3X2+4X+3=0的有根区间。

解：由于/(0)=3>0J(—1)=—5<0,于是，方程在(—1,0)内至少有一个实根，所以，

从x=-l,取步长〃=0.5向右进行根的搜索，即计算/(一0.5)得到/(—0.5)=J>0,从

8

而，原方程的有根区间缩小为(-I,-/)。

5.确定方程/+4/-10=0的有根区间。

解：由于函数/(x)=x3+4x2—10的定义域为(—8,+8),用逐步搜索法：由于

/(0)=-10<0,/(2)=14>0,于是，方程在(0,2)内至少有一个实根，所以，从x=0,

取步长。=0.5向右进行根的搜索，即计算/(0.5),/(1.0),/(1.5)的值得到

/(0.5)<0,/(1)<0,/(1.5)>0,从而原方程的有根区间缩小为(1,1.5)

6.二分发的基本思想是什么？

解：二分发的基本思想是将方程/(x)=0的有根区间逐步分半，通过判别/(x)在端

点的符号以及零点定理来缩小有根区间，使在足够小的区间内使方程/(x)=0有且仅有一

个根，并满足给定的精度要求为止。

7.以方程/")=0的有根区间为［a㈤为例(/⑷<0,/(/>)>0),简述二分法的具体作法。

解：第一步：将有根区间［a,"分半，用区间［a,4的中点审将［凡同分为两个相等

区间，计算中点的函数值/(巴9)。若/(审)=0,则x*=巴芋就是方程/(x)=0的

根；否则，若,/■(一)<()，由于/(x)在左半区间a,半内不变号，所以方程的有根

区间变为等■力。同理，若/(+)〉0,则方程的有根区间变为a,手,从而将

新的有根区间记为［外4］,且区间［四,仇］的长度仅为区间［。力］的一半，即仇-q=。

第二步：对压缩了的有根区间［q,6』又可施行同样的方法，即用中点生砂将区

间院,6』再分为两半，然后通过根的搜索判定所求的根位于哪半个区间，从而又确定一个

新的有根区间［%，&］，该区间的长度是区间［四,6』的一半。

如此反复可得出一系列有根区间且具有关系［a,b］=>［«,,/>,］=)­..=>瓦也］=)…，

其中后一个区间长是前一个区间长的一半，因此区间［知,“］的长度4-斯=空，当

左一00时，区间［%,4］的长度必趋于零，即这些区间最终收缩于一点显然x*就是方

程/(x)=0的根。

8.以方程/(x)=0的有根区间为［a,“，精度要求为£,试写出利用二分法求该方程的近

似根所需二分次数〃的计算公式。

1b-a

1_In----

解：若事先给定的精度要求为£>0,则只需卜*—〈舒<£,即:〉一］音

此时居就是满足给定精度要求的近似值，k为二分法的次数。

9.用二分法求下列方程在给定的有限区间及精度要求下的近似值及二分次数〃(编程)

(1)/(x)=xex-2[0.5,1]㈤=0.0001

解：4=0.852600k=12

(2)f(x)=x3-3x2+4x-3[1,1.5]J©=0.00001

解：4=1.499992左=15

(3)/(x)=x3+4x2-10[1,2]JD=0.0005

解：xk=1.364746k=10

(4)f(x)=x3-x-1[1,1.5]JD=0.00005

解：xk=1.324707人=13

10.若应用二分法求方程sin/=0在区间［0,1］上误差不超过!的近似值，应二分

多少次？

解：其近似根为0.437500,应分％=5次。

11.迭代法的基本思想是什么？

解：迭代法是一种逐次逼近法，首先给定方程/(x)=0的一个粗糙的初始近似根X。，

然后用一个固定公式反复校正这个根的近似值使之逐步精确化，直到满足预先给定的精度要

求为止。

12.迭代法的具体做法如何？

解：(1)将方程/(x)=0改写成等价形式》=奴x)，在根X*的附近任取一个初始近似根

X。°

(2)构造近似根序列：将为代入9(x)计算得到匹=°(%),一般%也勾,再把再作

为新的近似根代入夕(X)得到x2=火花)，重复上述步骤即可。

13.迭代法的几何意义是什么？

答：方程x=e(x)的求根问题在几何上就是确定曲线丁=夕G)与直线y=x交点p*的

横坐标x*。设迭代初值为演，曲线y=e(x)上以演为横坐标的点为20,Q(Xo)为Po点的

纵坐标，过Po点引平行于x轴的直线，并与直线_y=x相交于耳，其横坐标为否=员/),

然后过点以引平行线于y轴的直线，并与曲线y=9(x)的交点记作白，重复上述过程可得

点列0],02,…,p*,…,他们横坐标依次由迭代公式x«+i=8(xJ，左=0,1…所确定。如果

点列Pl，P2,…，小，…，逐步逼近P*，则迭代过程收敛，否则迭代过程发散。

14.叙述迭代过程收敛定理的内容。

解：假设迭代函数满足下列两个条件

(1)对任意的xe[a,b]有aVe(x)《b；

(2)存在正数£<1,使对任意“有

则(1)对任意初值X。e[a,可迭代过程=Q(XQ均收敛于方程x=0(x)的根x*,

即limxk=x*(kfoo)。

(2)误差事后估计公式为k*-居|K—!—,川-居|。

1—Z/

15.试构造收敛的迭代公式求解下列方程:

/、cosx+sinx

(1)x=-----------(2)x=4-2Ao

4

V2sin(x+—)

cosx+sinx

解：(1)将方程x=改写为x=-----------从而得到迭代公式

44

Xk+\k=0,1,2，…。

(2)将方程x=4-2x改写为x=ln(4-x)，从而得到迭代公式

xk+i=ln(4-xk)k=0,1,2,­­•<)

16.判断迭代法解方程/(x)=x-ln(x+2)=0在［0,2］内的根时所用的迭代过程的收敛性。

解：将方程x—ln(x+2)=0改写为x=ln(x+2),从而得到迭代公式

xi+I=ln(xt+2),A:=0,1,2,-o则夕(x)=ln(x+2)为迭代函数。由|"(刈=—<1,

由定理3.2可得该迭代法是收敛的。

17.用迭代法计算5=^6+也+46+^^的近似值。

19.牛顿法的基本思想是什么？具体做法如何？

解：基本思想：牛顿迭代法实质上是一种线性化的方法，其基本思想是将非线性方程

/(x)=0逐步归结为某种线性方程来求解的方法。

具体做法：设已知方程/(x)=0有近似根与，将/(x)在居作一阶泰勒展开，于是方

程/(%)=0可近似地表示为f(xk)+f\xk)(x-xJ=0是一个线性方程,设/'(居)H可

则》=居—华立,于是就有牛顿迭代公式=x.一半2,攵=0,1,2「-o

/(4)f(xk)

20.牛顿法的几何意义是什么？

解：牛顿迭代法实质上是用过点(%，/(4))的切线与x轴交点的横坐标居+1来逐步逼

近曲线y=/(x)与X轴交点的横坐标X*,所以牛顿法又叫切线法。

22.试证：用牛顿法求方程(x-2)2(x+3)=0在1,3］内的根r=2是线性收敛的。

证明：由牛顿迭代公式=x*-华星■，左=0,1,2,…，可得，

/(4)

(p(x)=x-^-=——二显然，d(2)w0,所以该迭代过程是线性收敛的。

/W3X+4

23.用牛顿法求方程1一。=0,导出求立方根〃■的迭代公式，并讨论其收敛性。

3

r-a

解：设/(x)=/-a=0,得牛顿迭代公式为x«+I=x*——J—=，左=0,1,…，牛顿

3Z

迭代函数8(X)=2；。.,"(X)=2t2A|^(3/^)|=0<1,所以该迭代公式收敛。

26.正割迭代法的基本思想是什么？具体做法如何？几何意义是什么？

解：基本思想：用过两点（4,/（4）），的直线的斜率这个差商来代替

牛堆迭代公式中的倒数厂（与）。

具体做法：对方程/（x）=0经过左次迭代后得到近似根Xg,x,,从而取

■Ta*）"＜乜）―/—）），于是牛顿迭代公式变为

（符一居-1）

-，/—；（%——X"。，此公式为正割法迭代公式。

几何意义：正割迭代法是用过两点/（5.,/（乙）），6（X&T，/（X«T））的直线与X轴交点

的横坐标4M来逐步逼近曲线/（x）与x轴交点的横坐标X*,因此正割迭代法又叫割线法。

27.简述正割迭代法与牛顿迭代法的区别。

解：牛顿迭代法在计算时只需要一个初值在计算乙+1只用到前一步的值Xk,但要

计算/'（乙）；而正割法在计算时需要两个初值不，在计算时要用到前两次的迭代

值4T，Xk，但不用计算导数。

30.使迭代法加速的方法有哪些？并分别写出它们的迭代公式。

答：使迭代法加速的方法有艾特肯加速公式和斯蒂芬森方法:

艾特肯加速公式：

校正：葭二=9（XK）

再校正：二=9鼠二）（K=0,1,2,--）

(XK+2-XK+J

改进：XK+1=XK+2

XK-2XK.]+XK+2

斯蒂芬森方法：

迭代:YK=9(XK),ZK=9(YjK=0,l,2,…I

加速：XK+1=XK-I^t(K=0,1,2,-).

第三章线性方程组的数值解法

1.线性方程组的数值解法有哪两大类？并简述他们的概念。

答：线性方程组的数值解法有两大类：

(1)直接法：直接法就是在没有舍入误差的情况下，经过有限步算术运算可求得

方程组精确解的算法。

(2)迭代法：迭代法就是用某种极限过程去逐步逼近线性方程组精确解的方法，即

先给定一个初始解向量，然后按新的迭代公式逐步求出解的更准确值的方法。

2.高斯消去法的基本思想是什么？

答：高斯消去法的基本思想是用逐次消去未知量的方法把原来方程组AX=b化为与其

同解的三角形方程组，而求解三角形方程组就容易了。

3.高斯主元素消去法是在何种情况下提出来的？

答：用高斯消去法解线性方程组ZX=6的消元过程中，可能会出现以下两种情况：第

一是主元素全是0的情形，致使消元过程无法进行下去；第二即使主元素不为0,但其绝对

值很小时作除数可能会导致其他元素数量级的严重增长和舍入误差的传播，使计算结果不可

靠。所以对于一般矩阵来说，最好每一步选取系数矩阵中绝对值大的元素作为主元素。

4.用高斯顺序消去法，完全主元素消去法和列主元素消去法解下列方程组，并写出高斯顺

序消去法的程序。

2匹+x2+2X3-53X]-x2+4X3=7

(1)〈5X]-+X3=8(2)<-X]+2%2—2xj=11o

X]-3%2—4xj——42X1-3X2-2X3=0

解：(1)将方程组的增广矩阵进行初等变化，并利用高斯顺序消去法得:

2124、’212」、'212⑸=1

5-1158f0-7-8:-9f078:9-Vx2=-1

1-3-4、0-7-10113,、002；4;d3=2

（2）将方程组的增广矩阵进行初等变化，并利用高斯顺序消去法得:

3-14.：7、'3-14.：7、’3-14,；7、Xi=2

-12-2f05-2:4f05-2''Af,£=1

2-3-26、o121002

利用完全主元素消去法得:

37

1；5

14

利用列主元素消去法得:

'3-14.：7、'3-147)’3-14X]=2

-12-2—05-2-A05-2-A—<X2=1o

3-3-26、o122)、o02]_

：・x

3~2

5.用矩阵的三角分解法解下列方程组，并掌握三角分解法的编程思路。

解：（1）对系数矩阵/作如下的三角分解:

--248-'1oo-U

\\%U13

=UU

-418-1612110°2223

-62-20」3II321_00U33

根据矩阵的乘法可得:

1•%]=—2=%]=—2,1•%2=4==4,1•=8=>〃i3=8；

/21-wH=-4=>/21=2,/21-w12+1-w22=18=>u22—10,

/21-w13+1-u23=-16=>u23=-32；Z3I-wn=-6=>/31=3,

,31,W12+132.〃22=2=/32=一1，hl'U\3+’32'"23+1,比33=一2°=〃33=-76°

100-248

于是有力210010-32LU,则原方程组可表示为

-11

4100必5

10解方程组Ly=b,即2108,得

03-117

291

?90

5-2485

21

y=-2。解方程组Ur=y,即010-32X2=-2，得x

95

-10_00-76X-10

_3_5

38

(2)对系数矩阵力作如下的三角分解:

13

U22U23

033

根据矩阵的乘法可得:

1•i/H=1=>uH=1,1-w122=>w12=2,1-w13=3=>wI3=3；

21•〃]]=2=>/2I=2,/21•u]24-1-U22=5=u22=1»

=

21"13+1,》23=2="23=—4；/31•〃]]=3=4]3,

‘31'"12+132•U22=1=’32——59%•+/j2,〃23+1*〃33=5—"33=-24o

100123

于是有力=21001-4=LU,则原方程组可表示为

3-5100-24

10012314114

21001-418。解方程组功=6,即218,得

3-5100-24x320320

1412314

y-10。解方程组=即01-4X2-10

-7200-24-72

6.用追赶法解下列方程组。

'2-100'丁-2-100'-6

X

-12-10x20-13-1021

(1)；(2)—

0-12-1X300-24-3X3-2

00-1234.000-35.X4.1

解：⑴由4=得:

-2-100-000一-1A00

-12-10Yi0001夕20

0

0-12-10%000iA

00-1200/4%0001

于是有必

=20=-1=4=-p/2=-1,/2■/3{+a2=2=>%=-»

a

2Pi=T=夕2-=-1,73，+%=2nai--,74=-1,%•尸3=-1

3

n夕3=--，/4,A

一2

2

-14

3

从而切=/为0Ux=y为

1

4

05

-16-

1

1--00'2~5

2「v1

F412

-2

01——0x,5

32二'，解得x=。

18

AA,3“3

001--y45

A14—

4A

0001.5.

-3_

(2)由/=LU得：

-2-100L000Ti00

-13-10/2fz20001/?20

0-24-3073a30。。\区

00-3500九%0001

_

于是有a,=20=—lnP\=~5/2=T//]+%=3n%=5,

216

a

2%=_《，,3=-2,%,62+%=4=%=—,/4=-3,%・夕3=-3

A15A<35

=四=-石"四+%=5n%=n

-2000

5卜J8

-100-6'

215

y2

从而Ly=于为16=，解得y=30Ux=y为

0-20为-2

T8

bv1

3534

00-3

76._35.

~142~

-3-

10035

-28

74

2

010x5

*5235

=3，解得x=0

X39

158

0017

"16儿34

34

0001.35.

.35.

7.设x=（项…eR〃，掌握常用向量范数的定义式IMJWLWLIx"。

解：IM=周+同+…+k』=ZWI；

/=1

国e=max|xj=max{xj，k2|，…,同}；（又叫最大范数）

H2=卜:+4+…+X；=（£X；）2；

8.已知x=(l,—1,3,0)"计算|胤,|叱』刈2』或。

解：|x|]=|匹|+L|+…+|x/=tk，J=5；

/=1

Hoo=max|x,.|=max{x1|,|x2|,---,|x„|}=3；

IWL=X+X；+•••+》；=(£x；)2=拒；

卜IP=（EM|）P=（2+3"）'。

%设/=(%)小€瞪，掌握常用矩阵范数的定义式|4,ML，|乩,MI.。

解：I4=max£同；

llxL=maxEK-l；

7=1

|网=）―/团;

MIF=度4）2。

1。.已知"2-1o'］，计算MMK，IK，K。

解：ML=max£|%|=5；

卜L=max£同=3：

六1

/12=）九”（/4）=,7+2而；

MF=（WW）2=E。

12.解线性方程组的迭代法有哪三种方法?

答：⑴雅可比迭代法(Jacobi)

(2)高斯-赛德尔迭代法(G—S)

(3)超松弛迭代法(SOR)

5X1+2X2+£=-12

13.设有方程组＜-再+4%2+2/=20

2x,-3X2+10x3=3

(1)写出用Jacobi迭代法解此方程组迭代公式的分量形式和矩阵形式。

(2)Jacobi迭代法是否收敛？为什么？

—

X1=-0.4X2O.2X3—2.4

解：该方程组可化为：Jx2=0.25X,-0.5X3+5,从而得到Jacobi迭代法的公式：

x3=-0.2xj+0.3—2+03

xf+1)=_04娉)_0.2x『_2.4

(k+i)(k)

•X”)=0.25xf)-0.5x*+5，其矩阵形式为:X=-D-'(U+L)X+D-'b,

铲)=-0.2xf)+0,3xf)+0.3

’500、’021、’000'-2、

其中：D=040,U=002,L=-100,b二20

、0010；、°o<2-30；、3,

，52I、

(2)用Jacobi迭代法解此方程组是收敛的。因为系数矩阵/=-142是严格对角占

、2-310?

优阵，所以Jacobi迭代法收敛。

20Xj+2X2+3X3=24

14.设有方程组<X]+8工2+X3=12

2xl-3X2+15X3=30

(1)写出用G・S迭代法解此方程组迭代公式的分量形式。

(2)G・S迭代法是否收敛？为什么？

—

X1=0.1%2―0.15工3+L2

解：该方程组可化为：k2=0.125%,-0.125X3+1.5,从而得到G・S迭代法的公式：

x3=-0.133Xj+0.2%2+2

xf”)=—0.卜，一0。5力*)一1.2

<铲)=0.125一旬一0.125-)+1.5,

铲)=-0.133x”+0.2铲)+2

’2023、

(2)用G-S迭代法解此方程组是收敛的。因为系数矩阵4=181是严格对角占优

3-315,

阵，所以G-S迭代法收敛。

一匹+8X2+0•/=7

15.设有方程组v