chapter041.pptVIP

* * 西安邮电学院应用物理系 第3章 矩阵特征值运算 许多实际问题的求解涉及到矩阵运算, 常见的矩阵运算有矩阵的基本运算(如四则运算, 幂运算, 指数、对数、开方运算, 转置、求逆运算), 矩阵的特殊运算(如特征值运算, 范数运算, 条件数运算, 奇异值运算, 秩运算, 迹运算, 正交运算, 伪逆运算等)和矩阵的分解运算(如特征值分解, Choleskey分解, QR分解, QZ分解, LU分解, Schur分解, 奇异值分解等)等。 本章学习目标 掌握矩阵基本知识； 掌握矩阵特征值的计算方法。 西安邮电学院应用物理系 3.1 矩阵有关概念 (1) 对角阵、上三角阵和下三角阵 D 称为对角阵; L / U称为下 / 上三角阵, 当其对角元全为1时, 称为单位下/上三角阵, 当其对角元全为0时, 称为严格下/上三角阵. det D= Πi di det L= Πi lii det U= Πi uii 西安邮电学院应用物理系 (2) 顺序主子阵与顺序主子式 3 阶 顺 序 主 子 阵 D3=det A3 3阶顺序主子式 西安邮电学院应用物理系 (3) 排列阵 每行每列只有一个元素等于1而其余元素为零的方阵称为排列阵 排列阵有以下性质 PT= P -1, det P=±1; (4) 相似阵 设A和B均为n阶矩阵, 若存在可逆矩阵S使得 B=S-1AS 或 A=SBS-1 则说A相似于B, 记作A~B 西安邮电学院应用物理系 (5) 正交阵 若 A?Rn?n 为 n 阶实矩阵, 且 AAT=I, 则 A 为正交阵。 特点: det A=±1 AT=A-1 (AT=A-1也是正交阵) 两个正交阵的乘积仍然是正交阵 在物理上 A 称为幺正矩阵。 西安邮电学院应用物理系 (6) 矩阵的特征值与谱半径 设A为n阶矩阵, 若存在某一个非零数?和非零列向量x 使得Ax=?x, 则称? 为矩阵A的特征值, x为A相应于特征值?的特征向量; ? 是A的特征根的充要条件是行列式 det(?I -A)=0 而A对应于特征根?的特征向量是齐次线性方程组 (?I -A)x=0 的非零解。 设A?Rn?n,?i 为其特征值, 则称? (A)=max|?i|为矩阵A的谱半径。 西安邮电学院应用物理系 (7) 对角占优阵 设A?Rn?n, 若 则称A为严格对角占优阵。 设A?Rn?n是不可约的, 若 则称A为不可约弱对角占优阵。 特点: aii ? 0 (i=1,2,…,n) A为非奇异矩阵( 即det A ?0 ) 西安邮电学院应用物理系 设A?Rn?n, 若存在一个排列阵P, 使得 则称A为可约的，若不存在这样的排列阵，则称A是不可约的。 西安邮电学院应用物理系 (8) 对称正定阵 设A?Rn?n, AT=A则称A为实对称阵; 设A?Rn?n, AT=A且满足(Ax, x) 0, ? x ?0, 则称A正定阵; 若满足(Ax, x)? 0, ? x ?0, 则称A为半正定阵。 特点: (1)实对称矩阵A正定的充要条件是A的所有特征值 ?i 0 (i=1,2,…,n); (2) 实对称矩阵A正定的充要条件是A的所有顺序主 子式都是正的; (3) 对称正定矩阵A的对角元素是正的, 即aii0; (4) 半正定矩阵A的特征值?i ?0 (i=1,2,…,n)。 西安邮电学院应用物理系 (9) 范数与条件数 向量范数 如果向量x?Rn的某个实值函数 f (x)= || x || 满足 (1)正定性: || x ||≥0, 当且仅当x=0时取等号 (2)齐次性:|| ?x ||=|? |.||x||, ? ? ?R (3)三角不等式:||x+y||?||x||+||y||, ? x, y?Rn 则称||x||为Rn上的一个向量范数。 常用的向量范数有: 西安邮电学院应用物理系 矩阵范数 设A?Rn?n的某个实值函数 f (A)= || A || 满足 (1)正定性: || A ||≥0, 当且仅当A=0时取等号 (2)齐次性:|| ?A ||=|? |?||A||, ? ? ?R (3)三角不等式:||A+B||?||A||+||B||, ? A, B?Rn (4) 相容性: ||AB||?||A||?||B||, ? A, B?Rn 则称||A||为Rn?n上的一个矩阵范数。 常用的矩阵范数有: 西安邮电学院应用物理系 矩阵条件数 若A-1存在,