数值分析 第7讲幂法和反幂法.ppt

第七讲矩阵特征值与特征向量的计算（1） ----幂法和反幂法 应用背景 求系统的固有振动频率 求机械系统的振动频谱（模态） 系统的稳定性分析 物理学中的某些临界值的确定 常用的数值求特征值的方法 幂法与反幂法 Jacobi法 QR分解法 ……. 作业 教材第65页：习题1、3，6，7 课外阅读：《C数值算法》第11章 数值分析 朱立永 北京航空航天大学 数学与系统科学学院 Email: numerical_analysis@ Password:beihang 答疑时间：星期四下午2:30－5:30 答疑地点：主216 第三章 矩阵特征值与特征向量的计算 幂法是一种求实矩阵A的按模最大的特征值λ1及其对应的特征向量x1的方法。特别适合于大型稀疏矩阵。 QR方法是一种变换方法，是计算一般(中小型)矩阵全部特征值问题的最有效方法之一. Jacobi是一种求实对称矩阵A的全部特征值及其对应的特征向量的方法。特别适合于中小型对称矩阵。 幂法与反幂法 我们来看下面的一个过程： 任取一个非零的n维向量u0，从u0出发，按照如下的递推公式 uk＝Auk-1 产生的向量序列 {uk}，我们来分析这个序列的收敛情况。 假设实矩阵A具有n个线性无关的特征向量x1, x2, …, xn 相应的特征值λ1, λ2, …, λn满足| λ1|| λ2|≥|λ3|≥…≥|λn|，于是存在不全为零的数组 α1, α2, …,αn,使得 当α1 ≠0时有 即 uk 可以近似的作为属于λ1的特征向量！ 几点注记： 在上述假设条件下λ1必为实数； 如果α1＝0将会怎样？ 该方法的本质是 uk＝Aku0，因此得名幂方法； 为防止 uk过大或过小，我们在每一步迭代中，对迭代向量单位化。 实际使用的迭代公式为： 于是可得 不同范数选取下的特征值的计算 1. 取范数为2－范数时 对应的迭代公式 1. 取范数为∞－范数时 对应的迭代公式 设uk＝(h1, h2, …, hn)T，其最大的分量为hr，令 1. A的最大特征值是m重时，幂方法仍然有效 2. 幂方法的收敛速度与|λ1/λ2|有关 反幂法 假设实矩阵A具有n个线性无关的特征向量x1, x2, …, xn 相应的特征值 λ1, λ2, …, λn 满足 | λ1|≥| λ2| ≥|λ3|≥…≥|λn-1| |λn| 对于 A-1 应用幂方法 取范数为2－范数时的迭代公式 对应的迭代公式 同理，我们可以对(A-pI)-1用幂法－称为带原点平移的幂法 取范数为2－范数时的迭代公式 对应的迭代公式 类似地，我们可以对(A-pI)用幂法得到A的一个特征值和特征向量 对应的迭代公式