矩阵特征值和特征向量的数值解法.docVIP

PAGE 1 第七章 矩阵特征值和特征向量的数值解法 7.2典型例题精解 例7.2.1 用幂法计算矩阵 的主特征值和相应的特征向量。 解：取初值向量u(0)=(0,0,1)T,由式（7.1.1）得 对k=1,2,…n,计算结果如表7.2.1所示。 表7.2.1 k V(k) mk 0 0.0 0.0 1.0000 1.00000 1 0.5 1.0 0.2500 4.00000 2 0.5 1.0 0.8611 9.00000 3 0.5 1.0 0.7306 11.4400 4 0.5 1.0 0.7535 10.9224 5 0.5 1.0 0.7493 11.0140 6 0.5 1.0 0.7501 10.9972 7 0.5 1.0 0.7500 11.0004 8 0.5 1.0 0.7500 11.0000 从表7.2.1看出，当k=8时，即可达到精度，所以矩阵A的主特征值λ1=11，对应的特征向量x1=(0.5,1.0,0.7500)T。 例 7.2.2 设矩阵 试取平移量,用幂法求出矩阵A的主特征值及对应的特征向量。如果对平移法的结果做一次Rayleigh商加速，求A之特征值和对应的特征向量。 解： 取，则 取u(0)=v(0)=(1,1,1)T,迭代结果如表7.2.2所示。 表7.2.2 k V(k) mk 0 1 1. 1. 1. 1 1 1. 0. 2. 2 1 0. 0. 2. 3 1 0. 0. 2. 4 1 0. 0, 2. 5 1 0. 0. 2. 6 1 0. 0. 2. 7 1 0. 0. 2. 8 1 0. 0. 2. 从表7.2.2知 若对上述结果做一次Rayleigh商加速，则有 与真值相比，误差 例 7.2.3 试用逆幂法求矩阵 矩阵近似于—6.42的特征值和特征向量。 解：分解A+6.42I=LR，得 取y(0)=(1,1,1)T做半迭代，计算结果如表7.2.3所示。 表7.2.3 k v(k) mk 1 —0. —0. 1 —820. 2 —0. —0. 1 —937. 3 —0. —0. 1 —937. 故得 例 7.2.4 试用Jacobi方法解矩阵 的全部特征值和特征向量。 解：令故取 得，则 所以 如此继续，可得 将计算中的元素近似为零，则最后得 故得A的特征值 对应的特征向量 例 7.2.5 设矩阵A的n个特征值互异，对任意的非零向量x0和y0做迭代 （k=0,1,2,…） (1)证明 （2）用上述方法求矩阵 的按模最大特征值。 （1）证明 设A的特征值为,对应的特征向量为。因为 A与AT有相同的特征值，所以设AT的特征值对应的特征向量。 首先证明，当i≠j时，，事实上，由得 再由得 因此 所以当i≠j时,故必有。 其次，设 由迭代公式及幂法知， 因此有 于是 所以成立。 （2）解 取，计算结果为 于是， 与准确值相比，误差为。 例 7.2.6 设方阵A的特征值均为实数，且满足