第四章 矩阵的分解.ppt

定义:设实数c与s满足c2+s2=1，称为Givens矩阵(初等旋转矩阵)，Givens矩阵确定的线性变换称为Givens变换．当c2+s2=1时，存在角度?，使得c＝cos?，s＝sin?。Givens变换ij第28页，共61页，星期日，2025年，2月5日Givens矩阵的性质性质1Givens矩阵是正交矩阵，且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．第29页，共61页，星期日，2025年，2月5日推论：设非零列向量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))TT1n(1)…T12(1)是有限个Givens矩阵的乘积．第30页，共61页，星期日，2025年，2月5日Householder变换定义：设单位列向量u?Rn，称H=I-2uuT为Householder矩阵(初等反射矩阵)，由Householder矩阵确定的线性变换称为Householder变换(初等反射变换)．对单位向量u，及任意与u垂直的列向量w，有Hu=(I-2uuT)u=u-2uuTu=-uHw=(I-2uuT)w=w-2uuTw=w因此，变换H是关于与u垂直平面的镜像。基本性质：(1)HT=H(对称矩阵)，(2)HTH=I（正交矩阵），(3)H2=H(对合矩阵)，(4)H-1=H(自逆矩阵)，(5)detH=-1。uwxHx第31页，共61页，星期日，2025年，2月5日Householder变换的性质定理：任意给定非零列向量x?Rn(n1)及单位列向量z?Rn，则存在Houscholder矩阵H，使得Hx=||x||z。证明：令定义矩阵H=I-2uuT，则有ux|x|z第32页，共61页，星期日，2025年，2月5日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(?)=HvHu初等反射矩阵不能由若干个初等旋转矩阵的乘积表示．第33页，共61页，星期日，2025年，2月5日基于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)=第34页，共61页，星期日，2025年，2月5日步n-1.detA(n-2)?0使A(n-2)的第1列b(n