几何的角度解释矩阵出现病态2.docx
《几何的角度解释矩阵出现病态2.docx》由会员分享,可在线阅读,更多相关《几何的角度解释矩阵出现病态2.docx(9页珍藏版)》请在冰豆网上搜索。
几何的角度解释矩阵出现病态2
矩阵的范数
范数的定义
若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有微小扰动时,参数估值
X也有扰动:
(A+ΔA)(X+ΔX)=B+ΔB
(1)
AΔX=ΔA(X+ΔX)+ΔB
ΔX=-A-1ΔA(X+ΔX)+A-1ΔB将上式两端取范数,并应用向量范数的三角不等式║AB║│≤│A║B║;║A+B║≤║A║+║B║。
||ΔX||≤||A-1||||ΔA||(||X||+||ΔX||)+||A-1|||ΔB||
把含有ΔX的项移到式子的左边有:
(1-||A-1||||ΔA||)||ΔX||≤||A-1||||ΔA||||X||+||A-1|||ΔB||
由于有1/||X||≤||A||/||B||
将上式两端同时乘以1/||X||得:
(1-||A-1||||ΔA||)||ΔX||/||X||≤||A-1||||ΔA||+||A-1||||ΔB||A||/||B||;
设K=||A||||A-||将上式整理的:
(1-K||ΔA||/||A||)||ΔX||/||X||≤K(||ΔA||/||A||+||ΔB||/||B||);
即有:
||ΔX||/||X||≤k/(1-K||ΔA||/||A||)(||ΔA||/||A||+||ΔB||/||B||);
问题与实验1:
试从几何的角度解释矩阵出现病态的原因,并用‘有说服力’的例子来支持你的观点;
线性方程组解的敏感性的几何解释(2x2矩阵)线性方程组求解:
两直线求交点下面两图分别反映了良态问题和病态问题两种情况。
良态情况病态情况
由图像可以看出当两条直线的夹角较大时其精确解与近似解距离较小(其中两条直线的交点是其精确解)当两条直线的夹角较小时其精确解与近似解距离较大。
虚线为直线的系数产生微小的扰动时的直线。
下面我们用两个例子说明:
例设方程组
,即
,它的精确解为
.A=[2,6;2,6.00001;]
B=[8;8.00001]
C=double(A\B)
cond(A)
norm(A)
A=
2.00006.0000
2.00006.0000
B=
8.0000
8.0000
C=11ans=4.0000e+006
ans=8.9443
下面我们来求上面两条直线的夹角:
计算直线夹角
根据给定的节点,计算夹角信息,示意图如下。
已知点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为计算以度为单位的夹角信息。
由这个公式计算这两条直线的夹角:
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)
holdon
plot(X1,y,'b')
plot(X1,z,'g')
X2=[9.9:
0.0001:
10.1];
n=-1/3*X2+4/3;
m=-2/5.99999*X2+8.00002/5.99999;
subplot(1,2,2)
holdon
plot(X2,n,'r')
plot(X2,m,'g')
现在考虑系数矩阵和右端项的微小变化对方程组解的影响,即考察方程组
,
其解变为
.A=[2,6;2,5.99999;]
B=[8;8.00002]
C=double(A\B)
cond(A)
norm(A)
A=
2.00006.0000
2.00006.0000
B=
8.0000
8.0000
C=10-2ans=4.0000e+006
ans=8.9443.同样我们来计算这两条直线产生夹角:
theta1=acosd(dot([2,6],[2,5.99999])/(norm([2,6])*norm([2,5.99999])))
theta1=2.8649e-005
从上面的结果可以看出系数矩阵的条件数较大,病态较严重,还可以看出对于系数矩阵微小的变化会对两条直线的夹角产生微小的影响,由于夹角的微小变化,会造成解的巨大变化。
下面我们来看两条相互垂直的直线系数矩阵微小的变化对解的影响:
3X+4Y=12;4X-3Y=6;其解为:
A=[3,4;4,-3;]
B=[12;6]
C=double(A\B)
cond(A)
norm(A)
A=34
4-3
B=12
6
C=2.4000
1.2000
ans=1(其中这是这个方程组的条件数。
)
ans=5
对于系数矩阵产生微小的偏差:
A=[2.999,4.09;4,-3;]
B=[12;6]
C=double(A\B)
cond(A)
norm(A)
A=2.99904.0900
4.0000-3.0000
B=12
6
C=2.3875
1.1833
ans=1.0180
ans=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)
holdon
plot(X1,y,'b')
plot(X1,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)
holdon
plot(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,3;]
det(B)
。
A=
2.00006.00005.0000
2.00006.00015.0002
1.00002.00003.0000
ans=5.0000e-004
B=
2.00006.00005.0000
1.99976.00034.9991
1.00002.00003.0000
ans=9.0000e-004