矩阵特征值与特征向量计算方法.docVIP

第九章 矩阵特征值与特征向量计算方法 教学目的 1. 掌握求矩阵特征值与特征向量的幂法及反幂法；2. 掌握求矩阵特征值的QR方法。 教学重点及难点 重点是求矩阵特征值与特征向量的幂法及反幂法求矩阵特征值的QR方法；难点是求矩阵特征值的带原点位移的QR方法。 教学时数 12学时 教学过程 §2 幂法及反幂法 2．1幂法 在一些工程、物理力学部标题中，需要我们求矩阵的按模最大的特征值（称为A的主特征值）和对应的特征向量。 幂法是一种计算矩阵的主特征值的一种迭代法，它最大优点是方法简单，适合于计算大型稀疏矩阵的主特征值。 设，其特征值为，对应特征向量为即 且线性无关。设特征值满足：（即为强占优） （2.1） 幂法的基本思想，是任取一个非零初始向量，由矩阵的乘幂构造一向量序列 (2.2) 称为迭代向量。 下面来分折。 由设为中一???基本，于是，有展开式 （且设） 且有 (2.3) 由假设（2.1）式,则 即 且收敛速度由比值确定。且有 （2．4） 这说明，当充分大时，有，或越来越接近特征向量。 下面考虑主特征值的计算。 用表示的第个分量，考虑相邻迭代向量的分量的比值。 从而是 (2.5) 说明相邻迭代向量分量的比值收敛到主特征,且收敛速度由比值来度量,越小收敛越快,但越小收敛越快,但,而接近于1时,收敛可能很慢。 定理7 （1）设n个线性无关的特征向量： （2）设特征值满足 （3）幂法： ） 则 （1）； （2） 如果主特征值为实的重根，即有 又设A有个线性无关的特征向量，其中 对于任意初始向量 则由幂法有 且有 （设不全为零） 由此，当充分大时，接近于与对应的特征向量的某个线性组合。 应用幂法计算的主特征值及对应的特征向量时，如果），迭代向量的各个不等于零的分量将随而趋于无究（或趋于零），这样电算时就可能溢出。为此，就南非要将迭代向量加以规范化。 设有非零向量 其中表示向量绝对值最大的元素，即如果有草药则 其中为所有绝对值最大的分量中最小指标。 显然有下面性性质： 设，则 在定理7条件下幂法可改进为： 任取初始向量。 迭代： 规范化： ， （2．6） 于是，由上式产生迭代向量序列及规范化向量 且改进幂法计算公式为： 设 对于 （2.7） 下面考查与计算的关系。 由 且有 （2.8） 其中 考查规范化向量序列: 由(2.7)及(2.8)式,则有 考查迭代向量序列: 于是， 定理8 (改进幂法) 设有个线性无关特征向量; 设特征值满足 且 (3)由改进幂法得到((2.7)式),则有 (a) (b) 且收敛速度由比值确定。 实现幂法，每迭代一次主要是计算一次矩阵乘向量，可编一个子程序。 计算矩阵主特征值及主特征向量 解 特征值为。由改进幂 法计算公式计算得（见表9-1）。 表9-1  0 1 2 3 4 5 6 7 8 9[0，0，1] [2，4，1] [4.5,9,7,75] [5.7222222,11. ,8.3611111] [5.461165,10.92233,8.2535556] [5.500238,11.0004696,8.2501174] [5.4987088,10.9974176,8.2493544] [5.5002348,11.000469,8.2501174] [5.4999574,10.9999148,8.2499787] [5.5000076,11.0000152,8.2500038]4 9 11.444444 10.92233 11.0142224 10.9974176 11.0004694 10.9999148 11.0000152[0.5,1,0.25] [0.5,1,0.8611111] [0.5,1,0.7305825] [