《数值计算》课后习题汇总
数值计算第一二章答案
第一章数值计算中的误差习题一1.1 下列各近似数的绝对误差限是最末位的半个单位，试指出它们各有几位有效数字。
1x =-3。
105 , 2x =0.001, 3x =0。
100, 4x =253。
40， 5x =5000， 6x =5⨯310.答案：4，1，3，6，4，1。
1。
2 设100〉*x >10,x 是*x 的有五位有效数字的的近似数,求x 的绝对误差限。
答案:当10<x 〈100时，因为有5位有效数字,所以绝对误差限为0。
005。
1。
3 求下列各近似数的相对误差限和有效数字位数： 1） 123x x x ++,2） 124x x x 3） 24x x 答案：()10.0005e x ≤()20.0005e x ≤()30.0005e x ≤ ()40.005e x ≤ ()50.5e x ≤ ()60.5e x ≤1）()()()()123123e x x x e x e x e x ++=++≤()()()123e x e x e x ++3221.5100.15100.510---≤⨯=⨯≤⨯2123()0.1510x x x ε-++=⨯123123123()()0.0004993...0.0004994r x x x e x x x x x x ε++++==≤++123x x x ++=-3。
004 精确到小数点后两位,所以有三位有效数字。
2）()()()()()()12424112424114224()e x x x x x e x x e x x x x e x x x e x x e x =+=++ =()()()241142124)x x e x x x e x x x e x ++()()()241142124x x e x x x e x x x e x ≤++ =660.5100.31050.0005 3.1050.510--⨯+⨯+⨯⨯ 所以43124() 1.71275100.510x x x ε--=⨯≤⨯124x x x =43.105100.0003105--⨯=-41241244124() 1.7127510()0.5515...3.10510r x x x e x x x x x x ε--⨯===⨯3）()()2222424244444()()1x x e x x e x e e x e x x x x x x ⎛⎫≈-≤+⎪⎝⎭325105420.5100.5100.197316100.77868100.1997100.510253.40253.40------⨯⨯=+=⨯+⨯≈⨯<⨯ 又由24x x 50.3946310-≈⨯知有0位有效数字 ∴522440.1997100.5r x e x x x -⎛⎫⨯≤≈ ⎪⎝⎭1。
《数值计算方法》习题答案
《数值计算方法》课后题答案详解吉 林 大 学第一章 习 题 答 案1. 已知(1)2,(1)1,(2)1f f f −===，求()f x 的Lagrange 插值多项式。
解：由题意知：()01201212001020211012012202121,1,2;2,1,1()()(1)(2)()()6()()(1)(2)()()2()()(1)(1)()()3(1)(2)(1)(2)()2162nj j j x x x y y y x x x x x x l x x x x x x x x x x l x x x x x x x x x x l x x x x x x x x L x y l x ==−=====−−−−==−−−−+−==−−−−−+−==−−−−+−==×+×−∴∑()2(1)(1)131386x x x x +−+×=−+2. 取节点01210,1,,2x x x ===对x y e −=建立Lagrange 型二次插值函数，并估计差。
解11201201210,1,;1,,2x x x y y e y e −−======1)由题意知：则根据二次Lagrange插值公式得：02011201201021012202110.510.520.51()()()()()()()()()()()()()2(1)(0.5)2(0.5)4(1)(224)(43)1x x x x x x x x x x x x L x y y y x x x x x x x x x x x x x x x x e x x e e e x e e x −−−−−−−−−−−−=++−−−−−−=−−+−−−=+−+−−+22)Lagrange 根据余项定理，其误差为(3)2210122()1|()||()||(1)(0.5)|3!61max |(1)(0.5)|,(0,1)6()(1)(0.5),()330.5030.2113()61()0.2113(0.21131)(0.21130.5)0.008026x f R x x e x x x x x x t x x x x t x x x x t x R x ξξωξ−+≤≤==−−≤−−∈′=−−=−+=−==≤××−×−=∴取 并令 可知当时，有极大值3. 已知函数y =在4, 6.25,9x x x ===处的函数值，试通过一个二次插值函数求的近似值，并估计其误差。
《数值计算方法》 课后题 答案(曾金平)湖南大学
= 0.105 × 10 2 − 0.144 × 10 2 + 0.657 × 10 1 − 1
= 0.167 × 101
g ( 2.19) = ((−0.81) × 0.219 × 101 + 3) × 0.219 × 101 − 1
= 0.123 × 10 1 × 0.219 × 10 1 − 1 = 0.169 × 101 1 1 即 f ( x) = 0.167 × 10 ， g ( x) = 0.169 × 10
而当 x = 2.19 时 x 3 − 3 x 2 + 3 x − 1 的精确值为 1.6852，故 g ( x ) 的算法较正确。 8．按照公式计算下面的和值（取十进制三位浮点数计算） ： （1）
∑3
i =1
6
1
i 6
；(2)
∑3
i =6
1
1
i
。
解： （1）
∑3
i =1 1
1
i
1 1 1 1 1 1 = + 2 + 3 + 4 + 5 + 6 = 0.333 + 0.111 + 0.037 + 0.012 + 0.004 + 0.001 3 3 3 3 3 3
xk
1 1.5 1.25 1.125 1.0625 1.09375 1.109375 1.1171875 1.11328125 1.115234375 1.1142578125 1.11376953125
f ( xk )
-0.1585 0.4962 0.1862 0.015051 -0.0718 -0.02835 -0.00664 0.004208 -0.001216 0.001496 0.001398 -0.000538
《数值计算方法》课后题答案(湖南大学-曾金平)
习题一1．设x >0相对误差为2%，4x 的相对误差。
解：由自变量的误差对函数值引起误差的公式：(())(())'()()()()f x xf x f x x f x f x δδ∆=≈得（1）()f x =11()()*2%1%22x x δδδ≈===；（2）4()f x x =时444()()'()4()4*2%8%x x x x x xδδδ≈===2．设下面各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出他们各有几位有效数字。
（1）12.1x =；（2）12.10x =；（3）12.100x =。
解：由教材9P 关于1212.m nx a a a bb b =±型数的有效数字的结论，易得上面三个数的有效数字位数分别为：3，4，53．用十进制四位浮点数计算 （1）31.97+2.456+0.1352； （2）31.97+（2.456+0.1352）哪个较精确？解：（1）31.97+2.456+0.1352 ≈21((0.3197100.245610)0.1352)fl fl ⨯+⨯+ =2(0.3443100.1352)fl ⨯+=0.3457210⨯（2）31.97+（2.456+0.1352）21(0.319710(0.245610))fl fl ≈⨯+⨯ = 21(0.3197100.259110)fl ⨯+⨯ =0.3456210⨯易见31.97+2.456+0.1352=0.210⨯，故（2）的计算结果较精确。
4．计算正方形面积时，若要求面积的允许相对误差为1%，测量边长所允许的相对误差限为多少？解：设该正方形的边长为x ，面积为2()f x x =，由(())(())'()()()()f x xf x f x x f x f x δδ∆=≈解得(())()()'()f x f x x xf x δδ≈=2(())(())22f x x f x x xδδ==0.5%5．下面计算y 的公式哪个算得准确些？为什么？（1）已知1x <<，（A ）11121xy x x-=-++，（B ）22(12)(1)x y x x =++； （2）已知1x >>，（A ）y=，（B ）y =； （3）已知1x <<，（A ）22sin x y x =，（B ）1cos 2xy x-=；（4）（A）9y =（B ）y =解：当两个同（异）号相近数相减（加）时，相对误差可能很大，会严重丧失有效数字；当两个数相乘（除）时，大因子（小除数）可能使积（商）的绝对值误差增大许多。
数值计算方法与算法第三版课后习题答案
数值计算方法与算法第三版课后习题答案1. 矩阵乘法问题描述给定两个矩阵A和B，尺寸分别为n×m和m×p，求矩阵A 和矩阵B的乘积矩阵C，尺寸为n×p。
算法实现import numpy as npdef matrix_multiplication(A, B):n, m = A.shapem, p = B.shapeC = np.zeros((n, p))for i in range(n):for j in range(p):for k in range(m):C[i][j] += A[i][k] * B[k][j] return C示例A = np.array([[1, 2], [3, 4]])B = np.array([[5, 6], [7, 8]])C = matrix_multiplication(A, B)print(C)输出结果：[[19. 22.][43. 50.]]2. 数值积分问题描述给定一个函数f(x)，以及积分区间[a, b]，求函数f(x)在区间[a, b]上的定积分值∫abf(x)dx。
算法实现简单的数值积分算法是采用小梯形法，将区间[a, b]均分成n个子区间，然后计算每个子区间的面积，最后将这些子区间面积相加得到定积分值。
def numerical_integration(f, a, b, n):h = (b - a) / nintegral =0for i in range(n):x1 = a + i * hx2 = a + (i +1) * hintegral += (f(x1) + f(x2)) * h /2 return integral示例import mathf =lambda x: math.sin(x)a =0b = math.pin =100result = numerical_integration(f, a, b, n) print(result)输出结果：1.99983550388744363. 非线性方程求解问题描述给定一个非线性方程f(x) = 0，求方程的根x。
数值计算方法课后习题答案
习题一1．设x >0相对误差为2%4x 的相对误差。
解：由自变量的误差对函数值引起误差的公式：(())(())'()()()()f x xf x f x x f x f x δδ∆=≈得 （1）()f x =11()()*2%1%22x x δδδ≈===；（2）4()f x x =时444()()'()4()4*2%8%x x x x x x δδδ≈===4．计算正方形面积时，若要求面积的允许相对误差为1%，测量边长所允许的相对误差限为多少？解：设该正方形的边长为x ，面积为2()f x x =，由(())(())'()()()()f x xf x f x x f x f x δδ∆=≈ 解得(())()()'()f x f x x xf x δδ≈=2(())(())22f x x f x x xδδ= =0.5%5．下面计算y 的公式哪个算得准确些？为什么？（1）已知1x <<，（A ）11121x y x x -=-++，（B ）22(12)(1)x y x x =++； （2）已知1x >>，（A）y =，（B）y =（3）已知1x <<，（A ）22sin x y x=，（B ）1cos 2xy x -=；（4）（A）9y =（B）y =解：当两个同（异）号相近数相减（加）时，相对误差可能很大，会严重丧失有效数字；当两个数相乘（除）时，大因子（小除数）可能使积（商）的绝对值误差增大许多。
故在设计算法时应尽量避免上述情况发生。
（1）（A ）中两个相近数相减，而（B ）中避免了这种情况。
故（B ）算得准确些。
（2）（B ）中两个相近数相减，而（A ）中避免了这种情况。
故（A ）算得准确些。
（3）（A ）中2sin x 使得误差增大，而（B ）中避免了这种情况发生。
故（B ）算得准确些。
（4）（A ）中两个相近数相减，而（B ）中避免了这种情况。
数值课后题答案
习 题 一1、解：根据绝对误差限不超过末位数的半个单位，相对误差限为绝对误差限除以有效数字本身，有效数字的位数根据有效数字的定义来求.因此 49×10-2 ：E = 0.005； r E = 0.0102； 2位有效数字. 0.0490 ：E = 0.00005；r E = 0.00102； 3位有效数字. 490.00 ：E = 0.005； r E = 0.0000102；5位有效数字.2、解：722= 3.1428 …… ， π = 3.1415 …… ， 取它们的相同部分3.14，故有3位有效数字.E = 3.1428 - 3.1415 = 0.0013 ；r E = 14.3E = 14.30013.0 = 0.00041.3、解：101的近似值的首位非0数字1α = 1,因此有 |)(*x E r |)1(10121--⨯⨯=n < = 21× 10-4 , 解之得n > = 5，所以 n = 5 . 4、证：)()(1)()(1)(*11**11**x x x nx E x n x E n n n-=≈--)(11)()(1)()(*****11****x E nx x x n x x x x nx x E x E r nnnn n r =-=-≈=- 5、解：(1)因为=20 4.4721…… ，又=)(*x E |*x x -| = |47.420-| = 0.0021 < 0.01, 所以 =*x 4.47.(2)20的近似值的首位非0数字1α = 4,因此有 |)(*x E r |)1(10421--⨯⨯=n < = 0.01 , 解之得n > = 3 .所以，=*x 4.47. 6、解：设正方形的边长为x ,则其面积为2x y =，由题设知x 的近似值为*x = 10 cm .记*y 为y 的近似值，则)(20)(20)(2)(*****x E x x x x x y E =-=-= < = 0.1，所以)(*x E < = 0.005 cm .7、解：因为)()(*1x x nx x E n n -≈-,所以n x nE x x x n xx E x E r nn nr 01.0)()()(*==-≈=. 9、证：)()()(**t gtE t t gt S S S E =-≈-=t t E gt t t gt S S S S E r )(22/)()(2**=-≈-= 由上述两式易知，结论. 10、解：代入求解,经过计算可知第(3)个计算结果最好.11、解：基本原则为：因式分解，分母分子有理化、三角函数恒等变形…… （1）通分；（2）分子有理化；（3）三角函数恒等变形.12、解： 因为20=x ,41.1*0=x ,所以|*0x x -| < = δ=⨯-21021 于是有|*11x x -| = |110110*00+--x x | = 10|*00x x -| < =δ10 |*22x x -| = |110110*11+--x x | = 10|*11x x -| < =δ210 类推有 |*1010x x -| < =810102110⨯=δ 即计算到10x ，其误差限为δ1010，亦即若在0x 处有误差限为δ，则10x 的误差将扩大1010倍，可见这个计算过程是不稳定的.习 题 二1、 解：只用一种方法.（1）方程组的增广矩阵为：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----11114423243112 → ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡----1010411101110112 → ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡---11041001110112 → 31=x , 12=x , 13=x . （2）方程组的增广矩阵为：⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡------017232221413 → ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--247210250413 → ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡--147200250413 → 21=x , 12=x , 2/13=x . （3）适用于计算机编程计算.2、 解：第一步：计算U 的第一行，L 的第一列，得611=u 212=u 113=u 114-=u3/1/112121==u a l 6/1/113131==u a l6/1/114141-==u a l第二步：计算U 的第二行，L 的第二列，得3/1012212222=-=u l a u 3/213212323=-=u l a u 3/114212424=-=u l a u 5/1/)(2212313232=-=u u l a l10/1/)(2212414242=-=u u l a l第三步：计算U 的第三行，L 的第三列，得10/37233213313333=--=u l u l a u 10/9243214313434-=--=u l u l a u 37/9/)(33234213414343-=--=u u l u l a l第四步：计算U 的第四行，得370/9553443244214414444-=---=u l u l u l a u从而， ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡----3101141101421126 = ⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡--137/910/16/1015/16/10013/10001⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡---370/95500010/910/37003/13/23/1001126 由b LY = , 解得Y =（6，-3，23/5，－955/370）T . 由Y UX = , 解得X =（1，-1，1，－1）T . 3、（1）解：首先检验系数矩阵的对称正定性，这可以通过计算其各阶顺序主子式是否大于零来判断.11a ＝ 3 > 0, 2223= 2 > 0, 301022123 = 4 > 0,所以系数矩阵是对称正定的.记系数矩阵为A ，则平方根法可按如下三步进行：第一步 分解：A = L L T . 由公式计算出矩阵的各元素：311=l 33221=l 3622=l 3331=l 3632-=l 233=l因此， L ＝⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡-23633036332003. 第二步 求解方程组LY = b . 解得Y = （335，36，2）T . 第三步 求解方程组L T X = Y . 解得X =（0，2，1）T .（2）解：首先检验系数矩阵的对称正定性，这可以通过计算其各阶顺序主子式是否大于零来判断.11a ＝ 3 > 0,2223= 2 > 0, 1203022323 = 6 > 0,所以系数矩阵是对称正定的.记系数矩阵为A ，则平方根法可按如下三步进行：第一步 分解：A = L L T . 由公式计算出矩阵的各元素：311=l 33221=l 3622=l 331=l 632-=l 333=l因此， L ＝⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡-363036332003 . 第二步 求解方程组LY = b . 解得Y = （335，66-，33）T. 第三步 求解方程组L T X = Y . 解得X = （1，21，31）T . 4、解： 对1=i , 2111==a d ；对2=i , 121-=t , 2121-=l , 252-=d ；对3=i , 131=t , 2732=t ,2131=l , 5732-=l ,5273=d .所以数组A 的形式为： ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡---=527572102521002A 求解方程组LY = b . 解得Y = （4，7，569）T . 求解方程组DL T X = Y . 解得X = （910，97，923）T .5、解：（1）设A = LU = ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡1010000000000010010015432l l l l ⎥⎥⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎢⎢⎣⎡5432106000000000600006006u u u u u计算各元素得： 51=u , 512=l , 1952=u , 1953=l , 19653=u ,65194=l , 652114=u , 211655=l , 2116655=u .求解方程组LY = d . 解得Y = （1，51-，191，651-，211212）T.求解方程组UX = Y . 解得X = （6651509，6651145，665703，665395-，665212）T.（2）设A = LU = ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡100100132l l ⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡3211001u u u 计算各元素得：51=u ，512=l ，5242=u ，2453=l ，241153=u . 求解方程组LY = d . 解得Y = （17，553，24115）T.求解方程组UX = Y . 解得X = （3，2，1）T. 6、证：（1）（2）相同.因为此方程组的系数矩阵为严格对角占优矩阵，所以雅可比迭代法和相应的高斯－赛德尔迭代法都收敛. （1）雅可比迭代公式：7107271)(3)(2)1(1+--=+k k k x x x14141)(3)(1)1(2+--=+k k k x x x329292)(2)(1)1(3+--=+k k k x x x高斯－赛德尔迭代公式：7107271)(3)(2)1(1+--=+k k k x x x14141)(3)1(1)1(2+--=++k k k x x x329292)1(2)1(1)1(3+--=+++k k k x x x（2）雅可比迭代公式：545152)(3)(2)1(1+-=+k k k x x x525351)(3)(1)1(2++-=+k k k x x x5115152)(2)(1)1(3++=+k k k x x x高斯－赛德尔迭代公式：545152)(3)(2)1(1+-=+k k k x x x525351)(3)1(1)1(2++-=++k k k x x x5115152)1(2)1(1)1(3++=+++k k k x x x7、(1)证：因为此方程组的系数矩阵为严格对角占优矩阵，所以雅可比迭代法和相应的高斯－赛德尔迭代法都收敛。
数值计算课后答案5
习 题 五 解 答1、用矩形公式、梯形公式、抛物线公式计算下列积分，并比较结果。
(1)120(8)4xdx n x =+⎰，(2)20sin (8)x xdx n π=⎰(3)1(4)n =⎰，(4)1(4)x e dxn -=⎰1*、用矩形公式、梯形公式、抛物线公式计算下列积分，并比较结果。
(1)120(4)4x dx n x =+⎰解：解：将区间［0，1］4等分，5个分点上的被积函数值列表如下（取2位小数）(1)矩形法。
用矩形法公式计算（取2位小数）或者 (2)梯形法用梯形法公式计算（取2位小数）： (3)抛物线法用抛物线法公式计算（取2位小数）：2、用复化梯形公式计算积分841dx x ⎰，由此计算ln2（注：841ln 2dx x=⎰），精度要求为410-。
解：8418ln8ln 4ln ln 24dx x =-==⎰，要求精度为410-，即误差不超过41102ε-=⨯。
将积分区间［4，8］n 等份，则步长844h n n -==在本题中，复化梯形公式的余项为2228484416()()()()12123r h f f f n nηηη--''''''=-=-=-注意到231(),(),()2f x f x x f x x x--'''==-=，所以在［4，8］区间上3()24f x -''≤⨯，则32232161621283346r n n n-⨯≤⨯⨯==⨯， 要使42111062n -≤⨯，需有42421110310577.36757862n n n n n -≤⨯⇒≥⇒≥⇒≥⇒=。
3、用复合梯形公式计算积分()baf x dx ⎰，问将积分区间[a,b]分成多少等份，才能保证误差不超过ε（不计舍入误差）？解：对于复合梯形公式来说，如果()f x ''在积分区间上连续，则其余项为2(),[,]12b a r h f a b ηη-''=-∈，设max ()a x bM f x ≤≤''=，则322()()()1212b a b a Mr h f nη--''=≤ 令32()12b a Mn ε-≤，得n ≥即当1n =+时，能保证计算的精度要求。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
习 题 一 解 答1．取3.14，3.15，227，355113作为π的近似值，求各自的绝对误差，相对误差和有效数字的位数。
2、用四舍五入原则写出下列各数的具有五位有效数字的近似数。
346．7854，7．000009，0．0001324580，0．6003003、下列各数都是对准确数进行四舍五入后得到的近似数，试分别指出他们的绝对误差限和相对误差限和有效数字的位数。
4.0.1％。
5、在计算机数系F(10,4,-77,77)中，对31120.14281100.31415910x x =⨯=-⨯与，试求它们的机器浮点数()(1,2)i fl x i =及其相对误差。
6、在机器数系F(10,8,L,U)中，取三个数4220.2337125810,0.3367842910,0.3367781110x y z -=⨯=⨯=-⨯，试按(),()x y z x y z ++++两种算法计算x y z ++的值，并将结果与精确结果比较。
7、某计算机的机器数系为F(10,2,L,U)，用浮点运算分别从左到右计算及从右到左计算10.40.30.20.040.030.020.01+++++++试比较所得结果。
8、对于有效数1233.105,0.001,0.100x x x =-==，估计下列算式的相对误差限21123212333,,x y x x x y x x x y x =++==9、试改变下列表达式，使其计算结果比较精确（其中1x 表示x 充分接近0，1x表示x 充分大）。
(1)1212ln ln ,x x x x -≈； (2)11,111x x x x---+；(3)1x ；(4)1cos ,01xx x x -≠且; (5)1cot ,01x x xx-≠且。
10、用4位三角函数表，怎样算才能保证1cos 2-有较高的精度？11、利用27.982≈求方程25610x x -+=的两个根，使它们至少具有4位有效数字。
12、试给出一种计算积分11(0,1,2,3,...)n x n I ex e dx n -==⎰， 近似值的稳定算法。
习 题 二 解 答1．用二分法求方程x 3-2x 2-4x-7=0在区间[3，4]内的根，精确到10-3，即误差不超过31102-⨯。
2．证明1sin 0x x --=在[0,1]内有一个根，使用二分法求误差不大于41102-⨯的根，需要迭代多少次？ 3．试用迭代公式10220,1210k k k x x x x +==++，求方程32210200x x x ++-=的根，要求精确到510-。
4．将一元非线性方程20cos x x e -=写成收敛的迭代公式，并求其在005.x =附近的根，要求精确到210-。
5．为求方程3210x x --=在015.x =附近的一个根，设将方程改为下列等价形式，并建立相应的迭代公式：12213223121121111121111311(),;(),();(),.()k kk kk k x x x x x x x x x x x x +++=+=+=+=+==--迭代公式迭代公式迭代公式试分析每种迭代公式的收敛性，并取一种公式求出具有4位有效数字的近似值。
6．设23()()x x C x ϕ=+-，应如何取C 才能使迭代公式1()k k x x ϕ+=具有局部收敛性？7．用牛顿法求方程3310x x --=在初始值02x =邻近的一个正根，要求3110k k x x -+-<。
8．用牛顿法解方程10c x-=，导出计算数c 的倒数而不用除法的一种简单的迭代公式。
用此公式求0.324的倒数，设初始值03x =，要求计算有5位有效数字。
9．设alimk 。
10．用快速弦截求方程3310x x --=在初始值02x =邻近的实根（取1 1.9x =，要求精确到310-）。
11、分别用下列方法求方程4cos x x e =在04x π=邻近的根，要求有三位有效数字。
（1）用牛顿法，取04x π=； （2）用弦截法，取01,42x x ππ==；（3）用快速弦截法，取01,42x x ππ==。
习 题 三 解 答1、用高斯消元法解下列方程组。
（1）12312312231425427x x x x x x x x -+=⎧⎪++=⎨⎪+=⎩①②③2、将矩阵1020011120110011A ⎛⎫ ⎪⎪= ⎪- ⎪⎝⎭ 作LU 分解。
3、用LU 分解紧凑格式分解法解方程组。
123457910168109171087157651x x x x ⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪ ⎪ ⎪ ⎪= ⎪ ⎪ ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭ 4、 用列主元的三角分解法解线性方程组。
1231231232213472320x x x x x x x x x -+-=-⎧⎪-+=⎨⎪--=⎩ 5．用追赶法解方程组2111210,12101210120A b -⎛⎫⎛⎫ ⎪ ⎪-- ⎪ ⎪⎪ ⎪==-- ⎪ ⎪-- ⎪ ⎪⎪ ⎪-⎝⎭⎝⎭。
6、用改进的Cholesky 分解法解方程组1231231234241021710341097x x x x x x x x x --=⎧⎪-++=⎨⎪-++=-⎩7、用改进的Cholesky 分解法解方程组 123441107131081152400246x x x x -⎛⎫⎛⎫⎛⎫ ⎪ ⎪ ⎪- ⎪ ⎪ ⎪= ⎪ ⎪ ⎪--- ⎪ ⎪ ⎪⎝⎭⎝⎭⎝⎭8、设(1,2,3)T x =-，求12,,x x x ∞9、设110223541A ⎛⎫ ⎪=- ⎪ ⎪⎝⎭，求12,,A A A ∞。
10、设1101223,35412A x -⎛⎫⎛⎫ ⎪ ⎪=-= ⎪ ⎪ ⎪ ⎪⎝⎭⎝⎭，计算,,x A Ax ∞∞∞，并比较Ax ∞与A x∞∞的大小。
11、给定方程组12312212111022110x x x --⎛⎫⎛⎫⎛⎫⎪⎪ ⎪--= ⎪⎪ ⎪ ⎪⎪ ⎪--⎝⎭⎝⎭⎝⎭。
（1）写出雅可比迭代法和高斯－赛德尔迭代法的迭代格式；（2）证明雅可比迭代法收敛而高斯－赛德尔迭代法发散。
（3）给定(0)(0,0,0)T x =，用迭代法求出该方程组的解，精确到(1)()31102k k x x +-∞-≤⨯。
12、给定方程组123211*********x x x ⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪= ⎪⎪ ⎪ ⎪⎪ ⎪⎝⎭⎝⎭⎝⎭。
（1）写出Jacobi 和Gauss-Seidel 迭代格式。
（2）证明Jacobi 迭代法发散而Gauss-Seidel 迭代法收敛。
（3）给定(0)(0,0,0)T x =，用迭代法求出该方程的解，精确到(1)()31102k k x x +-∞-≤⨯ 13、已知1234123412341234541012581034x x x x x x x x x x x x x x x x ---=-⎧⎪-+--=⎪⎨--+-=⎪⎪---+=⎩，考察Jacobi 迭代格式和Gauss-Seidel 迭代格式的收敛性。
14、方程组Ax=b ，其中31410,,01a a A a x b R a ⎛⎫ ⎪=∈ ⎪ ⎪⎝⎭利用迭代收敛的充分必要条件确定雅可比迭代法和高斯－赛德尔迭代法均收敛的a 的取值范围。
15、设方程组123430243413001424x x x ⎛⎫⎛⎫⎛⎫ ⎪⎪ ⎪-= ⎪⎪ ⎪ ⎪⎪ ⎪--⎝⎭⎝⎭⎝⎭，分别用Gauss-Seidel 迭代法和ω＝1.25的SOR 法求解此方程组，准确到4位有效数字（取(0)(1,1,1)T x =）。
17、设100999998A ⎛⎫= ⎪⎝⎭，计算A 的条件数(),2,p cond A p =∞。
18、设A 是n 阶非奇异方阵，B 是n 阶奇异方阵，试证明1()A B cond A A-≤。
19、举例说明一个非奇异矩阵不一定存在LU 分解。
习 题 四 解 答1、设010,1x x ==，写出()x f x e -=的一次插值多项式1()L x ，并估计插值误差。
2、给定函数表i x-0.1 0.3 0.7 1.1 ()i f x0.9950.9950.7650.454选用合适的三次插值多项式来近似计算f(0.2)和f(0.8)。
3、设(0,1,2,,)i x i n =是 n+1个互异节点，证明：（1）0()(0,1,2,,)nk k i i i x l x x k n ===∑；（2）0()()0(0,1,2,,)nk i i i x x l x k n =-==∑。
式中l i(x)是n 次Lagrange 插值基函数。
4、给定数据（()f x =x 2.0 2.1 2.2 2.4 f(x)1.4142141.4491381.483201.54919(1)试用线性插值计算f(2.3)的近似值，并估计误差；(2)试用二次Newton 插值多项式计算f(2.15)的近似值，并估计误差。
5、给定函数表x 0 1 2 4 5 y164688试求各阶差商，并写出牛顿插值多项式和插值余项。
5*、给定函数表x 0 1 2 3 4 y164688试求各阶差分，并求等距节点插值。
6试用三次牛顿差分插值公式计算f(0.1581)及f(0.636)。
7、设f(x)在[-4,4]有连续的4阶导数，且(1)1,(0)2,(0)0,(3)1,(3)1f f f f f ''-===== (1)试构造一个次数最低的插值多项式p(x)，使其满足(1)(1)1,(0)(0)2,(0)(0)0,(3)(3)1,(3)(3)1p f p f p f p f p f ''''-=-=-======== (2)给出并证明余项f(x)-p(x)的表达式。
7*、设f(x)在[-4,4]有连续的4阶导数，且(0)2,(0)0,(3)1,(3)1f f f f ''====试用两种方法构造三次埃尔米特插值多项式H(x)，使其满足 (0)(0)2,(0)(0)0,(3)(3)1,(3)(3)1p f p f p f p f ''''========。
8、设()(01)x f x e x =≤≤，试作一个二次多项式p(x)，使其满足(0)(0),(0)(0),(1)(1)p f p f p f ''===，并导出余项估计式。
9、给出sinx 在[0,π]上的等距结点函数表，用线性插值计算sinx 的近似值，使其截断误差为41102-⨯，问该函数表的步长h 取多少才能满足要求？10、求4()f x x =在区间[a,b]上的分段埃尔米特插值，并估计误差。