计算方法答案
计算方法习题及答案
计算方法习题及答案在学习计算方法的过程中,习题的练习和答案的掌握是非常重要的。
下面将为大家提供一些计算方法习题及答案,希望能够帮助大家更好地巩固知识。
一、整数运算习题1. 计算以下整数的和:-5 + 8 + (-3) + (-2) + 10。
答案:-5 + 8 + (-3) + (-2) + 10 = 8。
2. 计算以下整数的差:15 - (-6) - 10 + 3。
答案:15 - (-6) - 10 + 3 = 24。
3. 将 -3 × (-4) - 2 × 5 的结果化简。
答案:-3 × (-4) - 2 × 5 = 12 - 10 = 2。
二、分数运算习题1. 计算以下分数的和:1/2 + 2/3 + 3/4 + 4/5。
答案:1/2 + 2/3 + 3/4 + 4/5 = 47/20。
2. 计算以下分数的差:2/3 - 1/4 - 5/6。
答案:2/3 - 1/4 - 5/6 = -1/12。
3. 计算以下分数的积:2/3 × 3/4 × 4/5。
答案:2/3 × 3/4 × 4/5 = 4/15。
4. 将以下分数的除法化简为整数:3/8 ÷ 1/4。
答案:3/8 ÷ 1/4 = (3/8) × (4/1) = 3/2 = 1 1/2。
三、百分数运算习题1. 计算60% × 80%的结果。
答案:60% × 80% = 0.6 × 0.8 = 0.48 = 48%。
2. 计算40%除以20%的结果。
答案:40% ÷ 20% = (40/100) ÷ (20/100) = 2。
3. 计算200中的20%是多少。
答案:200 × 20% = 200 × 0.2 = 40。
四、多项式运算习题1. 计算以下多项式的和:(3x^2 + 4x + 5) + (2x^2 + x + 3)。
计算方法各习题及参考答案
计算⽅法各习题及参考答案第⼆章数值分析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]-上的3次最佳⼀致逼近多项式()p x .答案:()f x 的最佳⼀致逼近多项式为323()74p x x x =++. 3.4 ⽤幂级数缩合⽅法,求() (11)x f x e x =-≤≤上的3次近似多项式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 确定下列求积公式中的待定参数,使得求积公式的代数精度尽量⾼,并指明所确定的求积公式具有的代数精度.(1)101()()(0)()hh f x dx A f h A f A f h --≈-++?(2)11211()[(1)2()3()]3f x dx f f x f x -≈-++? (3)20()[(0)()][(0)()]2h h f x dx f f h h f f h α''≈++-?答案:(1)具有三次代数精确度(2)具有⼆次代数精确度(3)具有三次代数精确度.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 =(1)运⽤复化梯形公式计算上述积分值,使其截断误差不超过31102-?.(2)取同样的求积节点,改⽤复化⾟浦⽣公式计算时,截断误差是多少?(3)要求的截断误差不超过610-,若⽤复化⾟浦⽣公式,应取多少个节点处的函数值?答案:(1)只需7.5n ≥,取9个节点,0.946I ≈(2)4(4)46111|[]||()|()0.271102880288045n b a R f h f η--=-≤=? (3)取7个节点处的函数值.4.6 ⽤变步长的复化梯形公式和变步长的复化⾟浦⽣公式计算积分10sin xI dx x =?.要求⽤事后误差估计法时,截断误不超过31102-?和61102-?.答案:使⽤复化梯形公式时,80.946I T ≈=满⾜精度要求;使⽤复化⾟浦⽣公式时,40.946 083I s ≈=满⾜精度要求.4.7(1)利⽤埃尔⽶特插值公式推导带有导数值的求积公式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 ηη-=∈.(2)利⽤上述公式推导带修正项的复化梯形求积公式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 =.第六章解线性代数⽅程组的迭代法6.1对⽅程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 =.6.2讨论松弛因⼦ 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 =-.6.3给定线性⽅程组Ax b =,其中111221112211122A ?? ? ?=,证明⽤雅可⽐迭代法解此⽅程组发散,⽽⾼斯-赛得尔迭代法收敛.6.4设有⽅程组112233302021212x b x b x b -?????? ??? ?= ??? ? ??? ?-??????,讨论⽤雅可⽐⽅法和⾼斯-赛得尔⽅法解此⽅程组的收敛性.如果收敛,⽐较哪种⽅法收敛较快.答案:雅可⽐⽅法收敛,⾼斯-赛得尔⽅法收敛,且较快.6.5设矩阵A ⾮奇异.求证:⽅程组Ax b =的解总能通过⾼斯-赛得尔⽅法得到.6.6设()ij n nA a ?=为对称正定矩阵,对⾓阵1122(,,,)nn D diag a a a = .求证:⾼斯-赛得尔⽅法求解⽅程组1122D AD x b --=时对任意初始向量都收敛.第七章⾮线性⽅程求根例7.4对⽅程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 ≈=.例7.6对于迭代函数2()(3)x x c x ?=+-,试讨论:(1)当c 为何值时,1()k k x x ?+=产⽣的序列{}k x(2)c 取何值时收敛最快?(3)取1,2c =-()x ?51||10k k x x -+-<.答案:(1)(c ∈时迭代收敛.(2)c =时收敛最快.(3)分别取1, 2c =--,并取0 1.5x =,计算结果如下表7.7所⽰表7.7例7.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 ''=' 例7.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 ?? ? ?= ? ? ???的最⼩特征值,使它⾄少具有2位有效数字.答案:取5 2.234375λ≈-即有2位有效数字.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 =,⼩数点后保留8位.答案: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 ≈=??≈=?。
计算方法-刘师少版第一章课后习题完整答案
分, 试给出此递推公式误差的传播规律, 计算 I 10 时误差被放大了多少倍?这个算法是数值稳定的 吗? 解: I =
∫x
0 1 0
1
n
e x −1 dx , n = 0,1,2,L,10 ,由分部积分法有
1 0
n −1 x −1 I n = ∫ x n e x −1 dx = x n e x −1 1 e dx 0 − n∫ x
er ( x n ) =
e( x n ) nx n −1 ( x − x * ) x − x* = = n = n ⋅ er ( x) = αn% x xn xn
x n 的相对误差为 an%
1.10 设 x>0,x 的相对误差为 δ ,求 ln x 的误差。 解: e(ln x) ≈
1 ( x − x * ) = er ( x) = δ x
N +1
N
1 dx = arctan( N + 1) − arctan N 1+ x2 1 = arctan 1 + N ( N + 1) 1 2 gt ,假定 g 是准确的,而对 t 的测量有±0.1s 的误差,证明当 t 增加时,s 的绝对误差 2
1.12 设 s =
增加,而相对误差减少。 解:由题意知, e( s ) = s − s = gt (t − t ) = gt ⋅ e(t ) = 0.1gt
5
计算方法
于是
* * * * e( I 10 ) = −10e( I 9 ) = 10 ⋅ 9e( I 8 ) = L = 10!e( I 0 )
计算 I 10 时的误差被扩大了 10 倍,显然算法是数值不稳定的 1.14 设 f ( x) = 8 x − 0.4 x + 4 x − 9 x + 1 ,用秦九韶算法求 f (3)
计算方法-习题第一、二章答案.doc
第一章 误差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之间的数字。
%5101211021|*||*||)(|1211*=⨯⨯≤⨯≤-=+-+-n ra 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 ra x x x x ε解不等式411101021-+-≤⨯n a 知取n=4即可满足要求。
计算方法课后习题答案
0
1
2
3
2
1.888889
1.879452
1.879385
解得
用弦截法求解
取
依迭代公式为 进行计算。
计算结果列于下表,并和 比较
0
1
2
3
4
2
1.9
1.881094
1.879411
1.879385
解得
用抛物线法求解
则
故 则根号前的符号为正。
迭代公式为
取 计算
10.设
(3)如果要求截断误差不超过 ,那么使用复化Simpson公式计算时,应将积分区间分成多少等分?
解:(1)
= ,
当误差 时, 25.6,所以取 =26。
(2)
7.推导下列三种矩形求积公式:
证明: 将 在 处Taylor展开,得
两边在 上积分,得
将 在 处Taylor展开,得
两边在 上积分,得
将 在 处Taylor展开,得
(1)
依Taylor公式有
代人式(1)右端,则有
另一方面,
故隐式Euler格式的局部截断误差为
可见隐式Euler格式 是一阶方法。
证明 :对于Euler两步格式 : ,考察局部截断误差
,仍设 则有
注意到
于是
而
因此有
即Euler两步格式 是二阶方法。且其主项系数是2。
特别地,当 时,有
而当 时有
5.依据下列函数表分别建立次数不超过3的 插值多项式和 插值多项式,并验证插值多项式的唯一性。
0
1
2
4
1
9
23
3
解:
智慧树答案计算方法知到课后答案章节测试2022年
绪论1.工科人认知世界的时候,需要认可误差的存在。
答案:对2.不需要掌握很多很好的计算方法,单单凭借计算机强大的能力就可以解决大部分实际问题。
答案:错3.计算方法是一门理论数学课,可以获得寻求数学问题的精确解析解的知识答案:错4.解决某些实际问题时,选择不合适的计算方法有可能无法得到满意的结果。
答案:对5.求解高阶线性方程组(比如,大于150阶),用克莱姆法则来直接求解也是可以接受的。
答案:错第一章1.计算机进行乘除运算时按照先舍入后运算的原则。
()答案:错2.相对误差是个无名数,没有量纲。
()答案:对3.两近似值之商的相对误差等于被除数的相对误差与除数的相对误差之差。
()答案:对4.用 1+x近似表示所产生的误差是( )答案:截断误差5.设某数x,那么x的有四位有效数字且绝对误差限是的近似值是()答案:0.6930第二章1.若f(a)f(b)<0 ,则f(x)在(a,b)内一定有根。
()答案:错2.如果迭代格式在根的附近导数值的模大于1,则迭代发散。
()答案:对3.若x*是f(x)=0的重根,则牛顿不收敛。
()答案:错4.非线性方程的求根方法中,正割法收敛速度比Newton迭代法快。
()答案:错5.用牛顿迭代法求方程f(x)=在附近的根,第一次迭代值()答案:26.用简单迭代法求方程f(x)=0的实根,把方程f(x)=0表示成x=g(x),则f(x)=0的根是()答案:y=x与y=g(x)交点的横坐标7.以下对非线性方程的求根方法中哪些是线性收敛的?()答案:牛顿迭代法在重根的情况下;不动点迭代法第三章1.上面式子,用高斯消去法计算,解为()。
答案:x1=9,x2=-1,x3=-62.用三角分解法计算上式子,则L矩阵中第二行第一列元素的值为()。
答案:03.设上面矩阵,计算的大小为()。
答案:104.设矩阵,计算的大小为()。
答案:85.设上面矩阵,计算矩阵A的谱半径为()。
答案:56.高斯消去法属于线性方程组数值解法中的迭代法。
计算方法课后习题集规范标准答案
习 题 一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.掌握绝对误差(限)和相对误差(限)的定义公式。
计算方法_课后习题答案
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.下列各数都是经过四舍五入得到的近似数,指出他们有几位有效数字,并写出绝对误差限。
9800107480.566.385031.01021.1*65*5*4*3*2*1=⨯=====x x x x x x解: ,有5位有效数字,绝对误差限为; ,有2位有效数字,绝对误差限为; ;有4位有效数字,绝对误差限为; ;有5位有效数字,绝对误差限为; ;有1位有效数字,绝对误差限为; ;有4位有效数字,绝对误差限为。
2.要使的近似值的相对误差限小于,要取几位有效数字?解:由于,设要取位有效数字,则根据定理1.1.1,有()()%1.010811021111<⨯=⨯≤----n n r x ε,解得4≥n ,即要取4位有效数字。
3.序列满足递推关系若,计算到时误差有多大?这个计算过程数值稳定吗?解:,由于有3位有效数字,且,所以的绝对误差限为,因此的绝对误差限为。
很明显这个计算过程不是数值稳定的。
作业中出现的问题:第一题:主要是第五个数5*5107⨯=x ,不知道它有几位有效数字,很多同学认为有5或者6位有效数字,这是不对的,进而算错绝对误差限。
另外有个别同学分不清有效数字的概念,六个数的有效数字都弄错了。
第二题:主要是算错n ,不知道该取3还是4。
第三题:没有什么大的问题。
有个别同学一个数一个数的算出来了,这是不可取的。
直接迭代误差就行了。
附:地物1301班和1302班有几个同学花名册上没有名单,我添加上去了。
第二次作业1*11011021.01021.1⨯==x 4-5-1105.0105.0⨯=⨯1-*21031.0031.0⨯==x 3-2-1-105.0105.0⨯=⨯3*3103856.06.385⨯==x -14-3105.0105.0⨯=⨯2*41056480.0480.56⨯==x 3-5-2105.0105.0⨯=⨯65*5107.0107⨯=⨯=x 51-6105.0105.0⨯=⨯4*6109800.09800⨯==x 5.0105.04-4=⨯20%1.0110447213595.047213595.420⨯⋯=⋯=n {}n y ,,2,1,1101⋯=-=-n y y n n 41.120≈=y 10y ()()()*00*222*11*101010y y y y y y y y n n n n n n n -=⋯=-=-=-----*0y 1*010141.041.1⨯==y *0y 2-105.0⨯*10y 72-10105105.010⨯=⨯⨯1.利用二分法求方程在[2,3]内根的近似值,并指出误差。
《计算方法》练习题及答案
《计算方法》练习题及答案1. 单选题1. 数值3.1416的有效位数为()A. 3B. 4C. 5D. 6正确答案:C2. 常用的阶梯函数是简单的()次样条函数。
A. 零B. 一C. 二D. 三正确答案:A3. 设求方程f(x)=0的根的牛顿法收敛,则它具有()敛速。
A. 超线性B. 平方C. 线性D. 三次正确答案:C4. 构造拟合曲线不可以采用下列哪种准则()A. 使残差的最大绝对值为最小B. 使残差的绝对值之和为最小C. 使残差的平方和为最小D. 是残差的绝对值之差为最小正确答案:D5. 欧拉法的局部截断误差阶为()。
A. AB. BC.CD. D正确答案:B6. 依据3个样点(0,1),(1,2)(2,3),其插值多项式p(x)为()A. xB. x+1C. x-1D. x+2正确答案:B7. 题面如下,正确的是()A. 2B. 3C. -2D. 1正确答案:B8. 题面如下图所示,正确的是()A. AB. BC. CD. D正确答案:D9. 用列主元消去法解线性方程组,A. 3B. 4C. -4D. 9正确答案:C10. 利用克莱姆法则求解行列式时,求解一个n阶方程组,需要()个n阶行列式。
A. nB. n+1C. n-1D. n*n正确答案:C11. 线性方程组的解法大致可以分为()A. 直接法和间接法B. 直接法和替代法C. 直接法和迭代法D. 间接法和迭代法正确答案:C12. ()的优点是收敛的速度快,缺点是需要提供导数值。
A. 牛顿法B. 下山法C. 弦截法D. 迭代法正确答案:A13. 设x* = 1.234是真值x = 1.23445的近似值,则x*有()位有效数字。
A. 1B. 2C. 3D. 4正确答案:D14. 若a=2.42315是2.42247的近似值,则a有( )位有效数字.A. 1B. 2C. 3D. 4正确答案:C15. 所谓松弛法,实质上是()的一种加速方法。
现代数值计算方法习题解答
现代数值计算方法习题答案习 题 一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)因为=204.4721…… , 又=)(*x E |*x x −| = |47.420−| = 0.0021 < 0.01, 所以 =*x4.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 c m .记*y 为y 的近似值,则)(20)(20)(2)(*****x E x x x x x y E =−=−= < = 0.1,所以)(*x E< = 0.005 c m . 7、解:因为)()(*1x x nx x E n n −≈−,所以n x nE x x x n xx E x E r nn nr 01.0)()()(*==−≈=. 8、解: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 ,所以|*00x 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)方程组的增广矩阵为:−−−−11114423243112M M M → −−−−1010411101110112M M M →−−−11041001110112M M M → 31=x , 12=x , 13=x . (2)方程组的增广矩阵为:−−−−−−017232221413M M M → −−247210250413M M M → −−147200250413M M M → 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 l 6/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 l 10/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 l5432106000000000600006006u 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 l3211001u 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 x 525351)(3)(1)1(2++−=+k k k x x x 5115152)(2)(1)1(3++=+k k k x x x 高斯-赛德尔迭代公式:545152)(3)(2)1(1+−=+k k k x x x 525351)(3)1(1)1(2++−=++k k k x x x5115152)1(2)1(1)1(3++=+++k k k x x x7、(1)证:因为此方程组的系数矩阵为严格对角占优矩阵,所以雅可比迭代法和相应的高斯-赛德尔迭代法都收敛。
计算方法习题集及答案
得:
当方法为零稳定时 ,从而 ,故方法是二阶收敛的。
6.给出题(6.5)题中 时的公式的绝对稳定域.
解:
6.5中当 时,即为方法
其相应的差分方程的多项式为
令 ,
即方法的绝对稳定域为
7.指出Heun方法
0
0
0
0
1/3
1/3
0
0
2/3
0
2/3
0
1/4
0
3/4
的相容阶,并给出由该方法以步长h计算初值问题(6.45)的步骤.
即
取 。即
满足上述条件的多步方法即为一类三步四阶显示方法,令 可得
方法即为
3.形如
的k阶方法称为Gear方法,试确定一个三步Gear方法,并给出其截断误差主项。
解:线性k步公式为
由Gear法的定义知,三步Gear法满足
方法为 阶,故有
得:
取 得
得三步Gear方法:
其中
4.试用显式Euler法及改进的Euler法
证明:
且
即 为 的二阶零点
设
令
易知
又
由微分中值定理(Rolle定理) ,使得
进而 有三个零点, 有两个零点, 有一个零点,
即 使得
得
8.设 是Lagrange基函数,则 。
9.求一个次数不超过4次的多项式 ,使它满足
,并写出其余项表达式。
10.求一个四次插值多项式 ,使 时, ;而 时, ,并写出插值余项的表达式。
练习
班级
学号
姓名
1.试构造迭代收敛的公式求解下列方程:
(1) ; (2) 。
解:
(1)迭代公式 , 公式收敛
k
奥鹏2024计算方法网考答案
奥鹏2024计算方法网考答案
一、单项选择
1.计算机的主要功能包括(﹒)。
A.记忆
B.输入
C.输出
D.运算
答案:D.运算
2.在Windows中,在文档的右上角可以看到一个小窗口,上面有最小化、最大化和关闭三个按钮,这个小窗口的名称是(﹒)。
A.工具栏
B.状态栏
C.菜单栏
D.标题栏
答案:D.标题栏
3.在离散数学中,可以使用(﹒)来表示“属于”的关系。
A.=
B.>
C.∈
D.≤
答案:C.∈
4.下面哪种类型的语言能够完成最为复杂的编程任务?(﹒)
A.汇编语言
B.C语言
C.高级语言
D.机器语言
答案:C.高级语言
5.下列哪种技术是面向对象程序设计的基础内容?(﹒)
A.函数
B.继承
C.多态
D.类
答案:D.类
二、多项选择
1.网络的基本元素包括(﹒)
A.节点
B.中继器
C.主机
D.链路
答案:A.节点B.中继器C.主机D.链路
2.以下哪些是关系型数据库技术的操作?(﹒)
A.插入
B.更新
C.删除
D.查找
答案:A.插入B.更新C.删除D.查找
3.以下哪些是算法分析的内容?(﹒)
A.时间复杂度
B.空间复杂度
C.动态规划
D.可视化答案:A.时间复杂度B.空间复杂度C.动态规划4.以下哪些是计算机病毒的特征?(﹒)。
计算方法课后习题答案
计算方法课后习题答案在计算方法课程中,学生通常会接触到各种数学问题的求解方法,包括但不限于数值分析、线性代数、微分方程等。
以下是一些课后习题的解答示例:习题一:求解线性方程组设线性方程组为:\[ \begin{align*}a_{11}x_1 + a_{12}x_2 + \cdots + a_{1n}x_n &= b_1, \\a_{21}x_1 + a_{22}x_2 + \cdots + a_{2n}x_n &= b_2, \\\vdots \quad \quad & \ \vdots \\a_{m1}x_1 + a_{m2}x_2 + \cdots + a_{mn}x_n &= b_m,\end{align*} \]解答:使用高斯消元法或矩阵分解法求解上述方程组。
首先将系数矩阵转换为行简化阶梯形式,然后回代求解未知数 \( x_1, x_2,\ldots, x_n \)。
习题二:数值积分给定函数 \( f(x) \),需要在区间 \( [a, b] \) 上进行数值积分。
解答:可以使用梯形法、辛普森法等数值积分方法。
例如,使用梯形法的公式为:\[ \int_a^b f(x)dx \approx \frac{h}{2} \left( f(a) + 2f(a+h) + 2f(a+2h) + \cdots + 2f(b-h) + f(b) \right), \]其中 \( h = \frac{b-a}{n} \) 是区间的等分宽度,\( n \) 是等分数。
习题三:常微分方程的数值解给定一个常微分方程 \( y' = f(x, y) \),初始条件为 \( y(x_0) = y_0 \)。
解答:使用欧拉法或龙格-库塔法求解。
以欧拉法为例,其迭代公式为:\[ y_{n+1} = y_n + h f(x_n, y_n), \]其中 \( h \) 是步长,\( x_{n+1} = x_n + h \)。
计算方法各章作业答案
计算方法作业第二章插值1.(1(2)用二次Lagrange插值多项式求当X=0.15时Y的近似值。
(3)写出余项R(x)=f(x)-Pn(x)的表达式。
解:(1)Pn (x) =knknkjj jkj yxxxx)(00∑∏=≠=--n=3P 3(x)=321321))()(())()((yxxxxxxxxxxxx------+13121132))()(())()((yxxxxxxxxxxxx------+23212231))()(())()((yxxxxxxxxxxxx------+32313321))()(())()((yxxxxxxxxxxxx------x 0=0.0 x1=0.1 x2=0.2 x3=0.3y 0=0.0000 y1=0.0998 y2=0.1987 y3=0.2955P 3(x)=0000.0)3.00.0)(2.00.0)(1.00.0()3.0)(2.0)(1.0(⨯------xxx+0998.0)3.01.0)(2.01.0)(0.01.0()3.0)(2.0)(0.0(⨯------xxx+1987.0)3.02.0)(1.02.0)(0.02.0()3.0)(1.0)(0.0(⨯------xxx+2955.0)2.03.0)(1.03.0)(0.03.0()2.0)(1.0)(0.0(⨯------xxx(2) y(0.15) = P2(0.15) = 0.1494(3)R(x) = f(x)-Pn (x)=)!1()()1(++nf nξnk0=∏(x - x k)=!4)(4ξf(x – 0.0) (x – 0.1)(x – 0.2)(x – 0.3)第三章 方程求根5.求解方程12-3x+2cosx=0的迭代法n n x x cos 3241+=+(1)证明对于任意的x 0€R 均有*lim x x n x =∞→ (x *为方程的根)(2)取x 0=4,用此迭代法求方程根的近似值,误差不超过10-3,列出各次的迭代值。
计算方法习题及答案
第一章 绪论一.填空题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、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
上述三种迭代法都收敛. 是 A 的特征值, 故当 0 < θ < 可得 0 < θ <
2 . λn 1 √ 2 2 时, λn
有|1 − θλi | < 1, 从而迭代收敛. 反之, 若迭代收敛, 则 |1 − θλi | < 1,
2.11 当 |a| <
时, 迭代格式收敛. −8x1 + x2 + x3 = −7 2.12 提示: 将原方程组调整为: 上述方程组的系数矩阵是严格对角占优的, x1 − 5x2 + x3 = 14 x1 + x2 − 4x3 = −13
sin2 1◦ , 1+cos 1◦
1.10 提示: (1) sin(x + y ) − sin x = 2 sin y cos(x + y ), (2) 1 − cos 1◦ = 2 2 √ √ = − ln( 1010 + 1 + 105 ). ln( 1010 + 1 − 105 ) = ln √ 10 1 5
10 +1+10
1.11 (1) (A) 比较准确; (2) (A) 比较准确. 1.12 算法 2 准确. 在算法 1 中, ε0 ≈ 0.2231 带有误差 0.5 × 10−4 , 而这个误差在以后的每次计算中 顺次以 41 , 42 , · · · 传播到 In 中. 而算法 2 中的误差是按
1 4 11 33 11 22 1 22 2 4
0.374 ± 0.868i, ρ(BJ ) = 0.945 < 1, 0 0 −1 BS = (D − L) U = 0 0 0 0
迭代法收敛. (2) 将原方程组的系数矩阵调整为上述矩阵后, 写出迭代矩阵: BJ ∥BJ ∥∞ < 1, 故迭代法收敛.
a > 1, a(a − 14) + 12 > 0,
2
(k+1) (k ) (k+1) 2 (k ) 2 (k ) 2 (k ) 7 x1 x2 − 2 x +7 , x x2 − 9 x3 + 9 , = −9 = −9 9 3 9 1 (k ) (k ) (k+1) (k) (k+1) (k+1) 4 2 4 2 2 2 2.9 Jacobi: x2 x2 − 5 x3 − 5 = − 5 x1 − 5 x3 − 5 , Gauss-Seidel: = − 5 x1 , (k+1) (k+1) (k ) 2 (k ) 2 2 (k+1) 4 (k+1) 2 x3 x1 − 4 x . x x x . = −5 − = − − − 3 5 2 5 5 1 5 2 5 (k+1) (k ) (k ) (k ) 2 (k ) − x1 − 2 x −9 x3 ), = x1 + 1.32( 7 x1 9 9 2 (k+1) (k ) (k+1) (k ) (k ) SOR 迭代 (ω = 1.32): x2 −2 x − x2 − 4 x ), = x2 + 1.32(− 2 5 5 1 5 3 (k+1) (k+1) (k+1) (k) (k ) −2 x x x3 = x3 + 1.32(− 2 −4 − x3 ). 5 5 1 5 2 2.10 提示: x(k+1) = (I − θA)x(k) + θb, 故迭代矩阵 B = I − θA 的特征值为 1 − θλi , 其中 λi > 0
《现代数值计算方法(MATLAB版) 》
习题参考答案及部分习题解答提示 第一章
1.1 (1) 0.5, 0.00217%, 5; (2) 0.5×10−5 , 0.217%, 3; (3) 0.5×10−2 , 0.000217%, 6; (4) 0.5×102 , 0.0217%, 3. 1.2 (1) 0.5, 0.014%, 4; (2) 0.5×10−4 , 0.11%, 3; (3) 0.5×10−3 , 0.0017%, 5; (4) 0.5×10−9 , 0.017%, 4. 1.3 (1) 3.146, 0.5×10−4 ; (2) 3.1416, 0.5×10−4 ; (3) 3.14159. 1.4 提示:
∥A
−1
2.5 提示: 因系数矩阵 A 是严格对角占优的, 故 Jacobi 迭代法, Gauss-Seidel 迭代法以及 SOR 迭代 法在 0 < ω < 2 时, 都是收敛的.
2.6 提示: (1) BJ = −1 −2 0
−2 0 −2
λ −1 , |λI − BJ |= 1 0 2 0 −2 2 0 1 0 0 2 2 λ 2
1 4n
减少的, 是稳定的计算公式.
第二章
2.1 提示: 因 B 奇异, 故 ∃x ̸= 0, 使得 Bx = 0. 于是, Ax = (A − B )x,x = A−1 (A − B )x,∥x∥ ≤
1 . ∥∥A − B ∥∥x∥, 1 ≤ ∥A−1 ∥ · ∥A − B ∥,即∥A−1 ∥ ≥ ∥A− B∥ √ √ 2.2 ∥x∥1 = 9, ∥x∥2 = 29, ∥x∥∞ = 4; ∥A∥1 = 8, ∥A∥2 = 4 2, ∥A∥∞ = 6. 0 −0.4 −0.4 λ 0.4 0.4 2 2.3 提示: 迭代矩阵 BJ = −0.4 0 −0.8 , |λI − BJ |= 0.4 λ 0.8 =(λ − 0.8)(λ + −0.4 −0.8 0 0.4 0.8 λ √ √ 0.8λ − 0.32) = 0, λ1 = 0.8, λ2,3 = 0.4(−1 ± 3). 因 |λ3 | = 0.4(1 + 3) > 1, 故用 Jacobi 迭代法不收敛. −1 1 0 0 0 −2 2 1 0 0 0 −2 2 2.4 提示: Bs = 1 1 0 0 0 −1 = −1 1 0 0 0 −1 = 2 2 1 0 0 0 0 −2 1 0 0 0 0 −2 2 0 2 −3 . 谱半径 ρ(Bs ) = max |λi | = 2 > 1, 故用 Gauss-Seidel 迭代法不收敛. 0 0 2
1 2(a1 +1)
× 10−(n−1) = 10−4 ⇒ n = 5 − lg 2 − lg(a1 + 1) ⇒ 4 − lg 2 ≤ n ≤ 5 − 2 lg 2 ⇒ × 10−2 ≤ 0.5 × 10−2 . + 1) ⇒ 3 − lg 6 ≤ n ≤ 3 − lg 1.2 ⇒ √ 783 = (3)
. 因为
0
5 − 33
0
第三章
3.1 (1) x = (0, −1, 1)T ;
(i)
(2) x = (1.2, 2, −1.4).
3.2
1 3 n 3
+ n2 − 1 n. 3
3.3 提示: 必要性. 设 aii ̸= 0, i = 1, 2, · · · , k , 则可进行消去法的 k − 1 步. 每步 A(m) 由 A 逐次实 施 (−lij Ej + Ei ) → (Ei ) 的运算得到, 这些运算不改变相应顺序主子式之值, 所以有, a11 ∆m =
−1 2
0 0
−1
0 0 0
−1 = 0
1 2
−1 2
−1 2 0
1 , 其特征 −2 1 −2
1 0 1
故 ρ(Bs ) = < 1, 从而 Gauss-Seidel 迭代收敛. −1 √ 2 5 5 −2 , |λI − BJ | = λ(λ + 4 ) = 0,λ1 = 0, λ2,3 = ± 2 i, ρ(BJ ) = 0 a 1 a 2 3
3.699 ≤ n ≤ 4.3976 ⇒ n = 3. 1.5 |εr (x)| ≤ 1.6 提示: 1.7 提示:
1 2a1 1 × 10−(n−1) < 10−3 ⇒ n > 4 − 2 lg 2 ⇒ n = 4. 2×2 1 × 10−(n−1) = 3 × 10−3 ⇒ n = 4 − lg 6 − lg(a1 2(a1 +1) 56+ 562 −4 2
2.2218 ≤ n ≤ 2.9208 ⇒ n √ = 2. 1.8 提示: x1,2 =
282 − √781 28+ 783
= 28 ±
√
783, x1 = 28 + 27.982 = 55.982 ≈ 55.98, x2 = 28 −
1−cos2 1◦ 1+cos 1◦
=
1 55.982
≈ 0.01786. =
√
5 2
> 1, 故
2.8 提示: (1) A = 1 3 a > 1, ⇒ a3 − 14a + 12 > 0, Seidel 迭代收敛.
a > 0, a 2 − 1 > 0, ⇒ 2 , 当 |a| > 5 时, Jacobi 迭代收敛. (2) a3 − 14a + 12 > 0, a 所以, 当 a ≥ √ 14 时, A 对称正定, 从而 Gauss-
2
故 Jacobi 迭代, Gauss-Seidel 迭代均收敛. 2.13 提示: ρ(J ) = 0.9 < 1, 故迭代法收敛. 1 0 . 5 0. 5 2.14 提示: 容易验证 A = 0.5 1 0.5 是对称正定的, 故 Gauss-Seidel 迭代收敛, 但 2D − A = 0.5 0.5 1 1 −0.5 −0.5 −0.5 1 −0.5 不正定, 故 Jacobi 迭代发散. −0.5 −0.5 1 0 0 −1 3 2.15 提示: BJ = 1 0 0 . 特征方程 3λ + λ + 2 = 0, 特征值 λ1 = −0.478, λ2,3 = 1 2 0 3 3 故 Jacobi 迭代收敛. −1 −1 , 因为 ρ(BS ) = 1, 故 Gauss-Seidel 迭代发散. −1 −22 11 1 2.16 提示: (1) 将原方程组的系数矩阵调整为: 1 −4 2 , 显然为严格对角占优矩阵, 故 11 −5 −33 = 0