41第一节 幂法和反幂法.pptVIP

例2 已知矩阵A的一个近似特征值??5，用原点平移法求特征值?及其对应的特征向量，列表计算三次，取x0=(1,0,0)T ，保留两位小数. 其中 对B应用反幂法有 解 取p=5平移有 对B 进行LU分解得 其中 取 x0=(1,0,0)T, 迭代3步的计算结果见下表 k mk 1 0.50, -3.00, 1.00 -3 -0.17, 1, -0.33 2 -1.92, 7.83, -2.50 7.83 -0.25, 1, -0.32 3 -1.94, 8.10, -2.56 8.10 -0.24, 1, -0.32 可得B的按模最大特征值为 ??1/m3?1/8.10?0.12 ?=?+5?-0.12+5=5.12 对应的特征向量为 　　 v?(-0.24, 1, -0.32)T 故得A的一个特征值为 得对应的特征向量为 v?(-0.24, 1, -0.32)T 数学学院 信息与计算科学系 数学学院 信息与计算科学系 第四章矩阵特征值及特征向量的计算 工程技术中有多种振动问题, 如桥梁或建筑物的振动, 机械零件、飞机机翼的振动, 及一些稳定性分析和相关分析在数学上都可转化为求矩阵特征值与特征向量的问题. 问题的来源 二、反幂法 一、幂法 第一节 幂法与反幂法 三、原点平移法 幂法——计算矩阵的按模最大的特征值和相应 特征向量的一种向量迭代法。 适用矩阵——大型稀疏矩阵. 一、幂法 特征值问题及性质 　　设矩阵A?Rn×n，特征值问题是求??C和非零向量x?Rn，使 Ax=?x 其中?是矩阵A的特征值，x是矩阵A属于特征值?的特征向量. 　　求A的特征值问题等价于求A的特征方程 的根；求A的属于特征值?的特征向量等价于求 非零解. 设?为A∈Rn×n的特征值, x 称为A的与特征值 ? 相对应的一个特征向量,即Ax= ?x, (x≠0)则有 (1) cx (c≠0为常数)也是 A的与特征值 ? 相对应的一个特征向量，即A(cx)=?(cx)； (2) ?k为Ak的特征值, 且Akx=?kx ; 设n 阶矩阵A的特征值和相应特征向量为 1. 幂法的基本思想 v1 ,v2,…, vn 其中要求v1 ,v2,…, vn 线性无关，且有 即满足 称?1为A的按摸最大特征值(也称主特征值). 任取非零向量 　　　　　　　　　　, 则 设 a1?0, 由A构造向量序列{xk} 其中 就得到了向量序列 即对 用公式 　　由于 xk=Axk-1=Akx0 ，所以这种生成向量序列的方法就称为幂法. 将向量 xk 整理改写一下为 由于 所以得 2. 收敛性及幂法迭代公式 前面已知 故有 即有 就是 ?1 的特征向量. 而用向量 xk+1 的第i个分量与向量 xk 的第i个分量之比的极限 就是按摸最大特征值?1 . 因此当k充分大时，就可以近似得到按摸最大的特征值和对应的特征向量为 这样就得到两个极限 从而我们就得到幂法的迭代公式 幂法的迭代公式为 当k充分大时，有 收敛速度取决于比值 ,比值越小,收敛越快. 　　问题是：当|?1|1时, xk 所有不为零的分量将随 k?? 而趋于? , 在计算机计算时会造成“溢出”;而当|?1|1时, xk 的分量将随 k?? 而趋于0 , 成为0向量. 幂法的迭代公式为 当k充分大时, 有 3. 误差分析 　　为了避免这个问题，对每一步的xk进行规范化，这样就得到了修改计算公式为 4. 实用计算公式 当k充分大时, 有 其中 mk 是向量 yk中绝对值最大的第一个分量.这时xk分量的模最大为1. 例1 求矩阵 的主特征值与其对应的特征向量。 解　迭代公式为 取 x0=(0,0,1)T , 则有 直到k=8 时的计算结果见下表 k mk 1 2, 4, 1, 4 0.5, 1, 0.25 2 4.5, 9, 7.75 9 0.5, 1, 0.8611 3 5.7222, 11.4444, 8.361 11.4444 0.5, 1, 0.7360 4 5.4621, 10.9223, 8.2306 10.9223 0.5, 1, 0.7536 5 5.