最新(完美版)第八章习题答案_数值分析

第八章习题解答

3、设方程()0f x =有根，且'0()m f x M <≤≤。试证明由迭代格式1()k k k x x f x λ+=- (0,1,2,)k =产生的迭代序列{}0k k x ∞=对任意的初值0(,)x ∈-∞+∞，当20M λ<<时，均收敛于方程的根。

证明：

设()()x x f x ?λ=-，可知()x ?在(,)-∞∞上可导

对于任意给定的λ值，满足条件'0()m f x M <≤≤时

（1）''()1()x f x ?λ=- 则1'()11M x m λ?λ-≤≤-< 又20M

λ<<，M>0 则02M λ<<时，11M λ-<- 所以11'()11M x m λ?λ-<-≤≤-< 若令max{1,1}L M m λλ=--，则可知'()1x L ?≤<

（2）由0()(0)'()(0)'()x

x x dx x ?????ε=+=+? 则()lim 1x x L x ?→∞??≤< ???

所以，存在一个数a ，当x a >时，()x x ?<

同时，()x ?在[,]a a -内有界，即存在0b >使得[,]x a a ?∈-，()x b ?<

我们选取 max{,}c a b =，则对任意x 有0()max{,}x c x ?<

则对给定的任意初值0x ，设0max{,}d c x =

则0[,]x d d ∈-，于是在区间[,]d d -上有()x d ?<

即满足映内性

有（1）、（2）可知，()x ?满足收敛定理

迭代序列0{}k k x ∞=收敛于方程的根

6. 给出计算...222+++=x 的迭代格式，讨论迭代格式的收敛性，并证明2=x

解：构造迭代格式10,1,2,k x k +==???

2k x ≤

令()x ?=x ?∈?时，()x ??∈?

'()

x ?=，当x ?∈?时，1

'()12x ?<<

所以，迭代格式收敛，且收敛于()x x

在??上的根，即x=x=2

数值分析试题及答案汇总

