计算方法一二章答案

第二章数值分析已知多项式p(x) X1 X3 X2 X 1通过下列点：p(x)试构造一多项式q(x)通过下列点:表中p2(X)的某一个函数值有错误，试找出并校正它•答案：函数值表中P2( 1)错误，应有P2(1)O •利用差分的性质证明12 22n2 n(n 1)(2n1)/6 ・当用等距节点的分段二次插值多项式在区间[1,1]近似函数e x时，使用多少个节点能够保证误差不超过1 1062答案：需要个插值节点・设被插值函数f(x)C4[a,b] 出(叫x)是f(x)矢于等距节点baa Xo X1 Xn b的分段三次艾尔米特插值多项式，步长h •试估计n22I I f (x) H3(h)(x) I I .答案：| |f(x) H3(h) (x) | | M4 hl384第三章函数逼近求f(x) sin x, x [0, 0. 1]在空间span{l, x, x2} ±最佳平方逼近多项式，并给岀平方误差.答案：f (x) sin x的二次最佳平方逼近多项式为-52 sin x p2(x) 0. 832 440 7 10-5 1.000 999 lx 0. 024 985 lx2,二次最佳平方逼近的平方误差为0. 12 2 -12_ (sin x) p2 (x)) dx 0. 989 310 7 10~12・确定参数a, b和c ,使得积分[ax2 bx c 1 ] dx取最小值.l(a,b,c)求多项式f (x) 2x' x3 5x2 1在[1, 1]上的3 次最佳一致逼近多项式p(x) •8 10 a , b 0, c 33答案：f(X)的最佳一致逼近多项式为P(X) ； 7；4用幕级数缩合方法，求f (x) e s ( 1 x 1)上的3次近似多项式p6,3 ( x),并估计I f (x) P6,3(X)I ・答案：23 pe,3 ( x) 0. 994 574 65 0. 997 39583x 0. 542 968 75x2 0. 177 083 33x3,:f (x) P6,3 (x) | | 0. 006 572 327 7J求f (x) e s ( 1 x 1)上的关于权函数(x)的三次最佳平方逼近多1 X"项式S3 ( X),并估计误差I f(X)S3(X)〔2 和I I f(X)S3 (x) I •咎23、口Ss(x) 0. 994 571 0. 997 308x 0. 542 99lx2 0. 177 347 x3,丨丨 f (x) Ss(x) | 12 0. 006 894 83 , | | f (x) Ss( x) | | 0. 006 442 575 ・第四章数值积分与数值微分用梯形公式、辛浦生公式和柯特斯公式分别计算积分x n dx (n 1, 2, 3, 4),并与精确值比较答案：计算结果如下表所示式具有的代数术青度.版权文档，请勿用做商业用途h(1 ) h f (x) dx Aif ( h) Ao f (0) Ai f (h)X1(2 ) if (x) dx [f ( 1) 2f (xi) 3f (x?)]乜11 h 2(3) o f (x)dx 2h[ f (0) f (h)] h2[ f (0) f (h)]答案：(1)具有三次代数精确度(2)具有二次代数精确度(3)具有三次代数精确度. a h xi xo ,确定求积公式X12 31 (x xo) f (x) dx h2EAf (xo) Bf (xi) ] h3[Cf (xo) Df (xi) ] R[f]X中的待定参数A, B, C, D ,使得该求积公式的代数精确度尽量高，并给出余项表达式.2/103)取7个节点处的函数值.用变步长的复化梯形公式和变步长的复化辛浦生公式计算 】山心砥•要求积分13 1610 3和10 6・版权文档，请勿用做商业用途 22 Ts 0. 946满足精度要求；使用复化辛浦生公式时,2 0J 田上述i 公武推导帶修忑项韵営化梯形求积公式K2 其中余域(x)dx= [占(xd 予 CxoH ， &b).为 T N h [f po) 2f (xi) 2f (X2) 2f (XN 1) f (XN )],Xi xo in, (i 0, 1, 2, , N), Nh XN XO •$ x 9、用龙贝格方法计算椭圆 / y 2 1的周长，使结果具有五位有效数字. o 4 答案：1 41 9. 6884 .验证高斯型求积公 e f (x) dx Ao f (xo) Ai f (xi)的系数及节点分别为式f<4)()h 6,其中答案：A 3 , B 7 , C 30 20 1440 P2(x)是以 0, h, 口2h •为插值上的二次插值多项式，用3h0 f ( x)dx 的数值积分公式Ih,并用台劳展开法证明：P2 (x)导岀计算积分h 4 f (0) 0(h 5) • 8Ih 0 P2(X )dx°4给定积分Ih[ f(0) 3f (2h)]'sin x dx(2) (3)答5运用复化梯形公式#算上述积分值，使其截断误差不 聲萝改用复化辛浦生公式计算时，截断误差是多少？亠 10 “ •2取同样的求积节 要求的截断误差不超过106,若用复化辛浦生公式，应取多少个节点处的函数值？ (1)只需n 7.5，取 9个节点，I 0. 946 ba 4 ⑷"41 6h 1 f ⑷()2) |Rn[f]| |2880 2880 4 5(V 0. 271 10 6 用事后误差估计法时，截断误不超过答案：使用复化梯形公式时，I S4 0. 946 083满足精度要求. f (1) (x) dx插值公式推导带有导数值的求积公式(b i2a )[f (b) f (a)] R[f],其中 确定高斯型求积公式0 xf (x) dx Aof (xo) Aif (xi)闻 xo , xi 及系数Ao,Ai.答案：xo 0. 289 949xi 0. 821 162 , Ao 0. 277 556, Ai 0. 389 111. 利用埃尔米特 b%ba[f(R f 山)]Ao 2： 2S Ao 2： 21x 0 2 2, Xi 2 2 . 第五章解线性方程组的直接法1 11用按列选主元的高斯若当消去法求矩 A 的逆矩阵’其中A21 01 1 0答案：用追赶法求解三对角方程组21 X11 131X22111X3221x4欣X4 2, X3L X2 1， XI 0 .第六章解线性代数方程组的迭代法X! 8X2 7X! 9X2 8作简单调整，使得用高斯一赛得尔迭代法求解时对任9x1 X2 X3 7 意初始向量都收敛，并取初始向量X （O ） [0 0 0]T使（k 1）k （）3||x （k bx k （）|| 10.3版权文档，请勿用做商业用途答案：近似解为X” [1.0000 1. 0000 1. 0000] T . 6 . 2讨论松弛因子1. 25时，用方法求解方程组1020X150101x231243x3170103x4答案: xi 2、X3 2X 21,Xi 1.411XI6 1 4. 25 2. 75X20. 512. 753. 5 X31. 25 答 xi 2X2X3用平方根法（分解法）求解方程组3用矩阵的直接三角分解法解方程组4x1 3x2 16 3xi 4x2 X3 20X2 4x312〔121，证明用雅可比迭代法解此方程组发散，而高斯-赛得尔迭代法收敛・12 1 123 0 2 X1bi6・4 设有方程组0 21X2b 2讨论用雅可比方法和咼斯一赛得尔方21 2 X3b3法解此方程组的收敛性•如果收敛，比较哪种方法收敛较 版权文档，请勿用做商业用途为6 . 3给定线性方程组Ax b,其中答案：雅可比方法收敛，高斯一赛得尔方法收敛，且较快.6. 5设矩阵A 非奇异.求证：方程组Ax b 的解总能通过高斯一赛得尔方法得到. …Aaij n n 为对称正定矩阵，对角阵D diag (an, a22 , , ann)・求证：高斯u 一赛得尔方法求解方程组D 2 AD 2x b 时对任意初始向量都收 敛.第七章非线性方程求根例7. 4对方程3x 2 e s 0确定迭代函数(x)及区间［a, b ］,使对xo ［a, b ］,迭代过程 XR i (x), k 0, 1, 2,均收敛，并求解.要求 xk 1 xk | 10x X? 0.458960903 •在［3, 4］上，将原方程改写为e x 3 x 2 ,取对数得性条件，则迭代序列xki In(3 xk 2 ), k 0, 1,2,在［3, 4］中有惟一解.取x 0 3.5 , x xie 3.733067511 •例7 . 6对于迭代函数(x) x c(x 2 3)，试讨论：的收敛性・若收敛，则取 x (0)[0 0 0]T迭代求解，使 ||x (I )x (k)1104-X1 1.50001，X2 答案：方程组的近似解3.33333，X32.16667 •答案：若取(X )e 2 ,则在［1,0］中满足收敛性条件,因此迭代法e 2k , k 0,1,2,在（1,0）中有惟一解•取 X0 0. 5， 3取(X )9 e"i,在［0 ,上1满足收敛性条件， 迭代序列1Xk 1 k 1 03k 0, 1, 2,在[0,1]中有惟一解.取 xo 0. 5,X X140.910001967x 2 ) (x)・满足收敛x In （3(1)当c为何值时，x kl (x k)产生的序列｛x k｝收敛于3；(2)c取何值时收敛最快？顿法收敛，证明牛顿迭代序列｛Xk ｝有下列极限矢系：l k im xk i 2xk xk i第八章矩阵特征值用乘幕法求矩阵A 的按模最大的特征值与对应的特征向量，已矢口 5 5 0 A 0 5. 5 1，要求 x (k)| 10 6,这里 严表示|的第k 次近似值.3 1答案：1 5 ,对应的特征向量为[5,0,0] T :2 5 ,对应的特征向量为[5, 10, T 5 ・]1 1 0＞彳 2的按模最小的特征12例7设不动点迭代xki (x)的迭代函数(x)具有二阶连续导数，/是(x)的不动1 1 5取C,力別If 鼻(X 丿旳个动点3 '妥吞| XkiXkl 1U- •3) 223(1 ) c (,0)时矗代收敛•答案： 31c 时收敛最快• O 、 233)分别取c 1,123，并取xo1.5，计算结果如下表7• 7所示yk点，且(X*) 1,证明迭代式(xk ) , Zk (xk )(yk x k )2 , k 0, 1, 2,二阶收敛于x"・版Xk 1 Xk Zk 2yk Xk权文档，请勿用做商业用途 例设(x) x p(x) f (x) q(x)f 2)，试确定函数p(x)和q(x),使求解f (x) 0且以(x)为迭代函数的迭代法至少三阶收 敛.案:p(x) f X (x )・ q(x) ；[f f (W]3例7设f (x)在[a, b]上有高阶导数,x* (a, b)是 f(x) 0的m(m 2)重根，且牛知A 的按模较大的特征 值用反幕法求矩阵A的近似值为15,用p 5的原点平移法计算1及其对应的特征向量.版权文档，请勿用做商业用途 答案：0 A 的按模最小的特征值为3 0. 238442812212第九章 微分方程初值问题的数值解法用反复迭代(反复校正)的欧拉预估一校正法求解初值问题y © 0] 0<x 0.2 5 ,要求取步长h 0. 1,每步迭代误差不超过10 5 .答案:Y y(0. 1) yi y 】⑷ 0. 904 762 , y(0. 2) y 2 y?⑷ 0.818 594267 一x y , 0<x 0. 4用二阶中点格式和二阶休恩格式求初值问题"“ “嗜厲汀⑹1长h 0.2,运算过程中保留五位小数). 计算得用平面旋转变换和反射变换将向量X [23 0 5] T 变为与 ei [1 0 0 0]T 平行的向量.2/ 38 3/ 385/ 38答案： T3/ 13 2/ 13 0 00 1 010/ 49415/ 4940 13/4940. 324 442 840 0. 486 664 262 0 0. 811107 1040. 486 664 2620.812 176 0480 0.298 039 922H10.811 107 104 0. 298 039 922 00.530 266 798然后用QR 方法求A 的全部特征值.4 4 5答案：取5 2. 234375即有2位有效数字. 532若A 6 4 4 ,试把A 化为相似的上阵， 值， 21n 0 时，Ki 1.000 00, K2 1. 200 00, y(0. 2) yi=l. 240 00n 1 时，Ki 1. 737 60, 用二阶休恩格式, K 2 2. 298 72, 取初值yo 1计算得y(0. 4) y 2 =1. 699 740 1 5. 1248854 ,对应的特征向量为(8) _设方阵A 的特征值都是实数，且满足 n)时， [0.242 4310, 1 , 0. 320 011 7],为求1而作原1 2 n,点平移'试证：当平移量P 2,(2幕法收敛最快•用二分法求三对角对〈方 A的最小特征 使它至少具有2位有212 答案：用二阶中点格式，取初值yo 1n 0 时，Ki 1.000 00, Ka 1.266 67, y(0.2) yi=1.240 00n 1时，Ki 1.737 60, Ka 2.499 18, y(0.4) y 2 =1.701 76用如下四步四阶阿达姆斯显格式 y n 1 y n h(55f n 59 fn 137fn2 9fn 3)/24求初值问题y x y, y(0) 1在［0,0.5］上的数值解•取步长h 0.1 小数点后保留8位•答 y(0.4) y 40.583 640 216 ‘ y(0.5) y 51.797 421 984 ・ 为使二阶中点公式ym yn hf(Xn h 2h,yn h f(Xn, yn))，求解初值问题2 n nh 的大小应受到的限制条件・hf (Xn,yn)用如下反复迭代的欧拉预估T&榴式 yn (k 11) yn h[f(Xn,y n ) f(Xn1,y n (k)1)］'k 0,1,2,; n 0,1,2,求解初值问题心讪•小时，如何选择步长h ，使上述格式矢于k 的迭y(0) 1代收敛•2答案：h 时上述格式尖于K 的迭代是收敛的・e求系数a,b,c,d ，使求解初值问题y f (x, y), y(xo) a 的如下隐式二步法 yn2aynh(bfn2Cfmdfn)的误差阶尽可能高，并指出其阶数•高'为五阶。

