计算方法_习题第一、二章答案..
计算方法各习题及参考答案

第二章数值分析已知多项式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)的误差阶尽可能高,并指出其阶数•高'为五阶。
计算方法一二章答案

4
x3=0.3466 x7=0.3572 ∴ x ≈ 0.3574
1) 2) 4) 3)
x=1+1/x2 x3=1+x2 x2=x3-1 x2=1/(x-1)
方程求根
解:1) x 1 1 x2
|1’(x)|= | -2 1 x3 |= 2
(x)
1 1.53 | x0=1.5 =0.59 <1(收敛)
2) x 3 1 x 2
| 2’(x)|= | 1 3
(1)单调区间:
令f’(x)=ex-4=0, x=ln4≈1.4,所以有两个单调区间: [- ∞,1.4](递减)和[1.4, ∞](递增)
(2)有根区间:∴ 存在两个有根区间为:[0,1] 和[2,3]
[- ∞,1.4]区间:f(0)=1>0,f(1)=e-4<0,所以有根区间为:[0,1] [1.4,+ ∞]区间:f(2)=e2-8<0,f(3)=e3-12>0,所以有根区间为:[2,3]
方程求根
3:用简单迭代法求方程ex-4x=0的根,并验 证收敛性,精确到4位有效数字。
解:2.在区间[0,1]上构造收敛的公式并计算
x=ln(4x)= φ2(x) (1)两种等价形式: x=ex/4=φ1(x); xk (2) x=ex/4=φ1(x): e |φ1’(x)|=ex/4<1 (收敛), 迭代公式为: xk 1
计算方法课后习题集规范标准答案

习 题 一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.什么是计算方法?(狭义解释)答:计算方法就是将所求的的数学问题简化为一系列的算术运算和逻辑运算,以便在计算机上编程上机,求出问题的数值解,并对算法的收敛性、稳定性和误差进行分析、计算。
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位浮点数计算。
计算方法_课后习题答案

L3 x 的最高次项系数是 6,试确定 y1 。
解: l0 (x)
x x1 x0 x1
x x2 x0 x2
x x3 x0 x3
x 0.5 0 0.5
x 1 0 1
x2 02
= x3
7 2
x2
7 2
x 1
l1 ( x)
x x0 x1 x0
(2 2e1 4e0.5 )x2 (4e0.5 e1 3)x 1
2)根据Lagrange余项定理,其误差为
| R2 (x) ||
f
(3) ( 3!
)
21
(
x)
||
1 6
e
x(
x
1)(
x
0.5)
|
1 max | x(x 1)(x 0.5) |, (0,1) 6 0x1
x2 02
x4= 04
x3
7x2 14x 8 8
l1 ( x)
x x0 x1 x0
x x2 x1 x2
x x3 x1 x3
x0 1 0
x2 1 2
x4 1 4
=
x3
6x2 3
8x
l2 (x)
x x0 x2 x0
i j
而当 k 1时有
n
x jl j
j0
x
n
n
j0 i0 i j
x xi x j xi
x
j
计算方法习题二答案

计算方法习题二答案习题二1、利用二分法求方程f(x)=x3-2x-5=0,在2,3内根的近似值,并指出误差。
解:f(2)=-1<0 f(3)=19>0 f(2).f(3)<0f’(x)=3x2-2 在x∈2,3f’(x) >0所以在1,2上必仅有一根x=2 f(2)=-1 -x=3 f(3)=16 +x=2.5 f(2.5)=5.625 +x=2.25 f(2.25)=1.890625 +x=2.125 f(2.125) +x=2.0625 f(2.0625) -x=2.09375 f(2.09375) -x=2.109375 f(2.109375) +x=2.1015625 f(2.1015625) +所以x=2.109375+2010156252=2.097656252、证明方程1-x-sinx=0在0,1内有一个根,使用二分法求误差不大于12×10?4的根。
解:令f(x)=1-x-sinxf(0)=1f(1)=-sin1f(0).f(1)<0f’(x)=-1-cosx<0在0,1恒成立所以1-x-sinx=0在0,1内恒有一个根n≥ln1?0?ln?(12×10?4)ln2-1≈13.289所以n=14n a n b n x n+1f(x n+1)符号0 0 1 0.5 +1 0.5 1 0.75 +2 0.875 1 0.9375 +..143、能不能用迭代法求解下列方程,若不能时,将方程改写成能用迭代法的形式。
(1、)x=(cosx+sinx)/4 (2)x=4-2x解:(1、)f(x)=x=(cosx+sinx)/4f’(x)=?sinx+cosx4<1对x任何数恒成立所以可用迭代法设x0=0,则x1=0.25x2=0.2511x 3=0.2511所以x=0.251(2、)f(x)=4-2xf’(x)=x.2x ?1<0在x 为任意数不恒成立所以不能用迭代法令x=log 2(4?x)x 0=0x 1=2x 2=1x 3= |φ‘(x)|=|-14?x 1ln 2|对x ∈(1,2)<124、为求方程x 3-x 2-1=0在x 0=1.5附近的一个根,设将方程改写成下列等价形式,并建立相应的迭代公式。
计算方法教程(第2版)习题答案

