第四章 矩阵的分解.pptVIP

定义: 设实数c与s满足c2+s2=1，称 为Givens矩阵(初等旋转矩阵)，Givens矩阵确定的线性变换称为Givens变换．当c2+s2=1时，存在角度?，使得c＝cos?，s＝sin? 。 Givens变换 i j 第二十八页，共六十一页，2022年，8月28日 Givens矩阵的性质 性质1 Givens矩阵是正交矩阵，且detTij=1和[Tij(c, s)]-1= [Tij(c, s)]T= Tij(c, -s)。 性质2 设x=[?1,?2,…,?n]T，y=Tijx=[?1, ?2,…,?n]T，则有?i=c?i+s?j，?j=-s?i+c?j，?k=?k (k?i, j).当?i2+?j2?0时，选取c= ?i /[?i2+?j2]0.5，s= ?j /[?i2+?j2]0.5，就可使?i=[?i2+?j2]0.50， ?j=0。 定理：设x=[?1,?2,…,?n]T，则存在有限个Givens矩阵的乘积，记作T，使得Tx=||x||e1． 第二十九页，共六十一页，2022年，8月28日 推论：设非零列向量x?Rn及单位列向量z?Rn，则存在有限个Givens矩阵的乘积，记作T，使得Tx=||x||z。 证：根据定理，存在T(1)=T1n(1)…T12(1)，使得T(1)x=||x||e1；存在T(2)=T1n(2)…T12(2)，使得T(2)z=||z||e1=e1；于是有 T(1)x=||x||e1 =||x||T(2)z 或者 [T(2)]-1T(1)x=||x||z 于是 T=[T(2)]-1T(1) =[T1n(2)…T12(2)]-1T1n(1)…T12(1) =(T12(2))T…(T1n(2))T T1n(1)…T12(1) 是有限个Givens矩阵的乘积． 第三十页，共六十一页，2022年，8月28日 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 第三十一页，共六十一页，2022年，8月28日 Householder变换的性质 定理：任意给定非零列向量x?Rn (n1)及单位列向量z?Rn，则存在Houscholder矩阵H，使得Hx=||x||z。 证明：令 定义矩阵H=I-2uuT，则有 u x |x|z 第三十二页，共六十一页，2022年，8月28日 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 初等反射矩阵不能由若干个初等旋转矩阵的乘积表示． 第三十三页，共六十一页，2022年，8月28日 基于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)= 第三十四页，共六十一页，2022年，8月28日 步n-1.detA(n-2)?0使A(n-2)的第1列b(n-1)=(an-