数值计算方法第一章答案
2.（1）3580
绝对误差限：；相对误差限：
经过四舍五入得到的近似值2580，其各位都是有效数字，故有四位有效数字。

（2）0.00476
绝对误差限：；相对误差限：
有三位有效数字。

（3）
绝对误差限：；相对误差限：
精确到小数点后两位，所以有四位有效数字。

或者
近似值有四位有效数字。

（4）也可理解为)
(对于14300000精确到小数点前四位)
有四位有效数字。

3.解：取为最好用
，
准确到小数点后第二位，有三位有效数字
准确到小数点后位，有单位有效数字。

6.解：
方法一：
，
取
（）
由此，取4位有效数字.
方法二（有问题，看看错在什么地方？）：
即
则取三位有效数字，
验证
满足要求，所以应取三位有效数字。

8.解：
又
或
11.解：由求根公式
解得
或。

习 题 一3.已知函数y =4, 6.25,9x x x ===处的函数值，试通过一个二次插值函解：0120124, 6.25,9;2, 2.5,3y x x x y y y =======由题意 （1） 采用Lagrange插值多项式220()()j j j y L x l x y ==≈=∑27020112012010*********()|()()()()()()()()()()()()(7 6.25)(79)(74)(79)(74)(7 6.25)2 2.532.255 2.25 2.75 2.7552.6484848x y L x x x x x x x x x x x x x y y y x x x x x x x x x x x x ==≈------=++------------=⨯+⨯+⨯⨯-⨯⨯= 其误差为(3)25(3)25(3)2[4,9]2()(7)(74)(7 6.25)(79)3!3()83max |()|40.0117281|(7)|(4.5)(0.01172)0.008796f R f x x f x R ξ--=---==<∴<=又则（2）采用Newton插值多项式2()y N x =≈ 根据题意作差商表：224(7)2(74)()(74)(7 6.25) 2.64848489495N =+⨯-+-⨯-⨯-≈4. 设()()0,1,...,k f x x k n ==，试列出()f x 关于互异节点()0,1,...,i x i n =的Lagrange 插值多项式。

