计算方法习题册2
计算方法各习题及参考答案

计算⽅法各习题及参考答案第⼆章数值分析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 ≈=??≈=?。
计算方法(孙志忠)习题 第二章方程求根

第二章 方程求根
一、填空题
1、3 ()2510f x x x =−−=用二分法求方程在区间[1,3]内的根,进行
一步后根所在区间为( ),进行两步后根所在区间为( )。
2、解方程()0f x =的简单迭代法的迭代函数()x ϕ满足在有根区间内
( ),则在有根区间内任意取一点作为初始值,迭代解都收敛。
3、设)(x f 可微,求方程)(x f x =根的牛顿迭代格式是( )。
二、问答题
1、用牛顿法求方程30()10 1.5f x x x x =−−==在附近的一个根,并精确到6位有效数字。
[]42,
11,2x x =−2.给定方程()试说明该方程在内有根;
(2)构造一个求此根的收敛的迭代法,并说明理由。
3、给定方程()(1)10x f x x e =−−= 1) 分析该方程存在几个根;
2) 用迭代法求出这些根,精确至5位有效数;
3) 说明所用的迭代格式是收敛的。
计算方法练习题与答案

练习题与答案练习题一练习题二练习题三练习题四练习题五练习题六练习题七练习题八练习题答案练习题一一、是非题1.–作为x的近似值一定具有6位有效数字,且其误差限。
()2.对两个不同数的近似数,误差越小,有效数位越多。
()3.一个近似数的有效数位愈多,其相对误差限愈小。
()4.用近似表示cos x产生舍入误差。
( )5.和作为的近似值有效数字位数相同。
( )二、填空题1.为了使计算的乘除法次数尽量少,应将该表达式改写为;2.–是x舍入得到的近似值,它有位有效数字,误差限为,相对误差限为;3.误差的来源是;4.截断误差为;5.设计算法应遵循的原则是。
三、选择题1.–作为x的近似值,它的有效数字位数为( ) 。
(A) 7; (B) 3;(C) 不能确定 (D) 5.2.舍入误差是( )产生的误差。
(A) 只取有限位数 (B) 模型准确值与用数值方法求得的准确值(C) 观察与测量 (D) 数学模型准确值与实际值3.用 1+x近似表示e x所产生的误差是( )误差。
(A). 模型 (B). 观测 (C). 截断 (D). 舍入4.用s*=g t2表示自由落体运动距离与时间的关系式 (g为重力加速度),s t是在时间t内的实际距离,则s t s*是()误差。
(A). 舍入 (B). 观测 (C). 模型 (D). 截断5.作为的近似值,有( )位有效数字。
(A) 3; (B) 4; (C) 5; (D) 6。
四、计算题1.,,分别作为的近似值,各有几位有效数字?2.设计算球体积允许的相对误差限为1%,问测量球直径的相对误差限最大为多少?3.利用等价变换使下列表达式的计算结果比较精确:(1), (2)(3) , (4)4.真空中自由落体运动距离s与时间t的关系式是s=g t2,g为重力加速度。
现设g是精确的,而对t有秒的测量误差,证明:当t增加时,距离的绝对误差增加,而相对误差却减少。
5*. 采用迭代法计算,取k=0,1,…,若是的具有n位有效数字的近似值,求证是的具有2n位有效数字的近似值。
(完整版)四年级数学下册四则运算与简便计算练习题(2)

