QR方法求矩阵全部特征值.pdfVIP

数 值 分 析 课 程 设 计 QR方法求矩阵全部特征值 问题复述 用 算法求矩阵特征值： （i） （ii） 要求： （1）根据 算法原理编制求（i）与（ii）中矩阵全部特征值的程序并输出 计算结果（要求误差 ） （2）直接用现有的数学软件求（i），（ii）的全部特征值，并与（1）的结果比较。 问题分析 QR方法是目前公认为计算矩阵全部特征值的最有效的方法。它适用于任一种实矩。QR方法的原理是利 用矩阵的正交分解产生一个与矩阵A相似的矩阵迭代序列，这个序列将收敛于一个上三角阵或拟上三角 阵，从而求得原矩阵A的全部特征值。 QR是一个迭代算法，每一步迭代都要进行QR分解，再作逆序的矩阵乘法。要是对矩阵A直接用QR方 法，计算量太大，效率不高。因此，一个实际的QR方法计算过程包括两个步骤，首先要对矩阵A用初等 相似变换约化为上Hessenberg矩阵，再用QR方法求上Hessenberg矩阵的全部特征值。 示意如下： 对B矩阵的约化只需将每列次对角线上的元素约化为0。因此常常用平面旋转阵（Givens变换阵）来进行 约化。 一、QR方法原理及收敛性 设 已实现了QR分解，记 其中 是正交阵， 是上三角阵。因