数值分析试题一、填空题（2 0×2′） 1. ?? ????-=? ?????-=32,1223X A 设x =是精确值x *=的近似值，则x 有 2 位有效数字。 2. 若f (x )=x 7－x 3＋1，则f [20,21,22,23,24,25,26,27]= 1 ， f [20,21,22,23,24,25,26,27,28]= 0 。 3. 设，‖A ‖∞＝___5 ____，‖X ‖∞＝__ 3_____， ‖AX ‖∞≤_15_ __。 4. 非线性方程f (x )=0的迭代函数x =?(x )在有解区间满足 |?’(x )| <1 ，则使用该迭代函数的迭代解法一定是局部收敛的。 5. 区间[a ,b ]上的三次样条插值函数S (x )在[a ,b ]上具有直到 2 阶的连续导数。 6. 当插值节点为等距分布时，若所求节点靠近首节点，应该选用等距节点下牛顿差商公式的前插公式，若所求节点靠近尾节点，应该选用等距节点下牛顿差商公式的后插公式；如果要估计结果的舍入误差，应该选用插值公式中的拉格朗日插值公式。 7. 拉格朗日插值公式中f (x i )的系数a i (x )的特点是：=∑=n i i x a 0)( 1 ；所以当系数a i (x )满足 a i (x )>1 ，计算时不会放大f (x i )的误差。 8. 要使 20的近似值的相对误差小于%，至少要取 4 位有效数字。 9. 对任意初始向量X (0)及任意向量g ，线性方程组的迭代公式x (k +1)=Bx (k )+g (k =0,1,…)收敛于方程组的精确解x *的充分必要条件是 ?(B)<1 。 10. 由下列数据所确定的插值多项式的次数最高是 5 。 11. 牛顿下山法的下山条件为 |f(xn+1)|<|f(xn)| 。 12. 线性方程组的松弛迭代法是通过逐渐减少残差r i (i =0,1,…,n )来实现的，其中的残差 r i ＝ (b i -a i1x 1-a i2x 2-…-a in x n )/a ii ，(i =0,1,…,n )。 13. 在非线性方程f (x )=0使用各种切线法迭代求解时，若在迭代区间存在唯一解，且f (x )

数值分析最佳习题(含答案)

第一章绪论姓名学号班级习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。 1 若误差限为5 105.0-?,那么近似数有几位有效数字？（有效数字的计算）解：2*103400.0-?=x ，325* 102 1 1021---?=?≤-x x 故具有3位有效数字。 2 14159.3=π具有4位有效数字的近似值是多少？（有效数字的计算）解：10314159.0?= π，欲使其近似值* π具有4位有效数字，必需 41*1021 -?≤-ππ，3*3102 11021--?+≤≤?-πππ，即14209.314109.3*≤≤π 3 已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ?有几位有效数字？（有效数字的计算）解：3* 1021-?≤ -a a ，2*102 1 -?≤-b b ，而1811.2=+b a ，1766.1=?b a 2123****102 1 10211021)()(---?≤?+?≤-+-≤+-+b b a a b a b a 故b a +至少具有2位有效数字。 2 123*****102 1 0065.01022031.1102978.0)()(---?≤=?+?≤-+-≤-b b a a a b b a ab 故b a ?至少具有2位有效数字。 4 设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差？（误差的计算）解：已知 δ=-* *x x x ，则误差为 δ=-= -* **ln ln x x x x x 则相对误差为 * * ** * * ln ln 1ln ln ln x x x x x x x x δ = -= - 5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5* =，已知 cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v 2π=的绝对误差限与相对误差限。（误差限的计算）解： * 2******2),(),(h h r r r h r r h v r h v -+-≤-ππ 绝对误差限为 π ππ252.051.02052)5,20(),(2=??+????≤-v r h v

数值分析第8章作业

第八章矩阵特征值问题计算 3．用幂法计算下列矩阵的主特征值及对应的特征向量 12732343()341;()463213331a A b A --???? ????=-=-???? ????--???? 当特征值有3位小数稳定代终止。解：套用幂法公式 010,,,1,2,.... max()k k k k k v u v Au u k v -≠== = 取0(1,1,1)0T u =≠,将A 1代入上式，计算结果见下表则1A 的主特征值19.605572λ≈，特征向量1(10.6050.394369)T x ≈- 将2A 代入幂法公式，取0(1,1,1)T u =，计算结果见下表则2A 主住特征值18.869699λ≈，特征向量1(0.604228,1,0.160881)T x ≈- 4．用反幂法求矩阵 621231111A ?? ??=?? ???? 的最接近于6的特征向量。解：本题按带原点平移的反幂法计算。平移向量p=6,则将

021231115B A pI ?? ??=-=-?? ??-?? 进行三分解：PB=LU ，其中 1 002310101511 001,10,02 221004 2701005 5P L U ? ??? ????-??? ??? ??????===-???????????? ?? ?? ??? ??? 然后1(1,1,1)T Uv =，解得 1 111,max()v v u v = 1,,,2,3,.... max()k k k k k k k v Ly PU Uv y U k v -=== = 计算结果如下：

数值分析试卷及答案

二 1 求A的LU分解，并利用分解结果求解由紧凑格式故从而故 2求证：非奇异矩阵不一定有LU分解证明设非奇异，要说明A不一定能做LU分解，只需举出一个反例即可。现考虑矩阵，显然A为非奇异矩阵。若A有LU分解，则故，而，显然不能同时成立。这矛盾说明A不能做LU分解，故只假定A非奇异并不能保证A能做LU分解，只有在A的前阶顺序主子式时才能保证A一定有LU分解。 3用追赶法求解如下的三对角方程组解设有分解由公式其中分别是系数矩阵的主对角线元素及其下边和上边的次对角线元素，故有从而有故，，，故，，， 4设A是任一阶对称正定矩阵，证明是一种向量范数证明（1）因A正定对称，故当时，，而当时，（2）对任何实数，有（3）因A正定，故有分解，则故对任意向量和，总有综上可知，是一种向量范数。 5 设，，已知方程组的精确解为（1）计算条件数；（2）若近似解，计算剩余；（3）利用事后误差估计式计算不等式右端，并与不等式左边比较，此结果说明了什么？解（1）（2）

（3）由事后误差估计式，右端为而左端这表明当A为病态矩阵时，尽管剩余很小，误差估计仍然较大。因此，当A病态时，用大小作为检验解的准确度是不可靠的。 6矩阵第一行乘以一数成为，证明当时，有最小值证明设，则又故从而当时，即时，有最小值，且 7 讨论用雅可比法和高斯-赛德尔法解方程组时的收敛性。如果收敛，比较哪一种方法收敛较快，其中解对雅可比方法，迭代矩阵，故雅可比法收敛。对高斯-赛德尔法，迭代矩阵，故高斯-赛德尔法收敛。因=故高斯-赛德尔法较雅可比法收敛快。 8设，求解方程组，求雅可比迭代法与高斯-赛德尔迭代法收敛的充要条件。解雅可比法的迭代矩阵，故雅可比法收敛的充要条件是。高斯-赛德尔法的迭代矩阵，故高斯-赛德尔法收敛的充要条件是。 9 设求解方程组的雅可比迭代格式为，其中，求证：若，则相应的高斯-赛德尔法收敛。证明由于是雅可比法的迭代矩阵，故又，故，即，故故系数矩阵A按行严格对角占优，从而高斯-赛德尔法收敛。 10设A为对称正定矩阵，考虑迭代格式求证：（1）对任意初始向量，收敛；（2）收敛到的解。证明（1）所给格式可化为这里存在是因为，由A对称正定，，故也对称正定。设迭代矩阵的特征值为，为相应的特征向量，则与做内积，有因正定，故，从而，格式收敛。

数值分析试卷及答案

3用追赶法求解如下的三对角方程组解设有分解由公式其中分别是系数矩阵的主对角线元素及其下边和上边的次对角线元素，故有从而有故，，，故，，，

4设A是任一阶对称正定矩阵，证明是一种向量范数证明（1）因A正定对称，故当时，，而当时，（2）对任何实数，有（3）因A正定，故有分解，则故对任意向量和，总有综上可知，是一种向量范数。 5 设，，已知方程组的精确解为（1）计算条件数；（2）若近似解，计算剩余；（3）利用事后误差估计式计算不等式右端，并与不等式左边比较，此结果说明了什么？解（1）（2）（3）由事后误差估计式，右端为而左端

这表明当A为病态矩阵时，尽管剩余很小，误差估计仍然较大。因此，当A病态时，用大小作为检验解的准确度是不可靠的。 6矩阵第一行乘以一数成为，证明当时，有最小值证明设，则又故从而当时，即时，有最小值，且 7讨论用雅可比法和高斯-赛德尔法解方程组时的收敛性。如果收敛，比较哪一种方法收敛较快，其中解对雅可比方法，迭代矩阵，故雅可比法收敛。对高斯-赛德尔法，迭代矩阵

，故高斯-赛德尔法收敛。因=故高斯-赛德尔法较雅可比法收敛快。 8设，求解方程组，求雅可比迭代法与高斯-赛德尔迭代法收敛的充要条件。解雅可比法的迭代矩阵，故雅可比法收敛的充要条件是。高斯-赛德尔法的迭代矩阵，

数值分析习题与答案

第一章绪论习题一?1．设x＞0，x＊的相对误差为δ,求f(ｘ)=ln x的误差限。解:求lnx的误差极限就是求f(x)=lnx的误差限，由公式(１.2.４）有已知x*的相对误差满足，而 ,故即 2.下列各数都是经过四舍五入得到的近似值,试指出它们有几位有效数字，并给出其误差限与相对误差限。解:直接根据定义和式（１.2.２)(1.２.3）则得?有５位有效数字，其误差限,相对误差限有２位有效数字，有5位有效数字， 3．下列公式如何才比较准确？ (１)?(2）解:要使计算较准确，主要是避免两相近数相减，故应变换所给公式。

（1)?(2） 4.近似数x*=０.０３１0,是 3 位有数数字。 5.计算取,利用 :式计算误差最小。四个选项: 第二、三章插值与函数逼近习题二、三 1. 给定的数值表用线性插值与二次插值计算ln０．5４的近似值并估计误差限. 解：仍可使用n＝1及n＝２的Laｇrange插值或Neｗｔo ｎ插值，并应用误差估计（5.8)。线性插值时,用0．５及0．6两点，用Newton插值??误差限，因,

故? 二次插值时，用0.5，0．6,０.7三点,作二次Ｎewton插值 ?误差限，故? 2. 在－4≤x≤4上给出的等距节点函数表，若用二次插值法求的近似值，要使误差不超过,函数表的步长ｈ应取多少? 解：用误差估计式（５．8), ?令因?得 3. 若,求和.

解：由均差与导数关系 ?于是 4. 若互异,求的值,这里p≤n+1. 解:,由均差对称性可知当有?而当P＝n ＋1时 ?于是得 5. 求证. 解:解:只要按差分定义直接展开得 ? 6. 已知的函数表

郑州大学数值分析重点考察内容及各章习题

《数值分析》重点考察内容及各章作业答案学院：学号：姓名：

重点考察内容基本概念（收敛阶，收敛条件，收敛区域等）, 简单欧拉法。第一章基础掌握：误差的种类，截断误差，舍入误差的来源，有效数字的判断。了解：误差限，算法及要注意的问题。第二章插值掌握：Hermite插值，牛顿插值，差商计算，插值误差估计。了解：Lagrange插值第三章数据拟合掌握：给出几个点求线性拟合曲线。了解：最小二乘原理第四章数值积分微分掌握：梯形公式，Simpson公式，代数精度，Gauss积分，带权Gauss积分公式推导，复化梯形公式推导及算法。了解：数值微分，积分余项第五章直接法掌握：LU分解求线性方程组，运算量了解：Gauss消去法，LDL，追赶法第六章迭代法掌握：Jacobi,Gauss-Seidel迭代格式构造，敛散性分析，向量、矩阵的范数、谱半径了解：SOR迭代第七章Nolinear迭代法掌握：牛顿迭代格式构造，简单迭代法构造、敛散性分析，收敛阶。了解：二分法，弦截法第八章ODE解法掌握：Euler公式构造、收敛阶。了解：梯形Euler公式、收敛阶，改进Euler公式题目类型：填空，计算，证明综合题

第一章误差 1. 科学计算中的误差来源有4个，分别是________，________，________，________。 2. 用Taylor 展开近似计算函数000()()'()()f x f x f x x x ≈+-，这里产生是什么误差？ 3. 0.7499作 3 4 的近似值，是______位有效数字，65.380是舍入得到的近似值，有____几位有效数字，相对误差限为_______. 0.0032581是四舍五入得到的近似值，有_______位有效数字. 4. 改变下列表达式，使计算结果比较精确：（1）11,||1121x x x x --++ （2 ||1x （3） 1cos ,0,|| 1.x x x x -≠ （4）sin sin ,αβαβ-≈ 5. 采用下列各式计算61)时，哪个计算效果最好？并说明理由。（1）（2 ）99-3 ）6 (3-（4 6. 已知近似数*x 有4位有效数字，求其相对误差限。上机实验题： 1、利用Taylor 展开公式计算0! k x k x e k ∞ ==∑，编一段小程序，上机用单精度计算x e 的函数值. 分别取x =1，5，10，20，-1，-5，-10，-15，-20，观察所得结果是否合理，如不合理请分析原因并给出解决方法. 2、已知定积分1 ,0,1,2,,206 n n x I dx n x ==+? ,有如下的递推关系 111 110 0(6)61666 n n n n n x x x x I dx dx I x x n ---+-===++-? ? 可建立两种等价的计算公式 (1) 1016,0.154n n I I I n -= -=取； (2) 12011),0.6n n I nI I n -=-=（取

数值分析试卷及其答案

1、（本题5分）试确定7 22 作为π的近似值具有几位有效数字，并确定其相对误差限。解因为 7 22 =3.142857…=1103142857 .0-? π=3.141592… 所以 312102 11021005.0001264.0722--?=?=<=- π （2分）这里，3,21,0=-=+-=n n m m 由有效数字的定义可知7 22 作为π的近似值具有3位有效数字。（1分）而相对误差限 3102 1 0005.00004138.0001264.07 22-?= <≈= -= π π πε r （2分） 2、（本题6分）用改进平方根法解方程组：??? ?? ??=????? ??????? ??--654131*********x x x ；解设???? ? ??????? ? ?????? ??===????? ??--11111 1 131321112323121 32 132 31 21 l l l d d d l l l LDL A T 由矩阵乘法得： 5 7,21,215 27 ,25,2323121321- ==-== -==l l l d d d （3分）由y D x L b Ly T 1 ,-==解得 T T x y )9 23 ,97,910(,)563, 7,4(== （3分） 3、（本题6分）给定线性方程组???????=++-=+-+=-+-=-+17 7222382311387 510432143213 21431x x x x x x x x x x x x x x 1）写出Jacoib 迭代格式和Gauss-Seidel 迭代格式； 2）考查Jacoib 迭代格式和Gauss-Seidel 迭代格式的敛散性；解 1）Jacoib 迭代格式为

数值分析复习题要答案

第一章 1、ln2＝0.69314718…，精确到 10－3 的近似值是多少？解精确到 10－3＝0.001，即绝对误差限是 e ＝0.05％，故至少要保留小数点后三位才可以。 ln2≈0.693。 2、设115.80,1025.621≈≈x x 均具有5位有效数字，试估计由这些数据计算21x x ， 21x x +的绝对误差限解：记126.1025, 80.115x x == 则有11232411 10, | 102|||2 x x x x --≤?-≤?- 所以 121212121212211122||||||||||||x x x x x x x x x x x x x x x x x x -=-+-+≤-- 3411 80.11610 6.10102522 0.007057-==??+≤?? 1212112243|()|||11 |10100.0005522 |x x x x x x x x --≤≤?+?=+-+-+- 3、一个园柱体的工件,直径d 为10.250.25mm,高h 为40.00 1.00mm,则它的体积V 的近似值、误差和相对误差为多少。解： ()() 22222222 4 314210254000000330064 221025400002510251002436444 3300624362436 0073873833006 , .....; ()()()......, ..().()..% .r d h V d h V mm d h V dh d d h V mm V V V πππππεεεεε= ≈=??===+=???+?==±====第二章： 1、分别利用下面四个点的Lagrange 插值多项式和Newton 插值多项式N 3(x )，计算L 3(0.5)及N 3(-0.5) x －2 －1 0 1 f (x ) －1 1 2

《数值分析》杨大地-标准答案(第八章)

数值分析第8章数值积分与数值微分 8.1 填空题（1）n+1个点的插值型数值积分公式∫f(x)dx b a ≈∑A j n j=0f(x j )的代数精度至少是 n ，最高不超过 2n+1 。【注：第1空，见定理8.1】（2）梯形公式有 1 次代数精度，Simpson 公司有 3 次代数精度。【注：分别见定理8.1,8.3】（3）求积公式∫f(x)dx h 0≈h 2[f (0)+f (h )]+ah 2[f ′(0)?f ′(h)]中的参数a= 1/12 时，才能保证该求积公式的代数精度达到最高，最高代数精度为 3 。解：令f(x)=1,x,x 2带入有， { h 2[1+1]+ah 2[0?0]=h h 2[0+h ]+ah 2[1?1]=12 (h 2)h 2[0+h 2]+ah 2[0?2h ]=13 (h 3) //注：x 的导数=1 解之得，a=1/12，此时求积公式至少具有2次代数精度。 ∴ 积分公式为：∫f(x)dx h 0≈h 2[f (0)+f (h )]+h 2 12[f ′(0)?f ′(h)] 令 f(x)= x 3带入求积公式有：h 2 [0 +h 3]+ h 212 [0?3h 2]=14 (h 4)，与f(x)= x 4的定积分计算值1 4 (h 4)相等，所以，此求积公式至少具有3次代数精度。令f(x)= x 4带入求积公式有，h 2[0+h 4]+h 2 12[0?4h 3]=1 6(h 5)，与f(x)= x 5的定积分计算值1 5(h 5)不相等，所以，此求积公式的最高代数精度为3次代数精度。 8.2 确定下列求积公式的求积系数和求积节点，使其代数精度尽量高，并指出其最高代数精度。解题思路：按照P149 中8.3式进行求解，根据求积公式中未知量n 的数量决定代入多少f(x)，当积分公式代入求积节点x n 的计算结果与定积分的计算结果一致，继续代入求积节点X n+1,，若计算结果与对应的定积分计算结果不一致时，求积公式拥有最高n 次的代数精度。（1）∫f(x)dx 2h 0≈A 0f (0)+A 1f (h )+A 2f(2h) 解：令f(x)=1,x,x 2代入有,【注：本例中需求解A 0、A 1、A 2共3个未知量，故需3个相异求积节点f(x)】 {A 0+A 1+A 2=2h A 1h +A 22h =1 2(2h )2A 1h 2+A 2(2h )2=1 3(2h )3 求解得A 0=13h ，A 1=43h ，A 2=1 3h ， ∴求积公式为：∫f(x)dx 2h 0≈13hf (0)+43hf (h )+1 3 hf(2h) ∵该求积公式对3个相异节点1,x,x 2均有余项E (f )=0， //注：参见P149定理8.1 ∴该求积公式至少具有2次代数精度。令f(x)= x 3，代入求积公式有：4 3hh 3+1 3h (2h )3=4h 4 ∵函数f(x) = x 3的定积分结果为：∫x 3dx 2h 0=1 4(2h )4=4h 4 ，与求积公式计算值相等， ∴该求积公式具有3次代数精度。

数值分析试卷及其答案2

1、（本题5分）试确定7 22作为π的近似值具有几位有效数字，并确定其相对误差限。解因为 7 22=3.142857…=1103142857.0-? π=3.141592… 所以 3 12 10 2 110 21005.0001264.07 22--?= ?= <=- π （2分）这里，3,21,0=-=+-=n n m m 由有效数字的定义可知7 22作为π的近似值具有3位有效数字。（1分）而相对误差限 3 10 2 10005.00004138.0001264.07 22-?= <≈= -= π π πε r （2分） 2、（本题6分）用改进平方根法解方程组：???? ? ??=????? ??????? ??--654131321 112321x x x ；解设???? ? ? ?????? ? ?????? ??===????? ? ?--11 1 11113 1321 11232312132 1 32 31 21 l l l d d d l l l LDL A T 由矩阵乘法得： 5 7,21,21527,25,2323121321- == - == -==l l l d d d （3分）由y D x L b Ly T 1 ,-==解得 T T x y )9 23,97,910( ,)5 63, 7,4(== （3分） 3、（本题6分）给定线性方程组??? ? ? ??=++-=+-+=-+-=-+17722238231138751043214321 321431x x x x x x x x x x x x x x 1）写出Jacoib 迭代格式和Gauss-Seidel 迭代格式； 2）考查Jacoib 迭代格式和Gauss-Seidel 迭代格式的敛散性；解 1）Jacoib 迭代格式为

数值分析习题集及答案[1].(优选)

数值分析习题集（适合课程《数值方法A 》和《数值方法B 》）长沙理工大学第一章绪论 1. 设x >0,x 的相对误差为δ,求ln x 的误差. 2. 设x 的相对误差为2％,求n x 的相对误差. 3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字: *****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====? 4. 利用公式(3.3)求下列各近似值的误差限: ********12412324(),(),()/,i x x x ii x x x iii x x ++其中**** 1234 ,,,x x x x 均为第3题所给的数. 5. 计算球体积要使相对误差限为1％,问度量半径R 时允许的相对误差限是多少? 6. 设028,Y =按递推公式 1n n Y Y -=( n=1,2,…) 计算到100Y .27.982(五位有效数字),试问计算100Y 将有多大误差? 7. 求方程2 5610x x -+=的两个根,使它至少具有四位有效数字27.982). 8. 当N 充分大时,怎样求2 1 1N dx x +∞+?? 9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2 ? 10. 设 212S gt = 假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加,而相对误差却减小. 11. 序列 {}n y 满足递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字), 计算到 10y 时误差有多大?这个计算过程稳定吗? 12. 计算6 1)f =, 1.4≈,利用下列等式计算,哪一个得到的结果最好? 3 -- 13. ()ln(f x x =,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若

数值分析试题及答案

一、单项选择题（每小题3分，共15分） 1. 3.142和3.141分别作为π的近似数具有（）和（）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．4和4 2. 已知求积公式 ()()2 1 121 1()(2)636f x dx f Af f ≈ ++? ，则A ＝（） A ． 16 B ．13 C ．12 D ．2 3 3. 通过点 ()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足（） A ． ()00l x ＝0， ()110l x = B ． ()00l x ＝0， ()111l x = C ．() 00l x ＝1，()111 l x = D ． () 00l x ＝1，()111 l x = 4. 设求方程 ()0 f x =的根的牛顿法收敛，则它具有（）敛速。 A ．超线性 B ．平方 C ．线性 D ．三次 5. 用列主元消元法解线性方程组 1231231 220223332 x x x x x x x x ++=?? ++=??--=? 作第一次消元后得到的第3个方程（）. A ． 232 x x -+= B ．232 1.5 3.5 x x -+= C ． 2323 x x -+= D ． 230.5 1.5 x x -=- 单项选择题答案 1.A 2.D 3.D 4.C 5.B 得分评卷人二、填空题（每小题3分，共15分）

1. 设T X )4,3,2(-=, 则=1||||X ，2||||X = . 2. 一阶均差 ()01,f x x = 3. 已知3n =时，科茨系数()()() 33301213,88C C C ===，那么 () 33C = 4. 因为方程()420 x f x x =-+=在区间 []1,2上满足，所以()0f x =在区间内有根。 5. 取步长0.1h =，用欧拉法解初值问题 ()211y y y x y ?'=+?? ?=? 的计算公式 . 填空题答案 1. 9和29 2. ()() 0101 f x f x x x -- 3. 1 8 4. ()()120 f f < 5. ()12 00.1 1.1,0,1,210.11k k y y k k y +???? ?=+? ?=+???? =??L 得分评卷人三、计算题（每题15分，共60分） 1. 已知函数 21 1y x = +的一组数据：求分段线性插值函数，并计算 () 1.5f 的近似值. 计算题1.答案 1. 解 []0,1x ∈， ()1010.510.50110x x L x x --=?+?=---% []1,2x ∈，()210.50.20.30.81221x x L x x --=?+?=-+--%

数值分析1-4习题及答案

1、 0.1％,要取几位有效数字？（ c ） (a) 2 (b) 3 (c) 4 (d) 5 2、若* 12.30x =是经过四舍五入得到的近似数，则它有几位有效数字？（ c ） (a) 2 (b) 3 (c) 4 (d) 5 3、已知n ＋1个互异节点(x 0,y 0), (x 1,y 1),…, (x n ,y n )和过这些点的拉格朗日插值基函数l k (x )(k =0,1,2,…,n )，且ω(x )=(x －x 0) (x －x 1)… (x －x n )．则n 阶差商f (x 0,x 1,…, x n )= ( ) (a) ∑=n k k k y x l 0 )( (b) ∑='n k k k k x l y 0)( (c) ∑=n k k k x y 0)(ω (d) ∑='n k k k x y 0)(ω 4、已知由数据（0，0），（0.5，y ），（1，3），（2，2）构造出的三次插值多项式 33()6 P x x y 的的系数是，则等于（） (a) -1.5 (b) 1 (c) 5.5 (d) 4.25 5、设(0,1,2,3,4)i x i =为互异结点，()i l x 为拉格朗日插值基函数,则 4 2 () ()i i i x x l x =-∑等于（ a ） (a) 0 (b) 1 (c) 2 (d) 4 4()[,],()()(),()(),( )(), ' () ' (),22 ()()_________________________f x C a b H x a b a b H a f a H b f b H f H a f a f x H x ∈++====-=设是满足下列插值条件的三次多项式：则插值余项 1、是以0，1，2为节点的三次样条函数，则b=-2,c=3 2、已知(1)0,(1)3,(2)4,f f f =-=-=写出()f x 的牛顿插值多项式 2()P x =___2537 623x x +-__,其余项表达式 R(x)=__() (1)(1)(4) [1,4]6 f x x x ξξ'''-+-∈-_______________________ 3、确定求积公式1 0121 ()(1)(0)'(1)f x dx A f A f A f -≈-++? 中的待定参数，使其代数精度尽量高,则A 0=_ 29__________, A 1=__169________, A 2=_29 _______,代数精度=__2_________。

数值计算方法试题集及答案要点

《数值计算方法》复习试题一、填空题： 1、 ?? ??? ?????----=410141014A ，则A 的LU 分解为 A ? ???????? ???=????????? ?? ?。答案： ?? ????????--??????????--=1556141501 4115401411A 2、已知3.1)3(,2.1)2(, 0.1)1(===f f f ，则用辛普生（辛卜生）公式计算求得?≈3 1 _________ )(dx x f ,用三点式求得≈')1(f 。答案：2.367，0.25 3、1)3(,2)2(, 1)1(==-=f f f ，则过这三点的二次插值多项式中2x 的系数为，拉格朗日插值多项式为。答案：-1， )2)(1(21 )3)(1(2)3)(2(21)(2--------= x x x x x x x L 4、近似值*0.231x =关于真值229.0=x 有( 2 )位有效数字； 5、设)(x f 可微,求方程)(x f x =的牛顿迭代格式是( )；答案 )(1)(1n n n n n x f x f x x x '--- =+ 6、对 1)(3++=x x x f ,差商=]3,2,1,0[f ( 1 ),=]4,3,2,1,0[f ( 0 )； 7、计算方法主要研究( 截断 )误差和( 舍入 )误差； 8、用二分法求非线性方程f (x )=0在区间(a ,b )内的根时，二分n 次后的误差限为( 1 2+-n a b )； 9、求解一阶常微分方程初值问题y '= f (x ,y )，y (x 0)=y 0的改进的欧拉公

数值分析整理版试题及答案

例1、已知函数表 x -1 1 2 ()f x -3 0 4 求()f x 的Lagrange 二次插值多项式和Newton 二次插值多项式。解：（1）k x -1 1 2 k y -3 0 4 插值基函数分别为 ()()()()()()()()()() 1200102121()1211126 x x x x x x l x x x x x x x ----= ==-------- ()()()()()()()() ()()021******* ()1211122x x x x x x l x x x x x x x --+-= ==-+---+- ()()()()()()()()()()0122021111 ()1121213 x x x x x x l x x x x x x x --+-= ==-+--+- 故所求二次拉格朗日插值多项式为 () ()()()()()()()()()()2 20 2()11131201241162314 121123537623k k k L x y l x x x x x x x x x x x x x ==?? =-? --+?-+-+?+-????=---++-=+-∑ （2）一阶均差、二阶均差分别为

[]()()[]()()[][][]010********* 011201202303 ,11204 ,412 3 4,,5 2,,126 f x f x f x x x x f x f x f x x x x f x x f x x f x x x x x ---===-----= = =----=== --- k x ()k f x 一阶二阶 -1 -3 1 0 3/ 2 2 4 4 5/6 故所求Newton 二次插值多项式为 ()()[]()[]()() ()()()20010012012,,,35 311126537623P x f x f x x x x f x x x x x x x x x x x x =+-+--=-+ +++-=+- 例2、设2 ()32f x x x =++，[0,1]x ∈，试求()f x 在[0, 1]上关于()1x ρ=，{} span 1,x Φ=的最佳平方逼近多项式。解：若{}span 1,x Φ=，则0()1x ?=，1()x x ?=，且()1x ρ=，这样，有

数值分析第四版习题和答案解析

第四版数值分析习题第一章绪论 1.设x>0,x的相对误差为δ,求的误差. 2.设x的相对误差为2％,求的相对误差. 3.下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指出它们是几位有效数字: 4.利用公式求下列各近似值的误差限: 其中均为第3题所给的数. 5.计算球体积要使相对误差限为1％,问度量半径R时允许的相对误差限是多少 6.设按递推公式 ( n=1,2,…) 计算到.若取≈(五位有效数字),试问计算将有多大误差 7.求方程的两个根,使它至少具有四位有效数字(≈. 8.当N充分大时,怎样求 9.正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝ 10.设假定g是准确的,而对t的测量有±秒的误差,证明当t增加时S的绝对误差增加,而相对误差却减小. 11.序列满足递推关系(n=1,2,…),若(三位有效数字),计算到时误差有多大这个计算过程稳定吗 12.计算,取,利用下列等式计算,哪一个得到的结果最好 13.,求f(30)的值.若开平方用六位函数表,问求对数时误差有多大若改用另一等价公式计算,求对数时误差有多大 14.试用消元法解方程组假定只用三位数计算,问结果是否可靠 15.已知三角形面积其中c为弧度,,且测量a ,b ,c的误差分别为证明面积的误差满足第二章插值法 1.根据定义的范德蒙行列式,令证明是n次多项式,它的根是,且 . 2.当x= 1 , -1 , 2 时, f(x)= 0 , -3 , 4 ,求f(x)的二次插值多项式. 3.

4.给出cos x,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界. 5.设,k=0,1,2,3,求. 6.设为互异节点(j=0,1,…,n),求证: i) ii) 7.设且,求证 8.在上给出的等距节点函数表,若用二次插值求的近似值,要使截断误差不超过,问使用函数表的步长应取多少 9.若,求及. 10.如果是次多项式,记,证明的阶差分是次多项式,并且为正整数). 11.证明. 12.证明 13.证明 14.若有个不同实根,证明 15.证明阶均差有下列性质: i)若,则; ii)若,则. 16.,求及. 17.证明两点三次埃尔米特插值余项是并由此求出分段三次埃尔米特插值的误差限. 18.求一个次数不高于4次的多项式,使它满足并由此求出分段三次埃尔米特插值的误差限. 19.试求出一个最高次数不高于4次的函数多项式,以便使它能够满足以下边界条件,,. 20.设,把分为等分,试构造一个台阶形的零次分段插值函数并证明当时,在上一致收敛到. 21.设,在上取,按等距节点求分段线性插值函数,计算各节点间中点处的与的值,并估计误差. 22.求在上的分段线性插值函数,并估计误差. 23.求在上的分段埃尔米特插值,并估计误差. i) ii) 25.若,是三次样条函数,证明 i); ii)若,式中为插值节点,且,则. 26.编出计算三次样条函数系数及其在插值节点中点的值的程序框图(可用式的表达式). 第三章函数逼近与计算 1.(a)利用区间变换推出区间为的伯恩斯坦多项式. (b)对在上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较. 2.求证: (a)当时,. (b)当时,. 3.在次数不超过6的多项式中,求在的最佳一致逼近多项式.

数值分析试卷及其答案1

1. 已知325413.0,325413*2*1==X X 都有6位有效数字，求绝对误差限。（4分）解：由已知可知6 5.0102 1 ,0,6,10325413.0016*1=?= =-=?=ε绝对误差限n k k X 2分 620*2102 1 ,6,0,10325413.0-?= -=-=?=ε绝对误差限n k k X 2分 2. 已知?? ???=0 01 A 220- ?????440求21,,A A A ∞ （6分）解： {}, 88,4,1max 1==A 1分 {}, 66,6,1max ==∞A 1分 () A A A T max 2λ= 1分 ?????=0 1 A A T 4 2 ???? ? -420?????0 01 2 20 - ???? ?440= ?????0 01 80 ???? ?3200 2分 {}32 32,8,1max )(max ==A A T λ

1分 24322==A 3. 设32)()(a x x f -= （6分） ① 写出f(x)=0解的迭代格式 ② 当a 为何值时，)(1k k x x ?=+ （0,1……）产生的序列{}k x 收敛于 2 解： ①迭代格式为： x a x x x a x a x x a x x x f x f x x k k k k k k k k k k 665)(665)(6)()(')(2 2 32 1 += +=---=-=+? 3 分 ②时迭代收敛即当222,112 10)2(',665)('2<<-<-=-= a a x a x ?? 3分 4. 给定线性方程组，其中：?? ?=13A ?? ?2 2，?? ? ???-=13b 用迭代公式 )()()()1(k k k Ax b x x -+=+α（0,1……）求解，问取什么实数α ，可使迭代收敛（8分）解：所给迭代公式的迭代矩阵为?? ? --???--=-=ααααα21231A I B 2分

数值分析复习题及答案

数值分析复习题一、选择题 1. 3.142和3.141分别作为π的近似数具有（）和（）位有效数字. A ．4和3 B ．3和2 C ．3和4 D ．4和4 2. 已知求积公式()()2 11211()(2)636f x dx f Af f ≈++?，则A ＝（） A ． 16 B ．13 C ．12 D ．2 3 3. 通过点()()0011,,,x y x y 的拉格朗日插值基函数()()01,l x l x 满足（） A ．() 00l x ＝0，()110l x = B ． ()00l x ＝0，()111l x = C ．()00l x ＝1，()111l x = D ． ()00l x ＝1，()111l x = 4. 设求方程()0f x =的根的牛顿法收敛，则它具有（）敛速。 A ．超线性 B ．平方 C ．线性 D ．三次 5. 用列主元消元法解线性方程组1231231220223332x x x x x x x x ++=??++=??--=? 作第一次消元后得到的第3个方程（）. A ．232x x -+= B ．232 1.5 3.5x x -+= C ．2323x x -+= D ．230.5 1.5x x -=- 二、填空 1. 设 2.3149541...x *=，取5位有效数字，则所得的近似值x= . 2.设一阶差商 ()()()21122114,321f x f x f x x x x --= ==---， ()()()322332615,422f x f x f x x x x --===--

则二阶差商 ()123,,______f x x x = 3. 设(2,3,1)T X =--, 则2||||X = ，=∞||||X 。 4．求方程 2 1.250x x --= 的近似根，用迭代公式 1.25x x =+，取初始值 01x =，那么 1______x =。 5．解初始值问题 00'(,)()y f x y y x y =??=?近似解的梯形公式是 1______k y +≈。 6、 1151A ??= ?-??，则A 的谱半径＝。 7、设 2()35, , 0,1,2,... , k f x x x kh k =+== ，则[]12,,n n n f x x x ++= 和[]123,,,n n n n f x x x x +++= 。 8、若线性代数方程组AX=b 的系数矩阵A 为严格对角占优阵，则雅可比迭代和高斯-塞德尔迭代都。 9、解常微分方程初值问题的欧拉（Euler ）方法的局部截断误差为。 10、为了使计算 23123101(1)(1)y x x x =+ +----的乘除法运算次数尽量的少,应将表达式改写成。 11. 设T X )4,3,2(-=, 则=1||||X ，2||||X = . 12. 一阶均差()01,f x x = 13. 已知3n =时，科茨系数()()()33301213,88C C C ===，那么 ()33C = 14. 因为方程()420x f x x =-+=在区间[]1,2上满足，所以()0f x =在区间内有根。 15. 取步长0.1h =，用欧拉法解初值问题()211y y y x y ?'=+???=?的计算公式 . 16.设 * 2.40315x =是真值 2.40194x =的近似值，则*x 有位有效数字。