[理学]矩阵特征值与特征向量的数值计算.ppt

数值计算方法 第三章 矩阵特征值与特征向量的计算 一些工程技术问题需要用数值方法求得矩阵的全部或部分特征值及相关的特征向量。 3.1 特征值的估计 较粗估计?(A) ? ||A|| 欲将复平面上的特征值一个个用圆盘围起来。 3.1.1 盖氏圆 定义3.1-1 设A = [aij]n?n，称由不等式 所确定的复区域为A的第i个盖氏圆，记为Gi： i = 1，2，…，n。 定理3.1-1 若?为A的特征值，则 定理3.1-1 若?为A的特征值，则 证明：设Ax = ?x (x ? 0)，若k使得 因为 ? ? ? 例1 估计方阵特征值的范围 解： G1 = {z：|z – 1|? 0.6}；G2 = {z：|z – 3|? 0.8}； G3 = {z：|z + 1|? 1.8}；G4 = {z：|z + 4|? 0.6}。 注：定理称A的n个特征值全落在n个盖氏圆上，但未说明每个圆盘内都有一个特征值。 3.1.2 盖氏圆的连通部分 称相交盖氏圆之并构成的连通部分为连通部分。 孤立的盖氏圆本身也为一个连通部分。 定理3.1-2 若由A的k个盖氏圆组成的连通部分，含且仅含A的k个特征值。 定理3.1-2 若由A的k个盖氏圆组成的连通部分，含且仅含A的k个特征值。 证明: 令D = diag(a11，a22，…，ann)，M = A – D，记 则显然有A(1) = A，A(0) = D， 易知A(?)的特征多项式的系数是?的多项式， 从而A(?)的特征值?1(?)，?2(?)，…，?n(?)为?的连续函数。 A(?)的盖氏圆为： 因为A(0) = D的n个特征值a11，a22，…，ann，恰为A的盖氏圆圆心，当?由0增大到1时，?i(?)画出一条以?i(0) = aii为始点，?i(1) = ?i为终点的连续曲线，且始终不会越过Gi； 不失一般性，设A开头的k个圆盘是连通的，其并集为S，它与后n – k个圆盘严格分离，显然，A(?)的前k个盖氏圆盘与后n – k个圆盘严格分离。 当? = 0时，A(0) = D的前k个特征值刚好落在前k个圆盘G1，…，Gk中，而另n – k个特征值则在区域S之外，?从0变到1时， 与 始终分离（严格）。连续曲线始终在S中，所以S中有且仅有A的k个特征值。 注：1) 每个孤立圆中恰有一个特征值。 2) 例1中G2，G4为仅由一个盖氏圆构成的连通部分，故它们各有一个特征值，而G1，G3构成的连通部分应含有两个特征值。 3) 因为例1中A为实方阵，所以若?为A的特征值，则 也是A的特征值，所以G2，G4中各有一个实特征值。 3.1.3 盖氏圆与相似变换 由于特征值是相似不变量，所以代数上常用相似变换将矩阵化简以得到特征向量，这里也可用相似变换将盖氏圆的半径变小，以得到更好的估计。 原理：取对角阵作相似变换阵：P = diag(b1，b2，…，bn)其中bi 0，i = 1，2，…，n 则 与A有相同特征值. 而B的第i个盖氏圆为： ， 而B的第i个盖氏圆为： ， 适当选取b1，b2，…，bn就有可能使B的某些盖氏圆的半径比A的相应盖氏圆的半径小。 1) 欲缩小Gi，可取bi最大。 2) 欲缩小除Gi外的圆，可取bi最小。 例2，估计 的特征值范围。 解：A的三个盖氏圆分别为： G1 = {z：| z – 0.9|? 0.13}； G2 = {z：| z – 0.8|? 0.14}； G3 = {z：| z – 0.4|? 0.03} ?3 ? G3，较好。 为了更好地估计另外两个特征值可取b3最小： 取b1 = b2 = 1，b3 = 0.1即 则 所以G1 = {z：| z – 0.9|? 0.022}； G2 = {z：| z – 0.8|? 0.023}； G3 = {z：| z – 0.4|? 0.3} 三个盖氏圆分离，故有?1 ? G1，?2 ? G2，?3 ? G3。 3.2 幂法与反幂法 幂法是求方阵的最大特征值及对应特征向量的一种迭代法。 3.2.1 幂法 设An有n个线性无关的特征向量v1，v2，…，vn，对应的特征值?1，?2，…，?n，满足 |?1| |?2| ? …? |?n|