几何的角度解释矩阵出现病态2

《几何的角度解释矩阵出现病态2》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、矩阵的范数范数的定义　　若X是数域K上的线性空间，泛函║·║:X->R满足：　　1.正定性：║x║≥0，且║x║=0<=>x=0；　　2.正齐次性：║AB║│≤│A║B║；　　3.次可加性(三角不等式)：║A+B║≤║A║+║B║。　　那么║·║称为X上的一个范数。　　（注意到║x+y║≤║x║+║y║中如令y=-x，再利用║-x║=║x║可以得到║x║≥0，即║x║≥0在定义中不是必要的。）在AX=B式中，当A、B有微小扰动时，参数估值Ｘ也有扰动：（A＋ΔA）（Ｘ＋ΔＸ）＝B＋ΔB（1）AΔＸ＝ΔA（Ｘ＋ΔＸ）+ΔBΔＸ＝－A

2、-1ΔA（Ｘ＋ΔＸ）＋A－１ΔB将上式两端取范数，并应用向量范数的三角不等式║AB║│≤│A║B║；║A+B║≤║A║+║B║。

3、

4、ΔＸ

5、

6、≤

7、

8、A-1

9、

10、

11、

12、ΔA

13、

14、（

15、

16、Ｘ

17、

18、＋

19、

20、ΔＸ

21、

22、）＋

23、

24、A-1

25、

26、

27、ΔB

28、

29、把含有ΔＸ的项移到式子的左边有：(1-

30、

31、A-1

32、

33、

34、

35、ΔA

36、

37、)

38、

39、ΔＸ

40、

41、≤

42、

43、A-1

44、

45、

46、

47、ΔA

48、

49、

50、

51、Ｘ

52、

53、+

54、

55、A-1

56、

57、

58、ΔB

59、

60、由于有1/

61、

62、Ｘ

63、

64、≤

65、

66、A

67、

68、/

69、

70、B

71、

72、将上式两端同时乘以1/

73、

74、Ｘ

75、

76、得：(1-

77、

78、A-1

79、

80、

81、

82、ΔA

83、

84、)

85、

86、ΔＸ

87、

88、/

89、

90、X

91、

92、≤

93、

94、A-1

95、

96、

97、

98、ΔA

99、

100、+

101、

102、

103、A-1

104、

105、

106、

107、ΔB

108、

109、A

110、

111、/

112、

113、B

114、

115、;设K=

116、

117、A

118、

119、

120、

121、A-

122、

123、将上式整理的：（1-K

124、

125、ΔA

126、

127、/

128、

129、A

130、

131、）

132、

133、ΔＸ

134、

135、/

136、

137、X

138、

139、≤K(

140、

141、ΔA

142、

143、/

144、

145、A

146、

147、+

148、

149、ΔB

150、

151、/

152、

153、B

154、

155、);即有：

156、

157、ΔＸ

158、

159、/

160、

161、X

162、

163、≤k/（1-K

164、

165、ΔA

166、

167、/

168、

169、A

170、

171、）(

172、

173、ΔA

174、

175、/

176、

177、A

178、

179、+

180、

181、ΔB

182、

183、/

184、

185、B

186、

187、);问题与实验1:试从几何的角度解释矩阵出现病态的原因，并用‘有说服力’的例子来支持你的观点；线性方程组解的敏感性的几何解释(2x2矩阵)线性方程组求解：两直线求交点下面两图分别反映了良态问题和病态问题两种情

188、况。良态情况病态情况由图像可以看出当两条直线的夹角较大时其精确解与近似解距离较小（其中两条直线的交点是其精确解）当两条直线的夹角较小时其精确解与近似解距离较大。虚线为直线的系数产生微小的扰动时的直线。下面我们用两个例子说明：例设方程组，即，它的精确解为.A=[2,6;2,6.00001;]B=[8;8.00001]C=double(AB)cond(A)norm(A)A=2.00006.00002.00006.0000B=8.00008.0000C=11ans=4.0000e+006ans=8.9443下面我们来求上面两条直线的

189、夹角：计算直线夹角根据给定的节点，计算夹角信息，示意图如下。已知点X1(x1,y1)、点X2(x2,y2)、点X3(x3,y3)，计算直线由点X2，X1与X2，X3的夹角信息。根据向量内积，得到计算公式为：theta1=acosd(dot([x1-x2,y1-y2],[x3-x2,y3-y2])/(norm([x1-x2,y1-y2])*norm([x1-x2,y3-y2])));其中，dot([x1-x2,y1-y2],[x3-x2,y3-y2])为计算内积，norm([x1-x2,y1-y2])为计算向量长度，acosd为计

190、算以度为单位的夹角信息。由这个公式计算这两条直线的夹角：theta1=acosd(dot([2,6],[2,6.00001])/(norm([2,6])*norm([2,6.00001])))theta1=2.8649e-005它的图像是X1=[0.9:0.0001:1.1];y=-1/3*X1+4/3;z=-2/6.00001*X1+8.00001/6.00001;subplot(1,2,1)holdonplot(X1,y,'b')plot(X1,z,'g')X2=[9.9:0.0001:10.1];n=-1/3*X2+4/3

191、;m=-2/5.99999*X2+8.00002/5.99999;subplot(1,2,2)holdonplot(X2,n,'r')plot(X2,m,'g')现在考虑系数矩阵和右端项的微小变化对方程组解的影响，即考察方程组，其解变为.A=[2,6;2,5.99999;]B=[8;8.00002]C=double(AB)cond(A)norm(A)A=2.00006.00002.00006.0000B=8.00008.0000C=10-2ans=4.0000e+006ans=8.9443.同样我们来计算这两条直线产生夹角：t

192、heta1=acosd(dot([2,6],[2,5.99999])/(norm([2,6])*norm([2,5.99999])))theta1=2.8649e-005从上面的结果可以看出系数矩阵的条件数较大，病态较严重，还可以看出对于系数矩阵微小的变化会对两