四则运算与简便计算练习题课堂讲解(前三页)一、填空1.在一个算式里,如果只有加减法,要()计算,如果只有乘除法,要()计算。
2.在一个算式里,如果含有加、减、乘、除四种运算,要先算(),再算()。
3.在一个算式里如果含有小括号,要先算()。
二、口算36 ÷ 3 100 - 62 24 - 8 + 10 75 ×30 371 - 3715 + 24 - 12 200 ÷ 40 84 ÷ 4 48 ÷ 8 × 993÷100= 159+61= 600÷20= 78+222= 405-60= 1000÷8=17×11= 7600÷400= 480÷120= 695-75= 25×17×4=225-99= 640÷40= 468+199= 620-340= 3200÷80=三、比一比,算一算49 + 17 - 25 240 ÷ 40 × 5 300 - 50 × 249 -(17 + 25)240 + 40 × 5 300 - 50 × 20 × 0四、把下面几个分步式改写成综合算式.(1)960÷15=64 64-28=36 综合算式___________________.(2)75×24=1800 9000-1800=7200 综合算式___________ ______(3)810-19=791 791×2=1582 1582+216=1798 综合算式_____________(4)96×5=480 480+20=500 500÷4=125 综合算式_____________五、计算下面各题121 - 111 ÷ 37 (121 - 111 ÷ 37)× 5 280 + 650 ÷ 1345 × 20 × 3 1000 -(280 + 650 ÷ 13)(95 - 19 × 5 )÷74 (120 - 103)× 50 760 ÷ 10 ÷ 38 (270 + 180)÷(30 - 15)707 - 35 × 20 (95 - 19 × 5 )÷7419×96-962÷74 10000-(59+66)×64 5940÷45×(798-616)(270 + 180)÷(30 - 15)(315×40-364)÷7 12520÷8×(121÷11)(2010-906)×(65+15) 707 -35 × 2050+160÷40(58+370)÷(64-45)120-144÷18+35 347+45×2-4160÷52(58+37)÷(64-9×5)95÷(64-45) 178-145÷5×6+42 420+580-64×21÷28 812-700÷(9+31×11)(136+64)×(65-345÷23)85+14×(14+208÷26)121 -111 ÷ 37(120 - 103)× 50(284+16)×(512-8208÷18) 120-36×4÷18+35 (58+37)÷(64-9×5)45 × 20 × 3(121 -111 ÷ 37)× 5 280 +650 ÷ 13 1000 -(280 +650 ÷ 13)760 ÷ 10 ÷ 38 9846-87×(360÷60) 508×345÷(1526-1521)(124-85)×12÷26 (59+21)×(96÷8)325÷13×(266-250) 140-90÷5+678六、面各题,怎样简便就怎样计算。
实用计算方法习题2解答