注意到：若1n +个节点()0,1,...,i x i n =互异，则对任意次数n ≤的多项式()f x ，它关于节点()0,1,...,i x i n =满足条件(),0,1,...,i i P x y i n ==的插值多项式()P x 就是它本身。

可见，当k n ≤时幂函数()(0,1,...,)kf x x k n ==关于1n +个节点()0,1,...,i x i n =的插值多项式就是它本身，故依Lagrange 公式有()00(),0,1,...,nn n k kk i j j j j j i j ii jx x x l x x x k n x x ===≠-=≡=-∑∑∏特别地，当0k =时，有()0001nn n ij j j i j ii jx x l x x x ===≠-=≡-∑∑∏而当1k =时有()000nnn ij j j j j i j ii jx x x l x x x x x ===≠⎛⎫- ⎪=≡ ⎪- ⎪⎝⎭∑∑∏ 5.依据下列函数表分别建立次数不超过3的Lagrange 插值多项式和Newton 插值多项式，并验证插值多项式的唯一性。

计算方法引论课后答案第一章误差1.什么是模型误差，什么是方法误差？例如，将地球近似看为一个标准球体，利用公式 $A=4\pi r$ 计算其表面积，这个近似看为球体的过程产生的误差即为模型误差。

在计算过程中，要用到 $\pi$，我们利用无穷乘积公式计算 $\pi$ 的值：pi=2\cdot\frac{2}{1}\cdot\frac{2}{3}\cdot\frac{4}{3}\cdot\f rac{4}{5}\cdot\frac{6}{5}\cdot\frac{6}{7}\cdot\frac{8}{7}\cdot\ frac{8}{9}\cdot\cdots我们取前9项的乘积作为 $\pi$ 的近似值，得$\pi\approx3.xxxxxxxx5$。

这个去掉 $\pi$ 的无穷乘积公式中第9项后的部分产生的误差就是方法误差，也称为截断误差。

2.按照四舍五入的原则，将下列各数舍成五位有效数字：816.956，76.000，.322，501.235，.182，130.015，236.23.解：816.96，76.000，.501.24，.130.02，236.23.3.下列各数是按照四舍五入原则得到的近似数，它们各有几位有效数字？81.897，0.008，136.320，050.180.解：五位，三位，六位，四位。

4.若 $1/4$ 用 0.25 表示，问有多少位有效数字？解：两位。

5.若 $a=1.1062$，$b=0.947$，是经过舍入后得到的近似值，问：$a+b$，$a\times b$ 各有几位有效数字？已知 $da<\frac{1}{2}\cdot10^{-4}$，$db<\frac{1}{2}\cdot10^{-3}$，又 $a+b=0.\times10$。

begin{aligned}d(a+b)&=da+db\leq da+db=\frac{1}{2}\cdot10^{-4}+\frac{1}{2}\cdot10^{-3}=0.55\times10^{-3}<\frac{1}{2}\cdot10^{-2}end{aligned}所以 $a+b$ 有三位有效数字；因为 $a\timesb=0.xxxxxxxx\times10$。

计算方法第二版课后练习题含答案前言本文将为大家提供计算方法第二版课后练习题的答案，旨在帮助读者更好地学习和掌握计算方法的知识。

本文全部内容均为作者整理，尽可能保证每一题的答案正确性。

读者可以借助本文的答案，检验自己的练习成果，加强对计算方法知识的理解和掌握程度。

同时，读者也应该注意切勿直接复制答案，本文的答案仅供参考，希望读者能够通过自己的思考和探索，获得更深层次的学习感悟。

第一章引论1.1 计算方法的基本概念和思想练习题 1写出计算方法的三要素，并分别简要解释。

