计算方法 教学课件 徐士良 第三章.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 | |，| | 或 | |≤ε | |＋| | 考虑到在实际应用时，并不需要解出θ，而只需要求出sin2θ，sinθ与cosθ就可以了。因此，令 则可以得到 最后将sin2θ，sinθ与cosθ代入，就可以得到矩阵 的各元素。 由 的各元素还可以得到如下等式： 其中i,j≠p,q 。由此可以得到矩阵 的非对角线元素的平方之和 以及矩阵 的对角线元素的平方之和 由此可以看出，对称矩阵 经变换后（该变换称为旋转变换），就将选定的非对角线元素（一般选绝对值最大的）消去了，且其对角线元素的平方之和增加了2 ，而非对角线元素的平方之和减少了2 ，矩阵总的元素平方之和不变。但经过这样的变换之后，非对角线上的其他零元素就往往不再是零了。 每经过一次旋转变换，其矩阵的非对角线元素的平方之和总是“向零接近了一步”。如果反复选取绝对值最大的 ,经过旋转变换后，就可以逐步将矩阵 变为对角矩阵。实际上，作为一个迭代过程，只要满足一定的精度要求就可以了。 综上所述，可以得到用雅可比方法计算对称矩阵 的特征值的步骤如下： 1 2 选取非对角线元素中绝对值最大者 。若| |＜ε，则 i＝0，1，…，n-1 即为特征值 以及 第i列即为与 对应的特征向量；否则继续做下一步。 3 按下列公式计算sin2θ，sinθ与cosθ： 4 按下列公式计算矩阵 的新元素： i,j＝0，1，…，n-1；i,j≠p,q 5 。转第 2 步。 前面提到，按雅可比方法每迭代一次，对称矩阵 的非对角线元素的平方之和就“向零接近一步”。而在每一次迭代时，其旋转变换矩阵 为正交矩阵，即当进行到第m次迭代时， 也为正交矩阵。 因此， 与 具有相同的特征值。由此说明了在一定条件下，从这个迭代过程收敛到 其中 i＝0，1，…，n-1 为对称矩阵 的特征值，且 的第i列为与 对应的特征向量。 2. 雅可比过关法 首先计算对称矩阵 的所有非对角线元素平方之和的平方根，即 然后设置第1道关口r1＝E/n，对 中非对角线元素进行逐行（或逐列）扫描，分别与r1进行比较。若| |＜r1，则让其过关，否则用旋转变换 将 化成零。 需要指出的是，在某次旋转变换中变为零的元素中，在以后的旋转变换中可能又变为非零元素，因此，要重复进行上述的扫描过程，直到约化到对于所有的非对角线元素均满足条件 | |＜r1，i≠j 为止。 矩阵 中所有非对角线元素都过了第1道关口后，再设置第2道关口r2＝r1/n＝E/n2，对 中的非对角线元素再进行逐行（或逐列）扫描，对于不满足条件| |＜r1的所有 用旋转变换 将它化成零。直到 中所有的非对角线元素都满足条件 |aij|＜r2，i≠j 为止。即 中所有的非对角线元素都过了第2道关口。 重复以上过程，经过一系列的关口r1，r2，…，直到对于某个关口满足条件 rk＝E/nk≤ε 为止。其中ε为预先给定的精度要求。 3.5 求对称矩阵特征值的豪斯荷尔德方法 3.5.1 用豪斯荷尔德变换将一般实对称矩阵约化成对称三对角矩阵 在豪斯荷尔德方法中