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

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

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

3、；试冋其近似值 x = 2-7 , x2 = 2.71 , x2=2.71 , X3 = 2.718各有几位有效数字？并给出它们的相对误差限。【解】有效数字：1 1因为|e-x1 | = 0.01828：0.0510 ，所以x1 =2.7有两位有效数字；21 因为|e-x2 | = 0.00828：0.0510 ，所以x2 =2.71亦有两位有效数字；21 3因为|e-x3 |=0.00028：0.000510 ，所以x2.718有四位有效数字;2|e X1 |X10.05<2.7-1.85%| e - X2 |0.05x22.71-1.85%0005X32.718= 0.0184%精品

4、文档评 （1）经四舍五入得到的近似数；其所有数字均为有效数字;（2）近似数的所有数字并非都是有效数字2. （ p.12，题9）设捲=2.72 ； x2 =2.71828 ； x3 =0.0718均为经过四舍五入得出的近似值，试指明它们的绝对误差（限）与相对误差（限）。【解】“=0.005 ；r < 010051.84 10"；x12.720.000005 ；20.000005<x22.718281.84 10 ；邑 0.00005/；3 =0.00005 ；r3 36.96 10；x30.0718评经四舍五入得到的近似数，其绝对误差限为其末位数字所在位的半个单位3. （

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

6、x，求得 f (x) = cosx ； f (x)= f (x) = sinx ； f(5)-sin x ；(x)二 cosx ； f(6)(x)二-si nx，所以 f(2)(X0)(2)f(3)(x) = - cosx ；(X -X。)2P5(x)二 f(x°) f (x°)(x - x°)2!f (0)x2 .f(5)(0)5!-f(0) f(1)(0)x(5)5-x2!f (X。)口 (X-X°)55!13*15=X XX3!5!| f(6)(巴)| 插值误差：R5(x)f ( )|(xx°)6 6!3 0.33673 p5(0.336

7、7) =0.3367 -0.33676R5(0.3367):二刨)l(x _X0)6 岂丄 X6，若 x 二 0.5，则 6! 6!50.336750.3303742887，而5!652.02 10< 0.5 10 ，精度到小数点后 5位,6!故取 p5(0.3367) =0.33037，与精确值 f (0.3367) = sin(0.3367) = 0.330374191 相比 较，在插值误差的精度内完全吻合！2、( p.55，题12)给定节点x°二- 1,X1 =1,X2 =3,X3 =4，试分别对下列函数导出拉格朗 日余项：(1) f (x) = 4x3 -3x 2 ；4

8、3(2) f (x) = x -2xf(4)(E)3【解】依题意，n =3,拉格朗日余项公式为R3(x)(x-Xi)4!1=0(1) f (4)(xH0 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)R3(x)二亠 i：(x_Xi)4!3、( p.55，题13)依据下列数据表，试用线性插值和抛物线插值分别计算sin(0.3367)的近似值并估计误差。i012Xi0.320.340.36sin(xi)0.3145670.3334870.352274【解】依题意，n

9、 =3,拉格朗日余项公式为(1) 线性插值因为x =0.3367在节点Xo和Xi之间，先估计误差Ri(x)f”()普一 “严2!R(x)gsin(x°) 二 sin(xj =X1 _ X0X1 _X°Xo _Xi1(x _ x0 )sing )(捲 _ x)sin(x0) 1二丄104 ;须保留到小数点后 4为，计算过程多余两位。2(X1-Xo)2/4y=(X-X o)(X-Xi)Xi0二R(x)(0.3367 -0.32)sin(0.34)(0.34 - 0.3367)sin(0.32)0.020.0167 sin(0.34)0.0033 sin(0.32)10.02 0

