第三章 线性代数方程组的解法2.pptVIP

* §4 矩阵分解法 /* Matrix Factorization */ 我们知道对矩阵进行一次初等变换,就相当于用相应的初等矩阵去左乘原来的矩阵。因此我们这个观点来考察Gauss消元法用矩阵乘法来表示，即可得到求解线性方程组的另一种直接法：矩阵的三角分解。 Gauss消元法的矩阵形式 一、 三角分解法 本质 它是基本Gauss消元法的一种等价变形 ? Gauss消元法的矩阵形式 /* Matrix Form */ 为上三角矩阵 其中 为单位下三角 矩阵分解为单位下三角 和上三角矩阵的乘积 ?矩阵 分解的条件 如果矩阵 可以分解为一个单位下三角阵 和 一个上三角阵 的乘积，即 ，则称 此分解为LU分解，或三角分解。 定义1、 n阶方阵 的左上角 子矩阵的行列式 称为A的p阶顺序主子式。 定理1、 n阶方阵A有唯一LU分解的充分必要条件是： 矩阵A的各阶顺序主子式均不为零。 设 ?矩阵 分解的计算公式 根据矩阵的乘法，由 对应元素相等，得到 分解计算过程: Step1 Step2 Step3 Step4 Step5 Step6 Step2n-1 Step2(n-1) 先计算 的行 再计算 的列 依次交替进行 rk ?利用矩阵 分解求解方程组的思想 记 方程组可化为下面两个易求解的三角方程组 ?方程组求解的计算公式 先求解方程组 再求解方程组 工作量 例4：用三角分解法求解下列方程组 解： 系数矩阵 二、平方根法 当 为实对称正定矩阵时，三角分解法的变形。 ?实对称正定矩阵的几个重要性质 ? A?1 亦对称正定，且 aii 0 ? A 的顺序主子阵 Ak 亦对称正定 ? A 的特征值 ?i 0 ? A 的全部顺序主子式 det ( Ak ) 0（充要条件） 如果 是正定矩阵，则存在唯一的对角元素为正数的 下三角矩阵 ，使得 。 则 的所有顺序主子式为正 矩阵 存在LU分解 易证 其中 为 的主对角元素，且有 定理2、 证明： 设 正定, 单位上三角 记 由矩阵A 的LU分解的唯一性知： 其中 思想 ?平方根分解的计算公式 设 由 对应元素相等得 平方根分解的公式 回代公式见教材。 因对称性无需存储 Step1 Step2 Step3 Stepn 的计算过程: 逐 列 计 算 元素 仍然存放在矩阵 的相应位置上