数值分析第四版习题及答案

数值分析习题

第一章 绪 论

1. 设x>0,x的相对误差为δ,求lnx的误差. 2. 设x的相对误差为2％,求x的相对误差.

3. 下列各数都是经过四舍五入得到的近似数,即误差限不超过最后一位的半个单位,试指

出它们是几位有效数字:

*****x11.1021,x20.031,x3385.6,x456.430,x571.0.

n4. 利用公式求下列各近似值的误差限:

************(i)x1x2x4,(ii)x1x2x3,(iii)x2/x4,其中x1,x2,x3,x4均为第3题所给的数.

5. 计算球体积要使相对误差限为1％,问度量半径R时允许的相对误差限是多少 6. 设Y028,按递推公式

YnYn11783100 ( n=1,2,…)

Y计算到Y100.若取783≈(五位有效数字),试问计算100将有多大误差

27. 求方程x56x10的两个根,使它至少具有四位有效数字(783≈.

8. 当N充分大时,怎样求

N1dx1x2

29. 正方形的边长大约为100㎝,应怎样测量才能使其面积误差不超过1㎝

10. 设

S12gt2假定g是准确的,而对t的测量有±秒的误差,证明当t增加时S的绝对误

差增加,而相对误差却减小. 11. 序列

{yn}满足递推关系yn10yn11(n=1,2,…),若y021.41(三位有效数字),

计算到

y10时误差有多大这个计算过程稳定吗

6f(21)12. 计算,取21.4,利用下列等式计算,哪一个得到的结果最好

113,(322),,99702.63(21)(322) 2f(x)ln(xx1),求f(30)的值.若开平方用六位函数表,问求对数时误差有多大13.

若改用另一等价公式

计算,求对数时误差有多大

ln(xx21)ln(xx21)

14. 试用消元法解方程组

x11010x21010;x1x22.假定只用三位数计算,问结果是否可靠

15. 已知三角形面积

s1absinc,0c22,且测量a 其中c为弧度,,b ,c 的误差分别为

a,b,c.证明面积的误差s满足

sabc.sabc

第二章 插值法

1. 根据定义的范德蒙行列式,令

1Vn(x)Vn(x0,x1,,xn1,x)11x0xn1x2x0nx02xn1x2nxn1xn

证明Vn(x)是n次多项式,它的根是x0,,xn1,且

Vn(x)Vn1(x0,x1,,xn1)(xx0)(xxn1).

2. 当x= 1 , -1 , 2 时, f(x)= 0 , -3 , 4 ,求f(x)的二次插值多项式. 3. 给出f(x)=ln x 的数值表用线性插值及二次插值计算ln 的近似值.

x lnx

4. 给出cos x,0°≤x ≤90°的函数表,步长h =1′=(1/60)°,若函数表具有5位有效数

字,研究用线性插值求cos x 近似值时的总误差界.

maxl2(x)xxkhxk05. 设,k=0,1,2,3,求0xx3.

6. 设

xj为互异节点(j=0,1,…,n),求证:

i) ii) 7. 设

xl(x)xkjjj0nnk(k0,1,,n);

(xj0jx)klj(x)k1,2,2,n).

f(x)Ca,b且f(a)f(b)0,求证maxaxb1f(x)(ba)2maxf(x).8axb

xxf(x)e4x48. 在上给出的等距节点函数表,若用二次插值求e的近似值,要使截

断误差不超过10,问使用函数表的步长h应取多少

44nyy2nn9. 若,求及yn.

610. 如果f(x)是m次多项式,记f(x)f(xh)f(x),证明f(x)的k阶差分

kf(x)(0km)是mk次多项式,并且mlf(x)0(l为正整数).

11. 证明(fkgk)fkgkgk1fk.

12. 证明k0n1fgkn1kfngnf0g0gk1fk.k0n1

13. 证明

j02yjyny0.

14. 若f(x)a0a1xan1xn1anxn有n个不同实根x1,x2,,xn,证明

f(x)j1jnxkj0,0kn2;1an,kn1.15. 证明n阶均差有下列性质: i) 若F(x)cf(x),则

Fx0,x1,,xncfx0,x1,,xn;

,xngx0,x1,,28.

Fx,x,ii) 若F(x)f(x)g(x),则010174f2,2,f(x)xx3x116. ,求,xnfx0,x1,,2701f2,2,及,xn.

17. 证明两点三次埃尔米特插值余项是

R3(x)f(4)()(xxk)2(xxk1)2/4!,(xk,xk1)

并由此求出分段三次埃尔米特插值的误差限.

18. 求一个次数不高于4次的多项式P(x),使它满足P(0)P(k1)并由此求出分段三次

埃尔米特插值的误差限.

19. 试求出一个最高次数不高于4次的函数多项式P(x),以便使它能够满足以下边界条件

P(0)P(0)0,P(1)P(1)1,P(2)1.

20. 设

f(x)Ca,b,把a,b分为n等分,试构造一个台阶形的零次分段插值函数n(x)a,b上一致收敛到f(x).

并证明当n时,n(x)在2f(x)1/(1x),在5x5上取n10,按等距节点求分段线性插值函数Ih(x),21. 设

计算各节点间中点处的Ih(x)与f(x)的值,并估计误差.

a,b上的分段线性插值函数Ih(x),并估计误差.

22. 求f(x)x在2a,b上的分段埃尔米特插值,并估计误差.

23. 求f(x)x在424. 给定数据表如下:

xj yj试求三次样条插值S(x)并满足条件

S(0.25)1.0000,S(0.53)0.6868; ii) S(0.25)S(0.53)0.

2f(x)Ca,b,S(x)是三次样条函数,证明 25. 若

i) i)

f(x)dxS(x)dxf(x)S(x)dx2aaab2b2b2baS(x)f(x)S(x)dx;

ii) 若f(xi)S(xi)(i0,1,,n),式中xi为插值节点,且ax0x1xn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)可用式的

表达式).

第三章 函数逼近与计算

1. (a)利用区间变换推出区间为a,b的伯恩斯坦多项式.

0,/2上求1次和三次伯恩斯坦多项式并画出图形,并与相应的

(b)对f(x)sinx在马克劳林级数部分和误差做比较. 2. 求证:

(a)当mf(x)M时,mBn(f,x)M. (b)当f(x)x时,Bn(f,x)x.

0,2的最佳一致逼近多项式.

3. 在次数不超过6的多项式中,求f(x)sin4x在a,b上连续,求f(x)的零次最佳一致逼近多项式.

4. 假设f(x)在5. 选取常数a,使0x1maxx3ax达到极小,又问这个解是否唯一

0,/2上的最佳一次逼近多项式,并估计误差.

6. 求f(x)sinx在x0,1f(x)e7. 求在上的最佳一次逼近多项式.

2p(x)xr在1,1上与零偏差最小r是否唯一 r8. 如何选取,使

43f(x)x3x1,在0,1上求三次最佳逼近多项式. 9. 设

***T(x)T(2x1),x0,1T(x),T(x),T(x),T3(x). nn01210. 令,求

T(x)是在0,1上带权11. 试证

*n1xx2的正交多项式.

12. 在1,1上利用插值极小化求1f(x)tg1x的三次近似最佳逼近多项式.

x1,1上的插值极小化近似最佳逼近多项式为Ln(x),若fLnf(x)e13. 设在有界,

证明对任何n1,存在常数n、n,使

nTn1(x)f(x)Ln(x)nTn1(x)(1x1).

112331541655(x)1xxxxx1,128243843840,试将(x)降低到3次14. 设在上

多项式并估计误差.

15. 在1,1上利用幂级数项数求f(x)sinx的3次逼近多项式,使误差不超过.

a,a上的连续奇(偶)函数,证明不管n是奇数或偶数,f(x)的最佳逼近多项

16. f(x)是*Fn式(x)Hn也是奇(偶)函数.

b使axbsinxdx为最小.并与1题及6题的一次逼近多项式误差作比较. 17. 求a、

2021g(x)Ca,b,定义 f(x)18. 、

(a)(f,g)f(x)g(x)dx;(b)(f,g)f(x)g(x)dxf(a)g(a);aabb 问它们是否构成内积

1

x6dx01x19. 用许瓦兹不等式估计的上界,并用积分中值定理估计同一积分的上下界,并比

较其结果.

20. 选择a,使下列积分取得最小值:21. 设空间

11(xax)dx,xax2dx1221.

span1,x,2spanx100,x101,分别在1、2上求出一个元素,使得其

2xC0,1的最佳平方逼近,并比较其结果. 为

1span1,x2,x4f(x)x1,122. 在上,求在上的最佳平方逼近.

un(x)23.

sin(n1)arccosx1x2是第二类切比雪夫多项式,证明它有递推关系

un1x2xunxun1x.

1f(x)sinx2在1,1上按勒让德多项式及切比雪夫多项式展开,求三次最佳平方24. 将

逼近多项式并画出误差图形,再计算均方误差.

1,1上展成切比雪夫级数.

25. 把f(x)arccosx在2yabx26. 用最小二乘法求一个形如的经验公式,使它与下列数据拟合,并求均方误差.

xi yi 19 25 31 38 44 27. 观测物体的直线运动,得出以下数据:

时间t(秒) 0 10 30 50 80 110 距离s(米) 0 求运动方程. 28. 在某化学反应里,根据实验所得分解物的浓度与时间关系如下: 时间 浓度 0 5 0 10 15 20 25 30 35 40 45 50 55 用最小二乘拟合求yf(t).

29. 编出用正交多项式做最小二乘拟合的程序框图. 30. 编出改进FFT算法的程序框图. 31. 现给出一张记录xk4,3,2,1,0,1,2,3,试用改进FFT算法求出序列xk的离散频谱

Ck(k0,1,,7).

第四章 数值积分与数值微分

1. 确定下列求积公式中的待定参数,使其代数精度尽量高,并指明所构造出的求积公式所

具有的代数精度:

(1)

hh2hf(x)dxA1f(h)A0f(0)A1f(h);

