[特征值和特征向量的数值算法.pptVIP

第七章 特征值与特征向量的数值求法 7.4 QR算法 7.4.3 带原点位移的QR算法 7.4.2 QR算法及其收敛性 7.4.1 化矩阵为Hessenberg形 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 我们称这种分块上三交阵为矩阵A的Schur分块上三角阵，上三角阵 和对角阵是它的特殊情形。定理7.9并没有解决如何计算全部特征的问题。 7.4.1 化矩阵为Hessenberg形 对于实对称矩阵，可通过正交相似变换约化为对角矩阵。那么，对于一般的 实矩阵，通过正交相似变换可约化到什么程度呢？线性代数中有如下结果。 定理7.9 （实Schur分解定理） Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 定义7.2 (ij+1), 则称B为上Hessenberg矩阵，简称Hessenberg形，即B 的形状为 可以用平面旋转变换化矩阵为Hessenberg形,下面介绍另一种正交 变换。 为了节省运算工作量，实用的方法是先将矩阵约化为与Schur分块上三角 阵很近似的Hessenberg形。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 定义7.3 （7.4.1） 为（初等）镜面反射矩阵，或Householder变换矩阵。 Houholder矩阵H=H（w）有如下性质： (1) (2) (3) 记S为与w垂直的平面,则几何上x与y=Hx关于平面S对称。事实上,由 上式表明向量x-y与w平行，注意到y与x的长度相等，于是x经变换后的象y=Hx 是x关于s对称的向量，如图7-1所示。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. x w y x-y 图7-1 对应于性质（2），有下面的定理。 定理7.10 得Hx=y。 证 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 由此可得 定理得证。 （7.4.2） （7.4.3） Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 稳定性。（7.4.2）的意义是对向量作消元运算。与平面旋转不同的是，镜面反 射变换可成批的消去向量的非零元。 例7.4 解 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. （7.4.4） 定理7.11 为Hessenberg形。 证 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 变换，有 如此类推，经n-2步对称正交相似变换，得到Hessenberg形矩阵。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 推论7.1 对称三对角阵。 上述定理7.11的证明是构造的，即可以用镜面反射化矩阵Hessenberg形。 此定理可用平面旋转变换来证明，即也可用平面旋转变换化矩阵为Hessenberg 如此类推，最后得到的正交矩阵Q，是平面旋转矩阵的乘积。