10、.3304(2)抛物线插值插值误差：f'')_cosG)R2 (x)(X 'XoX Xi )(X X2 )(X Xo )(X1 X)(X X2 )363maX( x-Xo)(Xi - x)(X2 - x) 3 0.011 0_6- 6 、 6 _2yy=(x-xo)(x-xi)(x-x2)Max=3(x i-Xo)3/8XoXiX2x抛物线插值公式为:(X -Xi)(X -X2)P2(X)sin(xo)十(X Xo)(x-X2)sin(捲)+ 以一恥-冷)前(乂?)(X。Xi)(X°X2)(Xi- Xo)(Xi-X2)(X2Xi)(X2 Xo)1 (Xi _

11、X)(X2 _ x)0.022 IL 2sin(Xo) (x -x°)(X2x)sin(xj -(Xi -x)(x-x。)2sin (x2)F2 (0.3367)匹 2 3.8445 sin(0.32) 38.911 sin(0.34) - 2.7555 sin(0.36)1 0.023.8445 sin(0.32) 38.911 sin(0.34)-2.7555 sin(0.36)- 0.330374390.02经四舍五入后得：P2 (0.3367) = 0.330374，与 sin(0.3367)= 0.330374191 精确值相比较，在插值误差范围内完全吻合！1.3分段插值与

12、样条函数x3 +x21、( p.56，习题33)设分段多项式S(x)=°Nx3 +bx2 +cx1是以0,1,2为节点的三次样条函数，试确定系数b，c的值.【解】依题意，要求S(x)在x=1节点函数值连续：s_(1)=13 12=2 13 b 12 c 1一1 二 s .(1),即：b c =1(1)一阶导数连续：S_(1)=3 12 2 1=6 12 2 b 1 c=S(1),即: 2b c = 一1(2)解方程组(1)和(2),得b = _2,c = 3，即S(x)r 32x + x2x3 2x2 +3x1由于 S(1) =3 2 1*2=6 2 1 -22 =S".(

13、1)，所以 S(x)在 x=1 节点的二阶导数亦连续。12、已知函数 y2 的一组数据，x0 = 0, x1 = 1, x2 =2 和 y0 = 1, y1 = 0.5, y2 = 0.2，1 +x(1)求其分段线性插值函数；(2)计算f (1.5)的近似值，并根据余项表达式估计误差【解】(1)依题意，将x分为0,1和1,2两段，对应的插值函数为S1(x)和S2(x)，禾U用拉格朗日线性插值公式，求得X -X0y1XI _ X01 .0.5 二-0.5x 1 -0-1 1-0 'S2(x)X一紅 y1X1 X2亠y2X2 - X1x _ 2x _ 1鼻二 050.2 一0.3x0.81