(2)2h1f(x)dxA1f(h)A0f(0)A1f(h);

(3)

1hf(x)dxf(1)2f(x1)3f(x2)/3;

(4)

0f(x)dxhf(0)f(h)/1ah2f(0)f(h).

2. 分别用梯形公式和辛普森公式计算下列积分:

1(1e)xdx,n10dx,n8004x2x(1); (2);

11x2(3)91xdx,n4; (4)60sin2dx,n6.

3. 直接验证柯特斯公式具有5次代数精度.

e4. 用辛普森公式求积分01xdx并计算误差.

5. 推导下列三种矩形求积公式:

(1)

baf(x)dx(ba)f(a)f()(ba)22; f()(ba)22;

b(2)

af(x)dx(ba)f(b)b(3)

af(x)dx(ba)f(abf())(ba)3224.

6. 证明梯形公式和辛普森公式当n时收敛到积分

baf(x)dx.

7. 用复化梯形公式求积分

baf(x)dx,问要将积分区间a,b分成多少等分,才能保证误差

不超过(设不计舍入误差)

28. 用龙贝格方法计算积分10exdx,要求误差不超过10.

5cSa21()2sin2d0a9. 卫星轨道是一个椭圆,椭圆周长的计算公式是,这里a是椭

圆的半长轴,c是地球中心与轨道中心(椭圆中心)的距离,记h为近地点距离,H为远地点距离,R6371公里为地球半径,则a(2RHh)/2,c(Hh)/2.我国第一颗人造卫星近地点距离h439公里,远地点距离H2384公里,试求卫星轨道的周长.

10. 证明等式

nsinn33!n255!n4试依据nsin(/n)(n3,6,12)的值,用外推

算法求的近似值.

11. 用下列方法计算积分

(1) 龙贝格方法;

31dyy并比较结果.

(2) 三点及五点高斯公式;

(3) 将积分区间分为四等分,用复化两点高斯公式.

f(x)12. 用三点公式和五点公式分别求

差.f(x)的值由下表给出:

1(1x)2在x,和处的导数值,并估计误

x f(x) 第五章 常微分方程数值解法

1. 就初值问题yaxb,y(0)0分别导出尤拉方法和改进的尤拉方法的近似解的表达

式，并与准确解

y12axbx2相比较。

2. 用改进的尤拉方法解初值问题

yxy,0x1;y(0)1,

x取步长h=计算，并与准确解yx12e相比较。

3. 用改进的尤拉方法解

yx2xy;y(0)0,

x2取步长h=计算y(0.5)，并与准确解yexx1相比较。

4. 用梯形方法解初值问题

证明其近似解为

yy0;y(0)1,

n2hyn,2h

xyeh0并证明当时，它原初值问题的准确解。

5. 利用尤拉方法计算积分

在点x0.5,1,1.5,2的近似值。

x0edtt2

6. 取h=，用四阶经典的龙格－库塔方法求解下列初值问题：

yxy,0x1; 1）y(0)1, y3y/(1x),0x1; 2）y(0)1.

7. 证明对任意参数t，下列龙格－库塔公式是二阶的：

hyy(K2K3);nn12Kf(x,y);nn1K2f(xnth,ynthK1);K3f(xn(1t)h,yn(1t)hK1).8. 证明下列两种龙格－库塔方法是三阶的：

hyy(K13K3);nn14K1f(xn,yn);hhKf(x,yK1);nn233Kf(x2h,y2hK);nn23331）  hyy(2K13K24K3);nn19Kf(xn,yn);1hhKf(x,yK1);nn222Kf(x3h,y3hK).nn23442） 

9. 分别用二阶显式亚当姆斯方法和二阶隐式亚当姆斯方法解下列初值问题：

y1y,y(0)0,

xh0.2,y0,y0.181,y1ey(1.0)01取计算并与准确解相比较。

10. 证明解yf(x,y)的下列差分公式

yn1是二阶的，并求出截断误差的首项。 11. 导出具有下列形式的三阶方法：

1h1yn3yn1)(ynyn1)(4yn24

12. 将下列方程化为一阶方程组：

b1yn1b2yn2).yn1a0yna1yn1a2yn2h(b0yn

y3y2y0,1）y(0)1,y(0)1;

y0.1(1y2)yy0,2）y(0)1,y(0)0;

x(t)xy22,y(t),rxy,33rr

3）

x(0)0.4,x(0)0,y(0)0,y(0)2.

13. 取h=，用差分方法解边值问题

14. 对方程yf(x,y)可建立差分公式

yy0;y(0)0,y(1)1.68.

试用这一公式求解初值问题

yn12ynyn1h2f(xn,yn),

验证计算解恒等于准确解

y1;y(0)y(1)0,

x2xy(x).2

15. 取h=用差分方法解边值问题

(1x2)yxy3y6x3;y(0)y(0)1,y(1)2.

第六章 方程求根

21. 用二分法求方程xx10的正根，要求误差<。

2. 用比例求根法求f(x)1xsinx0在区间[0,1]内的一个根，直到近似根xk满足精度|f(xk)|0.005时终止计算。

32x1.5附近的一个根，设将方程改写成下列等价形式，并

3. 为求方程xx10在0建立相应的迭代公式。

22x11/xk1kx11/x1），迭代公式；

232xk131xkx1x2），迭代公式；

3）

x21x1，迭代公式xk11/xk1。

试分析每种迭代公式的收敛性，并选取一种公式求出具有四位有效数字的近似根。

x4. 比较求e10x20的根到三位小数所需的计算量；

1）在区间[0,1]内用二分法；

xkx(2e)/10，取初值x00。 k12) 用迭代法

5. 给定函数f(x)，设对一切x,f(x)存在且0mf(x)M，证明对于范围内

02/M的任意定数λ，迭代过程xk1xkf(xk)均收敛于f(x)的根x。

6. 已知x(x)在区间[a,b]内只有一根，而当a|(x)|k1，

试问如何将x(x)化为适于迭代的形式

将xtgx化为适于迭代的形式，并求x=（弧度）附近的根。

37. 用下列方法求f(x)x3x10在x02附近的根。根的准确值x＝1.…，要求

计算结果准确到四位有效数字。 1) 用牛顿法；

2）用弦截法，取x01,x11.9； 3）用抛物线法，取x01,x13,x22。 8. 用二分法和牛顿法求xtgx0的最小正根。 9. 研究求a的牛顿公式

xk1证明对一切k1,2,,xk1a(xk),x00,2xk

a且序列x1,x2,是递减的。

10. 对于f(x)0的牛顿公式xk1xkf(xk)/f(xk)，证明

Rk(xkxk1)/(xk1xk2)2

收敛到f(x)/(2f(x))，这里x为f(x)0的根。

11. 试就下列函数讨论牛顿法的收敛性和收敛速度：

x,x0;f(x)x,x0; 1)

23x,x0;f(x)23x,x0. 2)

3212. 应用牛顿法于方程xa0，导出求立方根a的迭代公式，并讨论其收敛性。

13. 应用牛顿法于方程值。

f(x)1a02x，导出求a的迭代公式，并用此公式求115的

14. 应用牛顿法于方程f(x)xa0和式，并求

nf(x)1a0nxn，分别导出求a的迭代公

lim(naxk1)/(naxk)2.k15. 证明迭代公式

xk1x(x3a)kk23xka

2是计算a的三阶方法。假定初值x0充分靠近根x，求

lim(axk1)/(axk)3.k

第七章 解线性方程组的直接方法

1. 考虑方程组：

0.4096x10.1234x20.2246x0.3872x120.3645x10.1920x20.1784x10.4002x20.3678x30.2943x40.4043;0.4015x30.1129x40.1550;0.3781x30.0643x40.4240;0.2786x30.3927x40.2557;

(a) 用高斯消去法解此方程组（用四位小数计算），

(b) 用列主元消去法解上述方程组并且与(a)比较结果。

2. (a) 设A是对称阵且a110，经过高斯消去法一步后，A约化为

a110Ta1A2

证明A2是对称矩阵。

(b)用高斯消去法解对称方程组：

4. 设A为n阶非奇异矩阵且有分解式A=LU，其中L为单位下三角阵，U为上三角阵，求证A的所有顺序主子式均不为零。

5. 由高斯消去法说明当i0(i1,2,,n1)时，则A=LU，其中L为单位下三角阵，U 为上三角阵。

0.6428x10.3475x20.8468x30.4127;0.3475x11.8423x20.4759x31.7321;0.8468x0.4759x1.2147x0.8621.123

|aii||aij|(i1,2,,n),6. 设A 为n阶矩阵，如果

j1jin称A为对角优势阵。证明：若A是

对角优势阵，经过高斯消去法一步后，A具有形式

a110Ta1A2。

7. 设A是对称正定矩阵，经过高斯消去法一步后，A约化为

a110Ta1A2，

A(aij)n,A2(aij其中

(2))n1;

证明 （1）A的对角元素aii0(i1,2,,n); （2）A2是对称正定矩阵；

(n)aaii,(i1,2,,n); n（3）

（4）A的绝对值最大的元素必在对角线上； （5）2i,jn(2)max|aij|max|aij|;2i,jn

（6）从（2），（3），（5）推出，如果

|aij|1，则对所有k

(k)|aij|1.

8. 设Lk为指标为k的初等下三角阵，即

11Lkmk1,kmnk11（除第k列对角元下元素外，和单位阵I相同）

~LIijLkIij也是一个指标为k的初等下三角阵，其中Iij为初等排

求证当i,jk时，k列阵。

9. 试推导矩阵A的Crout分解A=LU的计算公式，其中L为下三角阵，U为单位上三角阵。 10. 设Uxd，其中U为三角矩阵。

(a) 就U为上及下三角矩阵推导一般的求解公式，病写出算法。 (b) 计算解三角形方程组Uxd的乘除法次数。 (c) 设U为非奇异阵，试推导求U1的计算公式。

