数值计算课后答案
《数值计算方法》习题答案
《数值计算方法》课后题答案详解吉 林 大 学第一章 习 题 答 案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 ===处的函数值，试通过一个二次插值函数求的近似值，并估计其误差。
(湖南大学-曾金平《数值计算方法》课后题答案)
1习题一1．设x>0相对误差为2%，4x的相对误差。
解：由自变量的误差对函数值引起误差的公式：(())(())'()()()()f x xf x f x xf x f xδδ∆=≈得（1）()f x=11()()*2%1%22x xδδδ≈===；（2）4()f x x=时444()()'()4()4*2%8%xx x x xxδδδ≈===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.345612210⨯，故（2）的计算结果较精确。
4．计算正方形面积时，若要求面积的允许相对误差为1%，测量边长所允许的相对误差限为多少？2解：设该正方形的边长为x，面积为2()f x x=，由(())(())'()()()()f x xf x f x xf x f xδδ∆=≈解得(())()()'()f x f xxxf xδδ≈=2(())(())22f x x f xx xδδ==0.5%5．下面计算y的公式哪个算得准确些？为什么？（1）已知1x<<，（A）11121xyx x-=-++，（B）22(12)(1)xyx x=++；（2）已知1x>>，（A）y=，（B）y=；（3）已知1x<<，（A）22sin xyx=，（B）1cos2xyx-=；（4）（A）9y=（B）y=解：当两个同（异）号相近数相减（加）时，相对误差可能很大，会严重丧失有效数字；当两个数相乘（除）时，大因子（小除数）可能使积（商）的绝对值误差增大许多。
《数值计算方法》 课后题 答案(曾金平)湖南大学
= 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
数值计算方法第一章答案
数值计算方法第一章答案
2.（1）3580
绝对误差限：；相对误差限：
经过四舍五入得到的近似值2580，其各位都是有效数字，故有四位有效数字。
（2）0.00476
绝对误差限：；相对误差限：
有三位有效数字。
（3）
绝对误差限：；相对误差限：
精确到小数点后两位，所以有四位有效数字。
或者
近似值有四位有效数字。
（4）也可理解为)
(对于14300000精确到小数点前四位)
有四位有效数字。
3.解：取为最好用
，
准确到小数点后第二位，有三位有效数字
准确到小数点后位，有单位有效数字。
6.解：
方法一：
，
取
（）
由此，取4位有效数字.
方法二（有问题，看看错在什么地方？）：
即
则取三位有效数字，
验证
满足要求，所以应取三位有效数字。
8.解：
又
或
11.解：由求根公式
解得
或。
数值计算方法课后习题答案
第一章 绪论（12）1、设0>x ，x 的相对误差为δ，求x ln 的误差。
[解]设0*>x 为x 的近似值，则有相对误差为δε=)(*x r ，绝对误差为**)(x x δε=，从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x ， 相对误差为****ln ln )(ln )(ln x x x x rδεε==。
2、设x 的相对误差为2%，求n x 的相对误差。
[解]设*x 为x 的近似值，则有相对误差为%2)(*=x r ε，绝对误差为**%2)(x x =ε，从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε，相对误差为%2)()(ln )(ln ***n x x x nr==εε。
3、下列各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出它们是几位有效数字：1021.1*1=x ，031.0*2=x ，6.385*3=x ，430.56*4=x ，0.17*5⨯=x 。
[解]1021.1*1=x 有5位有效数字；0031.0*2=x 有2位有效数字；6.385*3=x 有4位有效数字；430.56*4=x 有5位有效数字；0.17*5⨯=x 有2位有效数字。
4、利用公式（3.3）求下列各近似值的误差限，其中*4*3*2*1,,,x x x x 均为第3题所给的数。
（1）*4*2*1x x x ++； [解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε；（2）*3*2*1x x x ；[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε；（3）*4*2/x x 。
数值课后题答案
习 题 一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)证：因为此方程组的系数矩阵为严格对角占优矩阵，所以雅可比迭代法和相应的高斯－赛德尔迭代法都收敛。
丁丽娟《数值计算方法》五章课后实验题答案(源程序很详细,且运行无误)
丁丽娟《数值计算方法》五章课后实验题答案（源程序都是自己写的，很详细，且保证运行无误）我做的五章数值实验作业题目如下：第二章：1、2、3、4题第三章：1、2题第四章：1、2题第六章：2、3题第八章：1、2题第二章1：(1)对A进行列主元素三角分解：function [l u]=myfun(A)n=size(A);for k=1:nfor i=k:nsum=0;m=k;for j=1:(k-1)sum=sum+A(i,j)*A(j,k);ends(i)=A(i,k)-sum;if abs(s(m))<abs(s(i))m=i;endendfor j=1:nc=A(m,j);A(m,j)=A(k,j);A(k,j)=c;endfor j=k:nsum=0;for r=1:(k-1)sum=sum+A(k,r)*A(r,j);endu(k,j)=A(k,j)-sum;A(k,j)=u(k,j);endfor i=1:nl(i,i)=1;endfor i=(k+1):nsum=0;for r=1:(k-1)sum=sum+A(i,r)*u(r,k);endl(i,k)=(A(i,k)-sum)/u(k,k);A(i,k)=l(i,k);endend求A的列主元素三角分解：>>A=[1 1 1 1 1;1 2 3 4 5;1 3 6 10 15;1 4 10 20 35;1 5 15 35 70]; >>[L,U]=myfun(A)结果：L =1.0000 0 0 0 01.0000 1.0000 0 0 01.0000 0.5000 1.0000 0 01.0000 0.7500 0.7500 1.0000 01.0000 0.2500 0.7500 -1.0000 1.0000U =1.0000 1.0000 1.0000 1.0000 1.00000 4.0000 14.0000 34.0000 69.00000 0 -2.0000 -8.0000 -20.50000 0 0 -0.5000 -2.37500 0 0 0 -0.2500(2)求矩阵的逆矩阵A-1:inv(A)结果为：ans =5 -10 10 -5 1-10 30 -35 19 -410 -35 46 -27 6-5 19 -27 17 -41 -4 6 -4 1（3）检验结果：E=diag([1 1 1 1 1])A\Eans =5 -10 10 -5 1-10 30 -35 19 -410 -35 46 -27 6-5 19 -27 17 -41 -4 6 -4 1 2:程序：function d=myfun(a,b,c,d,n)for i=2:nl(i)=a(i)/b(i-1);a(i)=l(i);u(i)=b(i)-c(i-1)*a(i);b(i)=u(i);y(i)=d(i)-a(i)*d(i-1);d(i)=y(i);endx(n)=d(n)/b(n);d(n)=x(n);for i=(n-1):-1:1x(i)=(d(i)-c(i)*d(i+1))/b(i);d(i)=x(i);end求各段电流量程序：for i=2:8endb=[2 5 5 5 5 5 5 5];c=[-2 -2 -2 -2 -2 -2 -2];V=220;R=27;d=[V/R 0 0 0 0 0 0 0];n=8;I=myfun(a,b,c,d,n)运行程序得：I =8.1478 4.0737 2.0365 1.0175 0.5073 0.2506 0.1194 0.04773：（1）求矩阵A和向量b的matlab程序：function [A b]=myfun(n)for i=1:nX(i)=1+0.1*i;endfor i=1:nfor j=1:nA(i,j)=X(i)^(j-1);endfor i=1:nb(i)=sum(A(i,:));end求n=5时A1，b1及A1的2-条件数程序运行结果如下：n=5；[A1,b1]=myfun(n)A1 =1.0000 1.1000 1.2100 1.3310 1.46411.0000 1.2000 1.4400 1.72802.07361.0000 1.3000 1.69002.1970 2.85611.0000 1.4000 1.96002.74403.84161.0000 1.50002.25003.3750 5.0625 b1 =6.10517.4416 9.0431 10.9456 13.1875cond2=cond(A1,2)cond2 =5.3615e+005求n=10时A2，b2及A2的2-条件数程序运行结果如下：n=10;[A2,b2]=myfun(n)A2 =1.0000 1.1000 1.2100 1.3310 1.4641 1.6105 1.7716 1.94872.1436 2.35791.0000 1.2000 1.4400 1.72802.0736 2.4883 2.98603.58324.29985.15981.0000 1.3000 1.69002.1970 2.85613.71294.8268 6.2749 8.1573 10.60451.0000 1.4000 1.96002.74403.8416 5.3782 7.5295 10.5414 14.7579 20.66101.0000 1.50002.25003.3750 5.0625 7.5938 11.3906 17.0859 25.6289 38.44341.0000 1.60002.5600 4.0960 6.5536 10.4858 16.7772 26.8435 42.9497 68.71951.0000 1.70002.8900 4.9130 8.3521 14.1986 24.1376 41.0339 69.7576 118.58791.0000 1.8000 3.2400 5.8320 10.4976 18.8957 34.0122 61.2220 110.1996 198.35931.0000 1.9000 3.6100 6.8590 13.0321 24.7610 47.0459 89.3872 169.8356 322.68771.00002.0000 4.0000 8.0000 16.0000 32.0000 64.0000 128.0000 256.0000 512.0000b2 =1.0e+003 *0.0159 0.0260 0.0426 0.0698 0.1133 0.1816 0.2866 0.4451 0.6801 1.0230cond2=cond(A2,2)cond2 =8.6823e+011求n=20时A3，b3及A3的2-条件数程序运行结果如下：n=20;[A3,b3]=myfun(n)A3 =1.0e+009 *Columns 1 through 100.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 11 through 200.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0002 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0005 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0006 0.0013 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0007 0.0015 0.0032 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0006 0.0014 0.0032 0.0075 0.0000 0.0000 0.0000 0.0001 0.0002 0.0005 0.0012 0.0029 0.0070 0.0167 0.0000 0.0000 0.0001 0.0001 0.0004 0.0009 0.0023 0.0058 0.0146 0.0364 0.0000 0.0000 0.0001 0.0002 0.0006 0.0017 0.0044 0.0113 0.0295 0.0766 0.0000 0.0001 0.0002 0.0004 0.0011 0.0030 0.0080 0.0215 0.0581 0.1570 0.0000 0.0001 0.0002 0.0007 0.0018 0.0051 0.0143 0.0400 0.1119 0.31330.0000 0.0001 0.0004 0.0010 0.0030 0.0086 0.0250 0.0726 0.2105 0.61030.0001 0.0002 0.0005 0.0016 0.0048 0.0143 0.0430 0.1291 0.3874 1.1623b3 =1.0e+009 *Columns 1 through 100.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0002 0.0004 0.0010Columns 11 through 200.0025 0.0059 0.0132 0.0287 0.0606 0.1246 0.2494 0.4874 0.9316 1.7434cond2=cond(A3,2)cond2 =3.2395e+022由上述运行结果可知：它们是病态的，而且随着n的增大，矩阵的病态变得严重。
数值计算方法答案
1数值计算方法习题一（2）习题二（6）习题三（15）习题四（29）习题五（37）习题六（62）习题七（70）2009．9，92习题一1．设x>0相对误差为2%，4x的相对误差。
解：由自变量的误差对函数值引起误差的公式：(())(())'()()()()f x xf x f x xf x f xδδ∆=≈得（1）()f x=11()()*2%1%22x xδδδ≈===；（2）4()f x x=时444()()'()4()4*2%8%xx x x xxδδδ≈===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.345612210⨯，故（2）的计算结果较精确。
4．计算正方形面积时，若要求面积的允许相对误差为1%，测量边长所允许的相对误差限为多少？3解：设该正方形的边长为x，面积为2()f x x=，由(())(())'()()()()f x xf x f x xf x f xδδ∆=≈解得(())()()'()f x f xxxf xδδ≈=2(())(())22f x x f xx xδδ==0.5%5．下面计算y的公式哪个算得准确些？为什么？（1）已知1x<<，（A）11121xyx x-=-++，（B）22(12)(1)xyx x=++；（2）已知1x>>，（A）y=，（B）y=；（3）已知1x<<，（A）22sin xyx=，（B）1cos2xyx-=；（4）（A）9y=（B）y=解：当两个同（异）号相近数相减（加）时，相对误差可能很大，会严重丧失有效数字；当两个数相乘（除）时，大因子（小除数）可能使积（商）的绝对值误差增大许多。
数值计算课后答案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
习 题 一 解 答1．取3.14，3.15，227，355113作为π的近似值，求各自的绝对误差，相对误差和有效数字的位数。
分析：求绝对误差的方法是按定义直接计算。
求相对误差的一般方法是先求出绝对误差再按定义式计算。
注意，不应先求相对误差再求绝对误差。
有效数字位数可以根据定义来求，即先由绝对误差确定近似数的绝对误差不超过那一位的半个单位，再确定有效数的末位是哪一位，进一步确定有效数字和有效数位。
有了定理2后，可以根据定理2更规范地解答。
根据定理2，首先要将数值转化为科学记数形式，然后解答。
解：（1）绝对误差:e(x)=π－3.14＝3.14159265…－3.14＝0.00159…≈0.0016。
相对误差：3()0.0016()0.51103.14r e x e x x -==≈⨯ 有效数字：因为π＝3.14159265…=0.314159265…×10，3.14＝0.314×10，m=1。
而π－3.14＝3.14159265…－3.14＝0.00159…所以│π－3.14│＝0.00159…≤0.005=0.5×10－2＝21311101022--⨯=⨯ 所以，3.14作为π的近似值有3个有效数字。
（2）绝对误差:e(x)=π－3.15＝3.14159265…－3.14＝－0.008407…≈－0.0085。
相对误差：2()0.0085()0.27103.15r e x e x x --==≈-⨯ 有效数字：因为π＝3.14159265…=0.314159265…×10，3.15＝0.315×10，m=1。
而π－3.15＝3.14159265…－3.15＝－0.008407…所以│π－3.15│＝0.008407……≤0.05=0.5×10－1＝11211101022--⨯=⨯ 所以，3.15作为π的近似值有2个有效数字。
（3）绝对误差:22() 3.14159265 3.1428571430.0012644930.00137e x π=-=-=-≈-L L 相对误差：3()0.0013()0.4110227r e x e x x--==≈-⨯ 有效数字：因为π＝3.14159265...=0.314159265 (10)22 3.1428571430.3142857143107==⨯，m=1。
数值计算课后习题答案
习 题 一 解 答1．取3.14，3.15，227，355113作为π的近似值，求各自的绝对误差，相对误差和有效数字的位数。
分析：求绝对误差的方法是按定义直接计算。
求相对误差的一般方法是先求出绝对误差再按定义式计算。
注意，不应先求相对误差再求绝对误差。
有效数字位数可以根据定义来求，即先由绝对误差确定近似数的绝对误差不超过那一位的半个单位，再确定有效数的末位是哪一位，进一步确定有效数字和有效数位。
有了定理2后，可以根据定理2更规范地解答。
根据定理2，首先要将数值转化为科学记数形式，然后解答。
解：（1）绝对误差:e(x)=π－3.14＝3.14159265…－3.14＝0.00159…≈0.0016。
相对误差：3()0.0016()0.51103.14r e x e x x-==≈⨯有效数字：因为π＝3.14159265…=0.314159265…×10，3.14＝0.314×10，m=1。
而π－3.14＝3.14159265…－3.14＝0.00159…所以│π－3.14│＝0.00159…≤0.005=0.5×10－2＝21311101022--⨯=⨯所以，3.14作为π的近似值有3个有效数字。
（2）绝对误差:e(x)=π－3.15＝3.14159265…－3.14＝－0.008407…≈－0.0085。
相对误差：2()0.0085()0.27103.15r e x e x x--==≈-⨯有效数字：因为π＝3.14159265…=0.314159265…×10，3.15＝0.315×10，m=1。
而π－3.15＝3.14159265…－3.15＝－0.008407…所以│π－3.15│＝0.008407……≤0.05=0.5×10－1＝11211101022--⨯=⨯所以，3.15作为π的近似值有2个有效数字。
（3）绝对误差:22() 3.14159265 3.1428571430.0012644930.00137e x π=-=-=-≈-相对误差：3()0.0013()0.4110227r e x e x x--==≈-⨯有效数字：因为π＝3.14159265...=0.314159265 (10)22 3.1428571430.3142857143107==⨯，m=1。
数值计算方法(宋岱才版)课后答案
第一章 绪论一 本章的学习要求（1）会求有效数字。
（2）会求函数的误差及误差限。
（3）能根据要求进行误差分析。
二 本章应掌握的重点公式（1）绝对误差：设x 为精确值，x *为x 的一个近似值，称e x x **=-为x *的绝对误差。
（2）相对误差：r e e x***=。
（3）绝对误差限：e x x ε***==-。
（4）相对误差限：r x x xxεε*****-==。
（5）一元函数的绝对误差限：设一元函数()()()0,df f x f x dx εε***⎛⎫==⋅ ⎪⎝⎭则。
（6）一元函数的相对误差限：()()1r df f x dx f εε****⎛⎫=⋅ ⎪⎝⎭。
（7）二元函数的绝对误差限：设一元函数()()(),0,f f x y f y y εε***⎛⎫∂==⋅ ⎪∂⎝⎭则。
（8）二元函数的相对误差限：()()()1r f f f x y x y f εεε******⎡⎤⎛⎫∂∂⎛⎫⎢⎥=⋅+⋅ ⎪ ⎪∂∂⎝⎭⎢⎥⎝⎭⎣⎦。
三 本章习题解析1. 下列各数都是经过四舍五入得到的近似值，（1）试指出它们有几位有效数字，（2）分别估计1123A X X X ***=及224X A X **=的相对误差限。
12341.1021,0.031,385.6,56.430x x x x ****====解：（1）1x *有5位有效数字，2x *有2位有效数字，3x *有4位有效数字，4x *有5位有效数字。
（2）1111123231312123,,,,A A AA x x x x x x x x x x x x ∂∂∂====∂∂∂由题可知：1A *为1A 的近似值，123,,x x x ***分别为123,,x x x 近似值。
所以()()111rA A Aεε***=()()()12311111123A A A x x x A X X X εεε*******⎡⎤⎢⎥=++⎢⎥⎢⎥⎣⎦⎛⎫⎛⎫⎛⎫∂∂∂ ⎪ ⎪ ⎪∂∂∂⎝⎭⎝⎭⎝⎭43123131212311111010100.215222x x x x x x x x x **-**-**-***⎡⎤=⨯⨯+⨯⨯+⨯⨯=⎢⎥⎣⎦()222222424441,,,X A Ax A X x x x x ∂∂===-∂∂则有同理有2A *为2A 的近似值，2x *，4x *为2x ，4x 的近似值，代入相对误差限公式：()()222rA A Aεε***=()()24212224A A X X A X X εε*****⎡⎤⎢⎥=+⎢⎥⎢⎥⎣⎦⎛⎫⎛⎫∂∂ ⎪ ⎪∂∂⎝⎭⎝⎭()33542224411*********X X X X X **--***⎡⎤⎢⎥=⨯⨯+⨯⨯=⎢⎥⎣⎦2. 正方形的边长大约为100cm ，怎样测量才能使其面积误差不超过21cm ? 解：设正方形的边长为x ，则面积为2S x =，2dsx dx=，在这里设x *为边长的近似值，S *为面积的近似值：由题可知：()()1ds s x dx εε***=≤⎛⎫ ⎪⎝⎭即：()21x x ε**⋅≤ 推出：()10.005200xcm ε*≤=。
数值计算课后规范标准答案2
习 题 二 解 答1．用二分法求方程x 3-2x 2-4x-7=0在区间[3，4]内的根，精确到10-3，即误差不超过31102-⨯。
分析：精确到10-3与误差不超过10-3不同。
解：因为f(3)＝-10＜0，f(4)=9＞0，所以，方程在区间［3，4］上有根。
由34311*1022222n n n n n n b a b a x x -----≤===<⨯ 有2n-1＞1000，又为210＝1024＞1000， 所以n ＝11，即只需要二分11次即可。
x *≈x 11=3.632。
指出：（1）注意精确度的不同表述。
精确到10-3和误差不超过10-3是不同的。
（2）在计算过程中按规定精度保留小数，最后两次计算结果相同。
如果计算过程中取4位小数，结果取3位，则如下表：（3）用秦九韶算法计算f(x n )比较简单。
1*．求方程x 3-2x 2-4x-7=0的隔根区间。
解：令32247y x x x =---， 则2344322()()y x x x x '=--=+-当23443220()()y x x x x '=--=+-=时，有12223,x x =-=。
函数单调区间列表分析如下：因为214902150327(),()y y -=-<=-<，所以方程在区间223(,)-上无根； 因为21490327()y -=-<，而函数在23(,)-∞-上单调增，函数值不可能变号，所以方程在该区间上无根；因为2150()y =-<，函数在(2,+∞)上单调增，所以方程在该区间上最多有一个根，而(3)=-10<0,y(4)=9>0，所以方程在区间(3,4)有一个根。
所以，该方程有一个根，隔根区间是(3.4)。
2．证明1sin 0x x --=在[0,1]内有一个根，使用二分法求误差不大于41102-⨯的根，需要迭代多少次？分析：证明方程在指定区间内有一个根，就是证明相应的函数在指定区间有至少一个零点。
丁丽娟《数值计算方法》五章课后实验题答案(源程序很详细,且运行无误)
丁丽娟《数值计算方法》五章课后实验题答案（源程序都是自己写的，很详细，且保证运行无误）我做的五章数值实验作业题目如下：第二章：1、2、3、4题第三章：1、2题第四章：1、2题第六章：2、3题第八章：1、2题第二章1：(1)对A进行列主元素三角分解：function [l u]=myfun(A)n=size(A);for k=1:nfor i=k:nsum=0;m=k;for j=1:(k-1)sum=sum+A(i,j)*A(j,k);ends(i)=A(i,k)-sum;if abs(s(m))<abs(s(i))m=i;endendfor j=1:nc=A(m,j);A(m,j)=A(k,j);A(k,j)=c;endfor j=k:nsum=0;for r=1:(k-1)sum=sum+A(k,r)*A(r,j);endu(k,j)=A(k,j)-sum;A(k,j)=u(k,j);endfor i=1:nl(i,i)=1;endfor i=(k+1):nsum=0;for r=1:(k-1)sum=sum+A(i,r)*u(r,k);endl(i,k)=(A(i,k)-sum)/u(k,k);A(i,k)=l(i,k);endend求A的列主元素三角分解：>>A=[1 1 1 1 1;1 2 3 4 5;1 3 6 10 15;1 4 10 20 35;1 5 15 35 70]; >>[L,U]=myfun(A)结果：L =1.0000 0 0 0 01.0000 1.0000 0 0 01.0000 0.5000 1.0000 0 01.0000 0.7500 0.7500 1.0000 01.0000 0.2500 0.7500 -1.0000 1.0000U =1.0000 1.0000 1.0000 1.0000 1.00000 4.0000 14.0000 34.0000 69.00000 0 -2.0000 -8.0000 -20.50000 0 0 -0.5000 -2.37500 0 0 0 -0.2500(2)求矩阵的逆矩阵A-1:inv(A)结果为：ans =5 -10 10 -5 1-10 30 -35 19 -410 -35 46 -27 6-5 19 -27 17 -41 -4 6 -4 1（3）检验结果：E=diag([1 1 1 1 1])A\Eans =5 -10 10 -5 1-10 30 -35 19 -410 -35 46 -27 6-5 19 -27 17 -41 -4 6 -4 1 2:程序：function d=myfun(a,b,c,d,n)for i=2:nl(i)=a(i)/b(i-1);a(i)=l(i);u(i)=b(i)-c(i-1)*a(i);b(i)=u(i);y(i)=d(i)-a(i)*d(i-1);d(i)=y(i);endx(n)=d(n)/b(n);d(n)=x(n);for i=(n-1):-1:1x(i)=(d(i)-c(i)*d(i+1))/b(i);d(i)=x(i);end求各段电流量程序：for i=2:8endb=[2 5 5 5 5 5 5 5];c=[-2 -2 -2 -2 -2 -2 -2];V=220;R=27;d=[V/R 0 0 0 0 0 0 0];n=8;I=myfun(a,b,c,d,n)运行程序得：I =8.1478 4.0737 2.0365 1.0175 0.5073 0.2506 0.1194 0.04773：（1）求矩阵A和向量b的matlab程序：function [A b]=myfun(n)for i=1:nX(i)=1+0.1*i;endfor i=1:nfor j=1:nA(i,j)=X(i)^(j-1);endfor i=1:nb(i)=sum(A(i,:));end求n=5时A1，b1及A1的2-条件数程序运行结果如下：n=5；[A1,b1]=myfun(n)A1 =1.0000 1.1000 1.2100 1.3310 1.46411.0000 1.2000 1.4400 1.72802.07361.0000 1.3000 1.69002.1970 2.85611.0000 1.4000 1.96002.74403.84161.0000 1.50002.25003.3750 5.0625 b1 =6.10517.4416 9.0431 10.9456 13.1875cond2=cond(A1,2)cond2 =5.3615e+005求n=10时A2，b2及A2的2-条件数程序运行结果如下：n=10;[A2,b2]=myfun(n)A2 =1.0000 1.1000 1.2100 1.3310 1.4641 1.6105 1.7716 1.94872.1436 2.35791.0000 1.2000 1.4400 1.72802.0736 2.4883 2.98603.58324.29985.15981.0000 1.3000 1.69002.1970 2.85613.71294.8268 6.2749 8.1573 10.60451.0000 1.4000 1.96002.74403.8416 5.3782 7.5295 10.5414 14.7579 20.66101.0000 1.50002.25003.3750 5.0625 7.5938 11.3906 17.0859 25.6289 38.44341.0000 1.60002.5600 4.0960 6.5536 10.4858 16.7772 26.8435 42.9497 68.71951.0000 1.70002.8900 4.9130 8.3521 14.1986 24.1376 41.0339 69.7576 118.58791.0000 1.8000 3.2400 5.8320 10.4976 18.8957 34.0122 61.2220 110.1996 198.35931.0000 1.9000 3.6100 6.8590 13.0321 24.7610 47.0459 89.3872 169.8356 322.68771.00002.0000 4.0000 8.0000 16.0000 32.0000 64.0000 128.0000 256.0000 512.0000b2 =1.0e+003 *0.0159 0.0260 0.0426 0.0698 0.1133 0.1816 0.2866 0.4451 0.6801 1.0230cond2=cond(A2,2)cond2 =8.6823e+011求n=20时A3，b3及A3的2-条件数程序运行结果如下：n=20;[A3,b3]=myfun(n)A3 =1.0e+009 *Columns 1 through 100.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000Columns 11 through 200.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0002 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0005 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0006 0.0013 0.0000 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0007 0.0015 0.0032 0.0000 0.0000 0.0000 0.0001 0.0001 0.0003 0.0006 0.0014 0.0032 0.0075 0.0000 0.0000 0.0000 0.0001 0.0002 0.0005 0.0012 0.0029 0.0070 0.0167 0.0000 0.0000 0.0001 0.0001 0.0004 0.0009 0.0023 0.0058 0.0146 0.0364 0.0000 0.0000 0.0001 0.0002 0.0006 0.0017 0.0044 0.0113 0.0295 0.0766 0.0000 0.0001 0.0002 0.0004 0.0011 0.0030 0.0080 0.0215 0.0581 0.1570 0.0000 0.0001 0.0002 0.0007 0.0018 0.0051 0.0143 0.0400 0.1119 0.31330.0000 0.0001 0.0004 0.0010 0.0030 0.0086 0.0250 0.0726 0.2105 0.61030.0001 0.0002 0.0005 0.0016 0.0048 0.0143 0.0430 0.1291 0.3874 1.1623b3 =1.0e+009 *Columns 1 through 100.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0001 0.0002 0.0004 0.0010Columns 11 through 200.0025 0.0059 0.0132 0.0287 0.0606 0.1246 0.2494 0.4874 0.9316 1.7434cond2=cond(A3,2)cond2 =3.2395e+022由上述运行结果可知：它们是病态的，而且随着n的增大，矩阵的病态变得严重。
数值计算方法教材答案
1习题一1．设x>0相对误差为2%，4x的相对误差。
解：由自变量的误差对函数值引起误差的公式：(())(())'()()()()f x xf x f x xf x f xδδ∆=≈得（1）()f x=11()()*2%1%22x xδδδ≈===；（2）4()f x x=时444()()'()4()4*2%8%xx x x xxδδδ≈===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.345612210⨯，故（2）的计算结果较精确。
4．计算正方形面积时，若要求面积的允许相对误差为1%，测量边长所允许的相对误差限为多少？2解：设该正方形的边长为x，面积为2()f x x=，由(())(())'()()()()f x xf x f x xf x f xδδ∆=≈解得(())()()'()f x f xxxf xδδ≈=2(())(())22f x x f xx xδδ==0.5%5．下面计算y的公式哪个算得准确些？为什么？（1）已知1x<<，（A）11121xyx x-=-++，（B）22(12)(1)xyx x=++；（2）已知1x>>，（A）y=，（B）y=；（3）已知1x<<，（A）22sin xyx=，（B）1cos2xyx-=；（4）（A）9y=（B）y=解：当两个同（异）号相近数相减（加）时，相对误差可能很大，会严重丧失有效数字；当两个数相乘（除）时，大因子（小除数）可能使积（商）的绝对值误差增大许多。
数值计算方法课后习题答案(李庆扬等)
第一章 绪论（12）1、设0>x ，x 的相对误差为δ，求x ln 的误差。
[解]设0*>x 为x 的近似值，则有相对误差为δε=)(*x r ，绝对误差为**)(x x δε=，从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x x x x x ， 相对误差为****ln ln )(ln )(ln x x x x rδεε==。
2、设x 的相对误差为2%，求n x 的相对误差。
[解]设*x 为x 的近似值，则有相对误差为%2)(*=x r ε，绝对误差为**%2)(x x =ε，从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε，相对误差为%2)()(ln )(ln ***n x x x nr==εε。
3、下列各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出它们是几位有效数字：1021.1*1=x ，031.0*2=x ，6.385*3=x ，430.56*4=x ，0.17*5⨯=x 。
[解]1021.1*1=x 有5位有效数字；0031.0*2=x 有2位有效数字；6.385*3=x 有4位有效数字；430.56*4=x 有5位有效数字；0.17*5⨯=x 有2位有效数字。
4、利用公式（3.3）求下列各近似值的误差限，其中*4*3*2*1,,,x x x x 均为第3题所给的数。
（1）*4*2*1x x x ++； [解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε；（2）*3*2*1x x x ；[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε；（3）*4*2/x x 。
数值计算方法》习题答案
《数值计算方法》课后题答案详解吉 林 大 学第一章 习 题 答 案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 ===处的函数值，试通过一个二次插值函数求的近似值，并估计其误差。
数值计算课后答案4
习 题 四 解 答1、设010,1x x ==，写出()x f x e -=的一次插值多项式1()L x ，并估计插值误差。
解:根据已知条件，有设插值函数为1()L x ax b =+，由插值条件，建立线性方程组为1011a b a b e -⨯+=⎧⎨⨯+=⎩ 解之得111a eb -⎧=-⎨=⎩则11()(1)1L x e x -=-+ 因为(),()x x y x e y x e --'''=-= 所以,插值余项为(1)(2)(2)011()()()()()(1)!1()()2!1()()()2!1(0)(1)((0,1))2n r x f x p x f x n f x f x x x x e x x ξξπξπξξ+-=-=+==--=--∈所以010101()max max (1)2111248x r x e x x e ξξ-≤≤≤≤-≤-=⨯⨯=。
2选用合适的三次插值多项式来近似计算f(0.2）和f(0.8)。
解：设三次插值多项式为230123()f x a a x a x a x =+++，由插值条件，建立方程组为230123230123230123230123(0.1)(0.1)(0.1)0.9950.30.30.30.9950.70.70.70.7651.1 1.1 1.10.454a a a a a a a a a a a a a a a a ⎧+⨯-+⨯-+⨯-=⎪+⨯+⨯+⨯=⎪⎨+⨯+⨯+⨯=⎪⎪+⨯+⨯+⨯=⎩即012301230123123012312301230.10.010.0010.9950.10.010.0010.9950.30.090.0270.9950.40.080.02800.70.490.3430.7650.80.480.344 1.761.1 1.21 1.3310.454a a a a a a a a a a a a a a a a a a a a a a a a a a -+-=-+-=⎧⎪+++=++=⎪⇒⎨+++=++=⎪⎪+++=⎩12301231232330.40.720.9880.3110.10.010.0010.9950.40.080.02800.320.288 1.760.384 3.831a a a a a a a a a a a a a ⎧⎪⎪⎨⎪⎪++=-⎩-+-=⎧⎪++=⎪⇒⎨+=⎪⎪-=-⎩解之得 01230.416.293.489.98a a a a =⎧⎪=-⎪⎨=-⎪⎪=⎩ 则所求的三次多项式为23()0.41 6.29 3.489.98f x x x x =--+. 所以2323(0.2)0.41 6.290.2 3.480.29.980.20.91(0.8)0.41 6.290.8 3.480.89.980.8 1.74f f =-⨯-⨯+⨯=-=-⨯-⨯+⨯=-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,,)n k i i i x x l x k n =-==∑。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
习 题 四 解 答1、设010,1x x ==，写出()x f x e -=的一次插值多项式1()L x ，并估计插值误差。
设插值函数为1()L x ax b =+，由插值条件，建立线性方程组为1011a b a b e -⨯+=⎧⎨⨯+=⎩ 解之得111a e b -⎧=-⎨=⎩则11()(1)1L x e x -=-+ 因为(),()x x y x e y x e --'''=-= 所以，插值余项为(1)(2)(2)011()()()()()(1)!1()()2!1()()()2!1(0)(1)((0,1))2n r x f x p x f x n f x f x x x x e x x ξξπξπξξ+-=-=+==--=--∈所以010101()max max (1)2111248x r x e x x e ξξ-≤≤≤≤-≤-=⨯⨯=。
2选用合适的三次插值多项式来近似计算f 和f 。
解：设三次插值多项式为230123()f x a a x a x a x =+++，由插值条件，建立方程组为230123230123230123230123(0.1)(0.1)(0.1)0.9950.30.30.30.9950.70.70.70.7651.1 1.1 1.10.454a a a a a a a a a a a a a a a a ⎧+⨯-+⨯-+⨯-=⎪+⨯+⨯+⨯=⎪⎨+⨯+⨯+⨯=⎪⎪+⨯+⨯+⨯=⎩即012301230123123012312301230.10.010.0010.9950.10.010.0010.9950.30.090.0270.9950.40.080.02800.70.490.3430.7650.80.480.344 1.761.1 1.21 1.3310.454a a a a a a a a a a a a a a a a a a a a a a a a a a -+-=-+-=⎧⎪+++=++=⎪⇒⎨+++=++=⎪⎪+++=⎩12301231232330.40.720.9880.3110.10.010.0010.9950.40.080.02800.320.288 1.760.384 3.831a a a a a a a a a a a a a ⎧⎪⎪⎨⎪⎪++=-⎩-+-=⎧⎪++=⎪⇒⎨+=⎪⎪-=-⎩解之得 01230.416.293.489.98a a a a =⎧⎪=-⎪⎨=-⎪⎪=⎩ 则所求的三次多项式为23()0.41 6.29 3.489.98f x x x x =--+。
所以2323(0.2)0.41 6.290.2 3.480.29.980.20.91(0.8)0.41 6.290.8 3.480.89.980.8 1.74f f =-⨯-⨯+⨯=-=-⨯-⨯+⨯=-3、设(0,1,2,,)i x i n =L 是 n+1个互异节点，证明： （1）0()(0,1,2,,)nk k i i i x l x x k n ===∑L ；（2）0()()0(0,1,2,,)n k i i i x x l x k n =-==∑L 。
证明： （1）由拉格朗日插值定理，以x 0,x 1,x 2,…x n 为插值节点，对y=f(x)=x k 作n 次插值，插值多项式为 0()()nn i i i p x l x y ==∑，而y i =x i k ,所以0()()()nnk n i i i i i i p x l x y l x x ====∑∑同时，插值余项(1)(1)11()()()()()()0(1)!(1)!n k n k n r x x p x f x x x n n ξξππ++=-===++所以0()nk k i i i l x x x ==∑结论得证。
（2）取函数()(),0,1,2,,k f x x t k n =-=L对此函数取节点(0,1,2,,)i x i n =L ，则对应的插值多项式为()()()nk n i i i p x x t l x ==-∑，由余项公式，得(1)(1)011()()()()()()()()0(1)!(1)!nn kk n ki i i r x x t x t l x f x x t x n n ξξππ++==---==-=++∑所以0()()()nkk i i i x t x t l x =-=-∑令t=x ，()()0nkiii x x l x =-=∑4、给定数据（()f x =(1)试用线性插值计算f 的近似值，并估计误差；(2)试用二次Newton 插值多项式计算f 的近似值，并估计误差。
解：用线性插值计算f ，取插值节点为和，则相应的线性插值多项式是1.54919 1.48320() 1.48320( 2.2)2.4 2.21.483200.32995(2.2)p x x x -=+--=+- 用x=代入，得(2.3) 1.483200.32995(2.3 2.2) 1.450205f ≈+⨯-= (2)根据定理2f(x)＝f(x 0)＋f[x 0，x 1](x-x 0)+f[x 0，x 1，x 2](x-x 0)(x-x 1)+…+f[x 0，x 1，…，x n ](x-x 0)(x-x 1)…(x-x n －1)+f[x 0，x 1，…，x n ，x]π(x) 。
以表中的上方一斜行中的数为系数，得 f ＝＋ × × × ＝ 指出：误差未讨论。
557()0167(1)(1)(2)(1)(2)(4)26p x x x x x x x x x x x =++--------。
指出： 余项未讨论。
解：由已知条件，显然，x 0=0，h=1，x=t 。
0(1)(1)(2)(1)(2)(3)()()01614(2)(140)2!3!4!(1)(2)35167(1)(1)(2)(3)36n n t t t t t t t t t p x th p t t t t t t t t t t t t ------+==+⨯+⨯+⨯-+⨯---=+------指出:在本题这种情况下,实际上()()n n p t p x =，也就是说,在这样的条件下,t 的多项式就是x 的多项式，可以直接转换。
一般情况下，把t 的关系转换为x 的关系需要根据x=x 0+th ，将t 用x 表示，即将0x x t h-=代入得到的多项式。
6解：所给节点是等距结点：000.125,0.125,,0,1,2,3,4,5i x h x x ih i ===+=。
令00()x x th t h=+=,根据等距结点插值公式，得 0(1)()()0.79618(0.02284)(0.00679)2!(1)(2)(1)(2)(3)(1)(2)(3)(4)(0.00316)0.00488(0.00460)3!4!5!n n t t p x th p t t t t t t t t t t t t t t -+==+⨯-+⨯----------+⨯-+⨯+⨯-则(0.1581)(0.1581)(0.1250.2648)0.790294822,(0.636)(0.6363)(0.125 4.088)0.651804826n n n n f p p h f p p h ≈=+=≈=+=。
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)的表达式。
解：(1)由7*可以求出满足(0)(0)2,(0)(0)0,(3)(3)1,(3)(3)1p f p f p f p f ''''======== 的三次埃尔米特插值多项式3252()2273H x x x =-+。
设22322252()()(3)2(3)273p x H x a x x x x a x x =+-=-++-，则p(x)满足(0)(0)2,(0)(0)0,(3)(3)1,(3)(3)1p f p f p f p f ''''========， 由(1)1f -=得 3222521(1)(1)2(13)(1)1273108a a ⨯--⨯-++---=⇒=-， 所以223222432521()()(3)2(3)27310811332108544p x H x a x x x x x x x x x =+-=-+--=-++-+。
(2)余项具有如下结构22()()()()(1)(3)r x f x p x k x x x x =-=+- 作辅助函数22()()()()(1)(3)t f t p t k x t t t ϕ=--+-则显然()t ϕ在点,1,0,3x -处有6个零点（其中0，3是二重零点），即 ()0,(1)0,(0)0,(0)0,(3)0,(3)0x ϕϕϕϕϕϕ''=-=====， 不妨假设(1,0)x ∈-。
由罗尔定理，存在123(1,),(,0),(0,3)x x ξξξ∈-∈∈， 使得123()0,()0,()0ϕξϕξϕξ'''===，再注意到(0)0,(3)0ϕϕ''==，即()t ϕ'有5个互异的零点12303ξξξ<<<<再次由罗尔定理得，存在111223343(,),(,0),(0,),(,3)ηξξηξηξηξ∈∈∈∈， 使得1234()0,()0,()0,()0ϕηϕηϕηϕη''''''''====第三次应用罗尔定理得，存在112223334(,),(,),(,)ξηηξηηξηη∈∈∈ 使得123()0,()0,()0ϕξϕξϕξ'''''''''===，第四次应用罗尔定理得，存在112223(,),(,)μξξμξξ∈∈ 使得(4)(4)12()0,()0ϕμϕμ==，第五次应用罗尔定理得，存在12(,)τμμ∈ 使得(5)()0ϕτ= 注意到(5)(5)(5)()()5!()()5!()t r t k x f t k x ϕ=-=-（()()()r t f t p t =-中p(t)是4次函数，其5次导数为0）。