数值分析课后习题与解答

合集下载

数值分析课后习题全解(可编辑优质文档)

数值分析课后习题全解（可编辑优质文档）(可以直接使用，可编辑完整版资料，欢迎下载）第5章数值分析课后习题全解第5章：解线性方程组的直接方法1．证明：由消元公式及A 的对称性得(2)211,,2,3,..........,1111111a a j i a a a a a a i j na a ijij j j iji =-=-== 故2A对称2．证明：（1）因A 对称正定，故,)0,1,2,......,e i ni >=aii=(Ae i其中i e =（0,…,0,1,0,...,0）T 为第i 个单位向量.(2)由A 的对称性及消元公式得111211122222n n nn n n u u u d d u u d d u d d ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(2)ij a =ija -1111i j a a a =ji a -111j a a 1i a =(2)ji a ，I,j=2,…,n故2A 也对称.又 11120Ta a A ⎡⎤⎢⎥⎣⎦=1L A 1TAL其中 1L =211111111.....1n a a a a ⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦显然1L 非其异,从而对任意的x ≠0,有1TL X ≠0,(x,1L A 1TL X)=(1TL x, A 1TL X)>0 (由A 的正定性) 故11T L AL 正定.又11T L AL =11200a A ⎡⎤⎢⎥⎣⎦,而11a >0,故2A 正定. 3.证明由矩阵乘法简单运算即得证.4.解设有分解4232125316⎡⎤⎢⎥-⎢⎥⎢⎥⎢⎥-⎣⎦=123431231αααα⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦1231111βββ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦由公式11,1111,2,3,,,2,3,.1i i i i ii i b a c b i n c i n αβαβααβ-==⎧⎪=+=⎨⎪==-⎩其中i b ,i a ,i c 分别是系数矩阵的主对角线元素及下边和上边的次对角线元素.故有112233414,272,27397,7138513αβαβαβα⎧==⎪⎪⎪=-=-⎪⎨⎪==⎪⎪⎪=⎩从而有4232125316⎡⎤⎢⎥-⎢⎥⎢⎥⎢⎥-⎣⎦=4732392785113⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎣⎦11221771131⎡⎤⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦故 1y =64=32, 2y =12372y --=573y =21022039137y -=, 4y =3518513y +=故4x =1,3x =420711313x -=,2x =352177x +=,1x =231122x -= 5. 解 (1)设U 为上三角阵1112112222n n nn n u u u x u u x u x ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦=12n d d d ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦因nn n u x =n d ，故n x =nnnd u . 因 10001010302171101⎡⎤⎢⎥-⎢⎥⎢⎥-⎢⎥-⎢⎦⎣ii i u x +1n ij j j i u x =+∑=id ,故i x =1ni ij ij i iid u xu =+-∑,i=n-1,n-2,,1当U 为下三角阵时11212212n n nn u u u u u u ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦12n x x x ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦= 12n d d d ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦得,1x =111d u , 1x =11i i ij jj iid u x u -=-∑,i=2,3,…,n.(2)除法次数为n,乘法次数为1+2+…+(n-1)=n(n-1)/2 故总的乘法次数为n+n(n-1)/2=n(n+1)/2. (3)设U 为上三角阵,1U-=S,侧S 也是上三角阵.由11121222n n nn u u u u s u ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦11121222n n nn s s s s s s ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦=111⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦得 1ii iis u =, i=1,2,…,nij s =-1jik kjk i iiusu =+∑,j=i+1,i+2,…,n; i=n-1,n-2,…,1当U 为下三角阵时,由11212212n n nn d d d d d d ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦11212212n n nn s s s s s s ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦= 111⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦得 1ii iis u =,i=1,2,…,n ij s =11i ik kjk iiusu -=-∑,i=2,3,…,n;j=1,2,…,i-16. 证明 (1)因A 是对称正定阵,故存在唯一的分解A=L TL ,其中L 是具有正对角元素的下三角阵.从而 1A -=(L TL )1-=(TL )1-L 1-=(L 1-)TL 1-(A 1-)T =11()TT L L --⎡⎤⎣⎦=11()T L L --=1A -故1A -是对称矩阵.又1L -非奇异,故对任意的 x ≠0,有1L -x ≠0,故1Tx A -X=11()T T x L L x --=11()()T L x L x -->0故1A -是对称正定矩阵,即1A -也对称正定.(2)由A 对称正盯,故A 的所有顺序主子式均不为零,从而A 有唯一的 Doolittle 分解A=L U.又U=1122nn u u u ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦1121111222111n n u u uu u u ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦=D 0U 其中D 为对三角阵, 0U 为单位上三角阵,于是 A=U L =D L 0U又 A=TA =TO U D TL由分解的唯一性即得T O U =L从而有 A=D L TL 又由A 的对称正定性知 1d =1D >0, i d =1ii D D ->0 (i=2,3,…,n) 故 D=12n d d d ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦=⎤⎥⎢⎥⎢⎥⎢⎢⎣⎤⎥⎢⎥⎢⎥⎢⎢⎣=12D 12D 故A=L D TL =L12D12D TL =(L 12D )(L 12D )T =LL T其中L=L 12D 为三角元为正的下三角矩阵.7. 解[A|I]=21311000310701001242001010150001⎡--⎤⎢⎥⎢⎥⎢⎥--⎢⎥-⎢⎦⎣-> 101500010138010302330011011111002⎡-⎤⎢⎥--⎢⎥⎢⎥⎢⎥---⎢⎦⎣-> 1015010138010300319021700431101⎡-⎤⎢⎥--⎢⎥⎢⎥--⎢⎥---⎢⎦⎣->421410003333010110114192170010333385542500013333⎡⎤---⎢⎥⎢⎥-⎢⎥⎢⎥---⎢⎥⎢⎥⎢⎥---⎢⎥⎣⎦-> 410231685178517100033641130100851785170010191538851785170001314585178517⎤--⎥⎥⎡⎥-⎢⎥⎢⎥⎢⎥---⎢⎥⎣⎥⎥--⎦->1A -=4102316851785173364113851785171953885178517314585178517⎡⎤--⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥---⎢⎥⎢⎥⎢⎥--⎣⎦=0.04705890.58823530.27058820.94117650.38823530.35294120.48235290.76470590.22352940.29411760.03529410.47058820.03529410.05882350.04705890.2941176--⎡⎤⎢⎥-⎢⎥⎢⎥---⎢⎥--⎣⎦8．解设有分解2112112112112-⎡⎤⎢⎥--⎢⎥⎢⎥--⎢⎥--⎢⎥⎢⎥-⎣⎦= 123451111ααααα⎡⎤⎢⎥-⎢⎥⎢⎥-⎢⎥-⎢⎥⎢⎥-⎣⎦123451111βββββ⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦由公式11,1111,(2,3,4,5),(2,3,4)i i i i i i i b c b i c i ααβαβααβ-==⎧⎪=+=⎨⎪==⎩其中i b ，i a ，i c 分别是系数矩阵的主角线元素及其下边和上边的次对角线元素，则有12α=， 232α=， 343α=， 454α=， 565α= 112β=-， 223β=-， 334β=-， 445β=-由12345231120410305140615y y y y y ⎡⎤⎢⎥⎢⎥-⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥-⎢⎥⎣⎦得1y =12，213y =，314y =，415y =，516y = 由112213314415⎡⎤-⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥-⎢⎥⎣⎦123451213141516x x x x x ⎡⎤⎢⎥⎢⎥⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎢⎥⎣⎦得5x =16，4x =13，3x =12，2x =23，1x =569．解设211123131-⎡⎤⎢⎥-⎢⎥⎢⎥⎣⎦=121312123231323111111d l l l d l l l d ⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦由矩阵乘法得1d =2， 2112l =-， 3112l = 252d =-，3275l =-3275d =由123141152617125y y y ⎡⎤⎢⎥⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥-=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦⎢⎥-⎣⎦得 14y =，27y =，3695y = 由12311122245717256927155x x x ⎡⎤⎡⎤-⎡⎤⎢⎥⎢⎥⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥--=⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎣⎦ 得 123111222247514175256927123559x x x ⎡⎤⎡⎤⎡⎤-⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥-=-=-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎣⎦⎣⎦故3x =239=2.555 555 6，279x ==0.777 777 8，1x =109=1.111 111 1 10．解 A 中2∆=0，故不能分解。

数值分析课后习题及答案

第一章绪论（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 ，并估计插值余项。

数值分析课后习题答案

0 1
0 10 1 1 0 0 0 1
0 0 12 1 1 2 0 0 0

1 2
0 0 0 1 1 0
1 2

1 2

1 2
1
0 0 0 1 0

1 2

1 2

0
1 2

1 2
0
0
0
341 1 1
2-5.对矩阵A进行LDLT分解和GGT分解，并求解方程组
Ax=b，其中
16 4 8
1
A 4 5 4 , b 2
8 4 22
3
解
16 A 4
4 5
84
44 11
2-3(1).对矩阵A进行LU分解,并求解方程组Ax=b,其中
2 1 1 A1 3 2
4 ,b6
1 2 2
5
解
2 A 1
1 3
1 2

2 11
22
1
5 2
1

3 21来自，所以 A12
1
2 1 1

5 3
2-2(1).用列主元Gauss消元法解方程组
3 2 6x1 4 10 7 0x2 7 5 1 5x3 6
解
3 2 6 4 10 7 0 7 10 7 0 7

r1r2
消元

10 7 0 7 3 2 6 4 0 0.1 6 6.1
r=0.5101-n/3.162…<0.5101-n/3<0.01% 因此只需n=5.即取101/2=3.1623

数值分析课后习题部分参考答案解析

数值分析课后习题部分参考答案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相⽐较。

电子科技大学-数值分析答案-钟尔杰

| x n +1 − 7 |=
而xn具有n位有效数，故
所以
| x n +1 − 7 |≤
由此得xn+1的误差限
1 2 7
| x n − 7 |2 ≤
1 × × 10 2− 2 n 2 7 4
1
| x n +1 − 7 |≤
1 × 10 1− 2 n 2
故，xn+1是 7 的具有 2n位有效数字的近似值。三、问题 1．假定 a0，b0是非负实数且a0≠b0，按如下递推公式
∑ [ai ∑ b j ]
i =1 j =1
n，仍为( n + 2 ) ( n – 1) / 2。，算法输出 11 试构造一个算法，对输入的数据 x0，x1，x2，……，xn，以及x（均为实数）为 ( x –x0) ( x –x1) ( x –x2)……( x –xn) 的计算结果。解算法如下：第一步：输入x；x0，x1，x2，……，xn，M Å (x – x0 )；k Å 0；第二步：M Å M×(x – x0 )；k Å k+1；第三步：判断，若 k ≤ n，则转第二步；否则输出 M，结束。 12 利用级数公式
4
π 1 dx = arctan 1 = 可以计算出无理数π 的值。将定积分表示为积分和 2 4 1+ x
R
H
∫
1
0
xn dx ( n = 1，2，…，20) 的递推 5+ x
关系，并研究递推算法的数值稳定性。 6．计算两个多项式Pn(x)和Qm(x)的乘积多项式Tn+m(x)的方法称为向量的卷积方法。设
第一章习题解答与问题
一、习题解答 1 设 x>0，x 的相对误差限为 δ，求 ln x 的误差。解：设 x的准确值为x*，则有 ( | x – x* | /|x*| ) ≤ δ 所以 e(ln x)=| ln x – ln x* | =| x – x* | ×| (ln x)’|x=ξ·≈ ( | x – x* | / | x*| ) ≤ δ 另解： e(ln x)=| ln x – ln x* | =| ln (x / x*) | = | ln (( x – x* + x*)/ x*) | = | ln (( x – x* )/ x* + 1) |≤( | x – x* | /|x*| ) ≤ δ 2 设 x = – 2.18 和 y = 2.1200 都是由准确值经四舍五入而得到的近似值。求绝对误差限 ε( x ) 和 ε( y ) 。解：| e(x) | = |e(– 2.18)|≤ 0.005，| e(y) | = |e( 2.1200)|≤ 0.00005，所以 ε( x )=0.005， ε( y ) = 0.00005。 3 下近似值的绝对误差限都是 0.005，问各近似值有几位有效数字 x1=1.38，x2= –0.0312，x3= 0.00086 解：根据有效数字定义，绝对误差限不超过末位数半个单位。由题设知，x1，x2， x3有效数末位数均为小数点后第二位。故x1具有三位有效数字，x2具有一位有效数字，x3具有零位有效数字。 4 已知近似数 x 有两位有效数字，试求其相对误差限。解：| er(x) | ≤ 5 × 10– 2 。 5 设 y0 = 28，按递推公式 yn = yn-1 –

数值分析课后习题答案

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.求以下矩阵的特征值和特征向量。

数值分析第五版课后习题答案

数值分析第五版课后习题答案数值分析是一门应用数学的分支学科，主要研究如何利用数值方法解决实际问题。

在学习这门课程的过程中，课后习题是不可或缺的一部分。

本文将对《数值分析第五版》的课后习题进行一些探讨和解答。

第一章是数值分析的导论，主要介绍了误差分析和计算方法的基本概念。

在课后习题中，有一道题目是关于误差传播的。

假设有一个函数f(x, y) = x^2 + y^2，其中x和y的测量误差分别为Δx和Δy，要求计算f(x, y)的误差。

解答：根据误差传播公式，可以得到f(x, y)的误差为Δf = √[(∂f/∂x)^2 *(Δx)^2 + (∂f/∂y)^2 * (Δy)^2]。

对于本题而言，∂f/∂x = 2x，∂f/∂y = 2y。

代入公式，得到Δf = √[(2x)^2 * (Δx)^2 + (2y)^2 * (Δy)^2] = 2√(x^2 * (Δx)^2+ y^2 * (Δy)^2)。

第二章是插值与多项式逼近的内容。

其中一道习题涉及到拉格朗日插值多项式。

给定n+1个数据点(x0, y0), (x1, y1), ..., (xn, yn)，要求构造一个n次多项式p(x)，使得p(xi) = yi (i = 0, 1, ..., n)。

解答：拉格朗日插值多项式的表达式为p(x) = Σ(yi * Li(x))，其中Li(x) = Π[(x - xj) / (xi - xj)]，j ≠ i。

将数据点代入表达式中，即可得到所求的多项式。

第三章是数值微积分的内容，其中一道习题是关于数值积分的。

给定一个函数f(x)，要求使用复化梯形公式计算定积分∫[a, b]f(x)dx。

解答：复化梯形公式的表达式为∫[a, b]f(x)dx ≈ h/2 * [f(a) + 2Σf(xi) + f(b)]，其中h = (b - a)/n，xi = a + i * h (i = 1, 2, ..., n-1)。

根据给定的函数f(x)，代入公式中的各个值，即可得到近似的定积分值。

数值分析课后习题答案

第一章习题解答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. 已知四个数：x 1=26.3，x 2=0.0250， x 3= 134.25，x 4=0.001。

试估计各近似数的有效位数和误差限，并估计运算μ1= x 1 x 2 x 3和μ1= x 3 x 4 /x 1的相对误差限。

解：x 1=26.3 ｎ=3 δx 1=0.05 δr x 1=δx 1/∣x 1∣=0.19011×10-2x 2=0.0250 ｎ=3 δx 2=0.00005 δr x 2=δx 2/∣x 2∣=0.2×10-2x 3= 134.25 ｎ=5 δx 3=0.005 δr x 3=δx 3/∣x 3∣=0.372×10-4x 4=0.001 ｎ=1 δx 4=0.0005 δr x 4=δx 4/∣x 4∣=0.5由公式：e r （μ）= e （μ）/∣μ∣≦1/∣μ∣Σni=1∣∂f/∂x i ∣δx ie r （μ1）≦1/∣μ1∣［x 2 x 3δx 1+ x 1 x 3δx 2 +x 1x 2δx 3］ =0.34468/88.269275 =0.0039049e r （μ2）≦1/∣μ2∣［-x 3 x 4/ x 21δx 1+ x 4/ x 1δx 3 + x 3/ x 1δx 4］ =0.497073. 设精确数ａ>0，x 是ａ的近似值，x 的相对误差限是0.2，求㏑x 的相对误差限。

数值分析课程第五版课后习题答案

数值分析课程第五版课后习题答案课后习题一：a) 求解非线性方程f(x) = x^3 - 2x - 5的根。