111. 证明（a）如果A是对称正定阵，则A也是正定阵；

（b）如果A是对称正定阵，则A可唯一写成ALL,其中L是具有正对角元的下三角阵。 12. 用高斯－约当方法求A的逆阵：

T21313107A12421015

13. 用追赶法解三对角方程组Axb，其中

14. 用改进的平方根法解方程组

210001121000A01210,b0001210000120

211x14123x5.2311x36

15. 下述矩阵能否分解为LU（其中L为单位下三角阵，U为上三角阵）若能分解，那么分解是否唯一

16. 试划出部分选主元素三角分解法框图，并且用此法解方程组

123111126,B221,C2515.A24146733161546

034x11111x22212x33.

17. 如果方阵A 有

aij0(|ij|t)，则称A为带宽2t+1的带状矩阵，设A满足三角分解

条件，试推导ALU的计算公式，对r1,2,,n.

uriari1）

rkkikmax(1,it)r1lr1u (ir,r1,,min(n,rt))；

lir(air2）18. 设

ikkrkmax(1,it)lu)/urr (ir1,,min(n,rt)).

0.60.5A0.10.3，

计算A的行范数，列范数，2-范数及F-范数。 19. 求证

(a) ||x||||x||1n||x||，

1(b)

n||A||F||A||2c2||A||Fnn。

n20. 设 PR且非奇异，又设||x||为R上一向量范数，定义

||x||p||Px||试证明

。

是R上的一种向量范数。 为对称正定阵，定义

n||x||p21. 设ARnn||x||A(Ax,x)1/2，

n试证明||x||A为R上向量的一种范数。 nTxR,x(xx,,x)12n22. 设，求证

lim(||xi||p)1/pmaxxi||x||yi11inn。

T23. 证明：当且尽当x和y线性相关且xy0时，才有

||xy||2||x||2||y||2。

24. 分别描述R中（画图）

2Sv{x|||x||v1,xR2},(v1,2,)。

25. 令

是Rn（或Cn）上的任意一种范数，而P是任意非奇异实（或复）矩阵，定义范

1||A||||PAP||。 ||x||||Px||数，证明

nn26. 设||A||s,||A||t为R上任意两种矩阵算子范数，证明存在常数c1,c20，使对一切

ARnn满足

c1||A||s||A||tc2||A||s

27. 设ARnnTT，求证AA与AA特征值相等，即求证(AA)(AA)。

TT28. 设A为非奇异矩阵，求证

1||A1||miny0||A||||y||。

11||A||||A||1(AA)29. 设A为非奇异矩阵，且，求证存在且有估计

30. 矩阵第一行乘以一数，成为

||A||||A1(AA)1||||A||.1||A||||A||1cond(A)||A||

cond(A)2A11。

证明当

23时，cond(A)有最小值。

TT1/2T31. 设A为对称正定矩阵，且其分解为ALDLWW，其中WDL，求证

2cond(A)[cond()]; 22(a)

(b) cond(A2)cond()2cond()2.

T32. 设

10099A9998

计算A的条件数。cond(A)v(v2,)

33. 证明：如果A是正交阵，则cond(A)21。 34. 设A,BRnn且

为上矩阵的算子范数，证明

cond(AB)cond(A)cond(B)。

第八章 解方程组的迭代法

1. 设方程组

(a) 考察用雅可比迭代法,高斯-塞德尔迭代法解此方程组的收敛性; (b) 用雅可比迭代法,高斯-塞德尔迭代法解此方程组,要求当||x终止．

(k1)5x12x2x312x14x22x3202x3x10x3231

x(k)||104时迭代

00A20, 证明:即使||A||1||A||1级数IAA2Ak也收敛． 2. 设

3. 证明对于任意选择的A, 序列

收敛于零． ４. 设方程组

111I,A,A2,A3,A4,23!4!

迭代公式为

a11x1a12x2b1;a21x1a22x2b2; (a11,a120);

1(k)(k1)x(bax);11122a11x(k)1(bax(k1));22211a22  (k1,2,).

(k){x}收敛的充要条件是 求证: 由上述迭代公式产生的向量序列

r5. 设方程组

a12a211.a11a22

x10.4x20.4x31x12x22x310.4x1x20.8x32x1x2x310.4x0.8xx32x2xx112323(a) (b) 1

试考察解此方程组的雅可比迭代法及高斯-塞德尔迭代法的收敛性。 6. 求证klimAkA的充要条件是对任何向量x，都有

7. 设Axb，其中A对称正定，问解此方程组的雅可比迭代法是否一定收敛试考察习题

klimAkxAx.5(a)方程组。 8. 设方程组

111xxx143442;x1x1x1;2434421x1xx1;341422111x1x2x4.42 4(a) 求解此方程组的雅可比迭代法的迭代矩阵B0的谱半径；

(b) 求解此方程组的高斯－塞德尔迭代法的迭代矩阵的谱半径； (c) 考察解此方程组的雅可比迭代法及高斯－塞德尔迭代法的收敛性。 9. 用SOR方法解方程组（分别取松弛因子1.03,1,1.1）

11x(,1,)T,(k)622要求当||xx||510时迭代终止，并且对每一个精确解

值确定迭代次数。

10. 用SOR方法解方程组（取＝）

4x1x21;x14x2x34;x4x3.32

5x12x2x312;x14x22x320;2x3x10x3.231

要求当||x(k1)x(k)||104时迭代终止。

11. 设有方程组Axb，其中A为对称正定阵，迭代公式

0试证明当

2x(k1)x(k)(bAx(k)), (k0,1,2,)

时上述迭代法收敛（其中0(A)）。

(k1)12. 用高斯－塞德尔方法解Axb，用xi记x(k1)的第i个分量，且

ri(k1)biaijxj1i1(k1)jaijxi(k)jin。

x(a) 证明 (b) 如果(k)(k1)ix(k)iri(k1)ai；

x(k)x，其中x是方程组的精确解，求证：

ri其中 (c) 设A是对称的，二次型

(k1)(k1)ii1(k)iri(k1)aii

aiji(k)jinaijj1(k1)j。

Q((k))(A(k),(k))

证明

Q((k1))Q((k))j1n(rj(k1))2ajj。

(d) 由此推出，如果A是具有正对角元素的非奇异矩阵，且高斯－塞德尔方法对任意初始向

量x(0)是收敛的，则A是正定阵。

13. 设A与B为n阶矩阵，A为非奇异，考虑解方程组

Az1Bz2b1,Bz1Az2b2,

nz,z,d,dR1212其中。

(a) 找出下列迭代方法收敛的充要条件

(m)(m1)Az1(m1)b1Bz2,Az2b2Bz1(m)(m0);

(b) 找出下列迭代方法收敛的充要条件

比较两个方法的收敛速度。 14. 证明矩阵

(m)(m1)Az1(m1)b1Bz2,Az2b2Bz1(m1)(m0);

1aaAa1aaa1

111a1a2是收敛的。 对于2是正定的，而雅可比迭代只对25102A310315. 设

3042107，试说明A为可约矩阵。 2(k1)(k)nnxCxgCR(k0,1,2,)，试证明：如果C的16. 给定迭代过程，，其中

特征值i(C)0(i1,2,)，则迭代过程最多迭代n次收敛于方程组的解。 17. 画出SOR迭代法的框图。

18. 设A为不可约弱对角优势阵且01，求证：解Axb的SOR方法收敛。 19. 设Axb，其中A为非奇异阵。 (a) 求证AA为对称正定阵；

T2cond(AA)(cond(A))22(b) 求证。

T第九章 矩阵的特征值与特征向量计算

1. 用幂法计算下列矩阵的主特征值及对应的特征向量：

327343463A1341A231213 , (b) 3, (a)

当特征值有3位小数稳定时迭代终止。 2. 方阵T分块形式为

T11T12T1nTT222nTTnn, 其中Tii(i1,2,,n)为方阵，T称为块上三角阵，如果对角块的阶数至多不超过2，则称T 为准三角形形式，用(T)记矩阵T的特征值集合，证明

(T)(Tii).i1n3. 利用反幂法求矩阵

的最接近于6的特征值及对应的特征向量。 4. 求矩阵

621231111

与特征值4对应的特征向量。 5. 用雅可比方法计算

400031013

的全部特征值及特征向量，用此计算结果给出例3的关于p的最优值。

6. (a)设A是对称矩阵，λ和x(||x||21)是A的一个特征值及相应的特征向量，又设P为一个正交阵，使

1.01.00.5A1.01.00.250.50.252.0

Pxe1(1,0,,0)T

证明BPAP的第一行和第一列除了λ外其余元素均为零。 (b)对于矩阵

T2102A10582811,

212x,,333是相应于9的特征向量，试求一初等反射阵P，使λ=9是其特征值，

TPxe1，并计算BPAPT。

7. 利用初等反射阵将

正交相似约化为对称三对角阵。 8. 设ARnn134A312421 Pij(2)a0的平面旋转阵，试推导计算PijA第i为使j1，且

ai1,aj1不全为零，

TAP行，第j行元素公式及ij第i列，第j列元素的计算公式。

9. 设An1是由豪斯荷尔德方法得到的矩阵，又设y是An1的一个特征向量。

1P2Pn2y； (a)证明矩阵A对应的特征向量是xP(b)对于给出的y应如何计算x 10. 用带位移的QR方法计算

120310121BA211011 013, (b) (a)

全部特征值。

11. 试用初等反射阵A分解为QR，其中Q为正交阵，R为上三角阵，

111A211245。

数值分析习题答案

第一章 绪论习题参考答案

(x*)1． ε（lnx）≈

x*r(x*)。

*n1r(x)2．

n(x)x*nnnx(x*)nx*n(x*)0.02n*x。

****xxxx12343． 有5位有效数字，有2位有效数字，有4位有效数字，有5位有效

*x5数字，有2位有效数字。

******4333(xxx)(x)(x)(x)0.5100.5100.5101.05101241244．

