第八章 特征值问题的计算方法1.pptVIP

幂法的加速 * 第八章 特征值问题计算 一、特征值和特征向量的基本概念与性质 §1 基本概念与性质 设 ，若存在向量 和复数 满足 ，则称 是矩阵 的特征值， 是特征值 相应的特征向量。 特征多项式 的根的集合：谱集 矩阵A的特征根模的最大值称为矩阵A的谱半径： 于是得到两个特征根分别为： 例8.1 求矩阵A的特征根与特征向量 特征多项式为： 对应的特征向量分别为： 于是得到两个特征根分别为： 例8.2 求矩阵A的特征根与特征向量 特征多项式为： 对应的特征向量分别为： 其中 称 为 的代数重数（简称重数）； 为 的几何重数。 设 ， 对于矩阵 的特征值 ，如果 ，则称该特征值 为 的一个半单特征值。 若 的所有特征值都是半单的，则称 是非亏损的。 是非亏损的等价条件是 有n个线性无关的特征向量 一般的，对矩阵A， 其特征多项式可表示为 设 ， 若存在矩阵 ，使得 则称 和 是相似的。 相似矩阵有相同的特征值 设 寻求已知矩阵 的相似矩阵 ，要求： 矩阵 的特征值和特征向量容易计算 本章QR算法的计算思想： 关于矩阵相似，有后面的结论 设 ，有 r个互不相同的特征值 ， 其重数分别为 ，则一定存在非奇异矩阵 使得 （Jordan分解） 其中 且除了 的排列 次序外， 是唯一的。 称作 的Jordan标准型 设 ，则存在酉矩阵 ，使得： （Schur分解） 其中 是上三角矩阵，且适当选择 ，可使 的元素 按任意指定的顺序排列。 设 ，令 （圆盘定理）/*Disc Theorem*/ 则 设 为对称矩阵，则存在正交矩阵 （谱分解定理）/*Spectral Decomposition*/ 其中 是 的n个特征值。 使得 设 为对称矩阵，且 的特征值为 （极大极小定理） 其中 表示 中所有k维子空间的全体。 则有 设 为对称矩阵，其特征值分别为 （Weyl定理） 则有 说明：对称矩阵的特征值总是良态的。 注意：实际问题中矩阵一般都是由计算或实验得到， 本身必然存在误差，不妨假设 例8.3求矩阵A的特征根与特征向量 其特征多项式为： 于是得到两个特征根分别为： 若取初始向量为： 先做xk+1 A*xk迭代，并计算|| xk+1 || / || xk || 可发现 对应的特征向量分别为： || xk ||表示的分量模长的最大值，即取无穷范数 1.座机电话号码000000 1.座机电话号码666667 1.座机电话号码000000 1.座机电话号码888889 1.座机电话号码058824 1.座机电话号码969697 1.座机电话号码461538 1.座机电话号码201550 1.座机电话号码941634 1.座机电话号码226121 1.9990 1.9995 1.9998 1.9999 1.9999 || xk+1 || / || xk || xk 0.5, 1.5 1.5, 2.5 3.5, 4.5 7.5, 8.5 15.5, 16.5 31.5,32.5 63.5,64.5 127.5, 128.5 255.5, 256.5 511.5, 512.5 1023.5, 1024.5 2047.5, 2048.5 4095.5, 4096.5 8191.5, 8192.5 1638.4, 1638.5 结果表明|| xk+1 || / || xk ||收敛到最大特征根， xk 收敛到对应的特征向量。 k 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 但xk 的绝对值越来越大， xk / || xk ||即为对应特征 向量[1;1] 考虑对每次迭代结果归一化， 若做如下迭代：xk+1 A*xk / || A*xk || 则有xk 收敛到对应的特征向量， || A*xk ||收敛到最大特征根。 1.座机电话号码000000 1.座机电话号码666667 1.座机电话号码000000 1.座机电话号码888889 1.座机电话号码058824 1.座机电话号码969697 1.座机电话号码461538 1.座机电话号码201550 1.座机电话号码941634 1.座机电话号码226121 1.9990 1.9995 1.9998 1.9999 1.9999 || A*xk || xk 0.3333,1.0000 0.6000,1.0000 0.7778,1.0000 0.8824,1.0000 0.9394,1.0000 0.9692,1.0000 0.9845,1.0000 0.9922,1.0000 0.9961,1.0000 0.998