答案计算方法的三要素为：模型、算法、误差分析。

•模型：计算方法所涉及的实际问题所对应的数学模型，是解决问题的基础；•算法：根据模型，构造相应的计算程序，即算法；•误差分析：计算结果与实际应用中所需的精度之间的差异，称为误差。

误差分析是对计算结果质量的保障。

1.2 算法的误差练习题 2写出二分法算法，并解释其误差。

答案算法：function binarySearch(a, target) {let low = 0;let high = a.length - 1;while (low <= high) {let midIndex = Math.floor((low + high) / 2);let midValue = a[midIndex];if (midValue === target) {return midIndex;} else if (midValue < target) {low = midIndex + 1;} else {high = midIndex - 1;}}return -1;}误差:二分法算法的误差上界为O(2−k)，其中k为迭代次数。

在二分法被成功应用时，k取决于与目标值x的距离，即 $k=\\log _{2}(\\frac{b-a}{\\epsilon})$，其中[a,b]是区间，$\\epsilon$ 是目标值的精度。

数值分析复习试题第一章绪论一.填空题1.为精确值的近似值；为一元函数的近似值；*xx ()**x f y =()x f y =1为二元函数的近似值，请写出下面的公式：：()**,*y x f y =()y x f y ,2=**e x x =-***r x xe x -=()()()*'1**y f x x εε≈⋅()()()()'***1**r r x f x y x f x εε≈⋅()()()()()**,**,*2**f x y f x y y x y x yεεε∂∂≈⋅+⋅∂∂()()()()()****,***,**222r f x y e x f x y e y y x y y y ε∂∂≈⋅+⋅∂∂2、计算方法实际计算时，对数据只能取有限位表示，这时所产生的误差叫舍入误差。

3、分别用2.718281，2.718282作数e 的近似值，则其有效数字分别有 6 位和7（三位有效数字）。

1.73≈-211.73 10 2-≤⨯4、设均具有3位有效数字，则的相对误差限为 0.0055 。

121.216, 3.654x x ==12x x 5、设均具有3位有效数字，则的误差限为 0.01 。

121.216, 3.654x x ==12x x +6、已知近似值是由真值经四舍五入得到,则相对误差限为0.0000204 .2.4560A x =T x 7、递推公式如果取作计算,则计算到时,误差为,⎧⎪⎨⎪⎩0n n-1y =y =10y -1,n =1,2,0 1.41y =≈10y ;这个计算公式数值稳定不稳定 不稳定 .8110 2⨯8、精确值，则近似值和分别有 3 位和14159265.3*=π141.3*1=π1415.3*2=π4 位有效数字。

9、若,则x 有 6 位有效数字，其绝对误差限为1/2*10-5 。

