正交变换与QR迭代矩阵特征值计算

《正交变换与QR迭代矩阵特征值计算》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第八章矩阵特征值计算计算方法——正交变换与QR迭代1本讲内容正交变换QR迭代Householder变换Givens变换QR分解Schur分解Hessenberg矩阵2Householder变换性质(1)对称：(2)正交：(3)对合：(4)保模：(5)定义：设且，称矩阵为Householder变换，或初等反射矩阵。3Householder变换定理：设x,yRn,xy且

2、

3、x

4、

5、2=

6、

7、y

8、

9、2，则存在n阶Householder变换H，使得y=Hx证：取4Householder变换定理：对任意的非零向量xRn，存在Householder变换H，使得Hx=e1其中=sgn

10、(x1)

11、

12、x

13、

14、2,e1=(1,0,...,0)T，的选取是为了防止在实际计算中与x1互相抵消若x1=0,则取=

15、

16、x

17、

18、25Givens变换定义：称矩阵为Givens变换，或旋转变换。ij6Givens变换性质(1)只有四个元素与单位矩阵不同(2)正交：(3)用G左乘一个矩阵时，只改变该矩阵中两行的值(4)用G右乘一个矩阵时，只改变该矩阵中两列的值7Givens变换定理：设x=(x1,...,xi,...,xj,...,xn)T，且xi,xj不全为零，则存在Givens变换G=G(i,j,)，使得8QR分解定理：（QR分解）设n阶实矩阵A非奇异，则存在正交分解A=

19、QR其中Q是正交矩阵，R是非奇异上三角矩阵。若限定R的对角线元素为正数，则此分解唯一。9QR分解算法设(j=1,...,n)(1)构造H1使得H1a1=1e1，令(2)构造使得，令算法（QR分解）10QR分解算法以此类推，经过n-1步，可得Householder矩阵H1,H2,...,Hn-1，使得令，即得11QR分解举例例：用Householder变换计算的QR分解解：(板书)12Schur分解定理：（Schur分解）设A为n阶实矩阵，则存在正交矩阵Q，使得其中Rii是一阶或二阶方阵。若Rii是一阶方阵，则它就是A的特征值；若Rii是二阶方阵，则其特征值为A的两个共轭复特

20、征值。拟上三角矩阵13QR迭代QR迭代算法计算矩阵的所有特征值和特征向量计算过程(1)令A1＝A(2)对k=1,2,...，计算Ak的QR分解计算直到Ak+1收敛到一个拟上三角阵14作业教材277页，习题1015