广西大学数值分析课后练习答案第3章
数值分析答案第三章习题
H0
=
125 125 = = 8.04 > u α 2 130.9 + 110.25 15.53
∴
拒绝
，有显著差异。
11.在十块田地上同时试种甲、乙两种品种作 物，根据产量计算得 x = 30.97， = 21.79 ， s x* = 26.7 y （ α = 1%）？假定两种品种作物产量分别服从 正态分布，且方差相等。 解： 1、（1）建立假设
10.48 − 10.5
接受 H 0 ，工作正常。
8.从某种实验物中取出24个样品，测量其发热 * 量，计算得 x = 11958 ，子样标准差 s = 323 ， 问以5%的显著水平是否可认为发热量的期望 值是12100（假定发热量服从正态分布）？ 解：（1）建立假设 H 0 ： µ = 12100 x = 11958 ， s * = 323 n=24， （2）在 H 0 成立的前提下，构造统计量
( x − µ1 ) 2
3.某批砂矿的5个样品中的镍含量经测定为 x(%) 3.25 3.27 3.24 3.26 3.24 设测定值服从正态分布。问在下α = 0.01 能 否接受假设：这批矿砂的(平均)镍含量为3.25。 解：设 ， 未知，计算 2 x ~ N (µ ,σ ) σ 2 x ＝３.252， =0.013。 （1）建立假设
在显著水平 α = 0.05 下，这两种枪弹(平均)速度 有无显著差异？
µ 解： （1）建立假设 H 0 ：1 = µ 2 （2）在 H 0 成立的前提下，构造统计量
u=
x1 − x2 s1 s2 + n1 n2
2 2
~ N (0,1)
数值分析课后习题及答案
第一章 绪论（12） 第二章 插值法（40-42）2、当2,1,1-=x 时，4,3,0)(-=x f ，求)(x f 的二次插值多项式。
[解]372365)1(34)23(21)12)(12()1)(1(4)21)(11()2)(1()3()21)(11()2)(1(0))(())(())(())(())(())(()(2221202102210120120102102-+=-++--=+-+-⨯+------⨯-+-+-+⨯=----+----+----=x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L 。
3、给出x x f ln )(=的数值表用线性插值及二次插值计算54.0ln 的近似值。
X 0.4 0.5 0.6 0.7 0.8 x ln -0.916291 -0.693147 -0.510826 -0.357765 -0.223144[解]若取5.00=x ，6.01=x ，则693147.0)5.0()(00-===f x f y ，510826.0)6.0()(11-===f x f y ，则604752.182321.1)5.0(10826.5)6.0(93147.65.06.05.0510826.06.05.06.0693147.0)(010110101-=---=--⨯---⨯-=--+--=x x x x x x x x x y x x x x y x L ，从而6202186.0604752.19845334.0604752.154.082321.1)54.0(1-=-=-⨯=L 。
若取4.00=x ，5.01=x ，6.02=x ，则916291.0)4.0()(00-===f x f y ，693147.0)5.0()(11-===f x f y ，510826.0)6.0()(22-===f x f y ，则 217097.2068475.404115.2)2.09.0(5413.25)24.0(3147.69)3.01.1(81455.45)5.06.0)(4.06.0()5.0)(4.0()510826.0()6.05.0)(4.05.0()6.0)(4.0()693147.0()6.04.0)(5.04.0()6.0)(5.0(916291.0))(())(())(())(())(())(()(22221202102210120120102102-+-=+--+-⨯++-⨯-=----⨯-+----⨯-+----⨯-=----+----+----=x x x x x x x x x x x x x x x x x x x x x x y x x x x x x x x y x x x x x x x x y x L ，从而61531984.0217097.21969765.259519934.0217097.254.0068475.454.004115.2)54.0(22-=-+-=-⨯+⨯-=L补充题：1、令00=x ，11=x ，写出x e x y -=)(的一次插值多项式)(1x L ，并估计插值余项。
数值分析课后习题部分参考答案解析
数值分析课后习题部分参考答案Chapter 1(P10) 5.求2的近似值x*，使其相对误差不超过0.1%。
解：...2 =1.4 …。
设X有n位有效数字，则|e(x)|乞0.5 10 10 。
*、，0.5 101』从而，| e r (x ) \<1故，若0.5 101』乞0.1%，则满足要求。
解之得，n丄4。
x =1.414。
(P10)7.正方形的边长约100cm，问测量边长时误差应多大，才能保证面积的误差不超过 1 cm2。
解：设边长为a，则a 100 cm。
设测量边长时的绝对误差为e，由误差在数值计算的传播，这时得到的面积的绝对误差有如下估计：：、2 100 e。
按测量要求，|2 100 e|_1解得，|eF0.5 10^。
Chapter 2(P47)5.用三角分解法求下列矩阵的逆矩阵：‘1 1 -1、A= 2 1 0 。
<1-10,解：设A—〔。
分别求如下线性方程组:■n■0 '0 'A。
= 0 ，A P= 1 ，AY = 04 ©先求A的LU分解(利用分解的紧凑格式)Ly = 1 和 U0 =y ，得，031 3Ly = 0 和 U = 丫，得，；丫1 32 31 31 3 所以，A ，=123 32 1 -1 — ——3 3丿广1 2 1 -3)/ 、 花『1'2 5 0 -5 X 22 1 0 14 1X 316 <_3_5115丿3<8>解：平方根法："(1)1 (1)1(2)2(1) -1 (0)2 .(1)1(-1)2(0)-3」「0 0^「1 1 —1'即，L =21 0 ，U = 0-1 2 。
2 h1°0 一3」经直接三角分解法的回代程，分别求解方程组,■1广0、Ly = 0 和 U° = y ，得，a =0 ；厂1」(P47) 6.分别用平方根法和改进平方根法求解方程组:先求系数矩阵 A 的Cholesky 分解(利用分解的紧凑格式)先求系数矩阵上=2」21 二 2d = 1；t 31改进平方根法:A 的形如A = L DL T的分解，其中L = (l j )4 4为单位下三角矩阵，D =diag{d 1,d 2,d 3,d 4}为对角矩阵。
数值分析第三版课本习题及答案
数值分析第三版课本习题及答案第⼀章绪论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.3. 下列各数都是经过四舍五⼊得到的近似数,即误差限不超过最后⼀位的半个单位,试指出它们是⼏位有效数字:*****123451.1021,0.031,385.6,56.430,7 1.0.x x x x x =====?4. 利⽤公式求下列各近似值的误差限:********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 -=…)计算到100Y .(五位有效数字),试问计算100Y 将有多⼤误差?7. 求⽅程25610x x -+=的两个根,使它⾄少具有四位有效数字.8. 当N 充分⼤时,怎样求211Ndx x +∞+?9. 正⽅形的边长⼤约为100㎝,应怎样测量才能使其⾯积误差不超过1㎝210. 设212S gt =假定g 是准确的,⽽对t 的测量有±秒的误差,证明当t 增加时S 的绝对误差增加,⽽相对误差却减⼩. 11. 序列{}n y 满⾜递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字),计算到10y 时误差有多⼤?这个计算过程稳定吗?12.计算61)f =,1.4≈,利⽤下列等式计算,哪⼀个得到的结果最好?13.()ln(f x x =,求f (30)的值.若开平⽅⽤六位函数表,问求对数时误差有多⼤?若改⽤另⼀等价公式ln(ln(x x =-计算,求对数时误差有多⼤?14. 试⽤消元法解⽅程组{101012121010;2.x x x x +=+=假定只⽤三位数计算,问结果是否可靠?15. 已知三⾓形⾯积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c 的误差分别为,,.a b c 证明⾯积的误差s ?满⾜.s a b cs a b c ≤++第⼆章插值法1. 根据定义的范德蒙⾏列式,令200011211121()(,,,,)11n n n n n n n n n x x x V x V x x x x x x x xx x ----==LLL L L L L L L证明()n V x 是n 次多项式,它的根是01,,n x x -L ,且101101()(,,,)()()n n n n V x V x x x x x x x ---=--L L .2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的⼆次插值多项式.3. 给出f (x )=ln x 的数值表⽤线性插值及⼆次插值计算ln 的近似值.x 为互异节点(j =0,1,…,n ),求证:i)0()(0,1,,);nkkj jj x l x x k n =≡=∑Lii)()()1,2,,).nk jj j xx l x k n =-≡0(=∑L7. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8max max a x ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()x f x e =的等距节点函数表,若⽤⼆次插值求x e 的近似值,要使截断误差不超过610-,问使⽤函数表的步长h 应取多少?9. 若2n n y =,求4n y ?及4n y δ.10. 如果()f x 是m 次多项式,记()()()f x f x h f x ?=+-,证明()f x 的k 阶差分()(0)kf x k m ?≤≤是m k -次多项式,并且()0(m lf x l +?=为正整数).11. 证明1()k k k k k k f g f g g f +?=?+?.12. 证明110010.n n kkn n k k k k f gf g f g g f --+==?=--?∑∑13. 证明12n j n j y y y -=?=?-?∑14. 若1011()n nn n f x a a x a x a x --=++++L 有n 个不同实根12,,,n x x x L ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质: i) 若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =L L ;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+L L L .16. 74()31f x x x x =+++,求0172,2,,2f L 及0182,2,,2f L . 17. 证明两点三次埃尔⽶特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔⽶特插值的误差限.18. 求⼀个次数不⾼于4次的多项式()P x ,使它满⾜(0)(1)P P k =-+并由此求出分段三次埃尔⽶特插值的误差限.19. 试求出⼀个最⾼次数不⾼于4次的函数多项式()P x ,以便使它能够满⾜以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造⼀个台阶形的零次分段插值函数()n x ?并证明当n →∞时,()n x ?在[],a b 上⼀致收敛到()f x .21. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()hI x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差.22. 求2()f x x =在[],a b 上的分段线性插值函数()h I x ,并估计误差. 23. 求4()f x x =在[],a b 上的分段埃尔⽶特插值,并估计误差.24. 给定数据表如下:(0.25)(0.53)0.S S "="=25. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明 i)[][][][]222()()()()2()()()bbbba a a a f x dx S x dx f x S x dx S x f x S x dx "-"="-"+""-";ii) 若()()(0,1,,)i i f x S x i n ==L ,式中i x 为插值节点,且01n a x x x b =<<<=L ,则[][][]()()()()()()()()()baS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'?.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可⽤式的表达式).第三章函数逼近与计算1. (a)利⽤区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做⽐较. 2. 求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =.3. 在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳⼀致逼近多项式.4. 假设()f x 在[],a b 上连续,求()f x 的零次最佳⼀致逼近多项式.5. 选取常数a ,使301max x x ax≤≤-达到极⼩,⼜问这个解是否唯⼀?6. 求()sin f x x =在[]0,/2π上的最佳⼀次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳⼀次逼近多项式.8. 如何选取r ,使2()p x x r =+在[]1,1-上与零偏差最⼩?r 是否唯⼀? 9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式.10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*的正交多项式.12. 在[]1,1-上利⽤插值极⼩化求11()f x tg x -=的三次近似最佳逼近多项式.13. 设()xf x e =在[]1,1-上的插值极⼩化近似最佳逼近多项式为()n L x ,若nf L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ?=-----,试将()x ?降低到3次多项式并估计误差. 15. 在[ ]1,1-上利⽤幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]220sin ax b x dxπ+-?为最⼩.并与1题及6题的⼀次逼近多项式误差作⽐较.18. ()f x 、[]1(),g x C a b ∈,定义 ()(,)()();()(,)()()()();b baaa f g f x g x dxb f g f x g x dx f a g a =''=''+??问它们是否构成内积?19. ⽤许⽡兹不等式估计6101x dx x +?的上界,并⽤积分中值定理估计同⼀积分的上下界,并⽐较其结果.20. 选择a ,使下列积分取得最⼩值:112221110010121,,,span x span x x 1?=?=,分别在1?、2?上求出⼀个元素,使得其为[]20,1x C ∈的最佳平⽅逼近,并⽐较其结果.22. ()f x x =在[]1,1-上,求在{}2411,,span x x ?=上的最佳平⽅逼近.23.sin (1)arccos ()n n x u x +=是第⼆类切⽐雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin 2f x x=在[]1,1-上按勒让德多项式及切⽐雪夫多项式展开,求三次最佳平⽅逼近多项式并画出误差图形,再计算均⽅误差.25. 把()arccos f x x =在[]1,1-上展成切⽐雪夫级数.26. ⽤最⼩⼆乘法求⼀个形如2y a bx =+的经验公式,使它与下列数据拟合,并求均⽅误差.28. 在某化学反应⾥,根据实验所得分解物的浓度与时间关系如下:⽤最⼩⼆乘拟合求.29. 编出⽤正交多项式做最⼩⼆乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出⼀张记录{}{}4,3,2,1,0,1,2,3k x =,试⽤改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =L第四章数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量⾼,并指明所构造出的求积公式所具有的代数精度: (1)101()()(0)()hh f x dx A f h A f A f h --≈-++?;(2)21012()()(0)()hh fx dx A f h A f A f h --≈-++?;(3)[]1121()(1)2()3()/3()(0)()/1(0)()hf x dx h f f h ah f f h ≈++'-'?.2. 分别⽤梯形公式和⾟普森公式计算下列积分: (1)120,84xdx n x =+?; (2)1210(1),10x e dx n x --=?;(3)1,4n =?; (4),6n =.3. 直接验证柯特斯公式具有5次代数精度.4. ⽤⾟普森公式求积分1xedx-?并计算误差.5. 推导下列三种矩形求积公式:(1)2()()()()()2baf f x dx b a f a b a 'η=-+-?; (2)2()()()()()2baf f x dx b a f b b a 'η=---?;(3)3()()()()()224baa b f f x dx b a f b a +"η=-+-?.6. 证明梯形公式和⾟普森公式当n →∞时收敛到积分7. ⽤复化梯形公式求积分()baf x dx,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍⼊误差)?8.1x e dx-,要求误差不超过510-.9. 卫星轨道是⼀个椭圆,椭圆周长的计算公式是S a =θ,这⾥a 是椭圆的半长轴,c 是地球中⼼与轨道中⼼(椭圆中⼼)的距离,记h 为近地点距离,H 为远地点距离,6371R =公⾥为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第⼀颗⼈造卫星近地点距离439h =公⾥,远地点距离2384H =公⾥,试求卫星轨道的周长.10. 证明等式3524sin3!5!n nn n ππππ=-+-L试依据sin(/)(3,6,12)n n n π=的值,⽤外推算法求π的近似值.11. ⽤下列⽅法计算积分31dyy ?并⽐较结果.(1) 龙贝格⽅法;(2) 三点及五点⾼斯公式;(3) 将积分区间分为四等分,⽤复化两点⾼斯公式.12. ⽤三点公式和五点公式分别求21()(1)f x x =第五章常微分⽅程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉⽅法和改进的尤拉⽅法的近似解的表达式，并与准确解bx ax y +=221相⽐较。
数值分析课后习题答案
7、计算的近似值，取。
利用以下四种计算格式，试问哪一种算法误差最小。
〔1〕〔2〕〔3〕〔4〕解：计算各项的条件数由计算知，第一种算法误差最小。
解：在计算机上计算该级数的是一个收敛的级数。
因为随着的增大，会出现大数吃小数的现象。
9、通过分析浮点数集合F=〔10，3，-2，2〕在数轴上的分布讨论一般浮点数集的分布情况。
10、试导出计算积分的递推计算公式，用此递推公式计算积分的近似值并分析计算误差，计算取三位有效数字。
解：此算法是数值稳定的。
第二章习题解答1.〔1〕 R n×n中的子集“上三角阵〞和“正交矩阵〞对矩阵乘法是封闭的。
〔2〕R n×n中的子集“正交矩阵〞，“非奇异的对称阵〞和“单位上〔下〕三角阵〞对矩阵求逆是封闭的。
设A是ｎ×ｎ的正交矩阵。
证明A-1也是ｎ×ｎ的正交矩阵。
证明：〔2〕A是ｎ×ｎ的正交矩阵∴A A-1 =A-1A=E 故〔A-1〕-1=A∴A-1〔A-1〕-1=〔A-1〕-1A-1 =E 故A-1也是ｎ×ｎ的正交矩阵。
设A是非奇异的对称阵，证A-1也是非奇异的对称阵。
A非奇异∴A可逆且A-1非奇异又A T=A ∴〔A-1〕T=〔A T〕-1=A-1故A-1也是非奇异的对称阵设A是单位上〔下〕三角阵。
证A-1也是单位上〔下〕三角阵。
证明：A是单位上三角阵，故|A|=1，∴A可逆，即A-1存在，记为〔b ij〕n×n由A A-1 =E，那么〔其中 j＞i时，〕故b nn=1, b ni=0 (n≠j)类似可得，b ii=1 (j=1…n) b jk=0 (k＞j)即A-1是单位上三角阵综上所述可得。
R n×n中的子集“正交矩阵〞，“非奇异的对称阵〞和“单位上〔下〕三角阵〞对矩阵求逆是封闭的。
2、试求齐次线行方程组Ax=0的根底解系。
A=解：A=～～～故齐次线行方程组Ax=0的根底解系为，3.求以下矩阵的特征值和特征向量。
数值分析习题(含标准答案)
]第一章 绪论姓名 学号 班级习题主要考察点：有效数字的计算、计算方法的比较选择、误差和误差限的计算。
1若误差限为5105.0-⨯,那么近似数有几位有效数字（有效数字的计算） 解：2*103400.0-⨯=x ，325*10211021---⨯=⨯≤-x x 故具有3位有效数字。
2 14159.3=π具有4位有效数字的近似值是多少（有效数字的计算） 解：10314159.0⨯= π，欲使其近似值*π具有4位有效数字，必需！41*1021-⨯≤-ππ，3*310211021--⨯+≤≤⨯-πππ，即14209.314109.3*≤≤π即取（ , ）之间的任意数，都具有4位有效数字。
3已知2031.1=a ，978.0=b 是经过四舍五入后得到的近似值，问b a +，b a ⨯有几位有效数字（有效数字的计算）解：3*1021-⨯≤-aa ，2*1021-⨯≤-b b ，而1811.2=+b a ，1766.1=⨯b a 2123****102110211021)()(---⨯≤⨯+⨯≤-+-≤+-+b b a a b a b a故b a +至少具有2位有效数字。
2123*****10210065.01022031.1102978.0)()(---⨯≤=⨯+⨯≤-+-≤-b b a a a b b a ab 故b a ⨯至少具有2位有效数字。
4设0>x ，x 的相对误差为δ，求x ln 的误差和相对误差（误差的计算）~解：已知δ=-**xx x ，则误差为 δ=-=-***ln ln xx x x x则相对误差为******ln ln 1ln ln ln xxx x xxx x δ=-=-5测得某圆柱体高度h 的值为cm h 20*=，底面半径r 的值为cm r 5*=，已知cm h h 2.0||*≤-，cm r r 1.0||*≤-，求圆柱体体积h r v2π=的绝对误差限与相对误差限。
（误差限的计算）解：*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相对误差限为%420120525)5,20()5,20(),(2==⋅⋅≤-ππv v r h v 6设x 的相对误差为%a ,求nx y =的相对误差。
数值分析第三版课本习题及答案
第一章 绪 论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.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 =按递推公式11783100n n Y Y -=-( n=1,2,…)计算到100Y .若取783≈27.982(五位有效数字),试问计算100Y 将有多大误差?7. 求方程25610x x -+=的两个根,使它至少具有四位有效数字(783≈27.982).8. 当N 充分大时,怎样求211Ndx x +∞+⎰?9. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝2?10. 设212S gt =假定g 是准确的,而对t 的测量有±0.1秒的误差,证明当t 增加时S 的绝对误差增加,而相对误差却减小. 11. 序列{}n y 满足递推关系1101n n y y -=-(n=1,2,…),若02 1.41y =≈(三位有效数字),计算到10y 时误差有多大?这个计算过程稳定吗?12. 计算6(21)f =-,取2 1.4≈,利用下列等式计算,哪一个得到的结果最好?36311,(322),,9970 2.(21)(322)--++13. 2()ln(1)f x x x =--,求f (30)的值.若开平方用六位函数表,问求对数时误差有多大?若改用另一等价公式22ln(1)ln(1)x x x x --=-++计算,求对数时误差有多大?14. 试用消元法解方程组{101012121010;2.x x x x +=+=假定只用三位数计算,问结果是否可靠?15. 已知三角形面积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c 的误差分别为,,.a b c ∆∆∆证明面积的误差s ∆满足.s a b cs a b c ∆∆∆∆≤++第二章 插值法1. 根据(2.2)定义的范德蒙行列式,令200011211121()(,,,,)11n n n n n n n n n x x x V x V x x x x x x x xx x ----==证明()n V x 是n 次多项式,它的根是01,,n x x - ,且101101()(,,,)()()n n n n V x V x x x x x x x ---=-- .2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的二次插值多项式.3. 给出f (x )=ln x 的数值表用线性插值及二次插值计算ln 0.54 的近似值.x 0.4 0.5 0.6 0.7 0.8 ln x -0.916291-0.693147-0.510826-0.357765-0.2231444. 给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究用线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6. 设jx 为互异节点(j =0,1,…,n ),求证:i) 0()(0,1,,);nk kj j j x l x x k n =≡=∑ii)()()1,2,,).nk jj j xx l x k n =-≡0(=∑7. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8max max a x ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()xf x e =的等距节点函数表,若用二次插值求xe 的近似值,要使截断误差不超过610-,问使用函数表的步长h 应取多少?9. 若2n n y =,求4n y ∆及4n y δ.10. 如果()f x 是m 次多项式,记()()()f x f x h f x ∆=+-,证明()f x 的k 阶差分()(0)kf x k m ∆≤≤是m k -次多项式,并且()0(m lf x l +∆=为正整数).11. 证明1()k k k k k k f g f g g f +∆=∆+∆.12. 证明110010.n n kkn n k k k k f gf g f g g f --+==∆=--∆∑∑13. 证明1200.n j n j y y y -=∆=∆-∆∑14. 若1011()n nn n f x a a x a x a x --=++++ 有n 个不同实根12,,,n x x x ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质: i)若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x = ;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+ .16. 74()31f x x x x =+++,求0172,2,,2f ⎡⎤⎣⎦ 及0182,2,,2f ⎡⎤⎣⎦ . 17. 证明两点三次埃尔米特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔米特插值的误差限.18. 求一个次数不高于4次的多项式()P x ,使它满足(0)(1)P P k =-+并由此求出分段三次埃尔米特插值的误差限.19. 试求出一个最高次数不高于4次的函数多项式()P x ,以便使它能够满足以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造一个台阶形的零次分段插值函数()n x ϕ并证明当n →∞时,()n x ϕ在[],a b 上一致收敛到()f x .21. 设2()1/(1)f x x =+,在55x -≤≤上取10n =,按等距节点求分段线性插值函数()h I x ,计算各节点间中点处的()h I x 与()f x 的值,并估计误差.22. 求2()f x x =在[],a b 上的分段线性插值函数()h I x ,并估计误差.23. 求4()f x x =在[],a b 上的分段埃尔米特插值,并估计误差. 24. 给定数据表如下:j x 0.25 0.30 0.39 0.45 0.53 j y0.50000.54770.62450.67080.7280试求三次样条插值()S x 并满足条件i) (0.25) 1.0000,(0.53)0.6868;S S '='= ii)(0.25)(0.53)0.S S "="=25. 若[]2(),f x C a b ∈,()S x 是三次样条函数,证明 i)[][][][]222()()()()2()()()bbbbaaaaf x dx S x dx f x S x dx S x f x S x dx"-"="-"+""-"⎰⎰⎰⎰;ii) 若()()(0,1,,)i i f x S x i n == ,式中i x 为插值节点,且01n a x x x b =<<<= ,则[][][]()()()()()()()()()baS x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'⎰.26. 编出计算三次样条函数()S x 系数及其在插值节点中点的值的程序框图(()S x 可用(8.7)式的表达式).第三章 函数逼近与计算1. (a)利用区间变换推出区间为[],a b 的伯恩斯坦多项式.(b)对()sin f x x =在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做比较. 2. 求证:(a)当()m f x M ≤≤时,(,)n m B f x M ≤≤. (b)当()f x x =时,(,)n B f x x =.3. 在次数不超过6的多项式中,求()sin 4f x x =在[]0,2π的最佳一致逼近多项式.4. 假设()f x 在[],a b 上连续,求()f x 的零次最佳一致逼近多项式.5. 选取常数a ,使301max x x ax≤≤-达到极小,又问这个解是否唯一?6. 求()sin f x x =在[]0,/2π上的最佳一次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳一次逼近多项式.8. 如何选取r ,使2()p x x r =+在[]1,1-上与零偏差最小?r 是否唯一? 9. 设43()31f x x x =+-,在[]0,1上求三次最佳逼近多项式.10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x .11. 试证{}*()nTx 是在[]0,1上带权21x x ρ=-的正交多项式.12. 在[]1,1-上利用插值极小化求11()f x tg x -=的三次近似最佳逼近多项式.13. 设()x f x e =在[]1,1-上的插值极小化近似最佳逼近多项式为()n L x ,若nf L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ϕ=-----,试将()x ϕ降低到3次多项式并估计误差. 15. 在[]1,1-上利用幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过0.005.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n nF x H ∈也是奇(偶)函数.17. 求a 、b 使[]220sin ax b x dx π+-⎰为最小.并与1题及6题的一次逼近多项式误差作比较.18. ()f x 、[]1(),g x C a b ∈,定义 ()(,)()();()(,)()()()();b baaa f g f x g x dxb f g f x g x dx f a g a =''=''+⎰⎰问它们是否构成内积?19. 用许瓦兹不等式(4.5)估计6101x dx x +⎰的上界,并用积分中值定理估计同一积分的上下界,并比较其结果.20. 选择a ,使下列积分取得最小值:1122211(),x ax dx x ax dx----⎰⎰.21. 设空间{}{}10010121,,,span x span x x 1ϕ=ϕ=,分别在1ϕ、2ϕ上求出一个元素,使得其为[]20,1x C ∈的最佳平方逼近,并比较其结果.22. ()f x x =在[]1,1-上,求在{}2411,,span x x ϕ=上的最佳平方逼近.23.[]2sin (1)arccos ()1n n x u x x +=-是第二类切比雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin2f x x =在[]1,1-上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方逼近多项式并画出误差图形,再计算均方误差.25. 把()arccos f x x =在[]1,1-上展成切比雪夫级数.26. 用最小二乘法求一个形如2y a bx =+的经验公式,使它与下列数据拟合,并求均方误差.i x 19 25 31 38 44 i y19.032.349.073.397.827. 观测物体的直线运动,得出以下数据:时间t (秒) 0 0.9 1.9 3.0 3.9 5.0 距离s (米) 010305080110求运动方程.28. 在某化学反应里,根据实验所得分解物的浓度与时间关系如下:时间 0 5 10 15 20 25 30 35 40 45 50 55 浓度0 1.272.162.863.443.874.154.374.514.584.624.64用最小二乘拟合求()y f t =.29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出一张记录{}{}4,3,2,1,0,1,2,3k x =,试用改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章 数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所具有的代数精度: (1)101()()(0)()hh f x dx A f h A f A f h --≈-++⎰; (2)21012()()(0)()hh f x dx A f h A f A f h --≈-++⎰;(3)[]1121()(1)2()3()/3f x dx f f x f x -≈-++⎰;(4)[][]20()(0)()/1(0)()hf x dx h f f h ah f f h ≈++'-'⎰.2. 分别用梯形公式和辛普森公式计算下列积分:(1)120,84xdx n x =+⎰; (2)1210(1),10x e dx n x --=⎰;(3)91,4xdx n =⎰; (4)260sin ,6dx n π-ϕ=⎰.3. 直接验证柯特斯公式(2.4)具有5次代数精度.4. 用辛普森公式求积分1x e dx-⎰并计算误差.5. 推导下列三种矩形求积公式:(1)2()()()()()2ba f f x dxb a f a b a 'η=-+-⎰; (2)2()()()()()2baf f x dx b a f b b a 'η=---⎰;(3)3()()()()()224baa b f f x dx b a f b a +"η=-+-⎰.6. 证明梯形公式(2.9)和辛普森公式(2.11)当n →∞时收敛到积分()baf x dx⎰.7. 用复化梯形公式求积分()baf x dx⎰,问要将积分区间[],a b 分成多少等分,才能保证误差不超过ε(设不计舍入误差)?8. 用龙贝格方法计算积分12x e dxπ-⎰,要求误差不超过510-.9. 卫星轨道是一个椭圆,椭圆周长的计算公式是22201()sin cS a d a π=-θθ⎰,这里a 是椭圆的半长轴,c是地球中心与轨道中心(椭圆中心)的距离,记h 为近地点距离,H 为远地点距离,6371R =公里为地球半径,则(2)/2,()/2a R H h c H h =++=-.我国第一颗人造卫星近地点距离439h =公里,远地点距离2384H =公里,试求卫星轨道的周长.10. 证明等式3524sin3!5!n nn n ππππ=-+-试依据sin(/)(3,6,12)n n n π=的值,用外推算法求π的近似值.11. 用下列方法计算积分31dyy ⎰并比较结果.(1) 龙贝格方法;(2) 三点及五点高斯公式;(3) 将积分区间分为四等分,用复化两点高斯公式.12. 用三点公式和五点公式分别求21()(1)f x x =+在x =1.0,1.1和1.2处的导数值,并估计误差.()f x 的值由下表给出:x1.0 1.1 1.2 1.3 1.4 ()f x0.25000.22680.20660.18900.1736第五章 常微分方程数值解法1. 就初值问题0)0(,=+='y b ax y 分别导出尤拉方法和改进的尤拉方法的近似解的表达式，并与准确解bx ax y +=221相比较。
数值分析第四版习题及答案
第四版数值分析习题第一章绪论设x>O,x 的相对误差为S ，求In x 的误差. 设x 的相对误差为2%,求x n 的相对误差. 下列各数都是经过四舍五入得到的近似数，即误差限不超过最后一位的半个单位 ，试指出它们是几位有效数字： x = 1.1021, x^ = 0.031, x^ = 385.6, x^ = 56.430, x^ = 7 1.0.利用公式(3.3)求下列各近似值的误差限：(i)x *+x ；+x 4,(ii)x *x ；x ；,(iii )x ；/x ；,其中 x ；,x ；,x 3,x ；均为第 3题所给的数.计算球体积要使相对误差限为 1%，问度量半径R 时允许的相对误差限是多少 ？设\)=28,按递推公式AY n =Y n d- _ .783100( n=1,2,…)计算到Y 00.若取7783衣27.982(五位有效数字)，试问计算^00将有多大误差？ 求方程X 2 -56X • 1 =0的两个根，使它至少具有四位有效数字 (■ 783沁27.982).\ ------ d x 当N 充分大时，怎样求N 1 x? 正方形的边长大约为 100 cm ,应怎样测量才能使其面积误差不超过 s *2设 2 假定g 是准确的，而对t 的测量有土 0.1秒的误差，证明当t 增加时s 的绝对 误差增加，而相对误差却减小. 序列｛yn｝满足递推关系y n _ 10y n _ 1(n=1,2,…)，若y0 _ X 2 1.41 (三位有效数字),计算到y 10时误差有多大？这个计算过程稳定吗？计算f = c- 2 一1)6,取' 2 : 1.4,利用下列等式计算，哪一个得到的结果最好？f (x) =1 n (x X -1)，求 f(30)的值.若开平方用六位函数表，问求对数时误差有多大改用另一等价公式ln(x_ Jx 2 T) = -ln(x +Jx 2 +1)计算,求对数时误差有多大？1. 2. 3. 4.5. 6.7.8.9.10.11.12.13.21 cm1 (、2 1)61 (32 . 2)3,99 -70、2.?若根据(2.2)定义的范德蒙行列式，令证明V n (x)是n 次多项式,它的根是X 0^L ，X nJ ,且当x= 1 , -1 , 2时，f(x)= 0 , -3,4 ,求f(x)的二次插值多项式.给出cos x,0 ° < x 90。
数值分析课后题
故
有
A B C 4 2 1dx 2 2 A( ) B 0 C (C A) 0 2xdx A( ) 2 B 0 2 C 2 ( A C ) 2 16 2 x 2 dx 3 2 2 A( ) 3 B 0 3 C 3 (C A) 3 0 x 3 dx 2 2 64 4 4 4 4 4 A( ) B 0 C ( A C ) 5 2x dx 2 A( ) 5 B 0 5 C 5 (C A) 5 0 x 6 dx 2
1 2 1 1 1 2 [ ] x 4 dx 3 4 4 3 5 1
得 数 值 积 分 公 式
2
2
f ( x)dx Af ( ) Bf (0) Cf ( ) 为 Gauss 型公式。
2 2
[解]要使数值积分公式 f ( x)dx Af ( ) Bf (0) Cf ( ) 为 Gauss 型公式，则 其具有 2n 1 5 次代数精度。依次将 f ( x) 1, x, x 2 , x 3 , x 4 , x 5 代入都应精确成
1
0
f ( x)dx Af (0) Bf ( x1 ) Cf (1) 具有尽可能高的代数精度。此时的代数精度是
多少？它是否是 Gauss 型公式？ [解]依次将 f ( x) 1, x, x 2 , x 3 代入求积公式，得到：
A B C 1 11dx A B C 1 0 1 1 Bx1 C 1 A 0 Bx C 1 Bx C xdx 1 1 2 2 0 ， 即 1 ， 解 得 2 1 1 Bx1 C 3 A 0 2 Bx12 C 12 Bx12 C x 2 dx 3 0 Bx3 C 1 1 1 3 3 3 3 3 1 A 0 Bx C 1 Bx C x dx 1 1 4 4 0 1 x1 2 B 2 3 ，从而求积公式为 1 f ( x)dx 1 f (0) 2 f ( 1 ) 1 f (1) ，令 f ( x) x 4 代 0 6 3 2 6 C 1 6 1 A 6 入得到：
数值分析课后习题答案
数值分析课后习题答案数值分析课后习题答案数值分析是一门应用数学的学科，主要研究用数值方法解决数学问题的理论和方法。
在学习数值分析课程时，习题是非常重要的一部分，通过解答习题可以加深对数值方法的理解和掌握。
下面将为大家提供一些数值分析课后习题的答案，希望能对大家的学习有所帮助。
1. 插值法是数值分析中常用的一种数值逼近方法。
给定一组已知数据点，我们希望通过插值方法找到一个函数，使得该函数在已知数据点上的取值与给定数据点的值尽可能接近。
常见的插值方法有拉格朗日插值法和牛顿插值法。
下面是一个使用拉格朗日插值法求解的习题：已知函数f(x)=sin(x)，求在区间[0, π/2]上的插值多项式P(x)，使得P(0)=0，P(π/2)=1。
解答：根据拉格朗日插值法的原理，我们需要构造一个满足条件的插值多项式。
首先，我们需要确定插值节点。
根据题目要求，我们取两个插值节点：x0=0，x1=π/2。
然后，我们需要确定插值多项式的系数。
设插值多项式为P(x)=a0+a1x，代入已知条件可得到两个方程：P(0)=a0=0P(π/2)=a0+a1(π/2)=1解方程组可得，a0=0，a1=2/π。
因此，插值多项式为P(x)=2x/π。
2. 数值积分是数值分析中的另一个重要内容。
它主要研究如何用数值方法计算函数的定积分。
常见的数值积分方法有梯形法则、辛普森法则等。
下面是一个使用梯形法则求解的习题：计算定积分∫[0, 1] e^(-x^2) dx。
解答：根据梯形法则的原理，我们可以将定积分转化为离散的求和问题。
首先，我们将积分区间[0, 1]等分为n个小区间，每个小区间的宽度为h=1/n。
然后，我们在每个小区间的两个端点上计算函数值，并将其加权求和。
根据梯形法则的公式，我们可以得到近似解为：∫[0, 1] e^(-x^2) dx ≈ h/2 * (f(0) + 2f(x1) + 2f(x2) + ... + 2f(x(n-1)) + f(1))其中，f(x)表示函数e^(-x^2)在点x处的取值。
数值分析课后习题答案
x2 6.6667x2 8.205
再解
1
15 56
x31.785,7得 x35.769
1 25069x4 0.47847x4 1.4872
1 x5 5.3718 x5 5.3718
2-10.证明下列不等式:
(1)x-yx-z+z-y; (2)|x-y|x-y;
证明 (1)x-y=(x-z)+(z-y)x-z+z-y
b.用Gauss消元法
102 x y 1 x y 2
回代得解: y=1, x=0.
102 x Байду номын сангаасy 1
100y 100
再用列主元Gauss消元法
102 x y 1 x y 2
回代得解: y=1, x=1.
x y
y 1
2
2-8.用追赶法求解方程组:
4 1
x1 100
1 4 1
x2 0
3-8.判定求解下列方程组的SOR方法的收敛性.
2 1 0 0 x1 1
1
0 0
2 1 0
1 2 1
0 12
x2 x3 x4
0 00
解 直接可验证系数矩阵A是负定矩阵,所以-A是对称
1-3.为了使101/2的相对误差小于0.01%,试问应取几位 有效数字?
解 因为101/2=3.162…=0.3162…10,若具有n位有效 数字,则其绝对误差限为0.5 101-n ,于是有
r=0.5101-n/3.162…<0.5101-n/3<0.01% 因此只需n=5.即取101/2=3.1623
1 2
0
12 1,
1 2
1 2
0
12
《数值分析》杨大地 答案(第三章)
数值分析第3章3.1 填空题（1）当A 具有严格对角优势或具有对角线优势且矩阵不可约时，线性方程组b A x =用Jacobi 迭代法和Gauss-Seidel 迭代法均收敛；（2）当线性方程组的系数矩阵A 对称正定时，Gauss-Seidel 迭代收敛；（3）线性方程组迭代法收敛的充分必要条件是迭代矩阵的谱半径小于1；SOR 法收敛的必要条件是0<w<2；（4）用迭代法求解线性方程组，若)1(),(≥=q B q ρ时不收敛，q 接近0时收敛较快，q 接近1时收敛较慢； （5）A=⎪⎪⎭⎫⎝⎛2111；J B =⎪⎪⎭⎫ ⎝⎛--02110；s B =⎪⎪⎭⎫⎝⎛-21010；)22()(=J B ρ；)21()(=s B ρ 解：∵ A=⎪⎪⎭⎫⎝⎛2111，∴ L=⎪⎪⎭⎫⎝⎛0100， D=⎪⎪⎭⎫ ⎝⎛2001， U=⎪⎪⎭⎫⎝⎛0010 //P49D −1=1ad−bc d −b−c a∴⎪⎪⎭⎫⎝⎛--=⎪⎪⎭⎫ ⎝⎛⨯⎪⎪⎭⎫⎝⎛-=+-=-02110011021001)(1U L D B J //P49 B J 指Jacobi 迭代法的矩阵形式 又 ∵ )(U LD +=⎪⎪⎭⎫⎝⎛00211001−−−→−=-=22/12)2122)1r r r r r ⎪⎪⎭⎫ ⎝⎛-210101001=))((1U L D I -+//P50 B S 指Seidel 迭代法的矩阵形式，U L D B S 1)(-+-= ∴ ⎪⎪⎭⎫⎝⎛-=21010s B ∴ 22)(=J B ρ21)(=s B ρ //谱半径ρ A =max 1≤i≤n λi3.2 用Jacobi 迭代法和Gauss-Seidel 迭代法求解方程组，各分量第三位稳定即可。
（1） ⎝⎛012121⎪⎪⎪⎭⎫210⎪⎪⎪⎭⎫ ⎝⎛321x x x =⎪⎪⎪⎭⎫⎝⎛-453 解：利用Jacobi 迭代法，有⎪⎪⎪⎩⎪⎪⎪⎨⎧-=---=-=+++k k k k k kk x x x x x x x 21331122112122521212123，令0x =()T 0,0,0，代入有：//P48k0 1 2 3 4 5 6 7 kx 1 0 1.5 2.75 3.625 4.25 4.6875 5 5.21875 k x 20 -2.5 -4.25 -5.5 -6.375 -7 -7.4375 -7.75 k x 323.254.1254.755.18755.55.71875k89101112131415k x 1 5.375 5.484375 5.5625 5.6171875 5.65625 5.68359375 5.703125 5.716796875 k x 2 -7.96875 -8.125 -8.234375 -8.3125 -8.3671875 -8.40625 -8.43359375 -8.453125 k x 35.8755.9843756.06256.11718756.156256.183593756.2031256.216796875k16 17 18 19 20 21 22 23 k x 1 5.726563 5.733398 5.738281 5.741699 5.744141 5.745850 5.747070 5.747925 k x 2 -8.466797 -8.476563 -8.483398 -8.488281 -8.491699 -8.494141 -8.495850 -8.497070 k x 36.2265636.2333986.2382816.2416996.2441416.2458506.2470706.247925由题，用Jacobi 迭代法进行第21次迭代后，前三位有效数字稳定，此时x=()T25.6,50.8,75.5-利用Gauss-Seidel 迭代法，得出⎪⎪⎪⎩⎪⎪⎪⎨⎧-=---=-=+++++1213211122112122521212123k k k k k kk x x x x x x x ，令0x =()T 0,0,0代入有//P50k0 1 2 3 4 5 6 7 k x 10 1.5 3.125 4.4375 5.09375 5.421875 5.585938 5.667969 kx 2 0 -3.25 -5.875 -7.1875 -7.84375 -8.171875 -8.335938 -8.417969 k x 33.6254.93755.593755.9218756.0859386.1679696.208984k8 9 10 11 12 13 kx 15.708984 5.729492 5.739746 5.744873 5.747437 5.748718 k x 2 -8.458984 -8.479492 -8.489746 -8.494873 -8.497437 -8.498718 k x 36.2294926.2397466.2448736.2474376.2487186.249359可见Gauss-Seidel 迭代法进行至12次迭代后稳定，x=()T 25.6,50.8,75.5-3.4 下面一些方程组的系数阵，试判断它们对Jacobi 迭代法，Gauss-Seidel 迭代法的收敛性。
(完整版)数值分析课后习题答案
第一章绪论习题一1.设x>0,x*的相对误差为δ，求f(x)=ln x的误差限。
解：求lnx的误差极限就是求f(x)=lnx的误差限，由公式(1.2.4)有已知x*的相对误差满足，而，故即2.下列各数都是经过四舍五入得到的近似值，试指出它们有几位有效数字，并给出其误差限与相对误差限。
解：直接根据定义和式(1.2.2)(1.2.3)则得有5位有效数字，其误差限，相对误差限有2位有效数字，有5位有效数字，3.下列公式如何才比较准确？（1）（2）解：要使计算较准确，主要是避免两相近数相减，故应变换所给公式。
（1）（2）4.近似数x*=0.0310,是 3 位有数数字。
5.计算取，利用：式计算误差最小。
四个选项：第二、三章插值与函数逼近习题二、三1. 给定的数值表用线性插值与二次插值计算ln0.54的近似值并估计误差限. 解：仍可使用n=1及n=2的Lagrange插值或Newton插值,并应用误差估计（5.8）。
线性插值时，用0.5及0.6两点，用Newton插值误差限，因，故二次插值时，用0.5，0.6，0.7三点，作二次Newton插值误差限，故2. 在-4≤x≤4上给出的等距节点函数表，若用二次插值法求的近似值，要使误差不超过，函数表的步长h 应取多少?解：用误差估计式（5.8），令因得3. 若，求和.解：由均差与导数关系于是4. 若互异，求的值，这里p≤n+1.解：，由均差对称性可知当有而当P＝n＋1时于是得5. 求证.解：解：只要按差分定义直接展开得6. 已知的函数表求出三次Newton均差插值多项式，计算f(0.23)的近似值并用均差的余项表达式估计误差.解：根据给定函数表构造均差表由式(5.14)当n=3时得Newton均差插值多项式N3(x)=1.0067x+0.08367x(x-0.2)+0.17400x(x-0.2)(x-0.3) 由此可得f(0.23) N3(0.23)=0.23203由余项表达式(5.15)可得由于7. 给定f(x)=cosx的函数表用Newton等距插值公式计算cos 0.048及cos 0.566的近似值并估计误差解：先构造差分表计算，用n=4得Newton前插公式误差估计由公式（5.17）得其中计算时用Newton后插公式（5.18)误差估计由公式（5.19）得这里仍为0.5658．求一个次数不高于四次的多项式p(x),使它满足解：这种题目可以有很多方法去做，但应以简单为宜。
数值分析第三版课本习题与答案
数值分析第三版课本习题与答案第⼀章绪论1. 设x >0,x 的相对误差为δ,求ln x 的误差.2. 设x 的相对误差为2％,求nx 的相对误差.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 -=…)计算到100Y .27.982(五位有效数字),试问计算100Y 将有多⼤误差?7. 求⽅程25610x x -+=的两个根,使它⾄少具有四位有效数字27.982).8. 当N 充分⼤时,怎样求211Ndx x +∞+?9. 正⽅形的边长⼤约为100㎝,应怎样测量才能使其⾯积误差不超过1㎝210. 设2⽽相对误差却减⼩.11. 序列{}n y 满⾜递推关系1101n n y y -=-(n=1,2,…),若0 1.41y =≈(三位有效数字),计算到10y 时误差有多⼤?这个计算过程稳定吗?12. 计算61)f =, 1.4≈,利⽤下列等式计算,哪⼀个得到的结果最好?3--13.()ln(f x x =,求f (30)的值.若开平⽅⽤六位函数表,问求对数时误差有多⼤?若改⽤另⼀等价公式ln(ln(x x =-计算,求对数时误差有多⼤?14. 试⽤消元法解⽅程组{101012121010;2.x x x x +=+=假定只⽤三位数计算,问结果是否可靠?15. 已知三⾓形⾯积1sin ,2s ab c =其中c 为弧度,02c π<<,且测量a ,b ,c 的误差分别为,,.a b c 证明⾯积的误差s ?满⾜.s a b cs a b c ≤++2.2)定义的德蒙⾏列式,令2000011211121()(,,,,)11n n n n n n n n n x x x V x V x x x x x x x xx x ----==证明()n V x 是n 次多项式,它的根是01,,n x x -,且101101()(,,,)()()n n n n V x V x x x x x x x ---=--.2. 当x = 1 , -1 , 2 时, f (x)= 0 , -3 , 4 ,求f (x )的⼆次插值多项式.3. 给出f (x )=ln x 的数值表⽤线性插值及⼆次插值计算ln 0.54 的近似值.4. 给出cos x ,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数字,研究⽤线性插值求cos x 近似值时的总误差界.5. 设0k x x kh =+,k =0,1,2,3,求032max ()x x x l x ≤≤.6. 设jx 为互异节点(j =0,1,…,n ),求证:i)()(0,1,,);nkkj jj x l x x k n =≡=∑ii)x l x k n =-≡0(=∑7. 设[]2(),f x C a b ∈且()()0f a f b ==,求证21()()().8maxmax a x ba xb f x b a f x ≤≤≤≤≤-"8. 在44x -≤≤上给出()x f x e =的等距节点函数表,若⽤⼆次插值求x e 的近似值,要使截断误差不超过610-,问使⽤函数表的步长h 应取多少?9. 若2n n y =,求4n y ?及4n y δ.10. 如果()f x 是m 次多项式,记()()()f x f x h f x ?=+-,证明()f x 的k 阶差分()(0)kf x k m ?≤≤是m k -次多项式,并且()0(m lf x l +?=为正整数).11. 证明1()k k k k k k f g f g g f +?=?+?.12. 证明110010.n n kkn n k k k k f gf g f g g f --+==?=--?∑∑13. 证明n j n j y y y -=?=?-?∑14. 若1011()n n n n f x a a x a x a x --=++++有n 个不同实根12,,,n x x x ,证明{10,02;, 1.1()n k njk n a k n j jx f x -≤≤-=-=='∑15. 证明n 阶均差有下列性质: i) 若()()F x cf x =,则[][]0101,,,,,,n n F x x x cf x x x =;ii) 若()()()F x f x g x =+,则[][][]010101,,,,,,,,,n n n F x x x f x x x g x x x =+.16. 74()31f x x x x =+++,求0172,2,,2f及0182,2,,2f.17. 证明两点三次埃尔⽶特插值余项是(4)22311()()()()/4!,(,)k k k k R x f x x x x x x ++=ξ--ξ∈并由此求出分段三次埃尔⽶特插值的误差限.19. 试求出⼀个最⾼次数不⾼于4次的函数多项式()P x ,以便使它能够满⾜以下边界条件(0)(0)0P P ='=,(1)(1)1P P ='=,(2)1P =.20. 设[](),f x C a b ∈,把[],a b 分为n 等分,试构造⼀个台阶形的零次分段插值函数()n x ?并证明当n →∞时,()n x ?在[],a b 上⼀致收敛到()f x .21.设2()1/(1)f x x=+,在55x-≤≤上取10n=,按等距节点求分段线性插值函数()hI x,计算各节点间中点处的()hI x与()f x的值,并估计误差.22.求2()f x x=在[],a b上的分段线性插值函数()hI x,并估计误差.()f x x=在[],a b上的分段埃尔⽶特插值,并估计误差.24.给定数据表如下:试求三次样条插值并满⾜条件i)(0.25) 1.0000,(0.53)0.6868; S S'='=ii)(0.25)(0.53)0. S S"="=25.若[]2(),f x C a b∈,()S x是三次样条函数,证明i)[][][][] 222()()()()2()()()b b b ba a a af x dx S x dx f x S x dx S x f x S x dx "-"="-"+""-";ii)若f x S x i n==,式中ix为插值节点,且01na x x x b=<<<=,则[][][]()()()()()()()()()b a S x f x S x dx S b f b S b S a f a S a ""-"="'-'-"'-'.26.编出计算三次样条函数()S x系数及其在插值节点中点的值的程序框图(()S x可⽤(8.7)式的表达式).第三章函数逼近与计算1.(a)利⽤区间变换推出区间为[],a b的伯恩斯坦多项式.(b)对()sinf x x=在[]0,/2π上求1次和三次伯恩斯坦多项式并画出图形,并与相应的马克劳林级数部分和误差做⽐较.2.求证:(a)当()m f x M≤≤时,(,)≤≤. (b)当()f x x=时,(,)nB f x x=.3.在次数不超过6的多项式中,求()sin4f x x=在[]0,2π的最佳⼀致逼近多项式.4.假设()f x在[],a b上连续,求()f x的零次最佳⼀致逼近多项式.5.选取常数a,使301maxxx ax≤≤-达到极⼩,⼜问这个解是否唯⼀?6. 求()sin f x x =在[]0,/2π上的最佳⼀次逼近多项式,并估计误差.7. 求()xf x e =在[]0,1上的最佳⼀次逼近多项式. 8. 如何选取r ,使()31f x x x =+-,在[]0,1上求三次最佳逼近多项式.10. 令[]()(21),0,1n n T x T x x =-∈,求***0123(),(),(),()T x T x T x T x . 11. 试证{}*()nT x 是在[]0,1上带权ρ=的正交多项式.12. 在[]1,1-上利⽤插值极⼩化求11()f x tg x -=的三次近似最佳逼近多项式.13. 设()xf x e =在[]1,1-上的插值极⼩化近似最佳逼近多项式为()n L x ,若n f L ∞-有界,证明对任何1n ≥,存在常数n α、n β,使11()()()()(11).n n n n n T x f x L x T x x ++α≤-≤β-≤≤14. 设在[]1,1-上234511315165()128243843840x x x x x x ?=-----,试将()x ?降低到3次多项式并估计误差. 15. 在[]1,1-上利⽤幂级数项数求()sin f x x =的3次逼近多项式,使误差不超过0.005.16. ()f x 是[],a a -上的连续奇(偶)函数,证明不管n 是奇数或偶数,()f x 的最佳逼近多项式*()n n F x H ∈也是奇(偶)函数.17. 求a 、b 使[]220sin ax b x dxπ+-?()(,)()();()(,)()()()();bbaaa f g f x g x dxb f g f x g x dx f a g a =''=''+??问它们是否构成积?19. ⽤许⽡兹不等式(4.5)估计6101x dx x +?的上界,并⽤积分中值定理估计同⼀积分的上下界,并⽐较其结果.20. 选择a ,使下列积分取得最⼩值:1122211(),x ax dx x ax dx----??.21. 设空间{}{}10010121,,,span x span x x 1?=?=,分别在1?、2?上求出⼀个元素,使得其为[]20,1x C ∈的最佳平⽅逼近,并⽐较其结果.22.()f x x =在[]1,1-上,求在{}2411,,span x x ?=上的最佳平⽅逼近.23.sin (1)arccos ()n n x u x +=是第⼆类切⽐雪夫多项式,证明它有递推关系()()()112n n n u x xu x u x +-=-.24. 将1()sin 2f x x=在[]1,1-上按勒让德多项式及切⽐雪夫多项式展开,求三次最佳平⽅逼近多项式并画出误差图形,再计算均⽅误差.25. 把()arccos f x x =在[]1,1-上展成切⽐雪夫级数.26. ⽤最⼩⼆乘法求⼀个形如2y a bx =+的经验公式,使它与下列数据拟合,并求均⽅误差.28. 在某化学反应⾥,根据实验所得分解物的浓度与时间关系如下:⽤最⼩⼆乘拟合求.29. 编出⽤正交多项式做最⼩⼆乘拟合的程序框图. 30. 编出改进FFT 算法的程序框图. 31. 现给出⼀记录{}{}4,3,2,1,0,1,2,3k x =,试⽤改进FFT 算法求出序列{}k x 的离散频谱{}k C (0,1,,7).k =第四章数值积分与数值微分1. 确定下列求积公式中的待定参数,使其代数精度尽量⾼,并指明所构造出的求积公式所具有的代数精度:。
数值分析分享
第一章习题解答1.在下列各对数中，X是精确值ａ的近似值（1）ａ=π，x=3.1 （2）ａ=1/7，x=0.143（3）ａ=π/1000，x=0.0031 （4）ａ=100/7，x=14.3试估计x的绝对误差和相对误差。
解：（1） e=∣3.1-π∣≈0.0416，δr= e/∣x∣≈0.0143（2） e=∣0.143-1/7∣≈0.0143 δr= e/∣x∣≈0.1（3） e=∣0.0031-π/1000∣≈0.0279 δr= e/∣x∣≈0.9(4) e=∣14.3-100/7∣≈0.0143 δr= e/∣x∣≈0.0012.已知四个数：x1=26.3，x2=0.0250， x3= 134.25，x4=0.001。
试估计各近似数的有效位数和误差限，并估计运算μ1= x1x2x3和μ1= x3x4/x1的相对误差限。
解：x1=26.3 ｎ=3 δx1=0.05 δrx1=δx1/∣x1∣=0.19011×10-2x2=0.0250 ｎ=3 δx2=0.00005 δrx2=δx2/∣x2∣=0.2×10-2x3= 134.25 ｎ=5 δx3=0.005 δrx3=δx3/∣x3∣=0.372×10-4x 4=0.001 ｎ=1 δx4=0.0005 δrx4=δx4/∣x4∣=0.5由公式：er （μ）= e（μ）/∣μ∣≦1/∣μ∣Σni=1∣∂f/∂xi∣δxie r （μ1）≦1/∣μ1∣［x2x3δx1+ x1x3δx2+x1x2δx3］=0.34468/88.269275 =0.0039049e r （μ2）≦1/∣μ2∣［-x3x4/ x21δx1+ x4/ x1δx3+ x3/ x1δx4］=0.497073.设精确数ａ>0，x是ａ的近似值，x的相对误差限是0.2，求㏑x的相对误差限。
解：δr ≦Σni=1∣∂f/∂xi∣δxi=1/㏑x·1/ x·δx=δr x/㏑x=0.2/㏑x即δr≦0.2/㏑x4.长方体的长宽高分别为50cm，20cm和10cm，试求测量误差满足什么条件时其表面积的误差不超过1cm2。
数值分析课后习题及答案
数值分析课后习题及答案第一章绪论（12）第二章插值法（40-42）2、当时，，求的二次插值多项式。
[解]。
3、给出的数值表用线性插值及二次插值计算的近似值。
X 0.4 0.5 0.6 0.7 0.8 -0.916291 -0.693147 -0.510826 -0.357765 -0.223144 [解]若取，，则，，则，从而。
若取，，，则，，，则，从而补充题：1、令，，写出的一次插值多项式，并估计插值余项。
[解]由，可知，，余项为，故。
2、设，试利用拉格朗日插值余项定理写出以为插值节点的三次插值多项式。
[解]由插值余项定理，有，从而。
5、给定数据表：，1 2 4 6 7 4 1 0 1 1 求4次牛顿插值多项式，并写出插值余项。
[解]一阶差商二阶差商三阶差商四阶差商 1 42 1 -34 0 6 17 1 0 由差商表可得4次牛顿插值多项式为：，插值余项为。
第三章函数逼近与计算（80-82）26、用最小二乘法求一个形如的经验公式，使它与下列数据相拟合，并求均方误差。
19 25 31 38 44 19.0 32.3 49.0 73.3 97.8[解]由。
又，，，故法方程为，解得。
均方误差为。
27、观测物体的直线运动，得出以下数据：时间t（秒）0 0.9 1.9 3.0 3.9 5.0 距离s（米）0 10 30 5080 110 [解]设直线运动为二次多项式，则由。
，。
又，，，故法方程为，解得。
故直线运动为。
补充题：1、现测得通过某电阻R的电流I及其两端的电压U如下表：I ……U ……试用最小二乘原理确定电阻R的大小。
[解]电流、电阻与电压之间满足如下关系：。
应用最小二乘原理，求R使得达到最小。
对求导得到：。
令，得到电阻R为。
2、对于某个长度测量了n次，得到n个近似值，通常取平均值作为所求长度，请说明理由。
[解]令，求x使得达到最小。
对求导得到：，令，得到，这说明取平均值在最小二乘意义下误差达到最小。