工程数值解线性方程组的直接法讲述.pptVIP

例7: 求矩阵A的三角分解 解：按紧凑格式 §2.2.2 Doolittle (杜利特尔)分解法 如果方程组 的系数矩阵 A 能分解成 其中 L 是下三角矩阵， U 是上三角矩阵， 这时方 程组就可化为两个容易求解的三角形方程组 先由 解出向量 y , 再由 解出向量 y 这就是原方程组 的解。 (2.4) 矩阵 A 分解成(2.4)的形式称为矩阵的三角分 解。 解方程组的直接三角分解法 定理 1 矩阵 有唯一的Doolittle分解的 充分必要条件是 A 的前 个顺序主子式 定理 2 矩阵 有唯一的Crout分解的充 分必要条件是 A 的前 个顺序主子式 若 L 是单位下三角阵， 则称相应的分解为 Doolittle分解(也称为LU分解)。 若 U 是单位上 三角阵， 则称相应的分解为Crout分解。 设矩阵 非奇异， 且 A 的前 个顺序 主子式都不为零。 则 A 有Doolittle分解 (2.5) 由分解式(2.5)及矩阵的乘法知 当 时有 §2.2.2 Doolittle (杜利特尔)分解法 Doolittle分解算法 对 对 Doolittle分解表上作业法 第1步 第2步 第n步 利用LU分解求解线性方程组的算法 先求解 即 所以 再求解 即 回代求解 例8： 利用Doolittle分解求解以下方程组 解： 回代求解x 设A为n 阶方阵, 若A的前(n?1)阶顺序主子式 均不为0, 则矩阵A可以唯一分解为 其中 L 为下三角阵, U 为单位上三角阵。 §2.2.3 Crout(克劳特)分解法 设n元线性方程组Ax=b的系数矩阵A为非奇异 的三对角矩阵 这类方程组具有许多明显的应用背景， 这种方程组称为三对角线性方程组。 分方程数值解、三次样条函数等问题中， 在求微 都会遇 到这样的线性方程组。 其中 §2.2.4追赶法求解三对角线性方程组 设n元线性方程组 的系数矩阵A为非奇异 的三对角矩阵。 定理: 三对角矩阵A= 计算公式: 计算顺序: 计算公式: 解三对角方程组的追赶法: 例: 用追赶法解方程组 解 ? 追赶法解三对角方程组 /* Crout Reduction for Tridiagonal Linear System */ Step 1: 对 A 作Crout 分解 直接比较等式两边的元素，可得到计算公式。 Step 2: 追——即解 ： Step 3: 赶——即解 ： 与G.E.类似，一旦 ?i = 0 则算法中断，故并非任何 三对角阵都可以用此方法分解。 总结： 定理 若 A 为对角占优 /* diagonally dominant */ 的三对角阵，且满足 ，则追赶法可解以 A 为系数矩阵的方程组。 Hey, what does diagonally dominant mean??? It means that the diagonal entries of the matrix are very LARGE. Well, how large is LARGE? They satisfy the following inequality: 注： ? 如果 A 是严格对角占优阵，则不要求三对角线上的所有元素非零。 ? 根据不等式 可知：分解过程中，矩阵元素不会过分增大，算法保证稳定。 ? 运算量为 O(6n)。 总结： 例9 利用追赶法求解以下方程组 解 * 第 k 步：在矩阵A(k)的第k列方框内选取绝对值最大的元素作为主元素；并通过行的互换将它换到第k行的位置，然后进行第k次消元，得矩阵 (A(k) | b(k)) → (A(k+1) | b(k+1)) 设用列主元消去法已完成 k-1 步，即有 2 .回代过程 则计算停止， 若 对 否则 注：列交换改变了 xk 的顺序，须记录交换次序，解完后再换回来。 例4: 全主元高斯消去法： 第 k 步消元时选 A(k) 中绝对值最大的元素为主元，即 ① 先选取全主元： if ik ? k then 交换第 k 行和第 ik 行； if jk ? k then 交换第 k 列和第 jk 列； ③ 消元 列交换改变了 xi 的顺序，须记录交换次序，解完后再换回来。 全主元高斯消去法具有很好的稳定性，但选全主元比较费时，故在实际计算中很少使用。 §2.1.3 全主元Gauss消去法 设增广矩阵为 第一步：首先在A