第五讲-矩阵特征值问题-自学纪注参考课件.pptxVIP

数 值 分 析 (5)Numerical Analysis计算机系 软件所 喻文健第五章 矩阵特征值问题矩阵特征值与特征向量的计算: ?数值线性代数, 矩阵计算(matrix computation)的重要内容矩阵特征值=特征多项式方程的解本章内容基本概念与特征值的分布幂法与反幂法矩阵的正交三角化, QR分解计算所有特征值的QR迭代算法一定是迭代解法Wenjian Yu矩阵特征值的基本概念Wenjian Yu基本概念矩阵的特征值与特征向量, (为方阵)特征值是特征方程的根, 复数域内有n个(含重根)特征值谱: 给定一特征值, 特征向量是方程的非零解对任一特征值, 特征向量都不唯一, 构成特征子空间假设, (1)特征值不一定是实数; (2)实特征值一定对应实特征向量; (3)非实特征值的共轭也是特征值, 其对应的特征向量一定不是实向量???Wenjian Yu特征值的有关性质?非奇异矩阵特征值均不为0; 0一定是奇异矩阵的特征值若为对角阵或上(下)三角阵, 则其特征值为其对角元若为分块对角阵或分块上(下)三角阵(对角块为方阵), 则 相似矩阵的特征值相等矩阵运算结果的特征值: 设为的特征值，则定义5.2 设有m个(m?n)不同的特征值, 若是特征方程的重根, 则称为的代数重数, 而的特征子空间的维数其几何重数. Th5.2?Th5.3Th5.4?Th5.5?Th5.9?;?;?;?;Wenjian Yu特征值的有关性质?设n阶方阵的m个不同的特征值为, 的代数重数为, 几何重数为, 则; ,不同特征值的特征向量线性无关, 所有特征子空间的个基形成一组线性无关向量若, 这种矩阵为非亏损阵, 否则为亏损阵非亏损阵有n个特征向量构成全空间的基 ()Th5.7: , 为对角阵 ? 为非亏损阵实对称阵非亏损, 且特征值为实数 (可正交对角化)Th5.8(Jordan分解): , Th5.6定义5.3(特征值分解)? ?的几何重数为, , 对应个约当块, 其阶数之和Wenjian Yu特征值的分布Motivation迭代法的收敛性, 看= 矩阵的2-范数、2-条件数: 是关于特征值上界的重要结论定义5.4 对于复矩阵, 设, 在复平面上以为圆心、为半径的圆称为的Gerschgorin圆盘Th5.10 (圆盘定理)的特征值必在某个圆盘上若n个圆盘中有m个连通, 且与其他分离, 则这m个圆盘恰包含m个特征值?思考: 如何证明????交互演示网站Wenjian Yu圆盘定理的应用Th5.11 严格对角占优阵, 正对角元, Re(?)0若A为对角占优的对称阵，则…例: 估计矩阵的特征值分布直接应用圆盘定理, D1与其他两个分离, 必含一个实特征值再对应用圆盘定理, 与其他圆盘分离, 必含一实特征值准确特征值为: 4.2030, -0.4429, -3.7601还可先做简单的相似变换, 再估计分布, 见pp.153, 例5.4??0D3D2D1-4 0 4???-4 0 4Wenjian Yu幂法与反幂法Wenjian Yu计算最大的特征值、特征向量定义5.5 模最大的特征值称为主特征值, 也叫”第一特征值”, 它对应的特征向量称为主特征向量 主特征值可能不唯一, 例如的模都是5幂法(power iteration): 取任意非零向量, 计算, 得到向量序列看一个例子: 取, 利用Matlab算向量序列通过实验看出: 相邻两个逐渐呈倍数关系, 倍数为3 趋近于3对应的特征值向量??(pp. 148, 例5.2, 其主特征值为3)Matlab demo更多例子：交互演示网站Wenjian Yu幂法幂法是否总能计算矩阵的主特征值?定理5.12 若矩阵有唯一的主特征值, 向量非零, 计算, 则一般有证明: 设非亏损阵, 且不是重特征值, ???条件: 的几何重数等于代数重数??的取值须保证有使为某个主特征向量, ???设的线性无关的单位特征向量为,???思考: 若=0呢?????对的取值要求Wenjian Yu幂法?特殊情况:不能为某特征向量!?关于幂法的说明为亏损矩阵的情况, 用矩阵的Jordan标准型进行证明?相邻项的第分量的比值?主特征值, 选为最大分量证明中假设了, 在实际应用时可随机选取关于幂法的问题实用的幂法: 用规格化向量的技术防止溢出定义5.6 记为向量的绝对值最大分量, 若不唯一则取最小编号的那个. 称为向量的规格化向量?? ,很大时, 可能出现上溢或下溢??? 收敛速度主要取决于Wenjian Yu实用的幂法规格化向量例: ,, 规格化向量为若为规格化向量, 则, 向量,的规格化向量分别为,, 若=则=在幂法的每步增加向量规格化操作?