************(x1x2x3)x2x3(x1)x1x3(x2)x1x2(x3)0.214790825**x2x21**(*)*(x2)*2(x4)8.855668106x4x4x4。

r(R)r(35．

3V)431(V)/236V33V1(V)1r(V)0.00333343V3。

6．

(Y100)10011110310310022。

7． x12878355.982，

x22878312878310.0178655.982。

1dxarctgNN1x228．

11(x)(S)S2(S)0.00529． 。

r(S)10．

(S)gt(t)0.1gt，

gt(t)2(t)0.212ttgt2，故t增加时S的

绝对误差增加，相对误差减小。

11． 12．

(y10)1010(y0)10812，计算过程不稳定。

6f(21)60.005051，如果令21.4，则f1(21)0.004096，

f2110.005233f0.0051254363f(322)0.008(21)(322)，3，，

f5997021，f4的结果最好。

13．

f(30)4.094622，开平方时用六位函数表计算所得的误差为

104中

12，分别代入等价公式f1(x)ln(xx1),f2(x)ln(xx1)22计算可得

(f1)ln(1xx21)1(xx21)6010431032xx21，

(f2)ln(1xx21)xx21111048.33107602。

14．

方程组的真解为

x110000000009999999981.000000,x21.000000999999999999999999，

而无论用方程一还是方程二代入消元均解得x11.00,x21.00，结果十分可靠。

tancc15．

第二章 插值法习题参考答案

n1i0sbsincaasincbabcosccabcsabsincabc

1.

Vn(x)(xxi)0jin1(xixj)；

Vn1(x0,x1,,xn1)0jin1(xixj).

2.

L2(x)0(x1)(x2)(x1)(x2)(x1)(x1)(3)4(11)(12)(11)(12)(21)(21)

5237xx623.

3. 线性插值：取x00.5,x10.6,y00.693147,y10.510826，则

ln0.54L1(0.54)y0y1y0(0.54x0)0.620219x1x0；

二次插值：取

x00.4,x10.5,x20.6,y00.916291,y10.693147,y20.510826，则

ln0.54L2(0.54)

(0.54x0)(0.54x2)(0.54x0)(0.54x1)(0.54x1)(0.54x2)y0y1y2(x0x1)(x0x2)(x1x0)(x1x2)(x2x0)(x2x1)

＝－ .

R1(x)f(x)L1(x)1f()(xx0)(xx1)2，其中[x0,x1].

4.

所以总误差界

|R1(x)|1max|cos(x)|max|(xx0)(xx1)|xx0xx120xx1

2(x1x0)2111811.061024860180 . l2(x)(xx0)(xx1)(xx3)(x2x0)(x2x1)(x2x3)

5.

xx0当

47h3 时，取得最大值 107727 .

x0xx3max|l2(x)|k6. i) 对f(x)x,(k0,1,,n)在x0,x1,,xn处进行n次拉格朗日插值，则有

xkPn(x)Rn(x)

由于f(n1)lj(x)xkji0n1f(n1)!j(n1)()(xx0)(xxn)

()0，故有i0lnk(x)xkjx.

kg(x)(xt),在x0,x1,,xn处进行n次拉格朗日插值，有 ii) 构造函数

nLn(x)(xjt)klj(x)i0.

g(n1)()n(xt)Ln(x)(xxj)(n1)!j0插值余项为 ，

k由于 g(n1)()0,(k1,2,,n).故有

(xt)Ln(x)(xjt)klj(x).ki0n

令tx,即得 i0(xnjt)klj(x)0.

7. 以a, b两点为插值节点作f(x)的一次插值多项式

L1(x)f(a)f(b)f(a)(xa)ba， f(x)L1(x)1f()(xa)(xb),[a,b]2，

据余项定理，

由于f(a)f(b)0,故

11max|f(x)|max|(xa)(xb)|(ba)2max|f(x)|.axbaxb2axb8

1R2(x)e(xx0)(xx1)(xx2),[4,4].68. 截断误差 |f(x)L1(x)||f(x)|其中 x0x1h,x2x1h, 则

xx13h3时取得最大值

4x4max|(xx0)(xx1)(xx2)|23h39 .

由题意，

|R2(x)|142e(3h3)106,69

所以，h0.006.

n1n2n2n1n1nny22,y(22)(22)2, 则可得 nn9.

yn2n1/22n1/2， 2yn(2n12n)(2n2n1)2n1，则可得

4yn2(2yn)2n.

10. 数学归纳法证

4yn2(2yn)2n2.

当k1时，f(x)f(xh)f(x)为m－1次多项式；

k假设 f(x)(0km)是m-k 次多项式，设为g(x)，则

k1f(x)g(xh)g(x)为m-(k+1)次多项式，得证。

11. 右fk(gk1gk)gk1(fk1fk)fk1gk1fkgk左

12. k0n1fn1kgkf0g1f0g0f1g2f1g1fn1gnfn1gn1,

gk0k1fkf1g1f0g1f2g2f1g2fngnfn1gn.2j

13.

yj0n1

(y2y1)(y1y0)(y3y2)(y2y1)(yn1yn)(ynyn1)

(yn1yn)(y1y0)yny0 .

14. 由于x1,x2,,xn是f(x)的n个互异的零点，所以

f(x)a0(xx1)(xx2)(xxn)

a0(xxi)a0(xxj)(xxi),i1i1ijnn对f(x)求导得

nnf(x)a0(xxi)(xxj)((xxi))i0i1ijij，

f(xj)a0(xjxi)则

i1ijn，

j1n1f(xj)a0xkjj1nxkj(xi1ijn.xi)

j0,0kn2,(n1)g(x)kg(x)x,(n1)!,kn1. 记k则 由以上两式得

j1n1f(xj)a0xkjj1ngk(xj)(xi1ijnjxi)1gk[x1,x2,,xn]a0

()0,0kn2,1gk1a0(n1)!a0,kn1.

nF(xj)F[x0,x1,,xn]j0(xjx0)(xjxj1)(xjxj1)(xjxn) 15. i)

(n1)

ncf(xj)j0(xjx0)(xjxj1)(xjxj1)(xjxn)cf[x0,x1,,xn].

ii) 证明同上。

f(7)()7!f[2,2,,2]1;7!7!16.

017f(8)()f[2,2,,2]0.7!

R3(xj)f(xj)p(xj)0,R3(xj)f(xj)p(xj)0,jk,k1.17.

018即xk,xk1均为R3(x)的二重零点。因而有形式：

R3(x)K(x)(xxk)2(xxk1)2.

22(t)f(t)p(t)K(x)(tx)(tx). kk1作辅助函数

则 (xk)0,(x)0,(xk1)0,(xk)0,(xk1)0. 由罗尔定理，存在1(xk,x),2(x,xk1),使得

(1)0,(2)0.

类似再用三次罗尔定理，存在(1,2)(xk,xk1),使得

(4)()0, 又 (4)(t)f(4)(t)4!K(x),

(4)K(x)f()4!, 可得

即 R3(x)f(4)()(xxk)2(xxk1)24!.,(xk,xk1).

18. 采用牛顿插值，作均差表：

xi f(xi) 一阶均差 二阶均差 0 1 2 0 1 1 1 0 -1/2 p(x)p(x0)(xx0)f[x0,x1](xx0)(xx1)f[x0,x1,x2]

(ABx)(xx0)(xx1)(xx2)

0xx(x1)(1/2)(ABx)x(x1)(x2)

31A,B,44 又由 p(0)0,p(1)1, 得

x2p(x)(x3)2.4所以

hba,xkakh.n 则

19. 记

xxi1xxif(xi1),x[xi,xi1].xixi1xi1xi

因为f(x)C[a,b]，所以f(x)在[a,b]上一致连续。

n(x)f(xi)当nN时，

axbhban，此时有

0in1xixxi1由定义知当n时，n(x)在[a,b]上一致收敛于f(x)。 20. Ih(x)在每个小区间[xk,xk1]上表示为

xxxximaxmax[f(x)f(xi)]i1[f(x)f(xi1)]0in1xixxi1xi1xixi1xi

xxxximaxmaxi1.0in1xixxi1xi1xixi1xi

xxxximaxmaxf(x)f(xi)i1f(xi1)0in1xixxi1xxxxi1ii1imax|f(x)n(x)|maxmax|f(x)n(x)|xxk1xxkfkfk1,(xkxxk1).xkxk1xk1xk

计算各值的C程序如下：

Ih(x)#include\"\" #include\"\" float f(float x) { return(1/(1+x*x)); }

float I(float x,float a,float b) {

return((x-b)/(a-b)*f(a)+(x-a)/(b-a)*f(b)); }

