§-2-------计算实对称矩阵特征值的同时迭代法.ppt

同时迭代法又称块乘幂法或子空间迭代法.它能同时求出求出一个实对称矩阵的几个特征值和特征向量.n阶实对称矩阵的特征值均为实数,且必存在一个标准特征向量系其中是与相应的特征向量.现设且.同时迭代法的计算步骤如下:取m个初始近似向量组成一个阶列直交矩阵V0,即有,对,执行以下各步:§2计算实对称矩阵特征值的同时迭代法

(1)计算(2.1)(2)计算阶实对称矩阵:(2.2)(3)计算Bk的特征值,设其排列次序为共计算其特征向量矩阵Wk,取它为一个直交矩阵;(4)计算UkWk;(5)对UkWk作直交分解:(2.3)其中Rk为上三角矩阵,Vk为阶列直交阵;(6)检验相邻两次迭代得到的和之间的差是否满足精度要求:若

其中TOL为预先给定的误差容限,则取作为的近似值,Vk+1的第j列向量作为与相应的近似特征向量;否则,对k增加1,转回到(1)继续进行迭代.下面讨论同时迭代法的收敛性.由于A为实对称矩阵,记,则记设(2.4)其中C为阶矩阵,并把C分块成

其中Ca为阶矩阵,其中Cb为阶矩阵,则设都非奇异,据(2.3)和(2.1)式有(2.6)其中从而,据(2.5)式就有

(2.7)令(2.8)若Ca非奇异,则