下载文档原格式
第一章数值计算中的误差习题一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。
计算⽅法各习题及参考答案第⼆章数值分析2.1 已知多项式432()1p x x x x x =-+-+通过下列点：试构造⼀多项式()q x 通过下列点：答案：54313()()()3122q x p x r x x x x x =-=-++-+． 2.2 观测得到⼆次多项式2()p x 的值：表中2()p x 的某⼀个函数值有错误，试找出并校正它．答案：函数值表中2(1)p -错误，应有2(1)0p -=．2.3 利⽤差分的性质证明22212(1)(21)/6n n n n +++=++ ．2.4 当⽤等距节点的分段⼆次插值多项式在区间[1,1]-近似函数xe 时，使⽤多少个节点能够保证误差不超过61102-?．答案：需要143个插值节点．2.5 设被插值函数4()[,]f x C a b ∈，()3()h H x 是()f x 关于等距节点01n a x x x b =<<<= 的分段三次艾尔⽶特插值多项式，步长b a h n-=．试估计()3||()()||h f x H x ∞-．答案：()443||()()||384h M f x H x h ∞-≤．第三章函数逼近3.1 求()sin ,[0,0.1]f x x x =∈在空间2{1,,}span x x Φ=上最佳平⽅逼近多项式，并给出平⽅误差．答案：()sin f x x =的⼆次最佳平⽅逼近多项式为-522sin ()0.832 440 710 1.000 999 10.024 985 1x p x x x ≈=-?+-，⼆次最佳平⽅逼近的平⽅误差为0.122-1220(sin )())0.989 310 710x p x dx δ=-=??．3.2 确定参数,a b c 和，使得积分2121(,,)[I a b c ax bx c -=++-?取最⼩值．答案：810, 0, 33a b c ππ=-== 3.3 求多项式432()251f x x x x =+++在[1,1]-上的３次最佳⼀致逼近多项式()p x ．答案：()f x 的最佳⼀致逼近多项式为323()74p x x x =++． 3.4 ⽤幂级数缩合⽅法，求() (11)x f x e x =-≤≤上的３次近似多项式6,3()p x ，并估计6,3||()()||f x p x ∞-．答案：236,3()0.994 574 650.997 395 830.542 968 750.177 083 33p x x x x =+++， 6,3||()()||0.006 572 327 7f x p x ∞-≤3.5 求() (11)xf x e x =-≤≤上的关于权函数()x ρ=的三次最佳平⽅逼近多项式3()S x ，并估计误差32||()()||f x S x -和3||()()||f x S x ∞-．答案：233()0.994 5710.997 3080.542 9910.177 347S x x x x =+++，32||()()||0.006 894 83f x S x -=，3||()()||0.006 442 575f x S x ∞-≤．第四章数值积分与数值微分4.1 ⽤梯形公式、⾟浦⽣公式和柯特斯公式分别计算积分1(1,2,3,4)n x dx n =?，并与精确值⽐较．答案：计算结果如下表所⽰4.2 确定下列求积公式中的待定参数，使得求积公式的代数精度尽量⾼，并指明所确定的求积公式具有的代数精度．（１）101()()(0)()hh f x dx A f h A f A f h --≈-++?（２）11211()[(1)2()3()]3f x dx f f x f x -≈-++? （３）20()[(0)()][(0)()]2h h f x dx f f h h f f h α''≈++-?答案：（１）具有三次代数精确度（２）具有⼆次代数精确度（３）具有三次代数精确度．4.3 设10h x x =-，确定求积公式12300101()()[()()][()()][]x x x x f x dx h Af x Bf x h Cf x Df x R f ''-=++++?中的待定参数,,,A B C D ，使得该求积公式的代数精确度尽量⾼，并给出余项表达式．答案：3711,,,20203020A B C D ====-，(4)6()[]1440f R f h η=，其中01(,)x x η∈．4.4 设2()P x 是以0,,2h h 为插值点的()f x 的⼆次插值多项式，⽤2()P x 导出计算积分30()hI f x dx =?的数值积分公式h I ，并⽤台劳展开法证明：453(0)()8h I I h f O h '''-=+．答案：3203()[(0)3(2)]4h h I p x dx h f f h ==+?．4.5 给定积分10sin xI dx x =（１）运⽤复化梯形公式计算上述积分值，使其截断误差不超过31102-?．（２）取同样的求积节点，改⽤复化⾟浦⽣公式计算时，截断误差是多少？（３）要求的截断误差不超过610-，若⽤复化⾟浦⽣公式，应取多少个节点处的函数值？答案：（１）只需7.5n ≥，取９个节点，0.946I ≈（２）4(4)46111|[]||()|()0.271102880288045n b a R f h f η--=-≤=? （３）取７个节点处的函数值．4.6 ⽤变步长的复化梯形公式和变步长的复化⾟浦⽣公式计算积分10sin xI dx x =?．要求⽤事后误差估计法时，截断误不超过31102-?和61102-?．答案：使⽤复化梯形公式时，80.946I T ≈=满⾜精度要求；使⽤复化⾟浦⽣公式时，40.946 083I s ≈=满⾜精度要求．4.7（１）利⽤埃尔⽶特插值公式推导带有导数值的求积公式2()()[()()][()()][]212ba b a b a f x dx f a f b f b f a R f --''=+--+?，其中余项为 5(4)()[](), (,)4!30b a R f f a b ηη-=∈．（２）利⽤上述公式推导带修正项的复化梯形求积公式020()[()()]12Nx N N x h f x dx T f x f x ''≈--?，其中 0121[()2()2()2()()]2N N N hT f x f x f x f x f x -=+++++ ，⽽ 00, (0,1,2,,), i N x x ih i N Nh x x =+==- ．4.8 ⽤龙贝格⽅法计算椭圆2214x y +=的周长，使结果具有五位有效数字．答案：49.6884l I =≈．4.9确定⾼斯型求积公式0011()()()x dx A f x A f x ≈+?的节点0x ，1x 及系数0A ，1A ．答案：00.289 949x =，10.821 162x =，00.277 556A =，10.389 111A =．4.10 验证⾼斯型求积公式00110()()()x e f x dx A f x A f x +∞-≈+?的系数及节点分别为0001 2 2A A x x ===-=+第五章解线性⽅程组的直接法5.1 ⽤按列选主元的⾼斯-若当消去法求矩阵A 的逆矩阵，其中11121 0110A -?? ?= ? ?-??．答案： 1110331203321133A -?? ? ?=---5.2 ⽤矩阵的直接三⾓分解法解⽅程组1234102050101312431701037x x x x= ? ? ? ? ? ? ? ? ??答案： 42x =，32x =，21x =，11x =．5.3 ⽤平⽅根法(Cholesky 分解法)求解⽅程组12341161 4.25 2.750.51 2.75 3.5 1.25x x x -?????? ??? ?-=- ??? ? ??? ???????答案： 12x =，21x =，31x =-．5.4 ⽤追赶法求解三对⾓⽅程组123421113121112210x x x x ?????? ? ? ? ? ? ?= ? ? ? ? ? ? ? ? ?????答案：42x =，31x =-，21x =，10x =．第六章解线性代数⽅程组的迭代法６.１对⽅程1212123879897x x x x x x x -+=??-+=??--=?作简单调整，使得⽤⾼斯－赛得尔迭代法求解时对任意初始向量都收敛，并取初始向量(0)[0 0 0]T x =，⽤该⽅法求近似解(1)k x+，使(1)()3||||10k k x x +-∞-≤．答案：近似解为(4)[1.0000 1.0000 1.0000]Tx =．６.２讨论松弛因⼦ 1.25ω=时，⽤SOR ⽅法求解⽅程组121232343163420412x x x x x x x +=??+-=??-+=-? 的收敛性．若收敛，则取(0)[0 0 0]T x=迭代求解，使(1)()41||||102k k x x +-∞-<．答案：⽅程组的近似解为*1 1.50001x =，*2 3.33333x =，*3 2.16667x =-．６.３给定线性⽅程组Ax b =，其中111221112211122A ?? ? ?=，证明⽤雅可⽐迭代法解此⽅程组发散，⽽⾼斯－赛得尔迭代法收敛．６.４设有⽅程组112233302021212x b x b x b -?????? ??? ?= ??? ? ??? ?-??????，讨论⽤雅可⽐⽅法和⾼斯－赛得尔⽅法解此⽅程组的收敛性．如果收敛，⽐较哪种⽅法收敛较快．答案：雅可⽐⽅法收敛，⾼斯－赛得尔⽅法收敛,且较快．６.５设矩阵A ⾮奇异．求证：⽅程组Ax b =的解总能通过⾼斯－赛得尔⽅法得到．６.６设()ij n nA a ?=为对称正定矩阵，对⾓阵1122(,,,)nn D diag a a a = ．求证：⾼斯－赛得尔⽅法求解⽅程组1122D AD x b --=时对任意初始向量都收敛．第七章⾮线性⽅程求根例７.４对⽅程230xx e -=确定迭代函数()x ?及区间[,]a b ，使对0[,]x a b ?∈，迭代过程1(), 0,1,2,k x x k ?+== 均收敛，并求解．要求51||10k k x x -+-<．答案：若取2()x x ?=，则在[1,0]-中满⾜收敛性条件，因此迭代法121, 0,1,2,k x k x k +== 在(1,0)-中有惟⼀解．取00.5x =-，*70.458960903x x ≈=-．取2()x x ?=，在[0,1上满⾜收敛性条件，迭代序列121, 0,1,2,k x k x k +== 在[0,1]中有惟⼀解．取00.5x =，*140.910001967x x ≈=- 在[3,4]上，将原⽅程改写为23xe x =，取对数得2ln(3)()x x x ?==．满⾜收敛性条件，则迭代序列21ln(3), 0,1,2,k k x x k +== 在[3,4]中有惟⼀解．取0 3.5x =， *16 3.733067511x x ≈=．例７.６对于迭代函数2()(3)x x c x ?=+-，试讨论：（１）当c 为何值时，1()k k x x ?+=产⽣的序列{}k x（２）c 取何值时收敛最快？（３）取1,2c =-()x ?51||10k k x x -+-<．答案：（１）(c ∈时迭代收敛．（２）c =时收敛最快．（３）分别取1, 2c =--，并取0 1.5x =，计算结果如下表７.７所⽰表７.７例７.13 设不动点迭代1()k x x ?+=的迭代函数()x ?具有⼆阶连续导数，*x 是()x ?的不动点，且*()1x ?'≠，证明Steffensen 迭代式21(), (), 0,1,2,()2k k k k k k k k k k k y x z x k y x x x z y x+===-?=-?-+?⼆阶收敛于*x ．例7.15 设2()()()()()x x p x f x q x f x ?=--，试确定函数()p x 和()q x ，使求解()0f x =且以()x ?为迭代函数的迭代法⾄少三阶收敛．答案：1()()p x f x ='，31()()2[()]f x q x f x ''=' 例７.19 设()f x 在[,]a b 上有⾼阶导数，*(,)x a b ∈是()0f x =的(2)m m ≥重根，且⽜顿法收敛，证明⽜顿迭代序列{}k x 有下列极限关系：111lim2k kk k k k x x m x x x -→∞-+-=-+．第⼋章矩阵特征值8.1 ⽤乘幂法求矩阵A 的按模最⼤的特征值与对应的特征向量，已知5500 5.51031A -?? ?=- ? ?-??，要求(1)()611||10k k λλ+--<，这⾥()1k λ表⽰1λ的第k 次近似值．答案：15λ≈，对应的特征向量为[5,0,0]T-；25λ≈-，对应的特征向量为[5,10,5]T --． 8.2 ⽤反幂法求矩阵110242012A -??=-- -的按模最⼩的特征值．知A 的按模较⼤的特征值的近似值为15λ=，⽤5p =的原点平移法计算1λ及其对应的特征向量．答案：(1) A 的按模最⼩的特征值为30.2384428λ≈(2) 1 5.1248854λ≈，对应的特征向量为(8)[0.242 4310, 1 ,0.320 011 7]T U =--．8.3 设⽅阵A 的特征值都是实数，且满⾜121, ||||n n λλλλλ>≥≥> ，为求1λ⽽作原点平移，试证：当平移量21()2n p λλ=+时，幂法收敛最快． 8.4 ⽤⼆分法求三对⾓对称⽅阵1221221221A ?? ? ?= ? ? ???的最⼩特征值，使它⾄少具有２位有效数字．答案：取5 2.234375λ≈-即有２位有效数字．8.5 ⽤平⾯旋转变换和反射变换将向量[2 3 0 5]T x =变为与1[1 0 0 0]Te =平⾏的向量．答案：203/2/00001010/0T ??- ?=--?0.324 442 8400.486 664 26200.811 107 1040.486 664 2620.812 176 04800.298 039 92200100.811 107 1040.298 039 92200.530 266 798H --??--= ? ?--8.6 若532644445A -??=- -，试把A 化为相似的上Hessenberg 阵，然后⽤QR ⽅法求A 的全部特征值．第九章微分⽅程初值问题的数值解法9.1 ⽤反复迭代（反复校正）的欧拉预估－校正法求解初值问题0, 0<0.2(0)1y y x y '+=≤??=?，要求取步长0.1h =，每步迭代误差不超过510-．答案： [4]11(0.1)0.904 762y y y ≈==，[4]22(0.2)0.818 594y y y ≈==9.2 ⽤⼆阶中点格式和⼆阶休恩格式求初值问题2, 0<0.4(0)1dy x y x dx y ?=+≤=?的数值解(取步长0.2h =，运算过程中保留五位⼩数)．答案：⽤⼆阶中点格式，取初值01y =计算得0n =时，1211.000 00, 1.200 00, (0.2)=1.240 00K K y y ==≈ 1n =时，1221.737 60, 2.298 72, (0.4)=1.699 74K K y y ==≈⽤⼆阶休恩格式，取初值01y =计算得0n =时，1211.000 00, 1.266 67, (0.2)=1.240 00K K y y ==≈ 1n =时，1221.737 60, 2.499 18, (0.4)=1.701 76K K y y ==≈9.3 ⽤如下四步四阶阿达姆斯显格式1123(5559379)/24n n n n n n y y h f f f f +---=+-+-求初值问题, (0)1y x y y '=+=在[0,0.5]上的数值解．取步长0.1h =，⼩数点后保留８位．答案：4(0.4)0.583 640 216y y ≈=，5(0.5) 1.797 421 984y y ≈=． 9.4 为使⼆阶中点公式1(,(,))22n n n n n n h hy y hf x y f x y +=+++，求解初值问题 , (0)y y y aλλ'=-??=?为实常数绝对稳定，试求步长h 的⼤⼩应受到的限制条件．答案：2h λ≤．9.5 ⽤如下反复迭代的欧拉预估－校正格式(0)1(1)()111(,)[(,)(,)]2 0,1,2,; 0,1,2,nn n n k k n n n n n n y y hf x y h y y f x y f x y k n +++++?=+??=++??==，求解初值问题sin(), 01(0)1x y e xy x y '?=<≤?=?时，如何选择步长h ，使上述格式关于k 的迭代收敛．答案：2h e<时上述格式关于k 的迭代是收敛的．9.6 求系数,,,a b c d ，使求解初值问题0(,), ()y f x y y x a '==的如下隐式⼆步法221()n n n n n y ay h bf cf df +++=+++的误差阶尽可能⾼，并指出其阶数．答案：系数为142,,33a b d c ====，此时⽅法的局部截断误差阶最⾼，为五阶5()O h ．9.7 试⽤欧拉预估－校正法求解初值问题, (0)=1, 0<0.2()/, (0)2dyxy z y dxx dz x y z z dx=-≤=+=，取步长0.1h =，⼩数点后⾄少保留六位．答案：由初值00(0)1, (0)2y y z z ====可计算得110.800 000z 2.050 000y =??=? ， 11(0.1)0.801 500(0.1) 2.046 951y y z z ≈=??≈=? 220.604 820z 2.090 992y =??=? ， 22 (0.2)0.604 659(0.2) 2.088 216y y z z ≈=??≈=?。
第一章 绪论（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 。
，。
第一章 绪论（12）1、设0>x ，x 的相对误差为δ，求x ln 的误差。
[解]设0*>x 为x 的近似值，则有相对误差为δε=)(*x r ，绝对误差为**)(x x δε=，从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x xx 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 问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、浮点数系),,,(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 问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位有效数字。
由π-=3.141 59 …-3.142 85…=-0.001 26…知722231021|722|1021--⨯≤-<⨯π因而x 3具有3位有效数字。
2 已知近似数x*有两位有效数字，试求其相对误差限。
分析 本题显然应利用有效数字与相对误差的关系。
解 利用有效数字与相对误差的关系。
这里n=2,a 1是1到9之间的数字。
%5101211021|*||*||)(|1211*=⨯⨯≤⨯≤-=+-+-n r a x x x x ε3 已知近似数的相对误差限为0.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 计算sin1.2，问要取几位有效数字才能保证相对误差限不大于0.01%。
分析 本题应利用有效数字与相对误差的关系。
解 设取n 位有效数字，由sin1.2=0.93…，故a 1=9。
411*10%01.01021|*||*||)(-+-=≤⨯≤-=n r a x x x x ε解不等式知取n=4即可满足要求。
411101021-+-≤⨯n a 5 计算，视已知数为精确值，用4位浮点数计算。
76017591-解0.131 8×10-2-0.131 6×10-2=0.2×10-5=-76017591结果只有一位有效数字，有效数字大量损失，造成相对误差的扩大，若通分后再计算：56101734.0105768.01760759176017591-⨯=⨯=⨯=-就得到4位有效数字的结果。
此例说明，在数值计算中，要特别注意两相近数作减法运算时，有效数字常会严重损失，遇到这种情况，一般采取两种办法：第一，应多留几位有效数字；第二，将算式恒等变形，然后再进行计算。
例如，当x 接近于0，计算时，应先把算式变形为xx sin cos 1-xx x x x x x cos 1sin )cos 1(sin cos 1sin cos 12+=+-=-再计算。
又例如，当x 充分大时，应作变换xx x x ++=-+111)1(1111+=+-x x x x 6 计算，取，采用下列算式计算：6)12(-=a 4.12≈（1）；6)12(1+（2）；27099-（3）；3)223(-（4）.3)223(1+问哪一个得到的结果最好？解 显然66666)12(1)12()12()12()12(+=++-=-=a []27099)223()12()12(3326-=-=-=-[]33266)223(1)12(1)12(1)12(+=+=+=-所以（1）≡（2）≡（3）≡（4），这4个算式是恒等的，但当取计算时，因为4.12≈（2），（3）都涉及到两个相近数相减，使有效数字损失，而（1）在分母算式上的乘幂数比算式（4）大，所以算式（4）最好，事实上，当取时，有|△x|<0.015，再由4.12≈的误差 也可直接估计出每个算式的误差，显然，算)(x f |||)4.1(||)()(x f x f x x f ∆'≈-∆+式（4）误差最小。
具体计算可行：（1）；36102.5)12(1-⨯≈+（2）0.127099≈-（3）；33100.8)223(-⨯≈-（4）.33101.5)223(1-⨯≈+比较可得用第（4）个算式所得的结果更接近于a 。
7 求二次方程x 2-(109+1)x+109=0的根。
解 由于x 2-(109+1)x+109=（x-109）(x-1)，所以方程的两个根分别为 x 1=109，x 2=1但如果应用一般二次方程ax 2+bx+c=0的求根公式：aac b b x 2422,1-±-=由于当遇到b 2>>4|ac|的情形时，有，则用上述公式求出的两个根中，ac b b 4||2-≈总有一个因用了两个相近的近似数相减而严重不可靠，如本例若在能将规格化的数表示到小数点后8位的计算机上进行计算，则-b=109+1=0.1×1010+0.000 000 0001×1010，由于第二项最后两位数“01”在机器上表示不出来，故它在上式的计算中不起作用，即在计算机运算时，-b=109.通过类似的分析可得9210||4=≈-b ac b 所以，求得的两个根分别为 99921102101024=+≈-+-=a ac b b x21010249922=+≈---=a ac b b x 显然，根x 2是严重失真的。
为了求得可靠的结果，可以利用根与系数的关系式：，在计算机上采用如下ac x x =21公式：aacb b b x 24)sgn(21---=12ax cx =其中，sgn （b ）是b 的符号函数，当b ≥0时sgn （b ）=1；当b<0时，sgn （b ）=-1。
显然，上述求根公式避免了相近数相减的可能性。
8 当N 充分大时，如何计算⎰++=111N Ndx xI 分析 函数的原函数已知，我们自然考虑用Newton-Leibniz 公式求这个定积分211x +的值。
由于N 很大，这样会遇到两个相近的数相减，因此，应采用一些变换公式来避免这种情况。
解 若用定积分的Newton-Leibniz 公式计算此题，有，则当N 充分大时，因为arctan （N+1）和arctanN⎰+-+=+12arctan )1arctan(11N NN N x 非常接近，两者相减会使有效数字严重损失，从而影响计算结果的精度，这在数值计算中是要尽量避免的，但是通过变换计算公式，例如：令tan θ1=N+1, tan θ2=N ，则由,得NN N N N N )1(11)1(11tan tan 1tan tan )tan(212121++=++-+=+-=-θθθθθθNN N N )1(11arctanarctan )1arctan(21++=-+=-θθ就可以避免两相近数相减引起的有效数字损失，从而得到较精确的结果。
所以，当N 充分大时，用计算积分的值较好。
⎰+++=+12211arctan 11N NN N x 9 计算积分.⎰==-11,2,1( n dx ex I x n n 分析 数值计算中应采用数值稳定的算法，因此在建立算法时，应首先考虑它的稳定性。
解 利用分部积分法，有⎰⎰⎰⎰--------=-==1011111110110111|dxe x n dx nx e e x de x dx e x x n n x x n x n x n 得递推公式：（1）1(1,2,)n n I I nI n -=-=⎰-==-1010011edx e x I x 利用公式（1）计算I n ，由于初值I 0有误差，不妨设求I 0的近似值时有大小为ε的误差，*0I 即ε+=*00I I 则由递推公式（1）得εε-=--=-=**1001I I I I I I εε!22222112+=+-=-=**I I I I I I εε!3!23333223+=⨯--=-=**I I I I I I εε!4!34444334+=⨯+-=-=**I I I I I I ┊ε!)1(n I I n n n-+=*显然初始数据的误差ε是按n!的倍数增长的，误差传播得快，例如当n=10时，10！≈3.629×106,，这表明I 10时已把初始误差ε扩大了很多倍，从而ε!10||1010=-*I I 的误差已把I 10的真值淹没掉了，计算结果完全失真。
*10I 但如果递推公式（1）改成)2,3,1,()(11 -=-=-k k n I I nI n n 于是，在从后往前计算时，I n 的误差减少为原来的，所以，若取n 足够大，误并逐步减n1小，显然，计算的结果是可靠的。
所以，在构造或选择一种算法时，必须考虑到它的数值稳定性问题，数值不稳定的算法是不能使用的。
10 为了使计算32)1(6)1(41310---+-+=x x x y 的乘除法运算次数尽量地少，应将表达式改写为怎样的形式？解 设.))64(3(10,11t t t y x t -++=-=在数值计算中，应注意简化运算步骤，减少运算次数，使计算量尽可能小。
11若x*=3587.64是x 的具有六位有效数字的近似值，求x 的绝对误差限。
12为使的近似值的相对误差小于0.1，问查开方表时，要取几位有效数字？7013利用四位数学用表求x=1-cos2°的近似值，采用下面等式计算：（1）1-cos2°（2）2sin 21°问哪一个结果较好？14求方程x 2-56x+1=0的两个根，使它至少具有四位有效数字（已知）。
982.27783≈15数列满足递推公式0}{=∞n x),2,1(1101 =-=-n x x n n 若取(三位有效数字)，问按上述递推公式，从x 0计算到x 10时误差有多大？这41.120≈x 个计算过程稳定吗？16如果近似值的相对误差限m n n a a a a x 10)101010(123121⨯⨯++⨯+⨯+±=+---* 小于，证明：这个数具有n 位有效数字。
1110)1(21+-⨯+n a第二章 插值法与数值微分1 已知，试利用插值法近似计算。
12144,11121,10100===115分析 由题中已知条件本题可利用三点二次Lagrange 插值，也可利用三点二次Newton 插值，它们所得结果相同。
解 利用三点二次Lagrange 插值。
记，则的二次12,11,10,144,121,100,)(210210=======y y y x x x x x f )(x f Lagrange 插值多项式为))(())(())(())(()(210120120102102x x x x x x x x y x x x x x x x x y x L ----+----= ))(())((1202102x x x x x x x x y ----+ )144121)(100121()144)(100(11)144100)(121100()144)(121(10----⨯+----⨯=x x x x )121144)(100144()121)(100(12----⨯+x x )115(115)115(2L f ≈=)144121)(100121()144115)(100115(11)144100)(121100()144115)(121115(10----⨯+----⨯=756722.10)121144)(100144()121115)(100115(12≈----⨯+因为，25232183)(,41)(,21)(---='''-=''='x x f x x f x x f )()()(22x L x f x R -= )144,100(),)()()((!31210∈---'''=ξξx x x x x x f 所以|)115()115(||)115(|22L f R -= |)144115)(121115)(100115(8361|25---⨯⨯=-ξ≤22510125163.029*******361--⨯=⨯⨯⨯⨯⨯2 已知的函数表)(x f y =x i 012y i8-7.5-18求函数在[0,2]之间的零点的近似值。