14、-22 -10.30769230769 ，而 S2(1.5) = -0.3 1.5 0.8 = 0.35 ,实际误差为：| f (1.5)S2(1.5) 1= 0.0423 乞 0.05。2 2 ，(1 X )f (XL2。®2)(1X2)32f(3)(x24x(1-x4),可(1x2)4|f (J|R(x"2!知M 2 = f (2)(1) =0.5，则余项表达式M 224| (x 1)(x -2) |2 0.5 = 0.5 = 0.0625 乞 0.51.4曲线拟合1、( p.57，习题35)用最小二乘法解下列超定方程组:"2x + 4y =113x -5y

15、=3x 2y = 62x y =7【解】构造残差平方和函数如下：2 2 2 2 Q（x, y） = （2x 4y11）（3x5y3）（x 2y6）（2x y - 7）,分别就Q对x和y求偏导数，并令其为零：Q(x,y) _0 :6xy = 17(1),dx：Q(x,y) 0:-3x 46y = 48(2),解方程组（1）和(2),得46 17486 483 17x =：3.04029,y1.241762732732、（ p.57，习题37）用最小二乘法求形如、二a bx2的多项式，使之与下列数据相拟合。【解】令X = X2，贝y y = a bX为线性拟合，根据公式（p.39,公式43），取m

16、=2, a仁0,N=5,求得55a +b迟 Xi5=5a + bZ x5i2=I： Vi(1)1i丝=1555555a送Xi2+ b瓦 X i = aZN2 + b工 Xi4=Xi）A =送 xyii=1i =1imi=1i=1依据上式中的求和项，列出下表xXi (=x i2)Xi2(=xi4)Xi yi (二Xi2y)191936113032168592532.362539062520187.53149961923521470893873.314442085136105845.24497.819363748096189340.8157271.453277277699369321.5将所求得的

17、系数代入方程组（1）和（2）,得5a0 十 5327b = 271.4(1)Q327a0 + 72776993 = 369321.5(2)271.4 7277699-369321.5 53277791878.1 小小 c ca0.97258 ；5 汉 7277699 - 5327 汉 53278011566u 5 369321.5 -5327 271.4400859.7b0.05004 ；5 7277699-5327 532780115662即: y =0.97258 0.05004x。精品文档1 ；精品文档2.1机械求积和插值求积1、( p.94，习题3)确定下列求积公式中的待定参数，使其代

18、数精度尽量高，并指明求积公 式所具有的代数精度：h(1) f(x)dx : A0f(-h) A1f(0) A2 f (h)；i113(2) 0f(x)dx ： Aof? AG?)人2七)；1 1 o f(x)dx 行 f(0) A0 f (x0) o【解】(1 )令f (x) =1, X,X1 2时等式精确成立，可列出如下方程组：Ao A A2 = 2h(1)-A A2 = 02AA2 二 一 hhhf (x)dx ： f (h) 4f (0)f(h)，可以3f(x) =x4不成立，故公式(1)具有3次代数精度。(2)令f (x) =1, x,x2时等式精确成立，可列出如下方程组：h4解得：A

19、 = A2, A1h ,33验证，对f(x) =x3公式亦成立，而对即:A1A2 =1代 +2A +3A2 =2(3)解得：2111113A0=A2S，A13，即：0f(x)dx：护(/-电)4求积公式，并指明该求积公式的代数精度。【解】依题意，先求插值求积系数： 七)，可以3A0 +12 A +27 A2 =16验证，对f(x) =x3公式亦成立，而对 f (x) = X4不成立，故公式(2)具有3次代数精度。(3)令f (x) =1,x时等式精确成立，可解得:A042x。-31132即： ° f (x)dx ” f (0)盲f (?)，可以验证，对f(x) =x3不成立，故公式(

20、3)具有2次代数精度。f(x) =x2公式亦成立，而对试构造计算积分l=ff(x)dx的插值型40dx =31 x -1 40 1 3dx 2 Qx223、；x)A°1 x _x°dxxi x°1 x-1-11 1 =1；左=右;2 2当f(x)=x，左边=(f(x)dx = 1x右边=1 11 ?2424i ；左 =右;当f (x) =x2,左边=f f(x)dx =3lx3右边=丄丄J鸟二2 16 2 1616 "右;插值求积公式:113f(x)dx 八 Akf(xkH-f(-) -f(；)1当f (x) = 1,左边=0 f (x)dx = 1 ；

21、右边=故该插值求积公式具有一次代数精度。2.2梯形公式和Simpson公式1、( p.95，习题9)设已给出f(x)=1 esin4x的数据表,x0.000.250.500.751.00f(x)1.000 001.655 341.551 521.066 660.721 591分别用复化梯形法与复化辛普生法求积分I f (X)dx的近似值。【解】 (1)用复化梯形法b a a = 0, b = 1, n = 5, h =nn 4 hT5 八；f(Xk) f(Xkd心210.254hn 4f(a) 2、f(xQ f(b)2k 40.25丄丄丄丄T5f(0.00) 2 f(0.25) f (0.50

22、) f (0.75) f (1.00)2T5 =0.125 1.000002 (1.655341.551521.06666)0.72159T5 = 1.28358(2 )用复化辛普生法b a 1 a = 0, b = 1, n = 2, h0.5n 2n 1 hhnn_1S2 八刃 f(xQ 4f(xkf(Xk.J =：f(a) 4 f(xk2、f(xQf(b)76k 267 k 2k=i0 5S2f (0.00) 4 f(0.25) f (0.75) 2 f (0.50) f (1.00)61S21.00000 10.888 3.10304 0.721591.309391 12、( p.95

23、，习题10)设用复化梯形法计算积分Iexdx，为使截断误差不超过 1 10问应当划分区间【0,1】为多少等分？如果改用复化辛普生法呢？【解】(1)用复化梯形法，a =0,b =1, f (x)二f'(x)二f''(x) =ex，设需划分n等分， 则其截断误差表达式为：|Rt HI-Tn卜(b - a)32 max12n2f''()(1 -0)312n31依题意，要求| Rt |10，即2-< 1 10* 二 n2 _: 212.849，可取 n = 213。(2)用复化辛普生法，12n226a =0,b =1, f (x) = f'(x)二

24、 f “''(x) =ex，截断误差表达式为:|Rs|I -Sn 卜(b-a)5180(2n)4maxf""()(1-0)5ee4 e4 ；2880n42880n41依题意，要求| RS |10*，即2e 1e 1010- n4 _ 370666，可取 n =4,划分 8 等分。2880n4214402.3数值微分1、( p.96，习题24)导出三点公式(51)、(52)和(53)的余项表达式精品文档f'(Xo)：丄_3f (xo) 4f(X1)_ f(X2)(51)2hf'(xi)：丄_f (xo) f(X2)(52)2hf'(X2

25、)：丄f(xo) -4f(xi) 3f(X2)(53)2h【解】如果只求节点上的导数值，利用插值型求导公式得到的余项表达式为 f (n (匕)nR(xjf'(Xk) - p'(Xk)二行Xk-Xj)(n +1)! y由三点公式(51)、(52)和(53)可知，n =2,h =x, -x0 = x2 -x,，贝U 、 fH(J)2 “、f''o)“、“、 f'"o)R(xorj(xo - Xj ) = T (X0 一 X1 )(Xo 一 X2)= T h2f (2 珂】)一 2f'''(,)R"爲肿"呼

26、(牛)h2j=1f(21)( 2)-_2R(X2)=2 ；区7)=亠化十)区7)亠人2(2 1)! j j 3!312、( p.96，习题25)设已给出f(x)2的数据表,(1+x)X1.0f(x)0.25001.11.20.22680.2066试用三点公式计算f'(1.0), f'(1.1), f'(1.2)的值，并估计误差【解】已知X。=1.0,为=1.1, X2=12 h =X1-X0=X2-X1 =0.1,用三点公式计算微商：1 1f'(1.0)-3f (1.0) 4f (1.1) - f (1.2)-3 0.25004 0.2268 - 0.2066

27、- -0.24702h2汇0.11f'(1.1)：f'(1.2)f(x)二1Hf (1.0) f (1.2)-0.25000.2066 =-0.21702h2 0.11 1f (1.0) -4f (1.1) 3f (1.2)0.2500 -4 0.2268 3 0.2066 = -0.18702h2 0.11f'(x)二(1 x)用余项表达式计算误f'''( 0)以h3f'''( ) 2h 3!R(1.2)二 f'''( 2)R(1.0)二R(1.1)二(1 x)主差2-24 0.1立5"

28、_2 3；= f''(x)=(1 x)6-24f'''(X"L，3 h23(1 1.0)224 0.13!(1 1.0)5 -24 0.12 3(1 1.1)-0.0025:0.00125-0.04967精品文档3、( p.96，习题 26)设f(x) =sinx，分别取步长h =0.1,0.01,0.001，用中点公式(52)计算f'(0.8)的值，令中间数据保留小数点后第6位。【解】中心差商公式：f'(a) ： f(a h) f (a，截断误差：R(h) = f (a) h2。可2h3!见步长h越小，截断误差亦越小。(1)

29、h = O1,Xo = 0.8h = 0.7公2 = 0.8 + h = 0.9,则1 1f'(0.8)sin(0.9) - sin(0.7)0.783327 - 0.644218 ： 0.695545 ；2h20.1 h = 0.01, x0 = 0.8h 二 0.79, x2 二 0.8 h = 0.81,则1 1f'(0.8)sin(0.81) -sin(0.79)0.724287 - 0.710353 ： 0.69672h2901(3) h = 0.001,x0 = 0.8 一 h 二 0.799,x2 = 0.8 h 0.801,则1 1f'(0.8)sin(

30、 0.801) -sin( 0.799)0.718052 - 0.716659 ： 0.69652h2 7.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 Euler 格式1、( p.124，题1)列出求解下列初值问题的欧拉格式2 2(1) y' = x2y2(0 乞 XE0.4)，y(0) =1，取 h=0.2 ；(2) y'=上 丫

31、 (1 乞 x iM.2)，y(0) =1，取 h=0.2 ；(X丿x【解】(1)yn 1 二 yn hy'n 二 y“ h(x； - y；) = y“ 0.2 (x； - y；)；2 2(2)Yn 1 二 ynh 台仏)=yn 0.2 (粤仏)。XnXnXnX.2、( p.124，题2)取h =0.2，用欧拉方法求解初值问题y'= -y -xy2(0 _ x_ 0.6)，y(0) =1。【解】欧拉格式：y 1= yn hy'yn h(-yn -Xny；) =yn °2 (-y. -Xny：);化 简后，yn 0.8yn -0.2xny2，计算结果见下表。n0

32、123Xn0.00.20.40.6Ya1.00.80.61440.46133、( p.124，题3)取h =0.1，用欧拉方法求解初值问题y_2y2(0乞x乞4),1 +x x1y(0) =0。并与精确解y2比较计算结果。1 +X21 2 1 2【解】欧拉格式：yn 1 =yn hy'n *n h( y - 2 y n) = y.，0.2 ( -2yn)；1+Xn1+Xn2 0 2化简后，yn 1二yn -0.4y22，计算结果见下表。1 +Xnyp =ynhf (Xn,yn) =ynh(-yn« yc = % +hf(Xn,yp) =yn +h(yp-Xny；) = 0.8

33、yn -0.2Xny；22-Xnyp)二 yn -0.2 (yp x.yp)。1、( p.124，题7)用改进的欧拉方法求解上述题2，并比较计算结果。【解】公式：因为 y' = f (x, y) - -y - xy2(0 _ x _ 0.6)，h =0.2，且 y(0) =1，则改进的欧拉(yp y2n0123Xn0.00.20.40.6yp1.00.67300.51470.3941yc0.760.70920.55640.4319yn0.880.69110.53560.413yn 1计算结果见下表。与原结果比较见下表n0123Xn0.00.20.40.6yn1.00.80.61440.

34、4613yn(改进)0.880.69110.53560.4133.3龙格-库塔方法1、( p.124，题11)用四阶经典的龙格-库塔方法求解初值问题y-8-3y，y(0) = 2，试取步长h =0.2计算y(0.4)的近似值，要求小数点后保留4位数字。【解】 四阶经典的龙格-库塔方法公式:yn .1 二 yn 6(Ki 2K2 2K3 K4)Ki = f (Xn, yn)h« K2 = f (x 1, yn + Ki) 乜 2hK3 = f (x 1 , yn + K2) 性 2iK 4 f(xn4yn + K3)列表求得y(0.4)如下:nXnyn00.02.00010.22.30

35、0420.42.46544.1迭代法及收敛定理(k =0,1,2,)，求201 ( p.153，题2、(p.153，题2)证明方程xcosx有且仅有一实根。试确定这样的区间a,b，使迭)试取Xo = 1，用迭代公式Xk 1 = -2Xk +2Xk +10方程X111 1【证明】设：g(x) cosx，则当xR时，g(x) cosx，,，且一阶导数2 22 21111g'(x) si nx 连续，|g'(x)|=|si nx|1，所以迭代过程Xk 1 cosxk 对22221x0 R均收敛。(压缩映像定理)，方程xcosx有且仅有一实根。 证毕2 2x1 代过程xk 1cosxk

36、对X。 a, b均收敛。2 10x - 20 = 0的根，要求准确到10 "。【解】 迭代计算结果列于下表kXkXk-Xk-11<0.001kXk|Xk-Xk-11<0.00111.538460.53846N61.365930.00937N21.295020.24344N71.370090.00416N31.401820.10680N81.368240.00185N41.354210.04761N91.369060.00082Y51.375300.02109N因为 | X9 - x8 | ： 0.00082 : 10，所以 x ” x9 二 1.36906。3、( p.1

37、53，题4)证明迭代过程二土 丄对任意初值2XkX0 1均收敛于 2。I证明】设:g(x|1，对于任意xxx 1，因为一21=2，所以 g(x) _2 。x1 _2xx1一阶导数g'(x)=1 1彳12x 2根据压缩映像定理，迭代公式初值x01均收敛。假设lim xk=X ，对迭代式xk 1+Xk11两边取极限，则有于占，则&讥2，解得宀朋，因宀2 不在x 1范围内，须舍去。-.2。证毕4.2牛顿迭代法1、(p.154，题17)试用牛顿迭代法求下列方程的根，要求计算结果有4位有效数字:(1)x3 -3x _1 = 0， x0 = 2(2)x2 -3x ex 2=0， x0 =1

38、【解】(1)Xk 1设 f (x) = x3 -3x -1，则f (Xk) _X3 -3XkXk f f'(Xk)3x2-3f'(x) =3x2 _3， -12x3+123(x： -1)牛顿迭代公式:(k =0,1,2,)，迭代计算过kXk|Xk-Xk-11<0.0001kXk|Xk-Xk-11<0.000111.888890.11111N31.879390.00006Y21.879450.00944N程见下列表。因为 | x3 - x2 | ： 0.00006 <10，所以 xx3 = 1.879。2 _3x(2)设 f (x)二 x2 -3x -ex2，则

39、 f' (x) = 2x -3 -ex，牛顿迭代公式：Xk _ e k (Xk -1)-2kk(k = 0,1,2,)2，则 f'(x) = 2x -3 -Xk 3Xk e * +2Xk 1 = XkXkx；f (Xk)2Xk -3- ef (Xk)kXk|Xk-Xk-11<0.0001kXk|Xk-Xk-11<0.00110.268940.73106N30.257530.00014N20.257390.01155N40.257530.00000Y因为 | X3 -X2 I： 0.00000 ： 10 '，所以 x : X4 二 0.2575。，迭代计算过程

40、见下列表。2、( p.154，题18)应用牛顿法于方程X3 - a =0，导出求立方根3 a(a0)的迭代公式,并证明该迭代公式具有二阶收敛性。33Xk -a 2Xk a3x23x2k =0,1,2,(2)由以上迭代公式，有:kXKHk3aa+3X22g(x ) =x ； g'(x ) =3(1x3)岛 2a=o； g'X)= x2/3a。Xk 1 -x 二 g(xj - g(x)二 g'(x )(Xk -x ) 2p(Xk -x )2Xk 1 -x = g''(x )2limk 诃-x)12!=3a，可见该迭代公式具有二阶收敛性。证毕【证明】(1)设：

41、f(x) = x3 a，则f'(x) =3x 10 “。 ，对任意x . 0,牛顿迭代公式、，、，f(Xk)、，xk 1 = xkxkf (Xk)5.1线性方程组迭代公式1、( p.17O，题1)用雅可比迭代与高斯-赛德尔迭代求解方程组:3x1 +x” 2，要求结N +2x2 =1X1(k d)1(k)3x2【解】 雅可比迭代公式:迭代计算结果列于下表。果有3位有效数字。kX1(k)x2k)|x1k)X1(2)|x2k)x2k)i< 0.0005?000-12/31/22/31/2N21/21/61/61/3N311/181/41/91/12N47/127/361/361/18N

42、50.601850.208330.018520.01389N60.597220.199080.00463 :0.00925N70.600310.201390.003090.00231N80.599540.199850.000770.00154N90.600050.200230.000510.00038N100.599920.199980.000030.00025Yx2k d)2十#)冷X1(k)x； : x!10)0.600;X2 ： x210) ： 0.200 ；由上表可见，所求根皆为小数点后第1位不为零的小数，要取3位有效数，则误差限为xr_ix2k1(x2k)高斯-赛德尔迭代公式：333

43、，迭代计算结果列于下表。xT 一 1x1(“)+1=1(1 + x2k)、一 2 2 6k(k)X(k)X2 (k)(J) |X1-X1|x2k)-x；j|<0.0005?000-12/31/62/31/6N20.61110.1944N30.60190.19910.00920.0047N40.60030.19990.00160.0008N50.60000.19990.00030.0000YX: x!5) : 0.600;X2 x25) 0.200 ；12、( p.171，题7)取川=1.25，用松弛法求解下列方程组，要求精度为10* 。2|4x1 - 3x2 = 163x14x2 - X

44、3 = 20-X2 4X3 二-12【解】欧先写出高斯-赛德尔迭代:(k 卅) 入1(k 1)X2(k 1)X33X(k)蔦X23 (k)_4X11 (k)x3k)5 二 943x2k).64-x2k)16-x3k)16(1)引入松弛因子，得(k 1)X1(k 1)X2(k 1)X35-4 5-4 5-4+ + +X7 X7 X7k k k对 站 必I - 4 丄 4 1-4- - - - -kk lr<k将方程组(1)代入(2)，并化简x；k d)(k 1) x2丿1 x(k)X14= x2k)6415x(k)_16X2-x3k)-162x3k 1)45(k)11 _(k)25256X

45、264计算结果见下表。k(k)(k)v(k)I (k)l (k)xz(2)|x3k)-xT |< eX1X2X3|X1-X1|X2-X2|?0000-152.5-3.12552.53.125N21.40625:2.65625-2.14844N32.158203.03223-2.28882N41.611733.15872-2.19860N51.635773.24423-2.19187N :61.549593.28508-2.17800N71.53284:3.30793-2.17320N81.515613.31978-2.17001N91.50880:3.32615-2.16847N :01

46、.504533.32951-2.16762N11.50245P 3.33130-2.16717N :21.501293.33225-2.16694N31.50069:3.33276-2.16672N41.500373.33306-2.16676N51.50016:3.33318-2.16670N61.500103.33325-2.16668N71.500053.33329-2.166680.000050.000040.00000Y迭代解：x1 -xf7)1.5001 ,X2 = x217) :- 3.3333, x3 二 x317) :- 2.1667.13x3石1667.3 10精确解：x11.5,x23.3333,235.1线性方程组迭代公式-赛德尔迭代公式，并考1、（ p.170，题2）试列出求解下列方程组的雅可比迭代公式与高斯 察迭代过程的收敛性。10洛 +x3 _5x4 = -7% 十 8x2 _ 3x3 =113x1 2x2 -8x3 X4 = 23% -2x2 2x3 7x4 =17【解】(1)雅可比迭代公式：X1(F-x3k)10710(k 1 x3)(k 1)X4_-x(k)8x1(k)二x1X2(k)11k)238(1)(k)X17(k)-X371710Gj