新计算方法 教学课件 徐士良 第三章.pptVIP

确定 的具体方法如下： 从矩阵 的右下角二阶子矩阵 中解出两个特征值 和 。若令 = + ， = ，则有 由于 是上H矩阵，可以验证矩阵 = 只有前三个元素非零，即矩阵的第1列为 = 其中 ， ， 的计算公式如下： 现在要构造一个正交对称 ，使 经线性变换 后能够消去其中的元素 和 。可以验证，如下形式的 可以满足这个要求： 其中 为3×3阶的矩阵，即 其中 最后，用 对上H矩阵 作相似变换，其结果为如下形式： 其中有上划线的为次对角线以下新增加的三个非零元素。 (2) 利用同样的方法，依次构造正交（且是对称的）矩阵 ， ，…， ，分别用它们对 ， ，…， 作相似变换： = ， i=1，2，…，n-2 最后得到上H矩阵 = 在这个过程中，一般有 若令它的第i列为 ，即 = 则需要构造的正交对称 应使 经线性变换 后能够消去其中的元素 和 。可以验证，如下形式的 可以满足这个要求： 其中 为3×3阶的矩阵，即 其中 。这就是 (i＝0，1，…，n－2)的一般形式。 其中参数 ， ， 由前面的计算公式给出，而一般的 ， ， (i＝1，2，…，n－2)是从相应的 的第i列中取得，即 , , 当i＝n－2时，有 ， ， 反复进行以上各步，直到将上H矩阵变为对角块全部是一阶块或二阶块为止，此时就可以直接从各对角块中解出全部特征值。在实际计算过程中，总是将矩阵分割成各不可约的上H矩阵，这样可以逐步降低主子矩阵的阶数，减少计算工作量。 QR方法是一种迭代的方法。在迭代过程中，如果次对角线元素的模小到一定程度，就可以把它们当作零看待。通常采用的判别准则有 | |≤εmin{| |，| |} 或 | |≤ε(| |＋| |) 因为 为对称矩阵， 为正交矩阵，所以 也为对称矩阵。若要求矩阵 的元素 = 0，则只需令 即 考虑到在实际应用时，并不需要解出θ，而只需要求出sin2θ，sinθ与cosθ就可以了。因此，令 则可以得到 最后将sin2θ，sinθ与cosθ代入，就可以得到矩阵 的各元素。 由 的各元素还可以得到如下等式： 其中i,j≠p,q 。由此可以得到矩阵 的非对角线元素的平方之和 以及矩阵 的对角线元素的平方之和 由此可以看出，对称矩阵 经变换后（该变换称为旋转变换），就将选定的非对角线元素（一般选绝对值最大的）消去了，且其对角线元素的平方之和增加了2 ，而非对角线元素的平方之和减少了2 ，矩阵总的元素平方之和不变。但经过这样的变换之后，非对角线上的其他零元素就往往不再是零了。 每经过一次旋转变换，其矩阵的非对角线元素的平方之和总是“向零接近了一步”。如果反复选取绝对值最大的 ,经过旋转变换后，就可以逐步将矩阵 变为对角矩阵。实际上，作为一个迭代过程，只要满足一定的精度要求就可以了。 综上所述，可以得到用雅可比方法计算对称矩阵 的特征值的步骤如下： (1) (2) 选取非对角线元素中绝对值最大者 。若| |＜ε，则 (i＝0，1，…，n-1)即为特征值 以及 第i列即为与 对应的特征向量；否则继续做下一步。 (3) 按下列公式计算sin2θ，sinθ与cosθ：