最新数值分析简明教程(第二版)课后习题答案

1、0.1算法1、(p.11，题1)用二分法求方程X -X-1 =0在1,2内的近似根，要求误差不超过10-3.【解】由二分法的误差估计式 &十黔=珀_；=10，得到2k1 _1000.两端取自然对数得k3-仁8.96，因此取k=9，即至少需In 2二分9次.求解过程见下表。kakbkXkf(xQ符号0121.5+1234567892、( p.11，题2)证明方程f(x)二ex，10x-2在区间0,1内有唯一个实根；使用二分法求这一实根，要求误差不超过-10。2【解】 由于f(x)二ex，10x-2，则f (x)在区间0,1上连续，且f(0) =e 10 0 - 2 = 1 ： 0，f (1) -

2、e1 10 1 2 = e 8 0 ,即卩 f(0) f(1) ： 0 , 由连续函数的介值定理知，f(x)在区间0,1上至少有一个零点.又f(x)二ex 10 0，即f (x)在区间0,1上是单调的，故f(x)在区间0,1内 有唯一实根.由二分法的误差估计式|x*-Xk匸唁10，得到2k -100.2 2 2两端取自然对数得k 一刀俱、2 3.3219二6.6438，因此取k = 7，即至少需二分In 27次.求解过程见下表。kakbkXkf(Xk)符号0010.512345670.2误差1. （p.12，题 8）已知 e=2.71828 ；试问其近似值 花=2.7 , x2 = 2.71

3、, X2=2.71 ,X3 =2.718各有几位有效数字？并给出它们的相对误差限。【解】有效数字：1 1因为| e-x1 | = 0.01828- 0.0510 ，所以 =2.7有两位有效数字；21 /因为|e-X2 |= 0.00828：0.0510 ，所以x 2.71亦有两位有效数字;，所以X3 =2.718有四位有效数字;1 3因为 |ex3 |=0.00028： 0.000510r1r2r3Xi=|e - X2 |X2|e-X3 |X3评 （1）（2）：竺85% ;2.7:005 =1.85% ；2.71=0.0184%。2.718经四舍五入得到的近似数，其所有数字均为有效数字;近似数

4、的所有数字并非都是有效数字2. （ p.12 ,题9）设 =2.72 , x2 =2.71828, x 0.0718均为经过四舍五入得出的近似值，试指明它们的绝对误差（限）与相对误差（限）。z1 0.005q.05, y Er1.8410 ；:：0000005 10$ ；2.718280.00005, -r 3 :Xq四眄 6.96 10,；0.0718评经四舍五入得到的近似数，其绝对误差限为其末位数字所在位的半个单位3. （ p.12，题 10）已知 x1 =1.42 , x -0.0184 , x3 =184 10,的绝对误差限均为0.5 10 一，问它们各有几位有效数字?【解】 由绝对误

5、差限均为0.5 10,知有效数字应从小数点后两位算起，故x 1.42，有三位；x2 二-0.0184 有一位；而 x3 =184 10，=0.0184，也是有一位。1.1泰勒插值和拉格朗日插值1、( p.54，习题1)求作f(x) =si nx在节点x0 =0的5次泰勒插值多项式 p5(x)，并计算P5 (0.3367)和估计插值误差，最后将P5 (0.5)有效数值与精确解进行比较。【解】由 f (x) =sin x，求得 f (x)二 cosx ； f (x)二-sinx ； f (x)二 si nx ； f(5)(x)二 cosx ； f(6)(x) - -sinx,所以P5(x)二 f(

6、Xo) f ()(Xo)(X-Xo)2! (X-Xo)f(3) (x) - - cosx ；f (X。)5!(x- x)5二 f(0) f (0)x f Qx2(0)X52!5!13*15=XXX3!5!插值误差：R5(x).|f(6)( )|(X0)|sin( )|(X0)1x6，若 x = 0.5,则6! 6! 6!0.33670.3367p5(0.3367) =0.33670.3303742887，而3!5!0.33676上二R5(0.3367)2.02 10=0.5 10 ，精度到小数点后 5位,6!故取 p5(0.3367) =0.33037,与精确值 f(0.3367) =sin(

7、0.3367) = 0.330374191 相比 较，在插值误差的精度内完全吻合!2、( p.55，题12)给定节点x0二-1,X1 =1,x2 =3,X3 =4，试分别对下列函数导出拉格朗日余项：(1) f(x) =4x3 -3x 2 ；(2) f(x)=x4-2x3f心3【解】依题意，n = 3，拉格朗日余项公式为&(x)(x-xj4!(1) f(4)(x0 t R3(x) =0 ；(2) 因为f(x) =4!，所以f(4)(料R3(x)(x 1)(x-1)(x-3)(x-4) =(x 1)(x-1)(x-3)(x-4)4!3、( p.55，题13)依据下列数据表，试用线性插值和抛物线插值

8、分别计算sin(0.3367)的近似值并估计误差。i012Xi0.320.340.36si n(N)0.3145670.3334870.352274【解】依题意，n =3,拉格朗日余项公式为Rj(x)二f)()3(1)线性插值因为x =0.3367在节点x和X1之间，先估计误差Ri(X)f”()=gx)iX0)(Xr0.012(x)(X -Xo)(X -Xi)(X -X2)(X - Xo)(Xi - x)(x - X2)3!63 0.013610max(x -XoXj -x)(x2 -x)6yy=(x-xo)(x-xi)(x-x2)Max=3(x i-Xo)3/8XoXiX2抛物线插值公式为:

9、P2(x)sin(X0)(X。-Xi)(X0 -X2)1灯&讥)sin(xj(Xi -X)(Xi -X2)(X-Xi)(X-X。)(X2 -Xi)(X2 -X。)sin(x2)i (Xi - x)(X2 - x)0.022 IL2sin(x0) (xx0)(x2x)sin(xi)-(% -x)(x-x。)2sin (x2)F2 (0.3367)匹 2 3.8445 sin(0.32) 38.9ii sin(0.34) - 2.7555 sin(0.36)1 0.02*1 n 3.8445 sin(0.32) 38.9ii sin(0.34) - 2.7555 sin(0.36)-0.33037

10、439 0.02经四舍五入后得：P2 (0.3367) =0.330374，与 sin(0.3367) = 0.330374i9i精确值相比较，在插值误差范围内完全吻合！i.3分段插值与样条函数i、( p.56，习题33)设分段多项式X3 X22x3 bx2 ex -1是以0,1,2为节点的三次样条函数，试确定系数b，c的值.【解】依题意，要求S(x)在x=1节点函数值连续：s_(1)=13 12=2 13 b 12c 1-1 二 S .(1),即：b c = 1(1)一阶导数连续：S (1) =3 12 2 1 = 6 122 b 1 c 二 S (1),即：2b c = 1(2)解方程组(

11、1)和(2),得b = 2,c = 3，即f 32S(x) = x +x2x3 2x2 +3x-1由于 S(1) = 3汇2x1+2 = 6x2x12x2 = S；(1),所以 S(x)在 x=1 节点的二阶 导数亦连续。12、已知函数y2 的一组数据，x = 0, x1=1,x2 =2 和y0= 1, yj = 0.5,y2= 0.2，1 +x(1) 求其分段线性插值函数；(2) 计算f (1.5)的近似值，并根据余项表达式估计误差。【解】(1)依题意，将x分为0,1和1,2两段，对应的插值函数为S(x)和S2(x)，禾U用拉格朗日线性插值公式，求得x-xXX0X1X0cl CL /S1 (

12、x) y0- y110.5=_0.5x 1;x X1X1 x0110_XX2X%X2X1cc CC CCS2(x)- y1 y20.50.2 = -0.3x 0.8X - X2X2 -x/1-22-11(2) f (1.5)2 ： 0.30769230769 ，而 S2(1.5)= -0.3 1.5 0.8 = 0.35,1+1.5f(1)(x)-2x2 2(1 X )f(2)(x) =-2(1 -3x2)(1x2)3f(3)(x)二24x(1-x2)(1 x2)4实际误差为：| f (1.5)-S2 (1.5)| = 0.0423乞 0.05。知M 2 = f (2) (1) =0.5，则余

13、项表达式|f(J|R(x2!M 224|(x_1)(x_2)|20.5 = 0.5470625 乞 0.51.4曲线拟合1、( p.57，习题35)用最小二乘法解下列超定方程组:”2x +4y =113x _5y =3x 2y = 62x y = 7【解】构造残差平方和函数如下：2 2 2 2Q(x, y)=(2x 4y-11)(3x_5y_3) (x 2y-6)(2x y - 7),分别就Q对x和：Q(x,y) =0 :2+y求偏导数，并令其为零:x：Q(x,y) =0 :解方程组(1)和46 1748x =2736x - y = 17- 3x 46y 二 48(2) , 得:3.04029

14、,(1)，6 483 176 48 3： 1.241762732、（ p.57，习题37）用最小二乘法求形如y = a bx2的多项式，使之与下列数据相拟合。【解】令X = X2，贝y y = a bX为线性拟合，根据公式（p.39,公式43），取m=2, a仁0,N=5,求得55a + bZ Xi5=5a + bZ x5i2 =送 yi(1)ii丝iT555555a迟Xi2+ b瓦 Xi = aZXi2 +b工 x4Xi)/i =瓦 xji=1i =1i吕i=1i=1依据上式中的求和项，列出下表xyiXi (=x i2)Xi2(=xi4)Xi yi (=xi2yi)1919361130321

15、68592532.362539062520187.53149961923521470893873.314442085136105845.24497.819363748096189340.8157271.453277277699369321.5将所求得的系数代入方程组（1）和（2）,得5a0 十 5327b = 271.4Q327a0 + 7277699b = 369321.5(2)271.4 7277699 -369321.5 53277791878.1a0.97258 ；5 乂 7277699 - 5327 汇 53278011566.5 369321.5 -5327 271.4400859

16、.7b0.05004 ；5 7277699-5327 53278011566即: y = 0.97258 0.05004x2。2.1机械求积和插值求积1、( p.94，习题3)确定下列求积公式中的待定参数，使其代数精度尽量高，并指明求积公 式所具有的代数精度：h(1) f (x)dx：Ao f(_h)A f (0)人 f (h)；i113(2) 0 f(x)dx ： Aof(；)几口；)A2fH)；042411(3) f(x)dx ： f (0) A0f(x0) o04【解】(1 )令f(x) =1,X,X2时等式精确成立，可列出如下方程组：A0 +A +A2 =2h* A0 + A2 = 0

17、+ A2= 2 h 3h4解得：A0 = A2= , A1= h33(1)f (x)dx 茫一f (h) + 4f (0) + f(h)，可以3f(X)=X4不成立，故公式(1)具有3次代数精度。(2)令f(X)=1,X,X2时等式精确成立，可列出如下方程组：即:验证，对f(x) =x3公式亦成立，而对Ao+A2 =1代 +2A +3A2 =23A0 +12A +27A2 =1621(1)(2)(3)11132f ( )-f( )2f()，可以3424验证，对解得：A0 = A2, A1，即： f (x)dx33f(x) =X3公式亦成立，而对 f(x) =X4不成立，故公式(2)具有3次代数

18、精度。(3)令f (x) =1,x时等式精确成立，可解得:A42x。-3即:f(X)二 X1130 f(x)dx ： ： f(0)-3不成立，故公式(3) 2兰f ()，可以验证，对43具有2次代数精度。f(x) =x2公式亦成立，而对1X0, X14求积公式，并指明该求积公式的代数精度。【解】依题意，先求插值求积系数：31 x -1 4 。1_3442、( p.95 ,习题6)给定求积节点试构造计算积分41I = o f (x)dx的插值型13dx0 X0 _ X1dx = -2 (1x2213、 4x) 什 01 x _Xo1 x-1A10X1X0dxZ 12(2x插值求积公式:1o f

19、(x)dx 八k =01 1Akf(xkH2f(4)3f(4)当f(X)=1，左边=0 f (x)dx=1 ；右边=gx1+x1=1 ； 左 =右;当1f(x) =x，左边=f(x)dx =1 ；左=右;当2 1f (x) = x ,左边= f (x)dx =右边J丄 1 22 16 2 1616 右;故该插值求积公式具有一次代数精度。2.2梯形公式和Simpson公式1、(p.95，习题9)设已给出f(x)=1in4x的数据表,x0.000.250.500.751.00f(x)1.000 001.655 341.551 521.066 660.721 591分别用复化梯形法与复化辛普生法求积

20、分I二o f (x) dx的近似值。【解】(1 )用复化梯形法：a =0,b =1, n =5,h二 1 =0.25n4Thh丫T5f(xQ m f(a) 2、f(xQ f(b)心22心T5 f(0.00) 2 f (0.25) f (0.50) f (0.75) f (1.00)2T5 =0.125 1.000002 (1.655341.551521.06666)0.72159T5 =1.28358(2 )用复化辛普生法：a = 0, b = 1, n = 2,h = b = 1 = 0.5n 2hnJ)W)匚n -4f(xk2)2、 f(xQ f(b)k T0 5S2 f (0.00) 4

21、 f (0.25) f (0.75)2 f(0.50) f(1.00)6S21 1.00000 10.888 3.10304 0.72159 1.30939122、(p.95，习题10)设用复化梯形法计算积分| = fexdx，为使截断误差不超过,问应当划分区间【0,1】为多少等分？如果改用复化辛普生法呢？【解】(1)用复化梯形法，a = 0,b = 1, f (x) = f(x) = f (x) = ex，设需划分n等分,则其截断误差表达式为:|Rt|=|I -几| =(b -a)312n2(1-0)3 ；12n3 e ；1依题意，要求I Rt N 102n2 -6e 105a：212.84

22、9，可取 n =213。(2)用复化辛普生法a = 0,b = 1, f (x) = f (x) = f(x) = ex，截断误差表达式为:|Rs|I -Sn |=(b-a)5180(2n)45maxf( r益討歳；1 5依题意，要求I Rs I10 ，即25 1 10* 二 n4_e103.70666，可取 n=4，划分 8 等分。2880n4214402.3数值微分1、( p.96，习题24)导出三点公式(51)、(52)和(53)的余项表达式1f(Xo)-3f(Xo) 4f(X1)-f(X2)(51)2h1f(xj-f(x) fg)(52)2h1fX)f(x)-4f(X1)3f(X2)(

23、53)2h【解】如果只求节点上的导数值，利用插值型求导公式得到的余项表达式为严* ) nR(xQ 二 f(xQ - p(xj严丨丨 X -Xj)(n +1)!jTj我由三点公式(51)、(52)和(53)可知，n =2,h-x0 =x2 -花，贝yR(X0)=畀 K (X0-Xjt3a(X0-X1)(X0-xA%h2(2 1)! j3!3R(xJR(X2)f(2 d)( 1)-2(X1 - Xj )= (2 1)! jjjj/f (2 巧()2(x2- x j )(21)! j 卫jj老2、( p.96，习题f(1)3!(x x0)(x1 x2)3!3 - (x2 - x0 )(x2 - x1

24、)= 3!x1.01.11.2f(x)0.25000.22680.2066的数据表,试用三点公式计算 f(1.0), f(1.1), f(1.2)的值，并估计误差。25)设已给出f(x)tx)2f( 2)h23【解】已知 x0 = 1.0, x1 = 1.1, x2 = 1.2, h = xi - x()= x2 - xi = 0.1,用三点公式计算微商：1 1f(1.0)-3f (1.0) 4f (1.1) - f (1.2)-3 0.2500 4 0.2268 - 0.2066 = -0.24702h2x0.11 1f(1.1) ：_f (1.0)f(1.2) =-0.2500 0.206

25、6 =0.21702h2 7.11 1f(1.2) f (1.0)-4f(1.1) 3f (1.2)0.2500-4 0.2268 3 0.2066 - -0.18702h2 7.11-2f(x)2 = f(x)齐：f(x)(1 +x)2(1 +x)3用余项表达式计算误R(1.0) f( 0)-24(1 x)5 ，R(1.1)二h23_f( 1)屮3!fh2差2-24 0.12as &5-0.00253、( p.96，习题 26)3(1 1.0)24 0.1253!(1 1.0)-24 0.12r3(1 1.1):0.00125-0.04967设f(x)二sinx，分别取步长h =0.1,0.

26、01,0.001，用中点公式(52)计算f(0.8)的值，令中间数据保留小数点后第6位。【解】中心差商公式：f(a) ： f(a h) f (a ，截断误差：R(h)二2h。可3!见步长(1)h越小，截断误差亦越小。=0.1, x0 二 0.8 - h = 0.7, x2 二 0.8 h = 0.9,则1 1sin(0.9) -sin(0.7)0.783327 - 0.644218 ： 0.695545 ；2h2 0.1f(0.8)h = 0.01, x0 二 0.8 - h 二 0.79, x2 二 0.8 h = 0.81,则1 1sin( 0.81) -sin(0.79)0.724287

27、 -0.710353 : 0.69672h2 0.01f(0.8)h - O.OO1,x0 = 0.8 - h = 0.799, x2 =0.8 h = 0.801,贝y1 1f(0.8) ： sin(0.801) -sin(0.799) ：0.718052 - 0.716659 ： 0.69652h2x0.01而精确值 f(0.8) =cos(0.8) =0.6967067，可见当 h =0.01时得到的误差最小。在h = 0.001时反而误差增大的原因是 f (0.8 h)与f (0.8 - h)很接近，直接相减会造成有效 数字的严重损失。因此，从舍入误差的角度看，步长不宜太小。3.1 E

28、uler 格式1、( p.124，题1)列出求解下列初值问题的欧拉格式2 2(1)y= x - y (0_x_0.4)， y(0) =1，取 h = 0.2 ；y=必 | +必 (1 EXE1.2)，y(0) =1，取 h =0.2 ；IX丿X【解】(1)yn 1 二 hyn 二 y h(x； - y；) = y“0.2 (x； - y；)；2 2(2)yn 1 = ynh(芻 仏)=yn 0.2 (聖 上)。XnXnXnXn2、( p.124，题2)取h = 0.2，用欧拉方法求解初值问题y - - y - xy2 (0 _ x _ 0.6)，y(0) =1。【解】欧拉格式：yn 1 二 y

29、nhyn二ynh(-yn -Xny；) = yn 0.2(-yn-Xnyi)；化简后，yn1 =.8yn -0.2xny2，计算结果见下表。n0123Xn0.00.20.40.6yn1.00.80.61440.46131 23、( p.124，题3)取h =0.1，用欧拉方法求解初值问题y2-2y2(0_x_4)，1 + xy(0) =0。并与精确解y 程 比较计算结果。1 +x【解】欧拉格式：1 2 1 2hyr h,址-2%)-% 0.2 -着2*化简后，yn j = yn -0.4y； -0电，计算结果见下表。1 +Xn1、( p.124，题7)用改进的欧拉方法求解上述题2，并比较计算结

30、果。ypy=ynhf (Xn, yn) = Yn 皿-丫.二 yn hf (Xn, yp)二 yn h(-yp (yp yc)22 2-Xnyn) = 0.8yn -0.2Xnyn2 2-Xnyp)二 yn -0.2 (ypx.yp)。n0123Xn0.00.20.40.6yp1.00.67300.51470.3941yc0.760.70920.55640.4319yn0.880.69110.53560.413计算结果见下表。yn 1【解】公式：因为 y= f (x,y) - -y - xy2(0 乞 x 0.6)，h = 0.2，且 y(0) =1，则改进的欧拉与原结果比较见下表n0123X

31、n0.00.20.40.6yn1.00.80.61440.4613yn(改进)0.880.69110.53560.4133.3龙格-库塔方法1、( p.124，题11)用四阶经典的龙格-库塔方法求解初值问题y = 8-3y , y(0)=2，试取步长h =0.2计算y(0.4)的近似值，要求小数点后保留4位数字。【解】 四阶经典的龙格-库塔方法公式:hVnVn -(K1 2K2 2K3 K4)6K1 二二 f(Xny)k2=f (x1 , ynK1)n+22K3 =f (xn甲，yn22K2)K4二 f(Xn十，yn+ hKa)列表求得y(0.4)如下:nXnyn00.02.00010.22.

32、300420.42.46544.1迭代法及收敛定理1、( p.153，题1)试取Xo =1，用迭代公式Xk 1 =202x2 2xk 10(k =0,1,2,)，求方程x3 2x210x-20=0的根，要求准确到10。【解】 迭代计算结果列于下表kXk|Xk-Xk-110.001kXk|Xk-Xk-110.00111.538460.53846N61.365930.00937N21 1.295020.24344N :71.370090.00416N31.401820.10680N81.368240.00185N41.354210.04761N91.369060.00082Y51.375300.0

33、2109N因为 | X9 -x8 卜 0.00082 :： 10，所以 X ： x9 =1.36906。12、 (p.153，题2)证明方程xcosx有且仅有一实根。试确定这样的区间a,b，使迭21代过程xk彳cosxk对x0 a,b均收敛。2111 1【证明】设：g(x) cosx，则当 x R时，g(x) cosx,，且一阶导数222 21 1 1 1 g(x)二Sinx连续，|g(x)|=|in1，所以迭代过程Xk 1 二?cosxk对1x0 R均收敛。(压缩映像定理)，方程x cosx有且仅有一实根。证毕23、 ( p.153，题4)证明迭代过程xk 丁 =乞丄对任意初值x0 1均收敛

34、于-2。2 Xk【证明】设：g(xX 1，对于任意x 1，因为-2 X-2，所以g(x) ,2。2 x2 x V 2 x111x1一阶导数g(x)21，根据压缩映像定理，迭代公式Xk 1 k 对任意2x22xkx 1初值x01均收敛。假设lim xk = x ，对迭代式xk 1 k 两边取极限，则有十2XkX”亡f，则X” 2讥，解得- 2，因X”2不在X 1范围内，须舍去。故x 2。证毕4.2牛顿迭代法4位有效数字:1、( p.154，题17)试用牛顿迭代法求下列方程的根，要求计算结果有(1)x _3x_1=0， x0 = 2(2)x2 -3x _ex 2=0， x0 = 1【解】(1)Xk

35、 1二 Xk设 f (x)二 x3 -3x -1，则3f (Xk)Xk -3XkXkf(Xk)2f(x) = 3x -3，牛顿迭代公式：=2xk 1(k二0,1,2,)，迭代计算过3x： 一 33(x： -1)kXk:k-Xk-110.0001kXk|Xk-Xk-110.00110.268940.73106N30.257530.00014N20.257390.01155N40.257530.00000YXkf(Xk)、， Xk 3Xk e +2xk 1 = xkxk -f (xk)2Xk -3-e，迭代计算过程见下列表。2、( p.154，题18)应用牛顿法于方程x3 - a = 0，导出求立

36、方根3 a(a 0)的迭代公式,并证明该迭代公式具有二阶收敛性。【证明】(1)设：f (x) = x3 - a，2f(x) = 3x，对任意X 0，牛顿迭代公式33xk -a2xk akk= 0,1,2,3x：3xi(2)由以上迭代公式，有:kXmKHk/Vg3-a+ 23X2-g3(1 予xMXk ! -X =g(Xk) -g(x ) =g(x )(Xk-x )豊Xk -x )22!程见下列表。kXk:k-Xk-110.0001kXk|Xk-Xk-110.000111.888890.11111N31.879390.00006Y21.879450.00944N :因为 | x3 -x2 |：

37、0.00006 ： 10鼻，所以 xx3 =1.879。(2)设 f (x) = x2 -3x-ex 2，则 f(x) = 2x-3-ex，牛顿迭代公式：2Xk _3_eXk=xj/(Xk-0-2化=0,1,2,)所以 x” : x4 =0.2575。因为 |x3 x2 I： 0.00000 ： 10*km(:二g”2；J：，可见该迭代公式具有二阶收敛性。证毕5.1线性方程组迭代公式1、( p.17O，题1)用雅可比迭代与高斯-赛德尔迭代求解方程组:3xi + X2 = 2Xi + 2x2 = 1要求结果有3位有效数字。X1(k 1)-X2k)【解】 雅可比迭代公式:=1(2 _X2)迭代计算

38、结果列于下表。X2(k 1)X1(k)1(k)VX1)k(k)X1()x2k)i(k)(k)-1 X；丿- X；丿 |Ix2k)-x2k4| 0.0005?000-12/31/22/31/2N21/21/61/61/3N311/181/41/91/12N47/127/361/361/18N50.601850.208330.018520.01389N60.597220.199080.004630.00925N70.600310.201390.003090.00231N80.599540.199850.000770.00154N90.600050.200230.000510.00038N100.5

39、99920.199980.000030.00025YX； x；10)0.600; X2x210)0.200 ；由上表可见，所求根皆为小数点后第1位不为零的小数，要取3位有效数，则误差限为1 10 。2X1(宀=x2k)+2 J(2-x2k)高斯-赛德尔迭代公式：彳33 3，迭代计算结果列于下表。x2k+-1X1(k+1(x2k). 2 2 6kX1(k)x2k)|xr-X1(kj)|x2k)-x2Z| 0.0005?000-12/31/62/31/6N20.61110.1944N30.60190.19910.00920.0047N40.60030.19990.00160.0008N50.600

40、00.19990.00030.0000Yx/ xj0.600;x2x25)0.200 ；12、( p.171，题7)取用-1.25，用松弛法求解下列方程组，要求精度为10 *。24x1 3x2 = 163x1 4x2 -x3 二 20-x2 4x3 = -12【解】欧先写出高斯-赛德尔迭代:(k出)X13 (k)-X24(k 1)X2” 5存2(1) (k 1)X31 (k)2_3=2x2k)x3k)6416引入松弛因子,(k 1)X1=(1(k)-)X1(k1)1(k)X1(k 1)1x2k 1)=(1k)k)(k 1)X2x3k 1)=(1-jx3k)(k1)3k)(k 1)X3(k 1)

41、1(k)15(k)X1X1-X25416(k 1)X2二29 x2k)-53k)-6416 X2(k 1)45(k)11(k)25X3二X2X3-256648(2)，并化简将方程组(1)代入计算结果见下表。k(k) x；)x2k)x3k)i (k)(k4) |X1()-x；) |x2k)-x2k)i|x3k)-x3k)ie?0000-152.5-3.12552.53.125N21.406252.65625-2.14844N32.158203.03223-2.28882N41.611733.15872-2.19860N51.635773.24423-2.19187N61.549593.28508

42、-2.17800N71.532843.30793-2.17320N81.515613.31978-2.17001N91.508803.32615-2.16847N01.504533.32951-2.16762N11.502453.33130-2.16717N21.501293.33225-2.16694N31.500693.33276-2.16672N41.500373.33306-2.16676N51.500163.33318-2.16670N61.500103.33325-2.16668N71.500053.33329-2.166680.000050.000040.00000Y迭代解：xx1(17) 1.5001,x2 =x217) ： 3.3333,x x317) : 2.1667.精确解:3Xi1.5,2X2103.3333,