四矩阵的分解.PPTVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * Householder变换 定义 ：设单位列向量u?Rn，称H=I-2uuT为Householder矩阵(初等反射矩阵)，由Householder矩阵确定的线性变换称为Householder变换(初等反射变换)． 对单位向量u，及任意与u垂直的列向量w，有 Hu=(I-2uuT)u=u-2uuTu=-u Hw=(I-2uuT)w=w-2uuTw=w 因此，变换H是关于与u垂直平面的镜像。 基本性质：(1) HT=H(对称矩阵)，(2) HTH=I（正交矩阵），(3) H2=H(对合矩阵)，(4)H-1=H (自逆矩阵)，(5)detH=-1。 u w x Hx Householder变换的性质 定理：任意给定非零列向量x?Rn (n1)及单位列向量z?Rn，则存在Houscholder矩阵H，使得Hx=||x||z。 证明：令 定义矩阵H=I-2uuT，则有 u x |x|z Givens变换与Householder变换的联系 定理：初等旋转矩阵(Givens变换)是两个初等反射矩阵(Householder变换)的乘积． 证：对Tij(?)，取Hu=I-2uuT和Hv=I-2vvT，其中单位向量u=(0,…,0,sin(0.25?),0,…,0, cos(0.25?),0, …,0)和单位向量v=(0,…,0, sin(0.75?),0,…,0, cos(0.75?),0, …,0)，则有Tij(?)= Hv Hu 初等反射矩阵不能由若干个初等旋转矩阵的乘积表示． 基于Givens旋转的QR分解 定理：任何n阶实非奇异矩阵A＝(aij)可通过左连乘初等旋转矩阵化为上三角矩阵。 证 步l.detA?0使A的第1列b(1)=(a11,a21,…,an1)T ?0；存在有限个Givens矩阵的乘积T1，使得T1b(1)=||b(1)||e1和 T1A= 步2.detA(1)?0使A(1)的第1列b(2)=(a22(1),a32(1),…,an2(1))T ?0；存在有限个Givens矩阵的乘积T2，使得T2b(2)=||b(2)||e1和 T2A(1)= 步n-1.detA(n-2)?0使A(n-2)的第1列b(n-1)=(an-1,n-1(n-2), an, n-1(n-2))T?0；存在有限个Givens矩阵的乘积Tn-1使Tn-1b(n-1)=||b(n-1)||e1和 Tn-1A(n-2)= 步n.令 则有A=QR，其中上三角阵R =TA，Q=T-1。因为T是有限个Givens矩阵的乘积，而Givens矩阵都是正交矩阵，所以T是正交矩阵，Q=T-1=TT也是正交矩阵。 基于Householder变换的QR分解 定理：任何n阶实非奇异矩阵A＝(aij) 可通过左连乘Househoulder矩阵化为上三角矩阵。 算法：类同基于Givens旋转的QR分解，仅需用Hi取代相应的Ti(?i=1,2, …,n-2)，并把T换成 Hi =In+1-I -uuT是n+1-i阶Househoulder矩阵，使得 都是n阶Househoulder矩阵。因此，S是有限个Househoulder矩阵的乘积，必是正交矩阵。 定理：设 为任意矩阵，则存在 使得 其中 B 为列满秩矩阵，C为行满秩矩阵(我们称此分解为矩阵的满秩分解)。 证明：假设矩阵 A 的前 r 个列向量是线性无关的，对矩阵 A 只实施初等变行换可以将其化成 矩阵的满秩分解 即存在 使得 于是有 其中 如果 A 的前 r 列线性相关，那么只需对 A 作初等列变换(只需做交换两列的变换)使得前 r 个列是线性无关的。然后重复上面的过程即可。这样存在 A=[Ar, An-r]=[B, BD] 从而 且满足 其中 例 分别求下面三个矩阵的满秩分解 解 （1）对此矩阵只实施初等行变换可以得到 由此可知 rank(A)=2，且该矩阵第一列，第三列是线性无关的。选取 （2）对此矩阵只实施初等行变换可以得到 所以 rank(A)=1，且此矩阵的