解答：可使用牛顿迭代法来求解非线性方程的根。

牛顿迭代法的迭代公式为：x_(n+1) = x_n - f(x_n)/f'(x_n)，其中x_n为第n次迭代的近似解。

对于给定的方程f(x) = x^3 - 2x - 5，计算f'(x)的导数为f'(x) = 3x^2 - 2。

选择一个初始近似解x_0，并进行迭代。

迭代的终止条件可以选择两次迭代间的解的差值小于某个预设的精度。

b) 计算矩阵加法和乘法的运算结果。

解答：设A和B为两个矩阵，A = [a_ij]，B = [b_ij]，则A和B的加法定义为C = A + B，其中C的元素为c_ij = a_ij + b_ij。

矩阵乘法定义为C = A * B，其中C的元素为c_ij = ∑(a_ik * b_kj)，k的取值范围为1到矩阵的列数。

c) 使用插值方法求解函数的近似值。

解答：插值方法可用于求解函数在一组给定点处的近似值。

其中，拉格朗日插值法是一种常用的方法。

对于给定的函数f(x)和一组插值节点x_i，i的取值范围为1到n，利用拉格朗日插值多项式可以构建近似函数P(x)，P(x) = ∑(f(x_i) * l_i(x))，其中l_i(x)为拉格朗日基函数，具体表达式为l_i(x) = ∏(x - x_j)/(x_i - x_j)，j的取值范围为1到n并且j ≠ i。

