《数值分析简明教程》第二版(王能超编著)课后习题答案高等教育出版社

1、【思路岛答案网】 整理提供0.1算法1、(p.11，题1)用二分法求方程X -X-1 =0在1,2内的近似根，要求误差不超过10-3.【解】由二分法的误差估计式 &十黔=珀_；=10，得到2k1 _1000.两端取自然对数得空10-仁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°

2、10 0 - 2 = 1 ： 0，f (1) -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 一 如0、2 3.3219二6.6438，因此取k = 7，即至少需二分ln27次.求解过程见下表。kakbkXkf(Xk)符号0010.512345670.2误差1.（p.12，题 8）已知 e=2.7182

3、8 ；试问其近似值 X, = 2.7 , x2 = 2.71 , X2=2.71 , x 2.718各有几位有效数字？并给出它们的相对误差限。【解】有效数字：1 ,因为| e-x, | = 0.01828：0.0510 ，所以X, =2.7有两位有效数字；21 /因为|e-X2 |= 0.00828：0.0510 ，所以x 2.71亦有两位有效数字；1 3因为|e-x3 |=0.00028：0.000510 ，所以x 2.718有四位有效数字;2|e X1 |0.05<2.7=1.85%=| e - X2 |X20.052.71= 1.85% ； |e-X3 |X30.00052.718

4、= 0.0184%。【思路岛答案网】 整理提供【思路岛答案网】 整理提供评 （1）经四舍五入得到的近似数；其所有数字均为有效数字;（2）近似数的所有数字并非都是有效数字2. （ p.12，题9）设 =2.72 ； x2 =2.71828； x 0.0718均为经过四舍五入得出的近 似值，试指明它们的绝对误差（限）与相对误差（限）。10.005o【解】M -0.005 ； ；r1 二一11.84 10"x12.72；2 =0.000005； ；r2 2 : 0.000005 1.84 10（；x22.71828-0.00005 ；X30.000050.07186.96 10,评经四舍五

5、入得到的近似数，其绝对误差限为其末位数字所在位的半个单位3. （ p.12，题 10）已知 x1 =1.42 ； x -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) =si nx在节点x0 =0的5次泰勒插值多项式 p5(x)，并计算P5 (0.3367)和估计插值误差，最后将P5 (0.5)有效数值与精