void main() { int i;

float x[11],xc,xx; x[0]=-5;

printf(\"x[0]=%f\\n\ for(i=1;i<=10;i++) { x[i]=x[i-1]+1;

printf(\"x[%d]=%f\\n\ }

for(i=0;i<10;i++) { xc=(x[i]+x[i+1])/2; I(xc,x[i],x[i+1]);

printf(\"I[%d]=%f\\n\ }

for(i=0;i<10;i++) { xx=(x[i]+x[i+1])/2; f(xx);

printf(\"f[%d]=%f\\n\

} }

21. Ih(x)在每个小区间[xk,xk1]上为

Ih(x)xxk12xxk2xkxk1(xk1xk)xxkxk1.xkxk1xk1xk

2(xk1xk)2h2|R(x)||Ih(x)f(x)||x(xk1xk)xxkxk1|.44

322. f(x)4x, 则Ih(x)在每个小区间[xk,xk1]上表示为

xxk1Ih(x)xxk1k22xxk12xk1xk4xxkxkxxkk122xxk1412xk1xkxk1xxk1xxk334(xxk)xk4(xxk1)xk1.xxxxk1kkk1|R(x)||f(4)()(xxk)2(xxk1)2/4!|

xk1xkh42xk1xk24!(xk)(xk1)4!.2216

23. h1x2x10.05, h2x3x20.09,

h3x4x30.06, h4x5x40.08.

hi1ihihi1

6yi1yiyiyi1hhhhii1i1i

则三次样条插值函数表达式为

imi1mymym(xix)3i(xxi1)3(i1i1hi)(xix)(iihi)(xxi1)6hi6hihi6hi6i) 由S(0.25)1.0000,S(0.53)0.6868，得 Si(x)10.6429,20.4,30.5714

06(f[x0,x1]1)0.276，14.3157,23.264,32.43

关于m0,m1,m2,m3,m4的方程组为

224. i) 因为f(x)C[a,b],所以

46(0.6868f[x3,x4])0.1692

右＝abb[f(x)S(x)]2dx2S(x)[f(x)S(x)]dxab

{[f(x)S(x)]22S(x)[f(x)S(x)]}dxa

(f(x)2S(x)2)dxab＝左。

ii) 由于S(x)为三次函数，故S(x)为常数，又f(xi)S(xi)，则

xi1xi[f(x)S(x)]dx0b，所以

nxi1aaS(x)[f(x)S(x)]dxi0xiS(x)[f(x)S(x)]S(x)[f(x)S(x)]dx

[S(x)(f(x)S(x))S(x)(f(x)S(x))]dxS(b)[f(b)S(b)]S(a)[f(a)S(a)]。

b第三章 函数逼近与计算习题参考答案

1. (a) 区间变换公式为

x'(ba)xa,xnx'aba，代入原公式可得新区间里的伯

kxakkBn(f,x)f((ba)a)Pk(),Pk(x)Cnx(1x)nknbak0恩斯坦多项式为;

B1(f,x)2x,B3(f,x)3x(12x)263x2(b)

2(12x)8x33，相应的麦克劳

1313P(x)x,P(x)xxR(x)0.6459641316，48林级数分别为部分和误差则为，R3(x)150.07969263840，大于伯恩斯坦多项式的误差。

mmPk(x)k0k0nnn2. mf(x)M，故

nkf()Pk(x)Bn(f,x)MPk(x)Mnk0，

nkkknkk1k1Bn(f,x)Cnx(1x)xCn(1x)(n1)(k1)x1xk0nk1当f(x)x时，。

sin4x01，

3. 对任意不超过6次的多项式g(x)，在

x(2k1),k1,8,8时，

若有g(x)sin4x1，则g(x)在0,2上至少有7个零点，这与g(x)不超过6次矛盾，所以g(x)sin4x1，g(x)0就是所求最佳一致逼近多项式。

(f,g)max(Mc,mc),Mmaxf(x),mminf(x)g(x)caxbaxb4. 设所求为，，

由47页定理4可知g(x)在a,b上至少有两个正负交错的偏差点，恰好分别

1Mc(mc),c(Mm)2为f(x)的最大值和最小值处，故由可以解得1g(x)(Mm)2即为所求。

x3a3a33a,x1()a()a1333，故，

5. 原函数与零的偏差极大值点分别为

解方程可得出唯一解

a120.636620a34。

2x2arccos20.880689,f(x2)0.7711786.

，故

cosx，得，

a0f(x2)xa120.10525722，故所求最佳一次逼近多项式为

P1(x)0.636620x0.105257，又因为两个偏差点必在区间端点，故误差限为

0xmaxsinxP1(x)P1(0)0.1052572。

x27. a1e11.71828，故由ee1可以解得x20.541325，f(x2)1.71828，

则有

a01f(x2)xa120.89406722，故所求最佳一次逼近多项式为

P1(x)1.71828x0.894067。

8. 切比雪夫多项式在1,1上对零偏差最小，

所求函数必为切比雪夫多项式的常

111p(x)T2(x)x2r22，解得唯一解 2。 数倍，

9. 作变换

x11153119tg(t)t4t3t2t22代入f(x)得1682816，则g(t)在

1,1上的三次最佳逼近多项式为

S(t)g(t)15251173T3(t)t3t2t1288168128，作逆变换t2x1代入S(t)，则

5211293P(x)5xxxf(x)在0,1上的三次最佳逼近多项式为44128。

*T1*(x)T1(2x1)2x1，T2*(x)T2(2x1)8x28x1，10. T0(x)T0(2x1)1，

T3*(x)T3(2x1)32x348x218x1，其中x0,1。

11.

1*Tn*(x)Tm(x)0xx2dx1Tn(2x1)Tm(2x1)xx20dx12Tn(x)Tm(x)dx*212T(x)n1x，故正

1交。

12. 用T4(x)的4个零点

xkcos2k1(k1,2,3,4)8做插值节点可求得三次近似最

23L(x)0.05240690.855066x0.0848212x0.0306032x3佳逼近多项式为。

xnxf(x)ef(x)e13. ，则有，其中x1,1。由拉格朗日插值的余项表达公式

可得出

f(x)Ln(x)fn1()eTn1(x)(n1)!2n(n1)!2n，令

e1en,n(n1)!2n(n1)!2n，则待证不等式成立，得证。

14. 由

泰勒级数项数节约，在

1,1，

上有

(x)M5,3(x)M5,3(x)(x)1511651T4(x)T5(x)3848384016即

1511651183321219931101T4(x)T5(x)xxx3848384016102412840961096151165131max(x)M5,3(x)0.007568361x138483840164096其中误差限为。

115f(x)sinxxx3x612015. 115P5(x)xx3x6120，取为f(x)的近似，

误差限为1x1可

以

得

maxf(x)P5(x)10.0001984137!，再对幂级数的项数进行节约就

到原函数的3次逼近多项式

M5,3(x)P5(x)115383T5(x)x31201632384，

其误差限为

1x1maxf(x)M5,3(x)1110.0007192467!12016，即为所求

*a,a上的奇函数时，设Fn(x)为原函数的最佳逼近多项式，则16. 当f(x)为Fn*(x)f(x)En*Fn*(x)f(x)Fn*(x)f(x)EnF(x)n，对有，所

*F(x)也是最佳逼近多项式，由最佳逼近多项式的唯一性，n以

Fn*(x)Fn*(x)，即Fn*(x)是奇函数。同理可证，当f(x)为a,a上的偶

*F函数时，最佳逼近多项式n(x)也是偶函数。

axbsinx17.

202dx324a22b224ab2a2b4，为使均方误差最小，

,b8243则有1218. (a)

a24b20,b24ab20，解得

a96243ba2。

(f,g)(g,f)f'(x)g'(x)dxa，

(cf,g)c(f,g)cf'(x)g'(x)dxbb，c为常

数，

(f1f2,g)(f1,g)(f2,g)f'1(x)g'(x)dxf'2(x)g'(x)dxaa，

(f,f)0，但当f(x)c时，(f,f)0,不满足定义，所以

(f,g)f'(x)g'(x)dxab不构成内积。（b）(f,g)(g,f)，(cf,g)c(f,g)，

(f1f2,g)(f1,g)(f2,g)，(f,f)0且当且仅当f(x)0时(f,f)0，满

足定义，所以(f,g)构成内积。

61x1x1112dxdx(dx)(xdx)0.196116000(1x)2012619. 1x，1x161121121210x6dx，

61x110.0714286dx0.14285701x7其中01，则14，由此可知用积分

中值定理估计比许瓦兹不等式估计更精确。

20.

11(xax2)2dx222a35，a0时最小。

11xax2dx在a1时，值为

2121a2a33a，a1时，值为1，a1时，值为33a2，a1时最小。

21. 要使

(x012axb)dx12最小，由拉格朗日乘子法可解得

a1,b16，误差为

1180，要使

0(x2ax100bx101)2dx最小，由拉格朗日乘子法可解得

a20097992009294,b53565356，误差为0.164063，前者误差小。

1(xabx2cx4)2dxx22. 1上均为偶函数，也为偶函数，则0最小，由拉格朗日

a乘子法可解得

15105950.1171875,b1.640625,c0.820312512864128。

un1(x)un1(x)23.

sin(n2)arccosxsinnarccosx1x22xun(x)，和差化积得

证。

(f,P0)sin1124. 由积分区间的对称性及勒让德多项式的奇偶性可知

1xdx02，

1311(f,P2)(x2)sinxdx01222，将原函数在此积分区间上按勒让德多项式

111(f,P)xsinxdx8sin4cos0.32507411222三次展开就可以求得，

1153111(f,P3)(x3x)sinxdx236cos432sin0.00234807122222，代入可*3S(x)0.487611P(x)0.00821825P(x)0.499938x0.0205456x313得，均方误

差为

n21371224sin2xdx(f,P)(f,P)2.44871013222112。

21f(x)Tk(x)21***Ckdxcoskdf(x)C0CkTk(x)21021xk125. ，其中。

10a10654b542.8010654a14748998b738643.0026. a，b，解方

程得a4.00955,b0.0471846，均方误差13.0346。

227. 经验公式为satbt，最小二乘法解得a2.31346,b10.65759，运动方程

2为s2.31346t10.65759。

28. 经验公式为yt/(atb)，最小二乘法解得a0.160744,b3.17914，浓度与

时间的函数关系为yt/0.160744t3.17914。

x,x,,x29.输入初始节点01m，权函数(x)0及正交多项式次数n。 k0,Pk(x)1，计算k1,k,Pk1(x)。 判断kn 否 是 计算k1,k,Pk1(x),*k1。 令kk1 F(x)n*kPk(x)k0

p输入初始数组{xk}，等分点数N2。 q1，计算m,m0,1,,(N/2)1。 判断q0(mod2) 否 是 qA(k2j)，2计算qA(k2j)，1计算A2(k2qj2q1)，k0,1,,2pq1，令qq1 A1(k2qj2q1)，k0,1,,2pq1，j0,1,,2q1。 否 j0,1,,2q11。 判断qp 是 判断q0(mod2) 否 是 CjA2(j),j0,1,,N1 CjA1(j),j0,1,,N1 30.

Ckxjexp(ikj)01,12i2,2i,32i24，j02222，C016，31.

7C1422，C20，C3422， C40，C5422，C60，C7422。

第四章 数值积分与数值微分习题参考答案

21. 1) 公式可对f(x)1,x,x均准确成立，即

解得

A1A1h4,A0h33，具有3次代数精度。

A1A0A12hhA1hA102h2A1h2A1h33 84A1A1h,A0h33，具有3次代数精度。 2)