实用计算方法习题2解答习题22.1 已知下,节点x1,x2及节点处函数值f(x1),f(x2),构造线性插值多项式p1(x). 解:p1(x)?(x?x2)(x1?x2)f(x1)?(x?x1)(x2?x1)f(x2)2.2 设f(xi)=i(i=0,1,2),构造二次式,使p2(x)满足:p2(xi)=f(xi)(i=0,1,2) 解: p2(x)?(x?x1)(x-x2)(x0?x1)(x0?x2)f(x0)?(x?x0)(x?x2)(x1?x0)(x1?x2)f(x1)?(x?x0)(x?x1).(x2?x0)(x2?x1)f(x2)2.3 设节点xi=i(i=0,1,2,3),f(0)=1,f(1)=0,f(2)=-7,f(3)=26,构造次数不超过3次的多项式p3(x),满足p3(xi)=f(xi),i=0,1,2,3 解:(x?x1)(x-x2)(x?x3)(x0?x1)(x0?x2)(x0?x3)(x?x0)(x?x1)(x?x3).(x2?x0)(x2?x1)(x?x3)(x?x0)(x?x2)(x?x3)(x1?x0)(x1?x2)(x?x3)(x-x0)(x-x1)(x-x2)(x3-x0)(x-x1)(x-x2)33p3(x)?f(x0)?f(x1)f(x)3 ?f(x2)?2.4已知函数y=f(x)的观察数据为-2 0 4 5 5 1 -3 1 试构造f(x)的拉格朗日多项式Pn(x),并计算f(-1)。
解先构造基函数xk yk所求三次多项式为2.5已知函数y=f(x)的数据如下表中第1,2列。
计算它的各阶均差。
k 0 1 2 3 4xk 0.40 0.55 0.65 0.80 0.90 f(xk) 0.410 75 0.578 15 0.696 75 0.888 111.201 52 一阶均差二阶均差三阶均差四阶均差解依据均差计算公式,计算结果列下表中。
k 0 1 xk 0.40 0.55 f(xk) 一阶均差二阶均差三阶均差四阶均差 0.410 75 0.578 15 1.116 00 2 3 40.65 0.80 0.90 0.696 75 1.168 00 0.280 00 0.888 11 1.275 73 0.358 93 0.197 33 1.201 52 1.384 10 0.433 48 0.213 00 0.031 34 2.6 设节点xi=i(i=0,1,2,3),f(0)=1,f(1)=0,f(2)=-7,f(3)=26,构造Newton插值多项式。
计算方法-习题第一、二章答案.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即可满足要求。
计算方法习题集及答案第四版
解:
y次迭代公式
k
0
1
2
3
3.5
3.64
3.63
3.63
6. 试证用牛顿法求方程在[1,3]内的根是线性收敛的。 解:
令
y次迭代公式 故
从而 ,时, 故, 故牛顿迭代公式是线性收敛的 7. 应用牛顿法于方程, 导出求立方根的迭代公式,并讨论其收敛
性。
解:
相应的牛顿迭代公式为 迭代函数,, 则,
习题1.1
1. 什么叫数值方法?数值方法的基本思想及其优劣的评价标准如 何?
数值方法是利用计算机求解数学问题近似解的方法 2. 试证明 及
证明: (1)令
即 又 即 ⑵ 设,不妨设, 令 即对任意非零,有 下面证明存在向量,使得, 设,取向量。其中。 显然且任意分量为, 故有即证。 3. 古代数学家祖冲之曾以作为圆周率的近似值,问此近似值具有
解: (1)迭代公式,公式收敛
k
0
1
2
3
0
(2),, 局部收敛 k0 1 2 3
0.25
0.25098 0.25098
456789
1.5 1.322 1.421 1.367 1.397 1.380 1.390 1.384 1.387 1.386
2. 方程在附近有根,把方程写成三种不同的等价形式:
(1),对应迭代公式;
9
10
11
12
13
14
15
16
1.4650 1.46593 1.4653 1.46572 1.46548 1.46563 1.465534 1.465595
迭代公式(2):
k
0
1
2
3
四年级运算定律与简便计算练习题大全 (2)
运算定律与简便计算(一)加减法运算定律1.加法交换律定义:两个加数交换位置.和不变字母表示:a+=a+bb例如:16+23=23+16 546+78=78+5462.加法结合律定义:先把前两个数相加.或者先把后两个数相加.和不变。
字母表示:)+a+b++=b(c)(ca注意:加法结合律有着广泛的应用.如果其中有两个加数的和刚好是整十、整百、整千的话.那么就可以利用加法交换律将原式中的加数进行调换位置.再将这两个加数结合起来先运算。
例1.用简便方法计算下式:(1)63+16+84 (2)76+15+24 (3)140+639+860举一反三:(1)46+67+54 (2)680+485+120 (3)155+657+2453.减法交换律、结合律注:减法交换律、结合律是由加法交换律和结合律衍生出来的。
减法交换律:如果一个数连续减去两个数.那么后面两个减数的位置可以互换。
字母表示:b=---aca-cb例2.简便计算:198-75-98减法结合律:如果一个数连续减去两个数.那么相当于从这个数当中减去后面两个数的和。
字母表示:)-=--a+bb(cca例 3.简便计算:(1)369-45-155 (2)896-580-1204.拆分、凑整法简便计算拆分法:当一个数比整百、整千稍微大一些的时候.我们可以把这个数拆分成整百、整千与一个较小数的和.然后利用加减法的交换、结合律进行简便计算。
例如:103=100+3.1006=1000+6.…凑整法:当一个数比整百、整千稍微小一些的时候.我们可以把这个数写成一个整百、整千的数减去一个较小的数的形式.然后利用加减法的运算定律进行简便计算。
例如:97=100-3.998=1000-2.…注意:拆分凑整法在加、减法中的简便不是很明显.但和乘除法的运算定律结合起来就具有很大的简便了。
例4.计算下式.能简便的进行简便计算:(1)89+106 (2)56+98 (3)658+997随堂练习:计算下式.怎么简便怎么计算(1)730+895+170 (2)820-456+280(3)900-456-244(4)89+997 (5)103-60 (6)458+996(7)876-580+220 (8)997+840+260 (9)956-197-56(二)乘除法运算定律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.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$ 是目标值的精度。