[计算机软件及应用]并行计算-矩阵特征值计算 -.docVIP

9 矩阵特征值计算 在实际的工程计算中，经常会遇到求 n 阶方阵 A 的特征值(Eigenvalue)与特征向量(Eigenvector)的问题。对于一个方阵 A，如果数值 λ 使方程组 Ax=λx 即 (A-λIn )x=0 有非零解向量(Solution Vector)x，则称λ为方阵A的特征值，而非零向量x为特征 值λ所对应的特征向量，其中In 为n阶单位矩阵。 由于根据定义直接求矩阵特征值的过程比较复杂，因此在实际计算中，往往采取一些 数值方法。本章主要介绍求一般方阵绝对值最大的特征值的乘幂(Power)法、求对称方阵特 征值的雅可比法和单侧旋转(One-side Rotation)法以及求一般矩阵全部特征值的 QR 方法及 一些相关的并行算法。 1.1 求解矩阵最大特征值的乘幂法 1.1.1 乘幂法及其串行算法 在许多实际问题中，只需要计算绝对值最大的特征值，而并不需要求矩阵的全部特征值。 乘幂法是一种求矩阵绝对值最大的特征值的方法。记实方阵A的n个特征值为λi i=(1,2, …,n)， 且满足： │λ1 │≥│λ2 │≥│λ3 │≥…≥│λn │ 特征值λi 对应的特征向量为xi 。乘幂法的做法是：①取n维非零向量v0 作为初始向量；②对于k=1,2, …,做如下迭代：  直至 u  k ?1 ? ? uk uk =Avk-1 vk = uk /║