数值分析ex12-13《数值分析》习题课II

《数值分析ex12-13《数值分析》习题课II》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、高斯消元法矩阵的三角分解雅可比迭代与赛德尔迭代迭代法收敛定理最速下降法《数值分析》习题课II2/20一、高斯消元法==三角方程组解法、顺序消元法、列主元法、追赶法二、矩阵的三角分解矩阵的紧凑格式分解、改进平方根法三、向量范数和矩阵范数常用的三种向量范数、常用的三种矩阵范数、条件数四、迭代法及收敛性分析雅可比迭代法、高斯-赛德尔迭代、收敛定理、误差定理、初等变分原理定理3.1约化主元ak+1,k+1(k)≠0(k=0,1,···,n-1)的充分必要条件是矩阵A的各阶顺序主子式不为零.Ex1.如果A是严格主对角占优矩阵,则det(A)≠0.证:

2、用反证法。设det(A)=0,则齐次方程组Ax=0有非零解u=[u1,u2,···,un]T.设考虑Au=0的第k个等式3/204/20两边约去

3、uk

4、，得这与主对角占优矛盾,故det(A)≠0。Ex2.设A对称且a11≠0,经过高斯消元法一步后,A约化为证明A2也是对称矩阵。证明:设经高斯消元一步后,得5/20所以,A2=A2T思考:1.若A是对称正定矩阵,经高斯消元一步后,右下角子矩阵A2也是对称正定矩阵;2.若A为对角占优矩阵，经过高斯消元法一步后，右下角子矩阵A2也是对角占优矩阵。Ex3.对任何一种矩阵的算子范数,证明矩阵A的谱半

5、径与A的范数有关系:ρ(A)≤

6、

7、A

8、

9、证:设是矩阵A任一特征值,x是对应的特征向量,则Ex4.若矩阵A是n阶对称矩阵,则有证：设是A的任一特征值,由于A对称,故2是矩阵ATA的特征值,即6/207/20由2-范数计算公式Ex5.对任意x，y∈Rn，利用向量范数的三角形不等式证明：证:

10、

11、x

12、

13、=

14、

15、(x–y)+y

16、

17、≤

18、

19、x–y

20、

21、+

22、

23、y

24、

25、

26、

27、x

28、

29、–

30、

31、y

32、

33、≤

34、

35、x–y

36、

37、同理,

38、

39、y

40、

41、–

42、

43、x

44、

45、≤

46、

47、y–x

48、

49、=

50、

51、x–y

52、

53、

54、

55、x

56、

57、–

58、

59、y

60、

61、≥–

62、

63、x–y

64、

65、–

66、

67、x–y

68、

69、≤

70、

71、x

72、

73、–

74、

75、y

76、

77、≤

78、

79、

80、x–y

81、

82、Jacobi迭代法的迭代矩阵8/20Gauss-Seidel迭代法的矩阵:BG-S=(D–L)-1UAx=b,将矩阵分裂:A=D–U–LBJ=D-1(U+L)特征多项式与特征方程:

83、I–D-1(U+L)

84、=

85、D-1

86、·

87、D–(U+L)

88、

89、D–(U+L)

90、=0特征多项式与特征方程:

91、I–(D–L)-1U

92、=

93、(D–L)-1

94、·

95、(D–L)–U

96、

97、(D–L)–U

98、=09/20Ex6.若A是严格主对角占优矩阵，求证解方程组AX=b的高斯-赛德尔迭代法收敛。证：高斯-赛德尔迭代矩阵为(D–L)-1U,该矩阵的特征方程

99、为

100、(D–L)–U

101、=0行列式对应的矩阵为当

102、

103、>1时，利用A矩阵的主对角占优性质，得故C()也是严格主对角占优矩阵。由于严格主对角占优矩阵的行列式不为零，故不是特征方程C()=

104、(D–L)–U

105、=0的根。所以当A是严格主对角占优矩阵时,(D–L)-1U的特征值必然满足：

106、

107、<1，从而高斯-赛德尔迭代矩阵谱半径小于1，迭代法收敛。10/2011/20Ex7.证明，当

108、

109、<1时，二阶约当块的方幂Jm极限值为零。证：由于假设则有由数学归纳法知12/20而

110、

111、<1，故思考:三阶约当块的方幂Jm表达式结构13/20Ex8.设A是

112、一个可逆矩阵，矩阵序列满足Xk+1=Xk(2I–AXk)，（k=0，1，2，……）证明:当时证明：由Xk+1=Xk(2I–AXk)，得I–AXk+1=I–AXk(2I–AXk)=(I–AXk)2于是I–AXk=(I–AXk-1)2=(I–AXk-2)2×2=··········14/2015/20练习2.设A=(aij)n×n为可逆下三角矩阵，证明A-1仍为下三角矩阵。练习1.分析求解三对角方程组追赶法的计算工作量。练习3.设A=(aij)n×n为可逆上三角矩阵，证明A-1仍为上三角矩阵。练习4.用列主元法解方程组练习5:求矩阵的2-范数

113、,以及2-范数意义下的条件数16/20练习6.设A=(aij)n×n为实对称正定矩阵,x∈Rn,b∈Rn,如果u使二次函数取极小值,证明u是线性方程组Ax=b的解。练习8.有方程组Ax=b,其中A为对称正定阵,且有迭代公式讨论使迭代序列收敛的的取值范围.练习7.写出n维向量序列{X(k)}收敛于向量X*的定义；设,而B是n阶方阵证明17/20(1)A–1=B(I+R+R2+……)；(2)任意给定n阶矩阵X0,由迭代格式Xk+1=XkR+B(k=0，1，2，……)产生的矩阵序列{Xk}收敛到矩阵A-1；(3)对矩阵序列{Xk}，有误差估计式

114、18/20练习9:设A是n阶可逆矩阵,有A的一个近似逆B,令R=I–AB如果

115、

116、R

117、

118、≤q<1，试证明