数值分析8-2.pptVIP

§2 幂法及反幂法 * 一、幂法 幂法是一种求实矩阵A的按模最大的特征值λ1及其对应的特征向量x1的方法。特别适合于大型稀疏矩阵。 幂法：由已知非零向量v0及矩阵A的乘幂Ak构造向量序 列{vk}以计算A的主特征值λ1及相应的特征向 量的方法。 的收敛速度由比值 来确定， R越小收敛速度越快，但当r≈1时，收敛可能就很慢。 应用幂法计算A的主特征值λ1及相应的特征向量时，如果|λ1|1(或|λ1|1)，迭代向量vk的各个不等于零的分量将趋向于无穷（或趋于零），这样在计算机实现时就有可能溢出。 解：使用8位浮点数进行运算， 特征值真值及相应特征向量为： 二、加速方法 由幂法原理可知，幂法收敛速度主要由 的大小决定。当r接近1时，敛速就缓慢。 进入矩阵 其中p为选择参数。若 A的特征值为 则，B的相应特征值为 而且A,B的特征向量相同。 如果需要计算A的主特征值λ1就要适当选择p，使得λ1-p仍然是B的主特征值，且使 对B应用幂法，使得在计算B的主特征值λ1-p中得到加速。 这种方法称为原点平移法。 对于A的特征值的某种分布，它是十分有效的。 例3：设A∈R4×4，其特征值为 比值 做变换 则B的特征值为 应用幂法计算B的主特征值的收敛速度的比值为 如果需要计算A的主特征值λ1就要适当选择p，使得λ1-p或者λn-p是B的主特征值，应使p的选择满足 且收敛速度的比值 当 时，w最小，这时收敛速度的比值为 三、反幂法 反幂法可求非奇异实矩阵的按模最小特征值及特征向量。 反幂法计算公式: *