(完整word版)几何角度解释矩阵出现病态

1、矩阵的范数范数的定义若X是数域K上的线性空间，泛函丨丨：X->R满 足：1. 正定性：| x |> 0,且 I x | =0 <=> x=0 ;2. 正齐次性：| AB |<| A | B |;3. 次可加性(三角不等式)：A+B w A + B o那么| |称为X上的一个范数。(注意到| x+y | w | x | + | y |中如令 y=-x，再利 用| -x | = | x |可以得到| x | >0, 即卩| x | >0在定义 中不是必要的。)在AX=B式中，当A、B有微小扰动时，参数估值X也有扰动：(A + A) (X + X) = B

2、+ B (1 )A X = A (X+AX) + B X = - A-1 A (X+AX) + A-1A B将上式两端取范数，并应 用向量范数的三角不等式 | AB | 1<1 A | B | ；A+B w A + B 。X |w |A-1| A| ( |X II + MX |)+ |A-1| | B| 把含有AX的项移到式子的左边有：(1-|A-1| A|)| AX | W |A-1| A| | X |+|A-1| | B|由于有 IX |W|A|/|B|将上式两端同时乘以|X |得：(1-|A-1| A A|)| AX |/|X| W |A-1| A A| +|A-1| A B|A|

3、/|B|;设K=|A|A-|将上式整理的：(1-K| A A|/|A|)|AX |/|X| W K(| A A|/|A|+| A B|/|B|)；即有：|AX |/|X| W k/ (1-K| A A|/|A| ) (| A A|/|A|+| A B|/|B|)；问题与实验1：试从几何的角度解释矩阵出现病态的原因，并用有说服力的例子来支持你的观点；线性方程组解的敏感性的几何解释(2x2矩阵熾 性方程组求解：两直线求交点下面两图分别反映了良 态问题和病态问题两种情况。良态情况病态情况由图像可以看出当两条直线的夹角较大时其精确 解与近似解距离较小(其中两条直线的交点是其精确 解)当两条直线的夹角较

4、小时其精确解与近似解距离 较大。虚线为直线的系数产生微小的扰动时的直线。下面我们用两个例子说明：例 设方程组 Ax b ，即 2 6x18 ，它的精确解为2 6.00001 x28.00001(1,1)T. A=2,6;2,6.00001;B=8;8.00001C=double(AB)cond(A)norm(A)A =2.0000 6.00002.0000 6.0000B =8.00008.0000C = 1 1 ans = 4.0000e+006ans =8.9443下面我们来求上面两条直线的夹角：计算直线夹角 根据给定的节点，计算夹角信息，示意图如下。已知点 X1(x1,y1)、点 X2(

5、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为计算以度为单位的夹角信息。由这个公式计算这两条直线的夹角：thetal = acosd(dot(2,6,2,6.00001)/( norm(2,6)* norm(2,6.0000

6、1) thetal =2.8649e-005它的图像是0.911.11.21.3-1.95 i-1.96 -1.97 -1.98 -1.99 -2 “-2.01 -2.02 -2.03 -2.04 1cc9.99.951010.0510.1X1二0.9:0.0001:1.1;y=-1/3*X1+4/3;10 / 10z=-2/6.00001*X1+8.00001/6.00001;subplot(1,2,1)hold onplot(X1,y, 'b' )plot(X1,z, 'g' )X2=9.9:0.0001:10.1;n=-1/3*X2+4/3;m=-2/5.

7、99999*X2+8.00002/5.99999; subplot(1,2,2)hold onplot(X2,n, 'r' )plot(X2,m, 'g' )现在考虑系数矩阵和右端项的微小变化对方程组解的影响， 即考察方程组2 6 x1 8 ，2 5.99999 x2 8.00002 ，其解变为 (10, 2)T . A=2,6;2,5.99999;B=8;8.00002C=double(AB)cond(A)norm(A)2.0000 6.00002.0000 6.0000B =8.00008.0000C =10 -2 ans =4.0000e+006ans =

8、 8.9443. 同样我们来计算这两条直线产生夹角： theta1 = acosd(dot(2,6,2,5.99999)/(norm(2,6)*norm(2,5.99 999)theta1 =2.8649e-005从上面的结果可以看出系数矩阵的条件数较大，病态较严重，还可以看 出对于系数矩阵微小的变化会对两条直线的夹角产生微小的影响， 由于夹角的微 小变化，会造成解的巨大变化。下面我们来看两条相互垂直的直线系数矩阵微小的变化对解的影响： 3X+4Y=12;4X-3Y=6; 其解为： A=3,4;4,-3;B=12;6C=double(AB)cond(A)norm(A)A = 3 44 -3B=

9、12C = 2.40001.2000ans =1 (其中这是这个方程组的条件数。)ans =5对于系数矩阵产生微小的偏差： A=2.999,4.09;4,-3; B=12;6C=double(AB)cond(A)norm(A)A = 2.9990 4.09004.0000 -3.0000B = 126C = 2.38751.1833ans =1.0180ans =5.0808 画上面的两条直线图像：X1=2.2:0.0001:2.5;y二 3/4*X1+12/4;z=4/3*X1 -6/3;subplot(1,2,1)hold onplot(X1,y, 'b')plot(X1,

10、z, 'g')X2=2.2:0.0001:2.5;n二2.999/4.09*X2+12/4.09;m=4/3*X2-6/3;subplot(1,2,2) hold onplot(X2,n, 'r' ) plot(X2,m, 'g' )可以看出这两条直线是相互垂直的， 当系数矩阵产生微小的 变化时， 会对两直线的夹角产生微小的变化， 这种变化对其解不 会产生太大影响。2.3 从空间几何图形体积的角度来认识矩阵病态态性在线性模型的最小二乘估计中，从空间几何图形的角度来理解病态性就是法矩阵列向量之间的复共线性， 即病 态越严重则复共线的列向量在多维空间的夹角就越小， 由列向 量 形成的超平行多面体的体积就越小A=2,6,5;2,6.0001,5.0002;1,2,3;det(A)B=2,6,5;1.9997,6.0003,4.9991;1,2,