[理学]第5章 矩阵特征值和特征向量.pptVIP

§3 Householder法 求一般实矩阵的全部特征值 Def 方阵B若满足：当ij+1时，bij=0,则称B为上Hessenberg阵(或准上三角阵)，即 i=j+1 i j+1 §2 平面旋转矩阵 ? 雅可比法的基本思想: 设法用一系列简单的正角阵Rk , 逐步地将 A 化为近似对角阵(非对角元近似化为0)。即选择Rk , 令 A的全部特征值 问题的关键：如何构造正交阵Rk ？ ? 平面旋转变换 ?雅可比算法: 设Ak-1 (k?1, A0 =A)未对角化，即非对角元中有较大的元素，设非对角元中按模最大的元素是 引入平面旋转矩阵 利用Rk(p,q)对Ak-1作旋转变换，使 中的非对角元 ? 应满足 常将?限制在 对Jacobi算法有几点说明： ? 1. 构造旋转矩阵时只需计算sin?,cos?,为了防止舍入误差扩大，sin?,cos?按下面公式计算： 否则， ? 2. 由于Ak 是对称阵，因此只要计算上三角（或下三角）元素即可，既节省计算量，有能保证Ak 严格对称。 ? 3. 的计算过程如下： ? 4. Ak 中经旋转变换化为零的元素，可能在Ak+1中又成为非零元素，因此不能期望通过有限次旋转变换将原矩阵A对角化，但可证 证明 Jacobi法的收敛性 由前面推论知 ? 5. 实际计算时，当k充分大或者当 时迭代终止， A的全部近似特征值 ? 6. 特征向量的计算：设经过m次旋转变换迭代结束，则 说明Pm的第j列就是?j的标准正交特征向量的近似值。 实际计算时，并不是保留 到最后才形成Pm，而是逐步形成的。 令 每一步的计算公式为 ? 7. 对经典Jacobi法的改进-----避免每次在非对角元中选主元素花费太多时间：循环雅可比法和雅可比过关法。 ?雅可比过关法： 1. 设阈值T0(一般取为 )，在A的非对角 元中按行（或列）扫描（只需扫描上（或下）三角元素），即按如下顺序与阈值T0作比 较： 若 |aij| T0 , 则过关； 否则，作一次旋转变换使aij =0。 2.设阈值 T1=T0 / n,重复上述 过程，直到满足精度要求为止。 按行扫描 按列扫描 ?循环雅可比法：不选主元素，直接将非对角元素按列（或行）的次序扫描并依次化为0。 由于前面已化为0的元素在后面又可能成为非零元素，因此要反复多次扫描，直到达到精度要求为止。 缺点：对一些已经足够小的元素也要作化零处理。 ?雅可比过关法+循环雅可比法：在前几次循环中使用雅可比过关法，经几次循环后，矩阵非对角元素的绝对值的大小已相差不大，这时再使用几次循环雅可比法，效果更好。 例：用雅可比方法计算下面实对称阵的特征值 解： (1) A0=A, 选非对角元中的主元素a12= -1, 因为a11=a22, 取 (1) 选非对角元中的主元素a13= 0.707107 (3) 选非对角元中的主元素a23= -0.627963 (5) 选非对角元中的主元素a13= 0.169525 (4) 选非对角元中的主元素a12= -0.276837 A的近似特征值为 ?1?3.414 209, ?2?0.585 986, ?3?1.999 800 A的准确特征值为 第5章 矩阵的特征值与特征向量 /* Eigen-values and Eigen-vectors of matrix */ ? 待求解的问题：矩阵的特征值?和特征向量x ?0，满足 : Ax=?x or (?I-A)x =0 Eigen-value Eigen-vector 工程技术中的许多问题例如电磁振荡、桥梁振动、机械振动等，都归结为求矩阵的特征值 和特征向量问题----代数计算中的重要课题。 ②特征向量: 已知A的特征值?，求齐次线性方程组 的非零解x, （ ,所以有非零解。）为 A对应于?的特征向量。 如何求解? 特征值：已知A=(aij)n?n，求A的特征多项式的根 有n个零点（实或复，计重数）： 即求解代数方程 A的特征值 从理论上讲，可利用代数方程求根求出特征值，再利用线性方程组的解法，求出特征向量。 ?缺点：工作量大且特征向量对矩阵的依赖很高；当矩阵阶数较高时，高次代数方程求根的计算