25讲-9—2Jacobi法和QR方法.pptVIP

矩阵的旋转变换 基本QR方法： 内容小结： 第九章 矩阵特征值与特征向量的计算 矩阵特征值与特征向量的计算 幂方法及程序 幂法的加速—反幂法及程序 Jacobi方法及程序 QR方法及程序 雅克比方法 第二节 雅克比方法 雅克比法的基本思想 雅克比方法理论与程序 雅克比方法概述 小结 雅克比方法概述 雅克比方法是求实对称矩阵的全部特征值及相应 特征向量的一种方法。 预备知识 （1）任意实对称矩阵A 可通过正交相似变化成 雅克比方法概述 （2）在正交相似变换下，矩阵元素的平方和不变。 基本思想 雅克比法的基本思想 下面以二阶实对称矩阵为例分析 矩阵形式为： 雅克比方法的实例分析 过程如下： 雅克比方法的实例分析 则上式简写为： 容易验证R是一个正交矩阵。正交变换R把对称 矩阵A变成为对角矩阵， 正交矩阵R的两个列向量分别为对应于两个特征 值的单位特征向量。 上述结果可推广到一般情况。 雅克比方法的实例分析 雅克比方法的一般推广 雅克比方法的一般推广 雅克比方法的一般推广 可证，虽然这种变换不一定能使矩阵中非对角元中零 元素的个数增加，但是可以保证非对角元的平方和递减。 雅克比方法的一般推广 上式表面，在旋转变换下，非对角元的平方和严格单调减。 因而对角元的平方和电调增，利用此点，则导出了Jocobi方法。 Jacobi方法的理论基础 雅克比法的理论基础 相似变换的过程中，每一步都选绝对值最大的非对角元素 如果在对A做 定理 例3 雅克比法 例题 雅克比法 例题 雅克比法 例题 雅克比法 例题 第三节 QR方法 QR分解法及程序 QR方法分析 基本QR分解法 小结 QR方法 QR算法是计算中小型矩阵的全部特征值与特征向量的最有效方法。适用于求实、复非奇异矩阵的特征值，是一种变换迭代法。 理论依据：任一非奇异实矩阵都可分解成一个正交矩阵Q和一个上三角矩阵R的乘积，而且当R的对角元符号取定时，分解是唯一的。 基本QR方法 可证，在一定条件下，基本QR方法产生的矩阵序列｛A(k)｝ “基本”收敛于一个上三角阵（或分块上三角阵）。即主对角线（或主对角线子块）及其以下元素均收敛，主对角线（或主对角线子块）以上元素可以不收敛。特别的，如果A是实对称阵，则｛A(k)｝ “基本”收敛于对角矩阵。 基本QR方法 因为上三角阵的主对角元（或分块上三角阵中，主 对角线子块的特征值）即为该矩阵的特征值，故当k充 分大时， A(k)的主对角元（或主对角线子块的特征值） 就可以作为A的特征值的近似。 基本的QR方法的主要运算是对矩阵QR分解，分解的 方法有多种。介绍一种Schmit正交化方法为例。 基本QR方法 基本QR方法 基本QR方法 基本QR方法 例4 QR方法例题 QR方法例题