第6章线性代数方程组的直接法学案.ppt

线性方程组数值解法的分类 ?直接法（适用于中等规模的n阶线性方程组） ◆ Gauss消去法及其变形 ◆矩阵的三角分解法 求解线性方程组的（Doolittle ）三角分解法 通过比较法直接导出L 和 U 的计算公式。 思路 ? 矩阵三角分解的实现 Doolittle三角分解： ? 直接三角分解法解AX = b的计算公式 对于r = 2, 3, …, n （2）计算U的第r行元素 （3）计算L的第r 列元素 (r ? n) (1) (4) (5) 例 用矩阵的Doolittle 三角分解法解方程组 ? Doolittle分解法的变形 紧凑格式的Doolittle三角分解法 例 所以 紧凑格式的列主元Doolittle分解法 例 §4 平方根法 1．矩阵的LDR分解 定理3：如果n阶矩阵A的所有顺序主子式均不等于零， 则矩阵A存在唯一的分解式A = LDR，其中L和R分别是 n阶单位下三角阵和单位上三角阵，D是对角元素不为零的n阶对角阵，上述分解称为A的LDR分解。 2．平方根法 如果A为对称正定矩阵，则存在一个实的非奇异下三角矩阵，使A=LLT ，且当限定的对角元素为正时， 这种分解是唯一的，称为矩阵A的cholesky分解。 定理4：（对称正定矩阵的三角分解） 将对称 正定阵 A 做 LU 分解 U = uij = u11 ui