§3豪斯荷尔德方法.pptVIP

§3 豪斯荷尔德方法 用正交相似变换（即用Householder矩阵）约化一般矩阵为 正交相似变换约化一般矩阵为上Hessenberg阵 过程： 用初等反射阵U作正交相似变换 （1）设 并设 。 上Hessenberg阵。 用正交相似变换（即用Householder矩阵）约化对称矩阵为 对称三对角阵。 选择初等反射阵 令 使 ，于是 其中 。 （2）第k步： 设对A已完成第1步，…，第k-1步正交相似变换, 即 其中 ，为k阶上Hessenberg阵， 设 于是，可选择初等反射阵Rk使 ， 其中 设 于是，可选择初等反射阵Rk使 ， 令 则 其中 为k+1阶上Hessenberg阵（当A为对称阵时，只需计算 （3）继续上述过程，最后有 例5 设 解 （1）选取 ，使 其中 , 试用初等反射阵正交相似约 化A为上Hessnberg阵。 约化计算，计算 。 得 计算 从而 为上Hessenberg阵。 酉相似变换约化一般矩阵的结论。 定理4（Schur定理） （上三角阵） 则存在酉阵U使 其中 （i =1，…，n ）为A的特征值。 说明:（1）A是酉变换：若A满足 标准正交基下的矩阵叫酉阵。 （2）UH是U的共轭转置阵。 酉变换在 证明：对矩阵的阶数n用数学归纳法证明。 当n=1时，A=（a），取U=（1）即可。 证明：对矩阵的阶数n用数学归纳法证明。 当n =1时，A =（a），取U =（1）即可。 假设对阶数小于或等于n-1的矩阵定理成立，下证对n阶矩阵 为A的特征值及对应的特征向量，即 第一个分量 其中 定理成立。 于是，存在初等反射阵H使 考虑HAH H 第一列元素 即 其中A1为n-1阶矩阵。 设对矩阵A1存在酉阵P1使P1HA1P1=R1（上三角阵） 令 P 于是 P H(HAHH)P = 令 U=HHP，则 UHAU=R（上三角阵）。 # 即 正交相似变换约化一般实矩阵的结论。 定理5（实Schur定理）设A=（aij）n×n，则存在正交矩阵Q使 （准上三角阵） 其中Rii(i =1,…,m)为一阶或二阶方阵，且每个一阶Rii是A的实特 征值，每个两阶对角块Rii的两个特征值是A的一对共轭复特征值。 (3) 两种方法都是数值稳定的。 (1) 正交相似变换约化对称阵为对称三对角阵的思路与正交相 (2) 计算量 正交相似变换约化对称阵为对称三对角阵大约需要 次乘法 似变换约化一般矩阵为上Hessenberg阵类似。 说明 ： 乘法运算 。 运算 ； 正交相 似 变换约化一般矩阵为上Hessenberg阵约需要 次 理解用正交相似变换（即用Householder矩阵）约化一般矩阵 (或对称阵)为上Hessenberg阵(或对称三对角阵)的思想方法。