2007用—矩阵特征值与特征向量的计算3.1.pptVIP

计算方法 第三章 反幂法的求解步骤 分解计算 　 反幂法迭代 (1)　 (2) k=1,2,… 1) 2) 　　　　　　　 求 求 * 设 x=∑ (αi)*(vi) 则 (Ax,x)=∑λi(αi)**2 (x,x)=∑λi * * * 7.若矩阵A=(aij)满足 则成A为严格对角占优矩阵.试证严格对角占优矩阵A经过一步Gauss消元后,所得的A(2)仍为严格对占优矩阵. 证明: 第三章 矩阵特征值 与特征向量的计算 计算方法 第三章 §3.1 乘幂法及其变体 §3.2 子空间迭代法 §3.3 Jacobi旋转法 §3.4 Householder方法 §3.5 QR算法 计算方法 第三章 §3.1 乘幂法及其变体 3.1.1 乘幂法 3.1.2 反幂法 3.1.3 乘幂法的加速 计算方法 第三章 设A为n阶方阵， 若 有数?使得 则称?为A的特征值，x为相应于?的特征向量。 (3.1.1) 计算方法 第三章 特征问题的求解包括 求特征值?，满足 求特征向量 满足齐次方程组 (3.1.2) (3.1.3) 称为A的特征多项式 计算方法 第三章 设 为　　　的特征值且 ,其中x≠0，则 1) 为cA的特征值（c为常数c≠0 ） 2)　　 为　　 的特征值，即 3) 为 Ak的特征值 计算方法 第三章 4)设A为非奇异矩阵，那么 且 为A-1 特征值，即 5) 6) 计算方法 第三章 7)设　　　为对称矩阵(其特征值次序记为　　　　　　),则对于任意非零 计算方法 第三章 乘幂法 用于求大型稀疏矩阵的主特征值的迭代方法。 其特点为公式简单，易于在计算机上实现。 计算方法 第三章 乘幂法的计算公式 设 取初始向量 令 …… 由此可得向量序列 (3.1.4) 计算方法 第三章 由(3.1.4)，有 (3.1.5) 称此方法为乘幂法 (3.1.4)或(3.1.5)称为乘幂公式 称为乘幂序列 计算方法 第三章 定理3.1 设 有完全特征向量系， 若?1,?2, …,?n为A的n个特征值，满足 |?1| ? |?2| ? … ? |?n| 对任何初始向量 由乘幂公式 所确定的迭代序列 有下面的结论： 计算方法 第三章 定理3.1 (1) 当|?1| |?2|时，有 此极限过程的收敛速度取决于 的程度。 r越小收敛越快。 r接近于1时，收敛很慢。 且 （当k充分大时）可取作为相应于?1的近似特征向量。 计算方法 第三章 (2) 当|?1| = |?2| |?3|时，有 ① 若?1= ?2，则主特征值?1及相应特征向量的求法同(1)。 ②若?1= －?2，则有 此极限过程的收敛速度取决于 的程度。 向量 分别作为主特征值?1，?2相应的近似特征向量 计算方法 第三章 乘幂法可用于近似计算矩阵按模最大的一个（或几个）特征值以及相应的特征向量 当比值 时，收敛速度快 计算公式简便，便于在计算机上实现。 计算方法 第三章 规范化的乘幂法公式 令 表示向量x个分量绝对值最大者，即 如果有某个i0，使得 则令 规范化的乘幂法公式 计算方法 第三章 规范化乘幂法的步骤 取n维异于0的初始向量x(0) 　一般让x(0)满足　　　　　　即 对于k=0,1,2,……按如下式迭代 计算方法 第三章 终止条件 输出结果 计算方法 第三章 例题 例3.1 用规范化乘幂法计算矩阵A的主特征值及相应的特征向量 计算方法 第三章 乘幂法的加速 引进参数 ，用矩阵 来代替A进行乘幂迭代。 设 为矩阵B的特征值， 则B与A特征值之间应有关系式 设 为矩阵A相应于 的特征向量， 则 也是 的特征向量。 计算方法 第三章 对任何 关于矩阵B的乘幂公式（3.1.5）可表示为 计算方法 第三章 为加速收敛速度，应如此选择参数 ，使 达到最小。 (3.1.29) 计算方法 第三章 原点移位法是一个矩阵