课后习题二：a) 解决数值积分问题。

解答：数值积分是求解定积分的数值近似值的方法。

常用的数值积分方法包括矩形法、梯形法和辛普森法。

矩形法采用矩形面积的和来近似曲边梯形的面积，梯形法采用等距离子区间上梯形面积的和来近似曲边梯形的面积，而辛普森法则利用等距离子区间上梯形和抛物线面积的加权和来近似曲边梯形的面积。

b) 使用迭代方法求解线性方程组。

解答：线性方程组的求解可以通过迭代方法来进行。

数值分析第三版课本习题及答案

第一章绪论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相比较。

数值分析课后习题答案

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

数值分析课后部分习题答案

证明由差商的定义 (a) 如果 F ( x ) = cf ( x ) ，则
F [ x0 , x1 ,⋯ , xn ] =
=
F [ x1 , x2 ,⋯ , xn ]-F [ x0 , x1 ,⋯ , xn− 1 ] x n − x0
cf [ x1 , x2 , ⋯ , xn ]-cf [ x0 , x1 ,⋯ , xn −1 ] x n − x0 f [ x1 , x2 , ⋯ , xn ]-f [ x0 , x1 ,⋯ , xn−1 ] = cf [ x0 , x1 , ⋯ , xn ] . x n − x0
1 1 1 1 |e( x*)| ≤ × 10m − n = × 10−2 ， |e( y*)| ≤ × 10m − n = × 10 −2 ， 2 2 2 2 1 1 |e( z*)| ≤ × 10 m − n = × 10 −2 ， 2 2 | e( y * z*) |≈| z * e ( y*) + y * e ( z *) |≤ z * | e ( y *) | + y * | e (z *) |
m − n = −3 ，所以， n = 4 ； z * = 0.00052 = 0.52 × 10−3 ，即 m = −3
1 1 × 10m − n = × 10−3 ， 2 2
由有效数字与绝对误差的关系得即
m − n = −3 ，所以， n = 0 .
1 1 × 10m − n = × 10−3 ， 2 2
1 1 ≤ 2.35 × × 10−2 + 1.84 × × 10−2 = 2.095 × 10−2 ， 2 2 1 | e( x * + y * z*) |≈| e( x*) + e( y * z*) |≤ × 10 −2 + 2.095 × 10−2 2 1 = 0.2595 × 10−1 ≤ × 10−1 ， 2