*2.71828x e x =≈=10、 设x*的相对误差为2％，求(x*)n 的相对误差0.02n11、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字；12、计算方法主要研究( 截断 )误差和( 舍入 )误差；13、为了使计算 的乘除法次数尽量地少，应将该表达式()()2334610111y x x x =++----改写为11,))64(3(10-=-++=x t t t t y ，为了减少舍入误差，应将表达式19992001-改写为199920012+。

《计算方法》习题答案第一章 数值计算中的误差1．什么是计算方法？（狭义解释）答：计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算，以便在计算机上编程上机，求出问题的数值解，并对算法的收敛性、稳定性和误差进行分析、计算。

2．一个实际问题利用计算机解决所采取的五个步骤是什么？答：一个实际问题当利用计算机来解决时，应采取以下五个步骤： 实际问题→建立数学模型→构造数值算法→编程上机→获得近似结果 4．利用秦九韶算法计算多项式4)(53-+-=x x x x P 在3-=x 处的值，并编程获得解。

解：400)(2345-+⋅+-⋅+=x x x x x x P ，从而 1 0 -1 0 1 -4 -3 -3 9 -24 72 -2191-38-2473-223所以，多项式4)(53-+-=x x x x P 在3-=x 处的值223)3(-=-P 。

5．叙述误差的种类及来源。

答：误差的种类及来源有如下四个方面：（1）模型误差：数学模型是对实际问题进行抽象，忽略一些次要因素简化得到的，它是原始问题的近似，即使数学模型能求出准确解，也与实际问题的真解不同，我们把数学模型与实际问题之间存在的误差称为模型误差。

（2）观测误差：在建模和具体运算过程中所用的一些原始数据往往都是通过观测、实验得来的，由于仪器的精密性，实验手段的局限性，周围环境的变化以及人们的工作态度和能力等因素，而使数据必然带有误差，这种误差称为观测误差。

（3）截断误差：理论上的精确值往往要求用无限次的运算才能得到，而实际运算时只能用有限次运算的结果来近似，这样引起的误差称为截断误差（或方法误差）。

（4）舍入误差：在数值计算过程中还会用到一些无穷小数，而计算机受机器字长的限制，它所能表示的数据只能是一定的有限数位，需要把数据按四舍五入成一定位数的近似的有理数来代替。

这样引起的误差称为舍入误差。

6．掌握绝对误差（限）和相对误差（限）的定义公式。

第一章 误差1 问3.142，3.141，722分别作为π的近似值各具有几位有效数字？分析 利用有效数字的概念可直接得出。

解 π=3.141 592 65…记x 1=3.142，x 2=3.141，x 3=722.由π- x 1=3.141 59…-3.142=-0.000 40…知3411110||1022x π--⨯<-≤⨯ 因而x 1具有4位有效数字。

由π- x 2=3.141 59…-3.141=-0.000 59…知2231021||1021--⨯≤-<⨯x π因而x 2具有3位有效数字。

由π-722=3.141 59 …-3.142 85…=-0.001 26…知231021|722|1021--⨯≤-<⨯π因而x 3具有3位有效数字。

2 已知近似数x*有两位有效数字，试求其相对误差限。

分析 本题显然应利用有效数字与相对误差的关系。

解 利用有效数字与相对误差的关系。

这里n=2,a 1是1到9之间的数字。

%5101101|*||)(|1211*=⨯≤⨯≤-=+-+-n rx x x ε3 已知近似数的相对误差限为0.3%，问x*至少有几位有效数字？分析 本题利用有效数字与相对误差的关系。

解 a 1是1到9间的数字。

1112*10110113%3.0)(--⨯≤⨯=<=x r ε 设x*具有n 位有效数字，令-n+1=-1，则n=2，从而x*至少具有2位有效数字。

4 计算sin1.2，问要取几位有效数字才能保证相对误差限不大于0.01%。

分析 本题应利用有效数字与相对误差的关系。

解 设取n 位有效数字，由sin1.2=0.93…，故a 1=9。

411*10%01.01021|*||*||)(-+-=≤⨯≤-=n r a x x x x ε解不等式411101021-+-≤⨯n a 知取n=4即可满足要求。

5 计算76017591-，视已知数为精确值，用4位浮点数计算。

计算⽅法作业集及答案第⼀章数值计算基本常识⼀.填空题1. ⽤四舍五⼊得到的近似数0.628，有_____位有效数字，其绝对误差限是____________。

2. ⽤四舍五⼊得到的近似数0.586，有_____位有效数字，其绝对误差限是____________。

3. ⽤四舍五⼊得到的近似数0.69，其绝对误差是__________，由此计算出的相对误差限是__________。

4. ⽤四舍五⼊得到的近似数0.7960，其绝对误差是__________，由此计算出的相对误差限是__________。

5. 设0.484是0.4900的近似值，那么0.484具有____位有效数字。

6. 设x*=0.231是真值x=0.229的近似值，则x*有_____位有效数字。

7. 设x*=0.23是真值x=0.229的近似值，则x*有_____位有效数字。

8. 设x=2.3149541…，取5位有效数字，则所得的近似值x*=_____。

9. 设x=2.3149541…，取4位有效数字，则所得的近似值x*=_____。

10. 若近似数0.1100有4位有效数字，由有效数字计算出的相对误差是____________。

11. 若近似数76.82有4位有效数字，由有效数字计算出的相对误差是____________。

12. 若近似数576.00有5位有效数字，由有效数字计算出的相对误差是____________。

13. ⽤3.15作为π的近似值有_____位有效数字。

14. ⽤3.14作为π的近似值有_____位有效数字。

15. ⽤3.1416作为π的近似值有_____位有效数字。

解答:1. 3、0.5*10-32. 3、0.5*10-33. 0.5*10-2、0.725%4. 0.5*10-4、0.00628%5. 16. 27. 28. 2.31509. 2.31510. 0.05%11. 0.007%12. 0.001%13. 214. 315. 5⼆.选择题1. 3.141580 是π的近似值,有( )位有效数字。

《计算方法教程(第二版)》习题答案第一章 习题答案1、浮点数系),,,(U L t F β共有 1)1()1(21++---L U t ββ 个数。

3、a ．4097b ．62211101110.0,211101000.0⨯⨯c ．6211111101.0⨯4、设实数R x ∈，则按β进制可表达为：,1,,,3,2,011)11221(+=<≤<≤⨯++++++±=t t j jd d lt t d t t d dd x βββββββ按四舍五入的原则，当它进入浮点数系),,,(U L t F β时，若β211<+t d ，则 l tt d dd x fl ββββ⨯++±=)221()(若 β211≥+t d ，则 l tt d d d x fl ββββ⨯+++±=)1221()(对第一种情况：t l lt l t t d x fl x -++=⨯≤⨯+=-βββββ21)21(1)()(11对第二种情况：t l ltl t t d x fl x -++=⨯≤⨯--=-ββββββ21)21(1)(11 就是说总有： tl x fl x -≤-β21)( 另一方面，浮点数要求 β<≤11d ， 故有l x ββ1≥，将此两者相除，便得t x x fl x -≤-121)(β 5、a ． 5960.1 b ． 5962.1 后一种准确6、最后一个计算式：00025509.0原因：避免相近数相减，避免大数相乘，减少运算次数7、a ．]!3)2(!2)2(2[2132 +++=x x x yb ．)21)(1(22x x x y ++=c ．)11(222-++=x x x yd ． +-+-=!2)2(!6)2(!4)2(!2)2(2642x x x x ye ．222qp p q y ++=8、01786.098.5521==x x9、 m )10(m f - 1 233406.0- 3 20757.0- 5 8.07 710计算宜采用：])!42151()!32141()!22131[()(2432+⨯-+⨯-+⨯--=x x x f第二章 习题答案1、a ．T x )2,1,3(= b ．T x )1,2,1,2(--= c ．无法解2、a ．与 b ．同上， c ．T T x )2188.1,3125.0,2188.1,5312.0()39,10,39,17(321---≈---=7、a ．⎪⎪⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛---14112111473123247212122123211231321213122 b ． ⎪⎪⎪⎪⎪⎭⎫⎝⎛--⎪⎪⎪⎪⎪⎭⎫ ⎝⎛-=⎪⎪⎪⎪⎪⎭⎫⎝⎛----333211212110211221213231532223522121⎪⎪⎪⎪⎪⎭⎫ ⎝⎛--⎪⎪⎪⎪⎪⎭⎫⎝⎛-=111211212130213221219、T x )3415.46,3659.85,1220.95,1220.95,3659.85,3415.46(1= T x )8293.26,3171.7,4390.2,4390.2,3171.7,8293.26(2= 10、T LDL 分解：)015.0,579.3,9.1,10(diag D =⎪⎪⎪⎪⎪⎭⎫⎝⎛=16030.07895.05.018947.07.019.01L Cholesky 分解⎪⎪⎪⎪⎪⎭⎫⎝⎛=1225.01408.10833.15811.18918.12333.12136.23784.18460.21623.3G 解：)1,1,2,2(--=x 12、16,12,1612111===∞A A A611,4083.1,61122212===∞A A A2)(940)()(12111===∞A Cond A Cond A Cond524)(748)()(22221===∞A C o n d A C o n d A C o n d⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎭⎫ ⎝⎛=--180.0000180.0000- 30.0000 180.0000- 192.0000 36.0000- 30.0000 36.0000- 9.0000,0.0139 0.1111- 0.0694- 0.1111- 0.0556 0.1111- 0.0694- 0.1111- 0.0139 1211A A1151.372,1666.0212211==--A A15、 1A ：对应 Seidel Gauss - 迭代收敛，Jacobi 迭代不收敛； 2A ：对应 Seidel Gauss - 迭代收敛，Jacobi 迭代不收敛； 3A ：对应 Seidel Gauss - 迭代收敛，Jacobi 迭代收敛；第三章 习题答案1、Lagrange 插值多项式：)80.466.5)(20.366.5)(70.266.5)(00.166.5()80.4)(20.3)(70.2)(00.1(7.51)66.580.4)(20.380.4)(70.280.4)(00.180.4()66.5)(20.3)(70.2)(00.1(3.38)66.520.3)(80.420.3)(70.220.3)(00.120.3()66.5)(80.4)(70.2)(00.1(0.22)66.570.2)(80.470.2)(20.370.2)(00.170.2()66.5)(80.4)(20.3)(00.1(8.17)66.500.1)(80.400.1)(20.300.1)(70.200.1()66.5)(80.4)(20.3)(70.2(2.14)(4--------⨯+--------⨯+--------⨯+--------⨯+--------⨯=x x x x x x x x x x x x x x x x x x x x x L Newton 插值多项式：)80.4)(20.3)(70.2)(00.1(21444779.0)20.3)(70.2)(00.1(527480131.0)70.2)(00.1(855614973.2)00.1(117647059.22.14)(4----+------+-+=x x x x x x x x x x x N2、设)(x y y =，其反函数是以y 为自变量的函数)(y x x =，对)(y x 作插值多项式：)1744.0)(1081.0)(4016.0)(7001.0(01253.0)1081.0)(4016.0)(7001.0(01531.0)4016.0)(7001.0(009640.0)7001.0(3350.01000.0)(----+---+--+--=y y y y y y y y y y y N 3376.0)0(=N 是0)(=x y 在]4.0,3.0[中的近似根。

《计算方法教程（第二版）》习题答案第一章习题答案1、浮点数系F（0丄L、U）共有2（0-l）0i（U-厶+1） + 1个数。

3、a.4097b.0.11101000 x 22 , 0.11101110 x 25 6c.0.11111101x264、设实数xeR,则按0进制可表达为：1"1 V00 <> d j < p , J = 2,3,…+ 1,…按四舍五入得原则，当它进入浮点数系F（PJ,LM）时，若心V丄0,则2/心）"（第+2+…2“P pZ P1cK （1 +1/（□"（卡+样+…丄厂）〃P P L P l对第一种情况妝一."（x）| = （滸 + …）X0**G）X0‘ =^0 一对第二种情况:卜_/心卜爭巴一…"V *（£）x0詁旷就就是说总有：心）&丄0一2另一方面，浮点数要求1M/V0,故有|A-|^1/7\将此两者相除，便得r5 a. 1.5960 b. 1.5962后一种准确6最后一个计算式:0.00025509原因:避免相近数相减，避免大数相乘,减少运算次数2\I~X (Jx ，+1 + J 牙2 _])(2x)2 (2x)4 (2x)6(2x)2"^! 41- ~6!2!_3 -0.20757 5 0.8 7107计算宜采用:去)+G -親)x+G - 土用+…]第二章习题答案1. a.x = (3,1, 2)7b.x = (2, — 1, 2, — 1 )zc.无法解2、 a.与b.同上，c.x = —(-17, 39, -10,-39)7 « (-0.5312,1.218&一0.3125,-1.2188)7(2 -2 -1、/ 、 1、2 -2 -15p -17、a.3-12 =% 12%=3 21J 23,、％ %1J 3%丿1 1J<12 1 -2〕1-2 1 -2、 25 3 -22 11 12 -2 -23 5 -2 2 13 -3、132 3 >、1 2 0 1；3 ,(1 、 (\2 1 -2]2 1 1 1 2 -2 2 31 -1\ 1 2 0 3； 1 19、=(46.3415 , 85.3659 , 95.1220 , 95.1220 , 85.3659 , 46.3415)b. y =2x 2(l + x)(l +8、 X| =55.98 9、m 1x 2 = 0.01786 /(10-H,) -0.233406x 2 =(26.8293, 7.3171, 2.4390,2.4390, 7.3171,26.8293/ 10、厶£）厶了分解:D = diag( 10,1.9, 3.579,0.015)12、阀“16, ||州厂 12, ||州8 = 16h||2 =1.4083, ||A|L=1%Cond x (A|) = Cond n (A 】)=4% Cond 2 (^) = 2 Cond { (A 2) = Cond^(A 2) = 748Cond 2(A 2) = 524第三章习题答案1、Lagrange 插值多项式:'0.0139 -0.1111 ・0.0694、( 9.0000 -36.0000 30.0000、v-0.1111 0.0556 -0.1111，^2 = -36.0000 192.0000 -180.0000,・ 0.0694 ・ 0.1111 0.0139>(30.0000 -180.0000 180.0000,A ；'= 372.1151 -眉— 0.1666…,0.91L =0.7 0.89471.0.5 0.7895 0.6030 Cholesky 分解、H1623 2.8460G =2.2136 1.2333 1.8918J.5811 1.0833 1.1408 0.1225丿15. A 】 :对应 Gauss — Seidel 迭代收敛,Jacobi 迭代不收敛;:对应 Gauss — Seidel 迭代收Jacobi 迭代不收敛;:对应 Gauss — Seidel 迭代收Jacobi 迭代收敛;1丿 解：2(2, — 2,1, —1)(x - 2.70)(x- 3・20)(x - 4.80)(x 一 5.66)(1.00 - 2.70)(l .00 - 3.20)(1.00 - 4.80)(l .00 - 5.66)(x 一 1.00)(% 一 3.20)(x 一 4.80)(x 一 5.66)(2.70 -1.00)(2.70 - 3.20)(2.70 - 4.80)(2.70 - 5.66)(x -1.00)(x- 2.70)(x 一 4.80)(x- 5.66) (3.20 -1.00)(3.20 - 2.70)(3.20 - 4.80)(3.20 - 5.66)… (x-l ・00)(x-2・70)(x-3・20)(x-5・66) + 3 & 3 x -------------------------------------------(4.80 一 1.00)(4.80 一 2.70)(4.80 一 3.20)(4.80 一 5.66) (x-1.00)(x 一 2.70)(x 一 3.20)(x- 4.80)+ 51.7 x ---------------------- ---------------------(5.66 一 1.00)(5.66 一 2.70)(5.66 一 3.20)(5.66 - 4.80)Newton 插值多项式：^4(x) = 14.2 + 2.117647059(% -1.00)+ 2.855614973(x- 1.00)(x 一 2.70)一 0.527480131(x-1.00)(x 一 2.70)(x- 3.20)+ 0.21444779(“ 一 1.00)(x- 2.70)(x - 3.20)(x 一 4.80)差商表:2、设y = y(x),其反函数就是以y 为自变量得函数x = x(y)^x(j)作插值多项式: N(y)= 0.1000-0・3350(y — 0.7001)+ 0.009640( y-0.700 l)(y - 0.4016)+ 0.0153 l(y - 0.700 l)(y - 0.4016)(y - 0.1081) + 0.01253(0.7001)( V - 0.4016)(y -0.108 l)(y - 0.1744)N(0) = 0.3376 就是 y(x) = 0在［0.3, 0.4 ］中得近似根。

计算⽅法引论课后答案.第⼀章误差1. 试举例,说明什么是模型误差,什么是⽅法误差.解: 例如,把地球近似看为⼀个标准球体,利⽤公式24A r π=计算其表⾯积,这个近似看为球体的过程产⽣的误差即为模型误差.在计算过程中,要⽤到π,我们利⽤⽆穷乘积公式计算π的值:12222...q q π=?其中112,3,...n q q n +?=??==?? 我们取前9项的乘积作为π的近似值,得3.141587725...π≈这个去掉π的⽆穷乘积公式中第9项后的部分产⽣的误差就是⽅法误差,也成为截断误差.2. 按照四舍五⼊的原则,将下列各数舍成五位有效数字:816.956 7 6.000 015 17.322 50 1.235 651 93.182 13 0.015 236 23 解: 816.96 6.000 0 17.323 1.235 7 93.182 0.015 236 3. 下列各数是按照四舍五⼊原则得到的近似数,它们各有⼏位有效数字? 81.897 0.008 13 6.320 05 0.180 0 解: 五位三位六位四位4. 若1/4⽤0.25表⽰,问有多少位有效数字? 解: 两位5. 若 1.1062,0.947a b ==,是经过舍⼊后得到的近似值,问:,a b a b +?各有⼏位有效数字?解: 已知4311d 10,d 1022a b --()433211110100.551010222d a b da db da db ----+=+≤+=?+?=?所以a b +有三位有效数字;因为0.1047571410a b ?=?,()43321110.94710 1.1062100.600451010222所以a b ?有三位有效数字.6. 设120.9863,0.0062y y ==,是经过舍⼊后作为12,x x 的近似值.求1211,y y 的计算值与真值的相对误差限及12y y ?与真值的相对误差限. 解: 已知-4-41112221211d ,d ,d =10,d 1022x y x x y x x x =+=+?=?, ()44111111110d d 12dr dr 0.50100.9863x xx x x y --==≈=≈? ???;()42222222110d d 12dr dr 0.81100.0062x xx x x y --==≈=≈? ???;()()()4221212dr dr dr 0.50100.81100.8210x x x x ---?=+≈?+?≈?.7. 正⽅形的边长约为100cm,应该怎样测量,才能使其⾯积的误差不超过1cm 2.解: 设正⽅形⾯积为S,边长为a,则S=a 2.所以要使:2d d 2d 1s a a a ==≤,则要求211d 0.5102200a a -≤==?.所以边长的误差不能超过20.510-?cm.8. ⽤观测恒星的⽅法求得某地维度为4502'''o(读到秒),试问:计算sin ?将有多⼤误差?解: ()()1d sin cos d cos 45022*''?'''==o.9 . 真空中⾃由落体运动距离s 与时间的关系由公式212s gt =确定,g 是重⼒加速度.现在假设g 是准确的,⽽对t 的测量有0.1s ±的误差,证明t 增加时,距离的绝对误差增加⽽相对误差却减⼩.证明: 因为:221d d d d d d d ;2.122s gt t gt t t s gt gt t s s t gt ??d s 与t 成正⽐,d s s与t 成反⽐,所以当d t 固定的时候, t 增加时,距离的绝对误差增加⽽相对误差却减⼩.10. 设0x >,x 的相对误差为δ,求ln x 的绝对误差. 解: 已知d x x δ=,所以ln x 的绝对误差()d d ln x x x δ==.11. 设x 的相对误差为%α,求nx 的相对误差.解: 1d d d %n n n n x nx x n xn x x xα-===.12. 计算球的体积,为了使相对误差限为1%,问度量半径R 时允许的相对误差限如何? 解: 已知34 3V R π=,设()d dr R R a R ==,则要使得 ()()3d dr dln d ln 3d ln 3d ln 3dr 31%V V V R R R R a V ========,则11%3a =?.第⼆章插值法与数值微分1.设y =在100,121,144x =三处的值是很容易求得的,试以这三个点建⽴y =的⼆次插值多项式,,且给出误差估计.⽤其中的任意两点,构造线性插值函数,⽤得到的三个线性插值函数,,并分析其结果不同的原因.解: 已知012012100,121,144;10,11,12x x x y y y ======,建⽴⼆次Lagrange 插值函数可得:()()()()21211441001441011100121100144121100121144121100 12144121144100x x x x L x x x ----= +------+--()211510.7228L ≈=.误差()()()()()()2012012,,,,3!f R x x x x x x x x x x ξξξ'''=---∈,所以20.00065550.001631R <<利⽤前两个节点建⽴线性插值函数可得:()()()()()11211001011100121121100x x L x --=+--()111510.7143L ≈=.利⽤后两个节点建⽴线性插值可得:()()()()()11441211112121144144121x x L x --=+--()111510.7391L ≈=.利⽤前后两个节点建⽴线性插值可得:()()()()()21441001012100144144100x x L x --=+()111510.6818L ≈=.,⼆次插值⽐线性插值效果好,利⽤前两个节点的线性插值⽐其他两个线性插值效果好.此说明,⼆次插值⽐线性插值效果好,插⽐外插效果好.2. 利⽤(2.9)式证明()()()0121001max ,8x x x x x R x f x x x x ≤≤-''≤≤≤证明: 由(2.9)式()()()()0101,2!f R x x x x x x x ξξ''=--<<当01x x x <<时,()()01max x x x f f x ξ≤≤''''≤,()()()01201101max 4x x x x x x x x x ≤≤--≤- 所以()()()0121001max ,8x x x x x R x f x x x x ≤≤-''≤≤≤3. 若()0,1,...,j x n 为互异节点,且有()()()()()()()()()011011............j j n j jj j j j j n x x x x x x x x l x xx x x x x x x -+-+----=证明()0,0,1,...,nk kj j j x l x xk n =≡=∑证明: 由于() 1 ;0 .j i ij i j l x i j δ=?==?≠? 且()0nk j j j x l x =∑和kx都为k 次多项式,⽽且在k+1个不同的节点处的函数值都相同0,1,...,k n =, 所以马上有()0,0,1,...,nk kj j j x l x xk n =≡=∑.4. 设给出sin x 在[],ππ-上的数值表,⽤⼆次插值进⾏计算,若希望截断误差⼩于5 10-,问函数表的步长最⼤能取多少? 解: 记插值函数为p(x),则()()()()()11sin sin 3!i i i x p x x x x x x x ξ-+'''-=--- 所以()()()()11cos max sin 3!i i i x x p x x x x x x ππξ-+-≤≤--=---()()()[]3112,0,2i g x th h t t t t -+=--∈⼜()()()[]12,0,2t t t t t ?=--∈的最⼤值为10.3849??= ?,所以有 350.3849max sin 106x x p h ππ--≤≤-≤< 所以 0.0538h ≤.5. ⽤拉格朗⽇插值和⽜顿插值找经过点()()()()3,1,0,2,3,2,6,10---的三次插值公式. 解: Lagrange 插值函数:()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()()12302330101020310121301301223202123303132 31033101622731033 .2781/5x x x x x x x x x x x x L x y y x x x x x x x x x x x x x x x x x x x x x x x x y y x x x x x x x x x x x x x x x x x x x x x x x x ------= +------------++--------+--=++-+-++⽜顿插值: ⾸先计算差商3 10 2 13 2 1.333 0.38896 104 0.8889 0.1420-----()()()()()3130.38893 1.142033.N x x x x x x x =-++-+++-也可以利⽤等距节点构造,⾸先计算差分。

第一章 绪论一.填空题1.*x 为精确值x 的近似值；()**x f y=为一元函数()x f y =1的近似值；()**,*y x f y =为二元函数()y x f y ,2=的近似值，请写出下面的公式：**e x x =-：***rx x e x -=()()()*'1**y f x x εε≈⋅ ()()()()'***1**r r x f x y x f x εε≈⋅ ()()()()()**,**,*2**f x y f x y y x y x yεεε∂∂≈⋅+⋅∂∂()()()()()****,***,**222r f x y e x f x y e y y x y y y ε∂∂≈⋅+⋅∂∂2、 计算方法实际计算时，对数据只能取有限位表示，这时所产生的误差叫 舍入误差 。

3、 分别用2.718281，2.718282作数e的近似值，则其有效数字分别有 6 位和 7 位；又取1.73≈（三位有效数字），则-211.73 10 2≤⨯。

4、设121.216, 3.654x x ==均具有3位有效数字，则12x x 的相对误差限为 0.0055 。

5、设121.216, 3.654x x ==均具有3位有效数字，则12x x +的误差限为 0.01 。

6、已知近似值 2.4560A x =是由真值T x 经四舍五入得到,则相对误差限为 0.000021 .7、递推公式,⎧⎪⎨⎪⎩0n n-1y =y =10y -1,n =1,2,如果取0 1.41y ≈作计算,则计算到10y 时,误差为8110 2⨯;这个计算公式数值稳定不稳定 不稳定 .8、精确值 14159265.3*=π，则近似值141.3*1=π和1415.3*2=π分别有 3 位和 4 位有效数字。

9、 若*2.71828x e x =≈=,则x 有 6 位有效数字，其绝对误差限为1/2*10-5。

10、 设x*的相对误差为2％，求(x*)n 的相对误差0.02n 二、计算题1. 有一个长方形水池,由测量知长为(50±0.01)米,宽为(25±0.01)米,深为(20±0.01)米,试按所给数据求出该水池的容积,并分析所得近似值的绝对误差和相对误差公式,并求出绝对误差限和相对误差限. 解：设长方形水池的长为L ，宽为W,深为H ，则该水池的面积为V=LWH当L=50,W=25,H=20时，有 V=50*25*20=25000(米3) 此时，该近似值的绝对误差可估计为()()()()()()()=V V VV L W H L W HWH L HL W LW H ∂∂∂∆≈∆+∆+∆∂∂∂∆+∆+∆ 相对误差可估计为：()()r V V V∆∆=而已知该水池的长、宽和高的数据的绝对误差满足()()()0.01,0.01,0.01L W H ∆≤∆≤∆≤故求得该水池容积的绝对误差限和相对误差限分别为()()()()()()325*20*0.0150*20*0.0150*25*0.0127.5027.501.1*1025000r V WH L HL W LW H V V V -∆≤∆+∆+∆≤++=∆∆=≤=2.已知测量某长方形场地的长a=110米,宽b=80米.若()()**0.1 0.1a a b b -≤-≤米，米试求其面积的绝对误差限和相对误差限. 解：设长方形的面积为s=ab当a=110,b=80时，有 s==110*80=8800(米2) 此时，该近似值的绝对误差可估计为()()()()()=b s ss a b a ba ab ∂∂∆≈∆+∆∂∂∆+∆ 相对误差可估计为：()()r s s s∆∆=而已知长方形长、宽的数据的绝对误差满足()()0.1,0.1a b ∆≤∆≤故求得该长方形的绝对误差限和相对误差限分别为()()()()() 80*0.1110*0.119.019.00.0021598800r s b a a b s s s ∆≤∆+∆≤+=∆∆=≤= 绝对误差限为19.0；相对误差限为0.002159。

第一章 误差1 问，，722分别作为π的近似值各具有几位有效数字分析 利用有效数字的概念可直接得出。

解 π= 592 65… 记x 1=，x 2=，x 3=722.由π- x 1= 59…= 40…知3411110||1022x π--⨯<-≤⨯ 因而x 1具有4位有效数字。

由π- x 2= 59…= 59…知2231021||1021--⨯≤-<⨯x π因而x 2具有3位有效数字。

由π-722= 59 … 85…= 26…知231021|722|1021--⨯≤-<⨯π因而x 3具有3位有效数字。

2 已知近似数x*有两位有效数字，试求其相对误差限。

分析 本题显然应利用有效数字与相对误差的关系。

解 利用有效数字与相对误差的关系。

这里n=2,a 1是1到9之间的数字。

%5101211021|*||*||)(|1211*=⨯⨯≤⨯≤-=+-+-n ra x x x x ε3 已知近似数的相对误差限为%，问x*至少有几位有效数字 分析 本题利用有效数字与相对误差的关系。

解 a 1是1到9间的数字。

1112*10)1(2110)19(21102110003%3.0)(--⨯+≤⨯+⨯=⨯<=a x r ε 设x*具有n 位有效数字，令-n+1=-1，则n=2，从而x*至少具有2位有效数字。

4 计算，问要取几位有效数字才能保证相对误差限不大于%。

分析 本题应利用有效数字与相对误差的关系。

解 设取n 位有效数字，由=…，故a 1=9。

411*10%01.01021|*||*||)(-+-=≤⨯≤-=n r a x x x x ε解不等式411101021-+-≤⨯n a知取n=4即可满足要求。

5 计算76017591-，视已知数为精确值，用4位浮点数计算。

解 =-76017591 8×10-2-0.131 6×10-2=×10-5结果只有一位有效数字，有效数字大量损失，造成相对误差的扩大，若通分后再计算：56101734.0105768.01760759176017591-⨯=⨯=⨯=- 就得到4位有效数字的结果。