6、确解进行比较。【解】由 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(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

7、) =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(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

8、)(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)依据下列数据表，试用线性插值和抛物线插值分别计算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之间，先估计误差R(x)f”()=gx)iX0)(Xr0.012<

9、;2二1 104；须保留到小数点后 4为，计算过程多余两位。2T y(X1-Xo)2/4y=(x-x 0)(x-X1)-Iij_XX1【思路岛答案网】 整理提供【思路岛答案网】 整理提供R(X)R(X)sin(x0) sin(xj =1 (x-x0)sin(x1) (% -x)sin(x0)lX。- X1X1 - X0X1 - X。1(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.3304(2)抛物线插值插值误差：f'''(©C

10、OS化)R>(x)(X -Xo)(X -Xi)(X -X2)(X - Xo)(Xi - x)(x - X2)3!63 0.013610max(x -XoXj -x)(x2 -x)6yy=(x-x°)(x-x1)(x-x2)Max=3(x 1-Xo)3/8XoX1X2【思路岛答案网】 整理提供【思路岛答案网】 整理提供抛物线插值公式为:P2(x)j) sin(Xo)(X。-XJ(Xo -X2)(x XoW-x?) sin(xj(% -X°)(X1 -X2)(X-Xj(X-Xo)(X2 -X1)(X2 - Xo)sin(x2)1(花 - X)(X2 - x)0.022 I

11、L2sin(x0) (xx0)(x2x)sin(x1)-(% -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.02*1 n 3.8445 sin(0.32) 38.911 sin(0.34) - 2.7555 sin(0.36)-0.33037439 0.02经四舍五入后得：P2 (0.3367) =0.330374，与 sin(0.3367) = 0.330374191精确值相比较，在插值误差范围内完全吻合！1.

12、3分段插值与样条函数1、( p.56，习题33)设分段多项式X3 X22x3 bx2 ex -1【思路岛答案网】 整理提供【思路岛答案网】 整理提供是以0,1,2为节点的三次样条函数，试确定系数b，c的值.【解】依题意，要求S(x)在x=1节点【思路岛答案网】 整理提供一阶导数连续：S' (1) =3 12 2 1 = 6 122 b 1 c 二 S (1),函数值连续:S_=13 12 =2 13b 12c 1 S (1),即：b c = 1(1)即: 2b c = 1(2)解方程组(1)和(2),得b = 2,c = 3，即f 32S(x) = «x +x2x3 2x2

13、+3x-1由于 S(1) = 3汇2x1+2 = 6x2x12x2 = S；(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两段，对应的插值函数为S(x)和S2(x)，禾U用拉格朗日线性插值公式，求得x-xXX0X1X0cl CL /Si(x) y0- y110.5=_0.5x 1;x° X1X1 x°01

14、10_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，则余项表达式|f

15、(J|R(x"2!M 224|(x_1)(x_2)|20.5 = 0.5470625 乞 0.5【思路岛答案网】 整理提供【思路岛答案网】 整理提供1.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-

16、3x 46y 二 48(2) , 得:3.04029,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

17、 (=xi2yi)191936113032168592532.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 369

18、321.5 -5327 271.4400859.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 ： A0f(；)几口；)A2fH)；042411(3) f(x)dx ： f (0) A0f(x0) o04【解】(1 )令f(x) =1,X,X2时等式精确成立，可列出如下方程组：A0

19、+A +A2 =2h* A0 + A2 = 0Ao + A2 = 2 h 3h 4解得：A0 = A2, A1h，即:33(1)(2)(3)hhf (x)dx ： f (-h) 4f (0)f(h)，可以验证，对f(x) =x3公式亦成立，而对f (x) =x4不成立，故公式1)具有3次代数精度。(2)令f(X)=1,X,X2时等式精确成立，可列出如下方程组:代十几+A2 =1(1)“代 +2入 +3A2 =2(2)3A0 +12A +27A2 =16(3)解得：2111113a-a?a-3，即："曲护叫比)2七)，可以验证，对f(x) =X3公式亦成立，而对 f(x) =X4不成立

20、，故公式(2)具有3次代数精度。A0(3)令f (x) =1,x时等式精确成立，可解得:1132即： 0 f (x)dx応：f (0)f (§)，可以验证，对f(x) =x3不成立，故公式(3)具有2次代数精度。12、( p.95 ,习题6)给定求积节点x0, x14求积公式，并指明该求积公式的代数精度。【解】依题意，先求插值求积系数：dx =31 x -1 40 1 342X0 -3f(x) =x2公式亦成立，而对试构造计算积分41I = o f (x)dx的插值型dx = -2 (x223、3X)3【思路岛答案网】 整理提供1 x _XoA10X1X0dx1x -41dx = 2

21、 (1 x22插值求积公式:10 f (x)dx 八k =01Akf(xkH2f(4)3f（4）当f（X）=1，左边=0 f (x)dx=1 ；右边=gx1+x1=1 ； 左 =右;当1f (x) = x，左边=o f (x) dx =i ；左=右;当2 1f (x) = x ,左边=° f (x)dx =右边=丄丄 1鸟2 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 5

22、21.066 660.721 591分别用复化梯形法与复化辛普生法求积分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)f(Xk1)V f

23、(a) jf(xk2n -4)2、 f(xQ f(b)k T0 5S2 f (0.00) 4 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)设用复化梯形法计算积分|二jldx，为使截断误差不超过i 10，问应当划分区间【0,1】为多少等分？如果改用复化辛普生法呢?【解】(1)用复化梯形法，a = 0,b = 1, f (x) = f'(x) = f ''(x) = ex，设需划分n等分,则其截断误差表达式为:|Rt|=|I -Tn|

24、 =(b -a)312n2(1-0)3 ；12n3 e ；1依题意，要求I Rt N 102n2 -6e 105a：212.849，可取 n =213。(2)用复化辛普生法a = 0,b = 1, f (x) = f' (x) = f''''(x) = ex，截断误差表达式为:|Rs|I -Sn |=(b-a)5180(2n)45maxf''''(厂翳”猛?；1 5依题意，要求| RS |10 ，即25 < 1 10* 二 n4_e103.70666，可取 n=4，划分 8 等分。2880n4214402.3数值微分

25、1、( p.96，习题24)导出三点公式(51)、(52)和(53)的余项表达式1f'(Xo)-3f(Xo) 4f(X1)-f(X2)(51)2h1f'(xj-f(x°) fg)(52)2h1f'X)f(x°)-4f(X1)3f(X2)(53)2h【解】如果只求节点上的导数值，利用插值型求导公式得到的余项表达式为f(n*)nR(xQ 二 f'(xQ - p'(xj7 H X W)(n +1)!jTj我由三点公式(51)、(52)和(53)可知，n =2,h-x0 =x2 -花，贝yR(X0)=有占 i2(X0-Xj"t3a(

26、X0-X1)(X0-xA%h2(2 1)! j3!3X1.01.11.2f(x)0.25000.22680.2066的数据表,设已给出f(X)'X)2R(xJ 二(21)( 1)丨丨(X1 Xj )二(2 1)! jj jj/2f'''(1)3!(x x0)(x1 x2) = 3!A "(X2_Xj)")(x2 - x0 )(x2 - Xj =3!f'''( 2)h22、( p.96，习题试用三点公式计算 f'(1.0), f'(1.1), f'(1.2)的值，并估计误差。f'(1.0)

27、f'(1.1)f'(1.2)f(X)5 x)21 1-3f (1.0) 4f (1.1) - f (1.2)-3 0.2500 4 0.2268 - 0.2066 = -0.24702h2 0.11 1_ f (1.0) f (1.2) =-0.2500 0.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.11 2；=f'(x) 卞；二 f''(x)(1+x)(1 x)4；=f'''(x)=-24(1 x)5

28、 '【解】已知 x0 = 1.0, Xr = 1.1, x2 = 1.2, h = X<| - x()= X2 - X<| = 0.1，用三点公式计算微商：【思路岛答案网】 整理提供用余项表达式计算误差R(1.0) = f ( 0)h2-24 0.123f'''( 1)屮3(1 1.0)5-0.0025R(1.1)二3!fTh224 0.123!(1 1.0)5-24 0.12:0.00125亍-0.049673(1 1.1)【思路岛答案网】 整理提供【思路岛答案网】 整理提供3、( p.96，习题 26)设f(x)二sinx，分别取步长h =0.1

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

30、则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,贝y1 1 f'(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)很接

31、近，直接相减会造成有效数字的严重损失。因此，从舍入误差的角度看，步长不宜太小。3.1 Euler 格式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 二 yn hy'n 二 £h(X； - 二 Yn 02 M - 丫7)；2 2(2)yn1 = ynh(磚仏)=yn0.2(聖上)。XnXnXnXn2、( p.124，题2)取h = 0.2，用欧拉方法求解初

32、值问题y' - - y - xy2 (0 _ x _ 0.6)，y(0) =1。【解】欧拉格式：yn 1 二 ynhy'n二ynh(-yn -Xny；) = yn 0.2(-yn-Xnyi)；化简后，yn1 =°.8yn -0.2xny；，计算结果见下表。n0123Xn0.00.20.40.6yn1.00.80.61440.46131 23、( p.124，题3)取h = 0.1，用欧拉方法求解初值问题y'2-2y2(0_x_4)，1 + Xy(0) =0。并与精确解y 程 比较计算结果。1 +x【解】欧拉格式：1 2 1 2hyr h, x2yn)-yn 0

33、.2 - x2yn化简后，yn j = yn -0.4y； -0电，计算结果见下表。1 +Xn1、( p.124，题7)用改进的欧拉方法求解上述题2，并比较计算结果。【思路岛答案网】 整理提供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.41

34、3计算结果见下表。yn 1【解】公式：因为 y'= f (x,y) - -y - xy2(0 乞 x < 0.6)，h = 0.2，且 y(0) =1，则改进的欧拉【思路岛答案网】 整理提供与原结果比较见下表n0123Xn0.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位数字。【解】 四阶经典的龙格-库塔方法公式：h%卅=

35、yn +：(K1 +2K2 +2K3 +K4)K1 二二 f(Xny)«k2=f (x1 ,ynK1)n+22K3 =f (xn甲，Yn22K2)K4二 f(Xn十，yn+ h©)列表求得y(0.4)如下:nXnyn00.02.00010.22.300420.42.46544.1迭代法及收敛定理(k =0,1,2,)，求201、( p.153，题1)试取Xo =1，用迭代公式Xk x： +2xk +10方程x3 2x210x-20=0的根，要求准确到10。【解】 迭代计算结果列于下表kXk|Xk-Xk-11<0.001kXk|Xk-Xk-11<0.00111.5

36、38460.53846N61.365930.00937N21 1.295020.24344N :71.370090.00416N31.401820.10680N81.368240.00185N41.354210.04761N91.369060.00082Y51.375300.02109N因为 | 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) c

37、osx,，且一阶导数2 22 21 1 1 1 g'(x) = - QSin x 连续，| g'(x)鬥 - in x |空？ 1，所以迭代过程Xk 1 二?cosxk 对1x0 R均收敛。(压缩映像定理)，方程xcosx有且仅有一实根。 证毕23、 ( p.153，题4)证明迭代过程xk 丁 =乞丄对任意初值x0 1均收敛于、2。2Xk【证明】设：g(x)=仝+1，对于任意x A1，因为2xi- 12 ，所以 g(x) .2。2 x2 x V 2 x111x1一阶导数g'(x)21，根据压缩映像定理，迭代公式Xk 1 k 对任意2x22xkx 1初值x0 1均收敛。假

38、设lim xk = x ，对迭代式xk 1 k 两边取极限，则有十2Xk宀吉，则22 = 2，解得x7，因宀出不在心范围内，须舍去。4.2牛顿迭代法1、(p.154，题17)试用牛顿迭代法求下列方程的根，要求计算结果有4位有效数字:(1)x3_3x_1=0， x0 = 2(2)x2 -3x _ex 2=0， x0 = 1【解】(1)Xk 1二 Xk设 f (x)二 x3 -3x -1，则3f (Xk)Xk -3XkXkf'(Xk)2f'(x) = 3x -3，牛顿迭代公式：=2xk1(k二0,1,2,)，迭代计算过3x： 一 33(x：-1)'【思路岛答案网】 整理提供

39、因为 | x3 -x2 |： 0.00006 ： 10鼻，所以 xx3 =1.879。(2)设 f (x) = x2 -3x-ex 2，则 f'(x) = 2x-3-ex，牛顿迭代公式：=xj/(Xk-0-2化=0,1,2,)f(Xk)、， Xk 3Xk e +2xk 1 = xkxk -'f (xk)2Xk -3-eXk2Xk _3_eXk程见下列表。kXk:k-Xk-11<0.0001kXk|Xk-Xk-1|<0.000111.888890.11111N31.879390.00006Y21.879450.00944N :【思路岛答案网】 整理提供，迭代计算过程见

40、下列表。kXk|Xk-Xk-11<0.0001kXk:k-Xk-11<0.00110.268940.73106N30.257530.00014N20.257390.01155N40.257530.00000Y所以 x” : x4 =0.2575。因为 |x3 x2 I： 0.00000 ： 10*2、( p.154，题18)应用牛顿法于方程x3 - a = 0，导出求立方根3 a(a 0)的迭代公式,并证明该迭代公式具有二阶收敛性。【证明】(1)设：f (x) = x3 - a，2f'(x) = 3x，对任意X 0，牛顿迭代公式33xk -a2xk akk= 0,1,2,3

41、x：3xi(2)由以上迭代公式，有:kXmKHk/Vg3-a+ 23X2-g"3(T)xMXk 1 -X ' =g(Xk) -g(x ) =g'(x )(Xk-x )豊Xk -x )22!【思路岛答案网】 整理提供【思路岛答案网】 整理提供km(:二；2= g”2；J：，可见该迭代公式具有二阶收敛性。证毕【思路岛答案网】 整理提供(1)【思路岛答案网】 整理提供(1)5.1线性方程组迭代公式1、（ p.17O，题1）用雅可比迭代与高斯-赛德尔迭代求解方程组:3X1 X 2，要求结X1 2x2 = 1果有3位有效数字。【解】 雅可比迭代公式:x；k d）1(k)21(k

42、)X2(2X2 )3331 (k)11(k)X1(1-花)2 22迭代计算结果列于下表。k(k)X1()x2k)i(k)(k)-| Xi- 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.2002