数值分析详细答案(全)

第二章插值法习题参考答案2.)12)(12()1)(1(4)21)(11()2)(1()3()21)(11()2)(1(0)(2+-+-⋅+------⋅-+-+-+⋅=x x x x x x x L3723652-+=x x . 3. 线性插值：取510826.0,693147.0,6.0,5.01010-=-===y y x x ，则620219.0)54.0()54.0(54.0ln 0010101-=-⋅--+=≈x x x y y y L ；二次插值：取510826.0,693147.0,916291.0,6.0,5.0,4.0210210-=-=-====y y y x x x ，则)54.0(54.0ln 2L ≈))(()54.0)(54.0())(()54.0)(54.0())(()54.0)(54.0(120210221012012010210x x x x x x y x x x x x x y x x x x x x y ----⋅+----⋅+----⋅=＝－0.616707 .6. i) 对),,1,0(,)(n k x x f k==在n x x x ,,,10 处进行n 次拉格朗日插值，则有)()(x R x P x n n k +=)())(()!1(1)(0)1(0n n ni k j j x x x x f n x x l --++=+=∑ ξ由于0)()1(=+ξn f，故有kni k j jxx x l≡∑=0)(.ii) 构造函数,)()(kt x x g -=在n x x x ,,,10 处进行n 次拉格朗日插值，有∑=-=ni j k j n x l t x x L 0)()()(.插值余项为 ∏=+-+=--nj j n n kx x n g x L t x 0)1()()!1()()()(ξ，由于).,,2,1(,0)()1(n k g n ==+ξ故有 .)()()()(0∑=-==-ni j k j n kx l t x x L t x令,x t =即得 ∑==-ni j k jx l t x)()(.8. 截断误差].4,4[),)()((61)(2102-∈---=ξξx x x x x x e x R其中 ,,1210h x x h x x +=-= 则hx x 331+=时取得最大值321044392|))()((|max h x x x x x x x ⋅=---≤≤- .由题意， ,10)392(61|)(|6342-=⋅⋅≤h e x R所以，.006.0≤h16. ;1!7!7!7)(]2,,2,2[)7(71===ξf f .0!7)(]2,,2,2[)8(810==ξf f19. 采用牛顿插值，作均差表：i x)(i x f一阶均差二阶均差0 1 20 1 11 0-1/2],,[))((],[)()()(210101000x x x f x x x x x x f x x x p x p --+-+=))()()((210x x x x x x Bx A ---++)2)(1()()2/1)(1(0--++--++=x x x Bx A x x x又由 ,1)1(,0)0(='='p p 得,41,43=-=B A 所以 .)3(4)(22-=x x x p第三章函数逼近与计算习题参考答案4.设所求为()g x c =，(,)max(,),max (),min ()a x ba x bf g M c m c M f x m f x ≤≤≤≤∆=--==，由47页定理4可知()g x 在[],a b 上至少有两个正负交错的偏差点，恰好分别为()f x 的最大值和最小值处，故由1(),()2M c m c c M m -=--=+可以解得1()()2g x M m =+即为所求。

数值分析习题(含标准答案)

数值分析习题（含标准答案）
一、选择题（每题5分，共20分）
1. 下列哪个选项不属于数值分析的研究范畴？
A. 数值微分
B. 数值积分
C. 数值逼近
D. 数据库管理
答案：D
2. 在数值分析中，求解线性方程组常用的方法有？
A. 高斯消元法
B. 迭代法
C. 拉格朗日乘数法
D. 上述所有方法
答案：D
3. 下列哪种方法适用于求解非线性方程组？
A. 牛顿法
B. 梯度下降法
C. 高斯消元法
D. 上述所有方法
答案：D
4. 在数值积分中，下列哪种方法具有最高的精度？
A. 梯形法则
B. 辛普森法则
C. 高斯求积法
D. 上述所有方法
答案：C
二、填空题（每题5分，共20分）
1. 数值分析的主要目的是通过有限步骤的运算，对数学问题进行近似求解。

2. 在数值微分中，常用的差分公式有前向差分、后向差分和中心差分。

3. 数值逼近的主要方法包括插值法和逼近法。

4. 在数值积分中，常用的方法有梯形法则、辛普森法则和高斯求积法。

三、解答题（每题10分，共30分）
1. 已知函数 f(x) = e^x，求其在 x = 0.5 处的导数。

答案：f'(0.5) ≈ 1.6487
2. 求解线性方程组 2x + 3y = 5，4x y = 1。

答案：x ≈ 0.625，y ≈ 1.25
3. 已知函数 f(x) = x^3 3x^2 + 4，求其在区间 [0, 2] 上的积分。

答案：f(x) 在区间 [0, 2] 上的积分≈ 3.6667。

(完整版)数值分析部分课后答案第二版朱晓临

数值分析第二版朱晓临第一章习题3.324.045≈324.0 60.0876≈60.090.00035167≈0.0003517 2.00043≈2.000 6.①**x x x-≤51441111111010100.005%222a a -+--⨯=⨯⨯≤⨯=（1≤1a ≤9）故它的相对误差限为0.005%②∵*12120....100....10n n n n x a a a a a a =±⨯=⨯＜()10.110na +⨯相对误差限＝0.03%***3311*n n n x x x x x x----=⨯⨯⨯⨯⨯⨯＜0.03%0.(a +1)10=0.3(0.a +1)10＜0.510 ∴至少有3位有效数字。

7.6*1), 1.4,0.004096A A =≈=则1.4≈时，⑴()610.005232781≈⑵(330.008-≈⑶()310.0051252613≈+⑷991-≈所以利用第三个得到的计算结果的绝对误差最小。

8.由函数的绝对误差公式：***(())'()()e f x f x e x ≈ ① 令2**2*(),()(),100f x x f x x x ===cm由题目得，*(())1e f x =，**'()2f x x = ②把②代入①，得： 1≈**2()x e x ⋅ 1≈*2100()e x ⨯⋅ *()e x 0.005cm ≈边长的测量误差不超过0.005cm 时，才能使其面积的误差不超过12cm 。

11.**()ln ,()ln f x x f x x ==令则由公式***(())'()()e f x f x e x ≈，得： ***1(())0.510e f x x x l x≈-<⨯- 又***()r x x x xε-≤，由此可知，*()0.510l r x ε-=⨯所以*x 的相对误差限为0.510l -⨯，有l 位有效数字。

数值分析课程第五版课后习题答案(李庆扬等)

第一章绪论（12）1、设0>x ，x 的相对误差为δ，求x ln 的误差。

[解]设0*>x 为x 的近似值，则有相对误差为δε=)(*x r ，绝对误差为**)(x x δε=，从而x ln 的误差为δδεε=='=*****1)()(ln )(ln x xx x x ，相对误差为****ln ln )(ln )(ln x x x x rδεε==。

2、设x 的相对误差为2%，求n x 的相对误差。

[解]设*x 为x 的近似值，则有相对误差为%2)(*=x r ε，绝对误差为**%2)(x x =ε，从而nx 的误差为nn x x nxn x x n x x x **1***%2%2)()()()(ln *⋅=='=-=εε，相对误差为%2)()(ln )(ln ***n x x x nr==εε。

3、下列各数都是经过四舍五入得到的近似数，即误差不超过最后一位的半个单位，试指出它们是几位有效数字：1021.1*1=x ，031.0*2=x ，6.385*3=x ，430.56*4=x ，0.17*5⨯=x 。

[解]1021.1*1=x 有5位有效数字；0031.0*2=x 有2位有效数字；6.385*3=x 有4位有效数字；430.56*4=x 有5位有效数字；0.17*5⨯=x 有2位有效数字。

4、利用公式（3.3）求下列各近似值的误差限，其中*4*3*2*1,,,x x x x 均为第3题所给的数。

（1）*4*2*1x x x ++；[解]3334*4*2*11***4*2*1*1005.1102110211021)()()()()(----=⨯=⨯+⨯+⨯=++=⎪⎪⎭⎫ ⎝⎛∂∂=++∑x x x x x f x x x e nk k k εεεε；（2）*3*2*1x x x ；[解]52130996425.010********.2131001708255.01048488.2121059768.01021)031.01021.1(1021)6.3851021.1(1021)6.385031.0()()()()()()()()(3333334*3*2*1*2*3*1*1*3*21***3*2*1*=⨯=⨯+⨯+⨯=⨯⨯+⨯⨯+⨯⨯=++=⎪⎪⎭⎫⎝⎛∂∂=-------=∑x x x x x x x x x x x f x x x e n k k kεεεε；（3）*4*2/x x 。

数值分析课后习题及答案

数值分析课后习题及答案第一章绪论（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使得达到最小。

对求导得到：，令，得到，这说明取平均值在最小二乘意义下误差达到最小。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

课后习题解答第一章绪论习题一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等距插值公式计算cos0.048及cos0.566的近似值并估计误差解：先构造差分表计算，用n=4得Newton前插公式误差估计由公式（5.17）得其中计算时用Newton后插公式（5.18)误差估计由公式（5.19）得这里仍为0.5658．求一个次数不高于四次的多项式p(x),使它满足解：这种题目可以有很多方法去做，但应以简单为宜。

此处可先造使它满足，显然，再令p(x)=x2(2-x)+Ax2(x-1)2由p(2)=1求出A＝，于是9.令称为第二类Chebyshev多项式，试求的表达式，并证明是［-1,1］上带权的正交多项式序列。

解：因10.用最小二乘法求一个形如的经验公式，使它拟合下列数据，并计算均方误差.解：本题给出拟合曲线，即，故法方程系数法方程为解得最小二乘拟合曲线为均方程为11.填空题(1)满足条件的插值多项式p(x)=().(2),则f［1,2,3,4］=()，f［1,2,3,4,5］=().(3)设为互异节点，为对应的四次插值基函数，则＝()，＝().(4)设是区间［0,1］上权函数为ρ(x)=x的最高项系数为1的正交多项式序列，其中，则＝()，＝()答：（1）（2）（3）（4）第4章数值积分与数值微分习题41.分别用复合梯形公式及复合Simpson公式计算下列积分.解本题只要根据复合梯形公式（6.11）及复合Simpson 公式（6.13）直接计算即可。

对，取n=8,在分点处计算f(x)的值构造函数表。

按式（6.11）求出，按式（6.13）求得，积分2.用Simpson公式求积分，并估计误差解：直接用Simpson公式（6.7）得由（6.8）式估计误差，因，故3.确定下列求积公式中的待定参数，使其代数精确度尽量高，并指明求积公式所具有的代数精确度.(1)(2)(3)解：本题直接利用求积公式精确度定义，则可突出求积公式的参数。

（1）令代入公式两端并使其相等，得解此方程组得，于是有再令，得故求积公式具有3次代数精确度。

（2）令代入公式两端使其相等，得解出得而对不准确成立，故求积公式具有3次代数精确度。

（3）令代入公式精确成立，得解得，得求积公式对故求积公式具有2次代数精确度。

4.计算积分，若用复合Simpson公式要使误差不超过，问区间要分为多少等分?若改用复合梯形公式达到同样精确度，区间应分为多少等分?解：由Simpson公式余项及得即，取n=6，即区间分为12等分可使误差不超过对梯形公式同样，由余项公式得即取n=255才更使复合梯形公式误差不超过5.用Romberg求积算法求积分，取解：本题只要对积分使用Romberg算法（6.20），计算到K＝3，结果如下表所示。

于是积分，积分准确值为0.7132726．用三点Gauss-Legendre求积公式计算积分.解：本题直接应用三点Gauss公式计算即可。

由于区间为，所以先做变换于是本题精确值7．用三点Gauss-Chebyshev求积公式计算积分解：本题直接用Gauss-Chebyshev求积公式计算即于是，因n=2,即为三点公式，于是，即故8.试确定常数A，B，C，及α，使求积公式有尽可能高的代数精确度，并指出所得求积公式的代数精确度是多少.它是否为Gauss型的求积公式?解：本题仍可根据代数精确度定义确定参数满足的方程，令对公式精确成立，得到由（2）（4）得A=C，这两个方程不独立。

故可令，得（5）由（3）（5）解得，代入（1）得则有求积公式令公式精确成立，故求积公式具有5次代数精确度。

三点求积公式最高代数精确度为5次，故它是Gauss型的。

第五章解线性方程组的直接法习题五1.用Gauss消去法求解下列方程组.解本题是Gauss消去法解具体方程组，只要直接用消元公式及回代公式直接计算即可。

故2.用列主元消去法求解方程组并求出系数矩阵A的行列式detA的值解：先选列主元，2行与1行交换得消元3行与2行交换消元回代得解行列式得3.用Doolittle分解法求的解.解：由矩阵乘法得再由求得由解得4.下述矩阵能否作Doolittle分解，若能分解，分解式是否唯一?解：A中，若A能分解，一步分解后，，相互矛盾，故A不能分解，但，若A中1行与2行交换，则可分解为LU对B，显然，但它仍可分解为分解不唯一，为一任意常数，且U奇异。

C可分解，且唯一。

5.用追赶法解三对角方程组Ax=b，其中解：用解对三角方程组的追赶法公式（3.1.2）和（3.1.3）计算得6.用平方根法解方程组解：用分解直接算得由及求得7.设，证明解：即，另一方面故8．设计算A的行范数，列范数及F-范数和2范数解：故9．设为上任一种范数，是非奇异的，定义，证明证明：根据矩阵算子定义和定义，得令，因P非奇异，故x与y为一对一，于是10.求下面两个方程组的解，并利用矩阵的条件数估计.，即，即解：记则的解，而的解故而由（3.12）的误差估计得表明估计略大，是符合实际的。

11.是非题（若"是"在末尾（）填+,"不是"填-）：题目中（1）若A对称正定,,则是上的一种向量范数（）（2）定义是一种范数矩阵（）（3）定义是一种范数矩阵（）（4）只要，则A总可分解为A=LU,其中L为单位下三角阵，U为非奇上三角阵（）（5）只要，则总可用列主元消去法求得方程组的解（）（6）若A对称正定,则A可分解为，其中L为对角元素为正的下三角阵（）（7）对任何都有（）（8）若A为正交矩阵，则（）答案：（1）（＋）（2）（－）（3）（＋）（4）（－）（5）（＋）（6）（＋）（7）（－）（8）（＋）第六章解线性方程组的迭代法习题六1.证明对于任意的矩阵A，序列收敛于零矩阵解：由于而故2.方程组(1)考查用Jacobi法和GS法解此方程组的收敛性.(2)写出用J法及GS法解此方程组的迭代公式并以计算到为止解：因为具有严格对角占优，故J法与GS法均收敛。

（2）J法得迭代公式是取，迭代到18次有GS迭代法计算公式为取3.设方程组证明解此方程的Jacobi迭代法与Gauss-Seidel迭代法同时收敛或发散解：Jacobi迭代为其迭代矩阵，谱半径为，而Gauss-Seide 迭代法为其迭代矩阵，其谱半径为由于，故Jacobi迭代法与Gauss-Seidel法同时收敛或同时发散。

4.下列两个方程组Ax=b，若分别用J法及GS法求解，是否收敛?解：Jacobi法的迭代矩阵是即，故，J法收敛、GS法的迭代矩阵为故，解此方程组的GS法不收敛。

5.设，detA≠0，用,b表示解方程组Ax=f 的J法及GS法收敛的充分必要条件.解J法迭代矩阵为，故J法收敛的充要条件是。

GS法迭代矩阵为由得GS法收敛得充要条件是6.用SOR方法解方程组(分别取ω=1.03,ω=1,ω=1.1)精确解，要求当时迭代终止，并对每一个ω值确定迭代次数解：用SOR方法解此方程组的迭代公式为取，当时，迭代5次达到要求若取，迭代6次得7.对上题求出SOR迭代法的最优松弛因子及渐近收敛速度，并求J法与GS法的渐近收敛速度.若要使那么J法GS法和SOR法各需迭代多少次?解：J法的迭代矩阵为，故，因A为对称正定三对角阵，最优松弛因子J法收敛速度由于，故若要求，于是迭代次数对于J法，取K＝15对于GS法，取K＝8对于SOR法，取K＝58.填空题(1)要使应满足（）.(2)已知方程组，则解此方程组的Jacobi迭代法是否收敛（）.它的渐近收敛速度R(B)=（）.(3)设方程组Ax=b,其中其J法的迭代矩阵是（）.GS法的迭代矩阵是（）.(4)用GS法解方程组，其中a为实数，方法收敛的充要条件是a满足（）.(5)给定方程组,a为实数.当a满足（），且0＜ω＜2时SOR迭代法收敛.答:(1)(2)J法是收敛的，(3)J法迭代矩阵是，GS法迭代矩阵(4)满足(5)满足第七章非线性方程求根习题七1.用二分法求方程的正根，使误差小于0.05解使用二分法先要确定有根区间。

本题f(x)=x2-x-1=0,因f(1)=-1,f(2)=1,故区间[1,2]为有根区间。

另一根在[-1,0]内，故正根在[1,2]内。

用二分法计算各次迭代值如表。

其误差2.求方程在=1.5附近的一个根，将方程改写成下列等价形式，并建立相应迭代公式.(1)，迭代公式.(2),迭代公式.(3)，迭代公式.试分析每种迭代公式的收敛性，并选取一种收敛最快的方法求具有4位有效数字的近似根解：（1）取区间且，在且，在中，则L<1，满足收敛定理条件，故迭代收敛。

（2），在中，且，在中有，故迭代收敛。

（3），在附近，故迭代法发散。

在迭代（1）及（2）中，因为（2）的迭代因子L较小，故它比（1）收敛快。