《计算方法教程(第二版)》习题答案第一章 习题答案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[中的近似根。
数值计算方法课后习题答案(李庆扬等)

第一章 绪论(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 。
计算方法-刘师少版第二章课后习题完整答案

0 < λf ′(x) < 2
− 2 < −λf ′(x) < 0
−1 < 1 − λf ′(x) < 1
1 − λf ′(x) < 1
即 ϕ ′(x) < 1 ,所以 xk+1 = ϕ (xk ) = xk − λf (xk ) 收敛于 f (x) =0 的根。
2.7 试用牛顿迭代法导出下列各式的迭代格式:
应的迭代公式:
(1)x
=1+
1 x2
,迭代公式
xk
+1
=1+
1
x
2 k
(2)x3 = 1 + x 2 ,迭代公式 xk+1 = 3 1 + xk2
(3) x 2
=
x
1 −
1
,迭代公式
xk
+1
=
1 xk −1
(4) x = x3 − 1 ,迭代公式
xk+1 = xk3 − 1
试分析每种迭代公式的收敛性,并选取一种收敛迭代公式求出具有四位有效数字的近似
x8 = 1.4656344
x9 = 1.4656000
x9
− x8
≤ 1 ×10−4 , 2
x9 = 1.4656000
2.5 对于迭代函数ϕ (x) = x + C(x 2 − 2) ,试讨论:
(1) 当 C 取何值时, xk+1 = ϕ (xk ), (k = 0,1,2,L) 产生的序列 {xk }收敛于 2 ;
6 6x2
63
ϕ ′(3 a ) == 5 − a (3 a )−3 = 5 − 1 = 1 ≠ 0
63
数值计算方法》习题答案

《数值计算方法》课后题答案详解吉 林 大 学第一章 习 题 答 案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.*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位有效数字的结果。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
解a1是1到9间的数字。
设x*具有n位有效数字,令-n+1=-1,则n=2,从而x*至少具有2位有效数字。
4计算sin1.2,问要取几位有效数字才能保证相对误差限不大于0.01%。
分析本题应利用有效数字与相对误差的关系。
解设取n位有效数字,由sin1.2=0.93…,故a1=9。
分析函数 的原函数已知,我们自然考虑用Newton-Leibniz公式求这个定积分的值。由于N很大,这样会遇到两个相近的数相减,因此,应采用一些变换公式来避免这种情况。
解若用定积分的Newton-Leibniz公式计算此题,有 ,则当N充分大时,因为arctan(N+1)和arctanN非常接近,两者相减会使有效数字严重损失,从而影响计算结果的精度,这在数值计算中是要尽量避免的,但是通过变换计算公式,例如:令tanθ1=N+1, tanθ2=N,则由 ,得
通过类似的分析可得
所以,求得的两个根分别为
显然,根x2是严重失真的。
为了求得可靠的结果,可以利用根与系数的关系式: ,在计算机上采用如下公式:
其中,sgn(b)是b的符号函数,当b≥0时sgn(b)=1;当b<0时,sgn(b)=-1。显然,上述求根公式避免了相近数相减的-cos2°
(2)2sin21°
问哪一个结果较好?
14求方程x2-56x+1=0的两个根,使它至少具有四位有效数字(已知 )。
15数列 满足递推公式
若取 (三位有效数字),问按上述递推公式,从x0计算到x10时误差有多大?这个计算过程稳定吗?
16如果近似值 的相对误差限小于 ,证明:这个数具有n位有效数字。
因而x2具有3位有效数字。
由π- =3.141 59…-3.142 85…=-0.001 26…知
因而x3具有3位有效数字。
2已知近似数x*有两位有效数字,试求其相对误差限。
分析本题显然应利用有效数字与相对误差的关系。
解利用有效数字与相对误差的关系。这里n=2,a1是1到9之间的数字。
3已知近似数的相对误差限为0.3%,问x*至少有几位有效数字?
具体计算可行:
(1) ;
(2)
(3) ;
(4) .
比较可得用第(4)个算式所得的结果更接近于a。
7求二次方程x2-(109+1)x+109=0的根。
解由于x2-(109+1)x+109=(x-109)(x-1),所以方程的两个根分别为
x1=109,x2=1
但如果应用一般二次方程ax2+bx+c=0的求根公式:
再计算。又例如,当x充分大时,应作变换
6计算 ,取 ,采用下列算式计算:
(1) ;
(2) ;
(3) ;
(4) .
问哪一个得到的结果最好?
解显然
所以(1)≡(2)≡(3)≡(4),这4个算式是恒等的,但当取 计算时,因为(2),(3)都涉及到两个相近数相减,使有效数字损失,而(1)在分母算式上的乘幂数比算式(4)大,所以算式(4)最好,事实上,当取 时,有|△x|<0.015,再由 的误差 也可直接估计出每个算式的误差,显然,算式(4)误差最小。
第二章插值法与数值微分
1已知 ,试利用插值法近似计算 。
分析由题中已知条件本题可利用三点二次Lagrange插值,也可利用三点二次Newton插值,它们所得结果相同。
解利用三点二次Lagrange插值。
就可以避免两相近数相减引起的有效数字损失,从而得到较精确的结果。所以,当N充分大时,用 计算积分的值较好。
9计算积分 .
分析数值计算中应采用数值稳定的算法,因此在建立算法时,应首先考虑它的稳定性。
解利用分部积分法,有
得递推公式:
(1)
利用公式(1)计算In,由于初值I0有误差,不妨设求I0的近似值 时有大小为ε的误差,即
解不等式 知取n=4即可满足要求。
5计算 ,视已知数为精确值,用4位浮点数计算。
解 0.131 8×10-2-0.131 6×10-2=0.2×10-5
结果只有一位有效数字,有效数字大量损失,造成相对误差的扩大,若通分后再计算:
就得到4位有效数字的结果。
此例说明,在数值计算中,要特别注意两相近数作减法运算时,有效数字常会严重损失,遇到这种情况,一般采取两种办法:第一,应多留几位有效数字;第二,将算式恒等变形,然后再进行计算。例如,当x接近于0,计算 时,应先把算式变形为
第一章误差
1问3.142,3.141, 分别作为π的近似值各具有几位有效数字?
分析利用有效数字的概念可直接得出。
解π=3.141 592 65…
记x1=3.142,x2=3.141,x3= .
由π-x1=3.141 59…-3.142=-0.000 40…知
因而x1具有4位有效数字。
由π-x2=3.141 59…-3.141=-0.00059…知
10为了使计算
的乘除法运算次数尽量地少,应将表达式改写为怎样的形式?
解设
在数值计算中,应注意简化运算步骤,减少运算次数,使计算量尽可能小。
11若x*=3587.64是x的具有六位有效数字的近似值,求x的绝对误差限。
12为使 的近似值的相对误差小于0.1,问查开方表时,要取几位有效数字?
13利用四位数学用表求x=1-cos2°的近似值,采用下面等式计算:
由于当遇到b2>>4|ac|的情形时,有 ,则用上述公式求出的两个根中,总有一个因用了两个相近的近似数相减而严重不可靠,如本例若在能将规格化的数表示到小数点后8位的计算机上进行计算,则-b=109+1=0.1×1010+0.000 000 0001×1010,由于第二项最后两位数“01”在机器上表示不出来,故它在上式的计算中不起作用,即在计算机运算时,-b=109.
则由递推公式(1)得
┊
显然初始数据的误差ε是按n!的倍数增长的,误差传播得快,例如当n=10时,10!≈3.629×106, ,这表明I10时已把初始误差ε扩大了很多倍,从而 的误差已把I10的真值淹没掉了,计算结果完全失真。
但如果递推公式(1)改成
于是,在从后往前计算时,In的误差减少为原来的 ,所以,若取n足够大,误并逐步减小,显然,计算的结果是可靠的。所以,在构造或选择一种算法时,必须考虑到它的数值稳定性问题,数值不稳定的算法是不能使用的。