非对称特征值问题6-3说课.pptVIP

第三节 加速方法 1、原点平移法 由前面讨论知道，应用幂法计算A的主特征值的收敛速度主要由比值 r=|λ2/λ1|来决定，但当r 接近于1时，收敛可能很慢. 这时，一个补救办法是采用加速收敛的方法. 其中p为参数，设A的特征值为?i，则对矩阵B的特征值为?i-p ，而且A, B的特征向量相同. 引进矩阵 B=A-pI . 如果要计算A的主特征值?1, 只要选择合适的数p，使?1-p为矩阵B=A-pI 的主特征值，且 那么，对矩阵B=A-pI应用幂法求其主特征值?1-p, 收敛速度将会加快. 这种通过求B=A-pI的主特征值和特征向量，而得到A的主特征值和特征向量的方法叫原点平移法. 对于A的特征值的某种分布，它是十分有效的. 　　例 设A∈R4×4有特征值 比值r=|λ2/λ1|≈0.9. 做变换 B=A-12I (p=12), 则B的特征值为 应用幂法计算B的主特征值μ1的收敛速度的比值为 且使收敛速度的比值 设A的特征值都是实数，且满足 则对实数p，使矩阵A-pI的主特征值为?1-p或?n-p时，当我们计算?1及x1时，首先应选取p使 显然，当?2-p=-(?n-p )时，即 P=(?2+?n)/2＝P* 时ω为最小值，这时收敛速度的比值为 当A的特征值都是实数，满足 且?2, ?n能初步估计出来，我们就能确定P*的近似值. 例 用原点平移加速法求矩阵A的主特征值与其对应的特征向量. 对B应用幂法，仍取 y0=(0,0,1)T , 则 解 取p=-2.5, 做平移变换B=A-pI，则 k 1 2, 4, 3.5 4 0.5, 1, 0.875 2 7, 14, 10.5625 14 0.5, 1, 0.7545 3 6.76, 13.5179, 10.1406 13.5179 0.5, 1, 0.7507 4 6.7503, 13.5007, 10.1256 13.5007 0.5, 1, 0.7500 5 6.7500, 13.5000, 10.1250 13.5000 0.5, 1, 0.7500 可得到B的主特征值为 ?1?13.5000, 主特征向量为 v1 ? (0.5 ,1.0, 0.7500)T , 因此，A的主特征值为 ?1 = ?1 +p ? 11.0000, 主特征向量仍为 x1 =(0.5,1,0.7500)T . 若知道某一特征根 ?i 的大致位置 p ，即对任意 j ? i 有| ?i ? p | | ?j ? p | ，并且如果 (A ? pI)?1存在，则可以用反幂法求(A ? pI)?1的主特征根 1/(?i ? p ) ，收敛将非常快。 设A∈Rn×n为对称矩阵，称 为向量x的瑞利商，其中(x, x)=xTx为内积. 由定理11知道，实对称矩阵A的特征值?1及?n可用瑞利商的极限值表示. 下面我们将瑞利商应用到用幂法计算实对称矩阵A的主特征值的加速上来. 2、对称矩阵的Rayleigh商加速 定理 设A∈Rn×n为对称矩阵，特征值满足 对应的特征向量vi满足(vi, vj)=δij (单位正交向量) ，应用幂法公式计算A的主特征值?1，则 由此可见，R(xk) 更快的收敛于?1. 证明 得 幂法的瑞利商加速算法： 其中A为n阶实对称矩阵. 对给定的误差限?，当|μ k–μk-1|?时，取近似值