3) x10.28990,x20.62660,或x10.68990,x20.12660. 具有2次代数精度。

4)

112，具有3次代数精度。

11234567[f(0)2(f()f()f()f()f()f()f())f(1)]288888888

2. 1)

T81[02(0.03110.06150.0906.011760.1644.1836)0.2]16 ＝

＝

S411357123[f(0)4(f()f()f()f())2(f()f()f())f(1)]648888444

1[04(0.03110.09060.14230.1836)2(0.06150.11760.1644)0.2]24 ＝  2) T101.3915

S51.4547 3) T417.2277 S217.3222 4) T61.0356 S31.0358 3. 柯特斯公式为

Cba[7f(x0)32f(x1)12f(x2)32f(x3)7f(x4)]90. xkakh,hba4.

b2345其中

6f(x)dxC验证对于f(x)1,x,x,x,x,x，a均成立，但f(x)x时不成立。

4.

S10(e4e1/2e1)6＝ baba4(4)()f()1802，

1111()4(e)()40.0003518021802。

RS所以

|RS|5. 1) 此差值型求积公式的余项为

由于xa在[a,b]上恒为正，故在[a,b]上存在一点，使

aRf()(xa)dxbRf()(xa)dxab所以有 2)

baf(x)dx(ba)f(a)bf()(ba)22。

f()(ba)2.2

Rf()(xb)dxab

f()(ba)2.2

3)

f()(xb)dxabaRf()ab2(x)dx22

f()bab2(x)dxa22

f()(ba)3.24

6. 梯形公式和辛甫森公式的余项分别为

ba2hf()12 bah4(4)RS()f()1802 Rr其中

[a,b],hban，

f(x)dx所以当n时，Rr0,RS0，即两公式均收敛到积分a，且分别为二阶和四阶收敛。

7. 设将积分区间分成n等分则应有

baba(ba)3|R|Mf()12n12n22b其中

Mmax|f(x)|axb

，

(ba)3Mn12解得。

8. 首先算出T1,T2,T4,T8，然后逐次应用3个加速公式

计算结果如下表

41S2kT2k1T2k,k0,1,233 161C2kS2k1S2k,k0,11515 641R2kC2k1C2k,k0,16363

k 0 1 T2k S2k C2k R2k 2 3

I2 所以，积分

0.632120.71327。

9. a(2RHh)/27782.5， c(Hh)/2972.5，

cS4a21()2sin2d0a所以



＝4×× ＝48728

47782.5210.01561sin2d0

（可任选一种数值积分方法，如柯特斯公式）。 10. 由泰勒展开式

x3x5sinxx3!5!

有

limnsinnsinn33!n255!n4

n由于n，用外推算法，令

T(h)1sinhh，则

111T()2.59808,T()3,T()3.10583,3 6 12 1411114111T1()T()T()3.13397T1()T()T()3.141113363331236， 6，

116111T2()T1()T1()3.141593156153，

即的近似值为。

I311.

11dyy

1) 计算结果如下表

k 0 1 2 3 T2k S2k C2k R2k

即积分I＝。

2)

311111dydtf(t)1t2yt2 ，令

三点高斯公式

5158515f()f(0)f()1.0980495995 五点高斯公式 II0.23693f(0.90618)0.47863f(0.53847)0.56889f(0)0.47863f(0.53847)0.23693f(0.90618) ＝。

3)

I311dyy

212.51311dydydydy1y1.5y22.5yy

111111111(dtdtdtdt)10.25t1.7510.25t2.2510.25t2.75410.25t1.25 133f(x)dxf()f()133对每个积分用高斯公式 ，得

3I＝。

此积分精确值为ln31.09861。 12. 三点公式：

1[3f(1.0)4f(1.1)f(1.2)]0.24720.1 1f(1.1)[f(1.0)f(1.2)]0.21720.1 f(1.0)f(1.2)1[f(1.1)f(1.3)]0.18920.1。

f(x)2(1x)3， f(x)6(1x)4， f(x)24(1x)5

h20.12|R1|f()24(11.2)51.5510333f(1.0)的误差 h20.12|R2|f()24(11.2)57.810466f(1.1)的误差

4f(1.2)的误差 |R3|6.210。

五点公式：

1[25f(1.0)48f(1.1)36f(1.2)16f(1.3)3f(1.4)]0.2483120.1 1f(1.1)[3f(1.0)10f(1.1)18f(1.2)6f(1.3)f(1.4)]0.2163120.1 1f(1.2)[f(1.0)8f(1.1)8f(1.3)f(1.4)]0.1883120.1。 f(1.0)误差分别为

R11.7103， R23.4104， R34.7104。

第五章 常微分方程数值解法习题参考答案

n(n1)2n2ahnbhah2，误差2，改进尤拉

1． 尤拉法表达式

yn1ynh(axnb)hn22yn1ynf(xn,yn)f(xnh,ynhf(xn,yn))ahnbh22法表达式，无误

差。

2． 近似解 准确解 近似解 准确解 3． 近似解 准确解 h2h2hnyn1yn(ynyn1)yn1ynyn1()y(0)122h，2h。4． ，即又由，则有

2nh2hn2h22hh22hhnlimynlim()lim(1)lime2hexnh02hh0h02h当h0时，h0。

5． 取步长h=，，f=，f(1)=，f=，f(2)=。 6． (1) 7． K1fn，

近似解 (2) ，

近似解 K2fnth(fxffy)no(h) K3fn(1t)h(fxffy)no(h)，

则

h2hynhy'ny\"nyn(fnth(fxffy)nfn(1t)h(fxffy)n)o(h2)22h2hhfn(fxffy)(2fnh(fxffy)n)o(h2)o(h2)22

D8． (1)

令

fxy，

泰勒展开可得

K1fn，

11K2fnhDfnh2D2fno(h2)318，

K3f(xn22212h,ynhK2)fnh(DhDfn)fnh2D2fno(h2)3333y9，

121hDfnh3(D2ffyD)fno(h3)26， 代入龙格-库塔

同理有

yn1ynhfn3o(h)。(2) 类似(1)展开可得K1fn，公式可得

11K2fnhDfnh2D2fno(h2)28，

K3f(xn33319h,ynhK2)fnh(DhDfn)fnh2D2fno(h2)4442y32，

121hDfnh3(D2ffyD)fno(h3)26， 代入龙格-库塔

同理有

yn1ynhfn3公式可得o(h)。

hyn1yn(2yn1yn)29． 二阶显式公式为，代入得y(1.0)0.626，二阶隐式hyn1yn(2ynyn1)2公式为，代入得y(1.0)0.633，真解为y(1.0)0.6321。

10．

h2(2)h3(3)yn1ynhyynyno(h3)26(1)n(1)n(2)n，

y'n1yhyh2(3)h2(2)h3(3)2(1)yno(h)yn1ynhynynyno(h3)226，，

y'n1yhy(1)n(2)nh2(3)5(3)yno(h2)h3yno(h3)o(h2)28，代入得，截断误

53(3)hyn8差首项为。

11．

h2(2)h3(3)yn1ynhyynyno(h3)26(1)n(1)n(2)n，

y'n1yhyh2(3)4h3(3)2(1)2(2)yno(h)yn2yn2hyn2hynyno(h3)23，，

(1)ny'n2h2(2)h3(3)yn1ynhyynyno(h3)(1)(2)2(3)2yn2hyn2hyno(h)，26，

代入待定系数的公式中可得系数之间的关系式为a0a1a21，

a12a2b0b1b21，a14a22b14b21，a18a23b112b21。

12．

(1)

y'z,z'3z2y，

，其其

中中

y(0)1,z(0)1。(2) 。

y'z,z'0.1(1y2)zyy(0)1,z(0)0x'a,a'x(x2y2)3,y'b,b'y(x2y2)3(3)，其中

x(0)0,a(0)0,y(0)0,b(0)2。

13．

用差商逼近导数的方法把原边值问题转化为等价差分法方程组可得

31y116y20,16y131y216y30,16y231y326.88，解此方程组可得y10.494380,y20.957862,y31.36148。

14．

h1,xnn，初值条件等于准确解，由数学归纳法代入差分公式中可得

22(n1)2h2(n1)h(n1)2h2(n1)hyn12ynyn11nhnh122，即

差分法求出的解恒等于准确解。

15．

差分方程6y05y11,26y054y125y21.8,29y159y227y30.6，

34y268y331y40.6,41y381y472.2，代入得y01.01487,y11.01785，

y21.07010,y31.21030,y41.51329。

第六章 方程求根

21. 令f(x)xx1，则f(0)10,f(2)10.

k 0 1 2 3 4 5 2.

ak bk xk f(xk)符号 0 1 2 2 2 1 － － ＋ ＋ － － 3. 1)

(x)112|(x)|122x，在x1.5附近，x，

迭代公式收敛。

|(x)|2x13(1x2)2/3

23 2) (x)1x，在x1.5附近，

迭代公式收敛，迭代得近似值。

(x) 3)

1x1，

|(x)|12(x1)3/2，|(1.5)|1.4141，

迭代公式发散。 4. 1) 二分14次得； 2) 迭代5次得。

5. 迭代函数(x)xf(x)， (x)1f(x)，

由已知

0f(x)M2，有0f(x)2，所以|(x)|1.

即迭代过程收敛。

1x(x)，此时 x(x)6. 将转化为

111.(x)k

在x4.5附近，xtgxtg(x)，所以迭代格式为xarctgx，迭代三次得。

|1(x)|xk13f(xk)2xk1xk2f(xk)3xk3，迭代三次得。

7. 1) 牛顿法迭代格式

2) 弦截法迭代格式

