第四章 第二节 矩阵三角分解法.ppt

* * 第二节 矩阵三角分解 法 一、 分解法 二、平方根法 1、 平方根法 2、 分解法 三、追赶法 设矩阵A存在LU分解，即有 比较两端的第一行，有 一、 分解法 比较两端的第一列，有 再比较两端第二行的其余元素得 再比较两端第二列的其余元素得 从上面的计算过程，不难归纳出一般的计算公式 再对 计算 同时，由 ，按上边顺序依次求得： 这个分解过程叫做 分解，用这种方 法求解方程组所需要的计算量和用 消元法 计算需要的计算量基本相同，但这种方法把 对系 数矩阵的计算和对右端项的计算分开了，这就使 我们在计算系数矩阵相同而右端项不同的一系列方程组时变得特别方便。 与此类似，如果A存在LU分解，则有 其中， 是以 为对角元素的对角矩阵， 是单位上三角矩阵， 是下三角矩阵，和 分解一样,我们可以依次得矩阵 和 的元素 和 ， 从而将矩阵A分解为 这样的分解称为 分解。 例2??????用 分解求解方程组 和 解 容易验证系数矩阵A的顺序主子式全不为零， 故A存在唯一的 分解，由 得 再由 可计算出 对第一个方程组，由 和 可求得 再由 和 ，求得解 对第二个方程组，由 和 得 由 和 得 二、平方根法 求解对称正定方程组，在许多工程实际计算问题中经常遇到，由于其本身的特点，使用前述的解法是不利的，应该采用适合其特点的解法。这里介绍的 分解法，其运算量和存储量较LU分解法均节省一半左右，是目前在计算机上解这类问题最有效的方法之一。设A为对称正定矩阵，则A的各阶主子式均大于零，于是有如下三角分解 其中 为下三角矩阵， 为单位上三角矩阵。 式中 为单位下三角矩阵， 为非奇异的对角 矩阵，记为 于是 成为 因为 ，于是有 由于 的对角元素皆为非零元素，则 也有分解式 这里， 和 同为下三角矩阵， 和 同为单位上三角矩阵，于是据分解式的唯一性，有 ，故 若用 表示对角矩阵 ， 则有 从而 式又可以表为 其中 已不是单位下三角矩阵了。 我们称 和 为对称正定矩阵A的 分解，下面据这两种分解式来讨论方程组 的计算方法，从分解式 出发而得到的计算方法称为平方根法； 从分解式 出发而得到的计算方法称为 分解法（改进的平方根法）， 这两种方法统称为 分解法 。 其中 和 分解一样，比较 式两端对应元素，可以依次算出 ，今设从第1列到第k-1 列的元素均已算出，现在来计算L的第k列元素 由关系式 ，有 1、 平方根法 设对称正定阵A具有形如 的分解： 于是当 时有 当 时，有 对 ，重复n次 的计算后， 便求得矩阵L。 这样，求解对称正定方程组 可化为解三角形方程组 及 ，它们的求解公式为： 由 式可以看出