第2章---解线性代数方程组的迭代法ppt课件.ppt

《第2章---解线性代数方程组的迭代法ppt课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第2章解线性代数方程组的迭代法求解线性代数方程组主要有直接法和迭代法两种常见方法。直接法一般适合小型的系数矩阵，为了求解现实当中常见的大型稀疏矩阵，下面我们将重点介绍迭代法。它是一种不断套用一个迭代公式，逐步逼近方程的解的方法。将讨论两类主要方法，一类是逐步逼近法，另一类是下降法，包括最速下降法和共轭梯度法。为了对线性方程组数值解的精确程度，以及方程组本身的性态进行分析，需要对向量和矩阵的“大小”引进某种度量，范数就是一种度量尺度，向量和矩阵的范数在线性方程组数值方法的研究中起着重要的作用。2.1向量、矩阵

2、范数与谱半径（2）齐次性：对任何实数和向量x

3、

4、kx

5、

6、＝

7、k

8、

9、

10、x

11、

12、（3）三角不等式：对任何向量x和y，都有

13、

14、x＋y

15、

16、≤

17、

18、x

19、

20、＋

21、

22、y

23、

24、（1）非负性：对任何向量x，

25、

26、x

27、

28、≥0，且

29、

30、x

31、

32、＝0当且仅当x＝02.2.1向量的范数定义2.1设是x的实值函数，且满足条件则称为Rn(或Cn)上的一个向量范数(或向量模)，

33、

34、x

35、

36、的值为向量x的范数。向量1-范数：向量2-范数：向量无穷范数：容易验证，以上三种范数都满足向量范数的三个条件。理论上存在多种多样的向量范数，但最常用的是如下三种。设向量

37、x＝(x1，x2，…，xn)T，定义解:对于向量x＝(1，-3，2，0)T，根据定义可以计算出：

38、

39、x

40、

41、1＝

42、1

43、＋

44、-3

45、＋

46、2

47、＋

48、0

49、＝6由此例可见，向量不同范数的值不一定相同，但这并不影响对向量大小做定性的描述，因为不同范数之间存在如下等价关系。例2.1设向量x＝(1，-3，2，0)T，求向量范数

50、

51、x

52、

53、p，P＝1，2，∞。定理2.1（范数的等价性）对于Rn上任何两种范数

54、

55、·

56、

57、α和

58、

59、·

60、

61、β，存在的正常数m，M，使得：范数的等价性表明，一个向量若按某种范数是一个小量，则它按任何一种范数也将

62、是一个小量。容易证明，常用的三种向量范数满足下述等价关系。

63、

64、x

65、

66、∞≤

67、

68、x

69、

70、1≤n

71、

72、x

73、

74、∞

75、

76、x

77、

78、∞≤

79、

80、x

81、

82、2≤

83、

84、x

85、

86、∞

87、

88、x

89、

90、∞≤

91、

92、x

93、

94、1≤

95、

96、x

97、

98、2定义2.2对于向量序列及向量如果则称向量序列x(k)收敛于向量x*。记作或2.1.2矩阵的范数矩阵范数是反映矩阵“大小”的一种度量，具体定义如下。定义2.3设为A的实值函数，且满足条件：（1）

99、

100、A

101、

102、≥0，且

103、

104、A

105、

106、＝0时，当且仅当A＝0（2）

107、

108、kA

109、

110、＝

111、k

112、

113、

114、A

115、

116、，k∈R（3）

117、

118、A＋B

119、

120、≤

121、

122、A

123、

124、＋

125、

126、B

127、

128、则

129、称上的一个矩阵范数，

130、

131、A

132、

133、的值为矩阵A的范数。设n阶矩阵A＝（aij），常用的矩阵范数有：矩阵1-范数：矩阵2-范数：矩阵无穷范数：列和行和以上三种范数都满足矩阵范数的条件，通常将这三种矩阵范数统一表示为

134、

135、A

136、

137、p，P＝1，2，∞。例2.2设矩阵求矩阵A的范数

138、

139、A

140、

141、p，P＝1，2，∞。解根据定义由于此方程的根为矩阵ATA的特征值，解得因此则它的特征方程为:在线性方程组的研究中，经常遇到矩阵与向量的乘积运算，若将矩阵范数与向量范数关联起来，将给问题的分析带来许多方便。设

142、

143、·

144、

145、是一种向量范数，由此

146、范数派生的矩阵范数定义为注意，此式左端

147、

148、A

149、

150、表示矩阵范数，而右端是向量Ax和x的范数，利用向量范数所具有的性质不难验证，由上式定义的矩阵范数满足矩阵范数的条件。通常将满足上式的矩阵范数称相容范数。由向量范数

151、

152、x

153、

154、p派生出的矩阵范数：通过向量范数定义的矩阵范数，满足不等式关系：称之为矩阵A的算子范数，其中p＝1，2或∞。定理2.2由上式所定义的矩阵范数为相容范数。证明：当x=0时，(1)式显然成立。2.1.3矩阵的谱半径矩阵范数同矩阵特征值之间有密切的联系，设λ是矩阵A相应于特征向量x的特征值，即Ax

155、＝λx,于是利用向量-矩阵范数的相容性，得到

156、λ

157、

158、

159、x

160、

161、=

162、

163、λx

164、

165、从而，对A的任何特征值λ均成立=

166、

167、Ax

168、

169、≤

170、

171、A

172、

173、

174、

175、x

176、

177、

178、λ

179、≤

180、

181、A

182、

183、（3）设n阶矩阵A的n个特征值为λ1,λ2，…λn，称为矩阵A的谱半径。从(3)式得知，对矩阵A的任何一种相容范数都有ρ(A)≤

184、

185、A

186、

187、。2.2迭代法的一般形式与收敛性定理2.2.1迭代法的一般形式已知线性代数方程组首先将方程组改写成等价的形式从而建立迭代式：称为迭代序列，并称H为迭代矩阵。则是线性方程组Ax=b的解。当给定初始向量后可得到迭代向量序

188、列，若在等式(2.2.3)两端取极限可得2.2.2迭代法的收敛性利用迭代公式(2.2.3)构造序列，以求得方程组(2.2.2)的近似解的算法称为解(2.2.2)式的简单迭代法。若迭代序列收敛，则称此迭代法是收敛的。两式相减，知误差向量满足下列迭代关系：由此递推：引理2.1：迭代法(2.2.3)式对任何初始近似均收敛的充分必要条件是引理2.2：的充要条件是定理2.4：迭代法(2.2.3)式对任何初始近似均收敛的充分