xk1xkf(xk)(xkxk1)f(xk)f(xk1)，迭代三次得。

3) 抛物线法 f(x0)3,f(x1)17,f(x2)1，

f[x1,x0]10,f[x2,x1]16,f[x2,x1,x0]6,

故 f[x2,x1]f[x2,x1,x0](x2x1)10，

xk1xk2f(xk)则

24f(xk)f[xk,xk1,xk2]，迭代三次得。

8. 最小正根为。

1a1(xk)2aaxk29. 2，即xka。

2xk11xka1a1a(1)(1)122xk2xk22axk，

即xk1xk，序列单调递减。

(x)x10. 迭代函数为

f(x)f(x)，且有

f(x)(x)x,(x)0,(x)f(x)， xkxf(x)lim2limxkxk(xx)2f(x). kk1，

1(xk1)[(xk1)(xk1)(xk2xk1)(k1)(xk2xk1)2]2

1(xk1)(xk2xk1)(k1)(xk2xk1)22

2其中k1介于xk1与xk2之间。将上式两边除以(xk1xk2)，并将(xk)在x处

xkxk1(xk1)(xk2)

泰勒展开得

xkxk1(xk1)1()k12xk1xk2 (xk1xk2)2(x)(k1)(xk1x)1(k1)2(xk1x)(xxk2)

xk1x(k1)(xk2x)21(k1)2xk1x12(xx)xxk2， k2

其中k1介于xk1与x之间。将上式两边取极限，及klimkx，klimkx，

得

xkxk11f(x)limlim(k1)2k(x2k2f(x)。 k1xk2)xk1xkf(xk)xkf(xk)，迭代格式发散。

11. 1)

2)

xk1xkf(xk)1xkf(xk)2，迭代格式收敛，且收敛到x0。

ek1xxlimk1plimkek(xx)kkk要使

lim1xk2C(C0为常数)pxk，

则p1，为一阶收敛。

312. 令f(x)xa，迭代公式为

xk133f(xk)xka2xkaxkxk22f(xk)3xk3xk。

(x)2x3ax，则(x)23a3(2)x332，所以(a)0，

又 (x)2ax4，所以

(a)2a1/30，因此迭代格式为线性收敛。 1axxf(x23k)xk3axkxkk1kf(xxkk)2a2a13.

x3k，

取a115,x010，迭代三次得11510.7238。

14. 求na的迭代公式分别为

xk(n1)xnkaxn11nxn1xk(11)xk，

k1annk xn1设迭代函数为 (x)an(11n1n1n)x，则(x)ax， n1lim(naxn2(na)knk1)/(axk)n12!2aa(1n)/2na.

15. 记迭代函数

(x)x(x23a)3x2a，则(a)a， 由上 (3x2a)(x)x33ax 两边求导得 6x(x)(3x2a)(x)3x23a

则可得 (a)0

对①式两边求二阶导数得 6(x)12x(x)(3x2a)(x)6x

则可得 (a)0

对①式两边求二阶导数得 36(x)3(6x)(x)(3x2a)(x)6

则可得 (a)32a

所以迭代公式是三阶方法，且

① lim(axk1)(axk)k3(a)3!14a.

第七章 解线性方程组的直接方法习题参考答案

(a)高斯消去法解得x10.1670,x21.6504,x32.1967,x40.4468；

列

主

元

消

去

法

解

得

1．

(b)

x10.181919,x21.66303,x32.21723,x40.446704。

2．

(a)

A2(ij)ai1,j1ai1,1a1,j1a11aj1,i1aj1,1a1,i1a11A2(ji)，故A2对称。(b) 高

斯消去法解得x14.58669,x20.631523,x32.73520。

a3．

(k)ija(a)

(k1)ijm(k1)i,k1k1,ja(t)amitatj(1)ijt1k1(r)ml,a，(b)由irirrjurj及(a)

r1的结论可得

4．

urja(r)rjarjlrkukjk1r1,

lirmira(r)ir/a(r)rr(airlikukr)/urrk1。

因为A非奇异，U的对角元uii不为零，又LU分解等价于高斯消去法，

(i)aiiuii0，由引理可知，矩阵A的顺序主子式均不为零。

5．

高斯消去法第k步等价于左乘单位下三角矩阵Lk，而顺序主子式均不为

L1A，

(n)UALnk1零保证所得矩阵对角元不为零，可进行第步消元，

1AL11LnULU。

A2iia6．

(2)iiaiia1iai1/a11aiia1iai1/a11j2,jinaijai1(a11a1i)/a11n(2)aij

j2,jinaijai1j2,jina1j/a11j2,jinaijai1a1j/a11j2,ji，则

A2是对角

优势阵，故高斯消去法与部分选主元高斯消去法对于对称的对角优势阵每一步均选取同样的主元，得出的是同样的结果。

7．

(2)(2)Taaaa/aaaa/aaaeAe0ijij1ji111ji1ij111ji，又有当iiii(1)，(2)

1eAeieb1i2时

TiTiT0a11a1Tei(ei)n1A2(ei)n10AIn10A2，故2是对称正

定矩阵，(3)

axymaxaij,xy(2)aiiaija1iai1/a11aiia1i2/a11aii，(4)若

，令

A'I1yAI1y，由于A'和A2'也是对称正定矩阵，代入

(2)2a'a'a'a'/a'a'a'/a11'0，矛盾，故A的绝对值最大的元xxxx1xx111xxx1得

素必在对角线上，(5)2i,jn称正定，

(k)(k1)aijaij(2)(2)maxaijaxxaxxmaxaij2i,jn，(6)对所有k均有Ak对

aij1。

8． 9．

11LkIijLkIijmk1,kmn,k对

1，其中mi,k与mj,k位置互换。 行

初

等

列

变

换

，

A施

(k)(k)(k1)(k)(k)mkjakj/akk,aijaijmkjaik,(i,jk1,(n)AL,mkjUkj即为所求。 令

,n)，进行n次初等列变换后，

10．

(a) 若U为n阶可逆下三角矩阵，Ukk0，则 当k1时x1d1/U11，而

xk(dkUk,ixi)/Ukki1k1当 k2,3,,n时，

，算法即从第一行开始顺序循环，

同理可知若U为n阶可逆上三角矩阵，则当k1时xndn/Unn，而当

xn1k(dn1kUn1k,nikxnik)/Un1k,n1ki2kk2,3,,n时，

，算法即从最后

一行开始逆序循环，(b)第k步循环进行k次乘除法，共进行n(n1)/2次乘除

U1ii法

1/Uii,U1ij1UikUkj/Uii,in1,n2,ki1j，

,1,ji1,i2,,n(c)。

11．

TT1TT11(a)AAIAAAA， 由此可知A也是对称矩阵，

x,(A1x)TA(A1x)xA1x0，由此可知A1也是对称正定矩阵，

(b)LnnAnn,LniAni/Lnn,Li,iAi,iL,Li,j(Ai,jLik,iLi,ik)/Li,i2i,kk1k1i1ni，

T得出唯一正对角元的下三角阵L使得ALL。

12． 13． 14． 15．

4/8510/1723/8516/1733/856/1741/8513/17A119/855/173/858/173/851/174/855/17。 x(5/6,2/3,1/2,1/3,1/6)T。 x(10/9,7/9,23/9)T。

按高斯消去法，A无法进行第二次消去，换行后可以分解，B第二次消

去可乘任意系数，分解不唯一，C可唯一分解。

002121I13A1/21003/201T711Ty(3,,2),x(,,)0210042632。 ，解得

16． 17． 18．

高斯消去法公式中去掉aij0(ijt)即可推出该公式。

A1.1,A10.8,A20.827853,AF0.842615。

x19．

(a)

1Anmaxxixix1nmaxxinx1ini11inn,(b)

i,j1Tan2ij/nT(Aii1nA)/nmax(AA)A2(Aii1nTA)i,j1an2ijA。

，，

20．

xPPx0,xP0Px0x0PxPPxPxxxyxPP(xy)PxPyxPyP，故xP是Rn上的向量范数。 (Ax,x)0,x1221．

，

AA0xAx0x0,xTA(Ax,x)x12AALTLLx22(Lx,Ly)Ly2Lx22Lx121222122故

xyA2Ly2Ly2212(Ax,x)(Ay,y)xAyA，故xA是Rn上的向量范数。

lim(xi)pi1np1/p22． 23．

maxxilim((xi/maxxi)p)1/pmaxxix1inpi11in1inn。

T充分性：若有x和y线性相关且xy0， 即xky(k0)，代入得

xy2(1k)y2x2y2；唯一性：若有xy2x2y2，由于 xy2xTx2xTyyTyxTxyTyx2y2，两边同时平方可得出

TTTxTy0，消去共同项可得xyxxyy，当且仅当x和y线性相关时等号成

立。

24． 2以上图像分别为x11，xA'max25．

y'01，xPAyPy1。

PAP1xxPAP1。

Ay'y'maxPy0maxx026． 由向量范数的相容性可知存在常数a1,a20，使得a1xsxta2xs，

nn于是令c1a1/a2>0，c2a2/a1>0，则对任意AR，均有不等式

c1Asmaxxs0a1Axa2xssmaxxt0AxxttAtmaxxs0a2Axa1xssc2As。

TTT27． 若(AA)，则x0,AAxx,就有AA(Ax)Ax，可推出

(AAT)即(ATA)(AAT)，同理可以推出(AAT)(ATA)，综合这两

TT点即可得(AA)(AA)。

28． 29． 在

1A11/maxx0A1xxmin1x1Ax0AxminyAyy0。

，故(AA)存

，

cond(A)1cond(A)1A1AA1A1，则

(IA1A)11/1A1AAAA1(AA)1A1A1(IA1A)1A1AA1A1A1A1AAA。

30．

cond(A)AA1，当2/3时，cond(A)36，当2/3时，cond(A)42/，当2/3时，cond(A)有最小值7。 31． (a)

2cond(A)2max/min(max(WTW)/min(WTW))2cond(W)2，

