特征值问题数值计算上机实验.docxVIP

特征值问题数值计算上机实验练习6.16：取初始向量为,用乘幂法计算主特征值及其相应的特征向量。请绘制主特征值误差下降曲线，以及特征子空间距离的下降曲线。然后采用Atiken加速技巧和Rayleigh 商技术分别对算法进行加速，并完成类似的工作。解：采用教科书上8.1算法编写程序。在n=100时，由于在主特征值对应的特征向量上投影为0，如果取，若无舍入误差的影响，就会导致乘幂法求解出的特征值为第二主特征值。故这里取初始向量为随机向量。为了便于比较，在每一步迭代中，将Atiken加速值和Rayleigh商加速值分别计算出来，并作出如下图像：（下图显示的为从10001步后开始）从上面两图可以看出，初始阶段乘幂法会剧烈震荡，在迭代一段时间后会突然特征值误差变大。在迭代6000次后，震荡现象不再出现。从第二张图中可以看出，采用Atiken和Rayleigh商加速得到的误差下降较快。由于Atiken加速和Rayleigh商加速均是对特征值进行的，对特征向量无影响，故特征子空间的距离的下降情形一致。下图为特征子空间距离的下降曲线。n=101时情况类似。练习6.17用反幂法求解最靠近2的特征值，及其对应的特征向量。观察是否有所谓的“一次收敛”的特性。解：当n=100时，的靠近2的特征值有1.9688963761593和2.0311036238407，可以看出特征值并不非常靠近2，预期反幂法的“一次