43、30.000510.00038N100.599920.199980.000030.00025Y3位有效数，则误差限为(k 1)X20.200 ；1位不为零的小数，要取,x110)0.600;x2、x210)由上表可见，所求根皆为小数点后第1 10 "。2X1(k 1)高斯-赛德尔迭代公式:1X23(k)沁（心），迭代计算结果列于下表。丄X1(k1)J(1 x2k)2 2 6kX1(k)x2k)|X1(k)-xI< 0.0005 ?000-12/31/62/31/6N20.61110.1944N30.60190.19910.00920.0047N40.60030.19990.00

44、160.0008N50.60000.19990.00030.0000Yk 1)x；x；5)0.600;x2x25)0.200 ；2、（ p.171，题17）取用-1.25，用松弛法求解下列方程组，要求精度为10 *。24x1 3x2 = 163x1 4x2 -x3 二 20-x2 ' 4x3 = -12【思路岛答案网】 整理提供(1)【解】欧先写出高斯-赛德尔迭代:(k出)X13 (k)-X24(k 1)X23 (k)_4X1x2k)16 (k 1)X31 (k)2-3 9x2k) 丄x3k)6416引入松弛因子,(k 1)X1=(1(k)-')X1(k1)1(k)X1（k 1

45、）x2k d)=(1k)2(k1)(k 1)X2(2)x3k d)=(1）x3k）(k1)3k)将方程组（1）代入（2），并化简(k 1)X31(k)4X115(k)(k 1)X29 (k)_ 64 X2J16(k 1)X345(k)11 (k)64 X325256【思路岛答案网】 整理提供(1)计算结果见下表。k(k) x；)x2k)x3k)1(k)(k4) |X1()-x；) |x2k)-x2k)i|x3k)-x3k)i<e?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-2.17800N71.532843.30793-2.17320N81.515613.31978-2.17001N91.508803.32615-2