TTcond(WT)2max(WTW)/min(WTW)cond(W)2(WW)(WW)(b),，

cond(A)2cond(WT)2cond(W)2。

32． 33． 34．

1cond(A)2max/min39206.0，cond(A)AA39601。

cond(A)2max(ATA)max(AAT)maxImaxI1。

。

cond(AB)ABB1A1ABB1A1cond(A)cond(B)第八章 解线性方程组的迭代法习题参考答案

1. (a) Jacobi迭代矩阵

0.40.20BD1(LU)0.2500.50.20.30

3特征方程为 |IB|0.210.0550

特征根均小于1，Jacobi迭代法收敛。 Gauss-Seidel迭代矩阵

00.40.21G(DL)U00.40.700.040.17特征方程为 |IG|0.570.0960 特征根均小于1，Gauss-Seidel迭代法收敛。 (b) Jacobi迭代格式为

X(k1)BX(k)f1

1TfDb(1.250.3)1其中B如上，，

32

迭代18次得

X3.99999642.99997391.9999999T，

Gauss-Seidel迭代格式为

X(k1)GX(k)f2

1T其中G如上，f2(DL)b(2.42.61.53)，

迭代8次得

X4.0000362.9999852.000003。

T00A20， 2. 证：

2A0, 故Ak0(k2,3,4,)， 则

10IAA2AkIA21， 因此

即级数收敛。

3. 证： 设||A||a，

Anlim0nn!一方面，，

An||An||||A||nanlimlimlimlim0nn!nnnn!n!n!另一方面，

Anlim0nn!因此，即序列收敛于零。 4. 证：由已知迭代公式得迭代矩阵

a12a110 aa|IG|212210a11a22则特征多项式为

0Ga21a22解得

a12a21a11a22，

ra12a211a11a22。

向量序列x收敛的充要条件是 1，即

(k)5. (a) 谱半径(B)1.0931，Jacobi迭代法不收敛； 矩阵A对称正定，故Gauss-Seidel迭代法收敛。 (b) 谱半径(B)01，Jacobi迭代法收敛； 谱半径(B)21，Gauss-Seidel迭代法不收敛； 6. 证：必要性

limAkAk，则 klimAkA0，

对任意向量x，有

limAkxAxlim(AkA)xlimAkAx0kkk

limAkxAx0limAxAx因而有 k，即kk。

limAkxAxx充分性 因对任何向量，都有k，令xiei，则

即当k时，Ak的任一列向量的极限为A的对应的列向量，因而有

klimAkeiAeilimAkAk。

7. A对称正定，Jacobi迭代法不一定收敛，如题5(a)。

8. (a) Jacobi迭代矩阵的谱半径

(B)12；

(b) Gauss-Seidel迭代矩阵的谱半径(B)0.25； (c) 两种方法的谱半径均小于1，所以两种方法均收敛。

事实上，对于方程组Axb，矩阵A为严格对角占优则Jacobi和Gauss-Seidel迭代法均收敛。 9. 取X(0)0，迭代公式为

使当

XX(k)(k1)(k)(k)(k)XX(14XX112)14(k1)(k)(k)X2(4X1(k1)4X2X3(k))X24(k1)(k)(k1)XX(3X24X3(k))334

5106时迭代终止，

取1.03时，迭代5次达到

X(5)0.50000431.00000010.4999999T；

取1时，迭代6次达到

X(6)0.50000381.00000020.4999995T；

取1.1时，迭代6次达到

X(6)0.50000350.99999890.5000003T。

10. 迭代公式为

取X(0)(k1)(k)(k)(k)XX(125X1(k)2X2X3)115(k1)(k)(k)XX(20X1(k1)4X22X3(k))224(k1)(k)(k1)(k)XX(32X1(k1)3X210X3)3310

0，0.9，迭代8次达到精度要求

2.999892.00003。

TX(8)4.0002711. 证：所给迭代公式的迭代矩阵为BIA，

其n个特征值分别为11,12,1n,(0i,i1,,n)，

02当

时，有11i1,(i1,2,,n)，

因而(B)1，迭代法收敛。

X(k1)i12. 证：(a)

X(k)iri(k1)iia即为Gauss-Seidel迭代格式。 ri(k1)ii(k)(k)XXiia及，可得

(b) 由

X(k1)iX(k)i(k1)i(k)iri(k1)iia；

其中，

i1ri(k1)biaijxj1ni1(k1)jaijxjin(k)jaijxaijxjj1j1i1nni1(k1)jaijx(jk)jin

aij(xxjj1(k1)j)aij(xxjji(k)j)(aijj1(k1)jaij(jk))ji。

(c) (d)

13. (a) 由已知，有z1(m1)(m)(m1)A1Bz2A1b1，及z2A1Bz1(m)A1b2，

(m1)12(m1)111z(AB)zABAbAb1， 12则 1(m1)(m1)12(m1)(m)1zz(AB)zz即由1到1的迭代矩阵为，所以由1到1的迭代矩阵为AB，

1则迭代方法收敛的充要条件为(AB)1。

(b) 由已知可推得z1(m1)(A1B)2z1(m)A1BA1b2A1b1，所以迭代矩阵为

(A1B)2，则迭代方法收敛的充要条件为[(A1B)2]1。

由迭代矩阵可以看出，(b)迭代法的收敛速度是(a)的2倍。

1a11a20,a114. 证：由于10，当2时，a1

|A|(12a)(1a)0，所以A正定。

11a2收敛。 Jacobi迭代矩阵谱半径为(B)2|a|，所以只对215. 取排列阵PI23，则

53PTAP003211024037 21A为可约矩阵。

16. 证： 迭代矩阵的特征方程为|i(C)IC|0，

n(0)若i(C)0,(i1,2,,n)，则|C|0，所以C0，即对任给向量X，迭代n

(n)n(0)n1XCXfffCgcgg，则 次后，，其中

X(n1)CngCn1gcggf

即最多迭代n次收敛于方程组的解f。

17. 用SOR方法解方程组AXb，其中A对称正定，数组x用来存放解向量，用

|p0|max|xi(k1)xi(k)|1in控制迭代终止，k表示迭代次数。

k=0, i=1 xi0(i1,,n) kk1,p00 否 i<=n 是 p(biaijxjaijxj)/aii;j1jii1nxixip;ii1; |P|>|P0| P0=P; 是 |P|>ε 否 输出x, k; 1L(DL)((1)DU)， 18. 证：方程组的SOR迭代矩阵为

1|(DL)||(1)DLU|0， 特征方程

即|(1)DLU|0， 记G(1)DLU

(1)a11a21an1a12a13(1)a22a23an2an3(1)ann

a1na2n只要当||1时，|G|0，则|G|0的根均满足||1。

A不可约则G也不可约，又A为弱对角优势阵，则当||1且01时，

|gii|(1)|aii|(1)(|aij||aij|)j1ji1i1n

即||1时，G为不可约弱对角占优，于是有|G|0，故(G)1，SOR方法收敛。

j1ji1jiTTTn19. 证：(a) (AA)AA，xR,x0,设Ax(a1,,an)，则 2xT(ATA)x(Ax)T(Ax)a12an0，ATA为对称正定阵。

T|aij|i1|anij||gii|T (b) 因为AA为对称阵，所以

max(ATA)cond(AA)2T(AA) min左

T(cond(A)2)(A22A1右

(AA))maxT2min(AA)左。

T2220. 证：A为严格对角占优，则A存在。

A1maxA1xx1

第九章 矩阵的特征值与特征向量计算习题参考答案

1．(a)取初始值(1,1,1)得

k 1 3 Tuk max(vk) 5 7 8 (b)取初始值(1,1,1)得 k 1 5 9 13 14 Tuk max(vk) n TITiiI0i12．(T),x0，使得(TI)x0，即

n，一定存在

ni使得TiiI0，则(Tii)，

n(T)i1(Tii)，反之

(Tii)(T)i1，故

(T)i1(Tii)。

14/2711/275/27(A6I)111/272/271/275/271/274/271(A6I)0.750741，原max3．，由幂法得

矩阵最接近6的特征值为maxA7.33202，对应的特征向量为

xT(1,0.485551,0.185986)。

T4．设特征向量为x(a,b,c)，则有4a4a,3bc4b,b3c4c，解得对应的TTx(1,0,0),x(0,1,1)。 12特征向量为

5．雅可比迭代进行五步可得12.53652,20.0166474,31.48023，对应的特征向量分别为

Tx1(0.531632,0.461334,0.710313)，

Tx2(0.721102,0.686454,0.0938701)，

Tx3(0.444293,0.562111,0.697592)，

最优值p1.26658。

TTTBPAPPAxPxe1，PxePx11，P正交，则P第一列16．(a) ，1又BPAP是对称矩阵，B的第一行和第一列除外均为零。(b) P为反射阵，

TPxe1Pe1x，解得

2/32/31/3900P1/32/32/3BPAPT01802/32/31/3009，。

7．由豪斯荷尔德方法得

005011U00.60.8BUAU52.920.5600.80.600.560.92，。

(2)a8．j1ai1sinaj1cos0，

sin解得代入得

aj12ai21aj1,cosai12ai21aj1，

(PijA)ikaikcosajksin9.(a)

AxAPP12Pn2yP1aikai1ajkaj1aa2i12j1,(PijA)jkajkcosaiksinajkai1aikaj1aa2i12j1,Pn2PnT2TP1AP1Pn2yP1Pn2An1yP1Pn2y，(b)

TAPAn1Pn2可求出初等反射阵Pn2，依次类推，xPPn212n2由

(k)sak3310．(a)令，带位移QR方法计算可得

Pn3Pn2y。

13.37215,21.99872,32.37079 ，

(k)(b) 令ska33，带位移QR方法计算可得

13.73169,22.00036,30.267949。

105R052505500325015250355253250550325555003330A0330031，

11．

1/3AQR2/3故有

2/3

2/32/331/32/331/302/0333303。