第3章解线性方程组的直接法(免费阅读).pptVIP

数 值 分 析 ——第3章 解线性方程组的直接法 §1 直接法与三角形方程组求解 1.1 线性方程组的分类 不同类型的方程组由不同的解法： 按系数矩阵零元素多少分：稠密、稀疏（占80℅）. 按系数矩阵阶数分：高阶4000阶以上、低阶 按系数矩阵形状分：对称正定、对角线、对角占优、三角形方程组 §1 直接法与三角形方程组求解 1.2 三角形线性方程组的解法 ① 下三角形：设 若lii ? 0，则 §1 直接法与三角形方程组求解 1.2 三角形线性方程组的解法 ② 上三角形：设 （Ux=b） 若uii ? 0，则 §1 直接法与三角形方程组求解 1.3 直接法概述 ① 高斯消去法 设方程组Ax = b，其中 A = (aij)n×n， 增广矩阵： 其中A(n)为上三角阵，则A(n)x = b(n)与原方程组同解。 此方法称为高斯消去法。 §1 直接法与三角形方程组求解 1.3 直接法概述 ② 三角分解法 若A = LU. 其中L，U为三角阵，则 通常为L下三角，U为上三角。 此方法称为三角分解法。 §2 Gauss消去法 2.1 消元与回代计算 1. 消元 例： §2 Gauss消去法 2.1 消元与回代计算 1. 消元 一般情况，设detA ? 0 ① 第1次消元 2.1 消去与回代计算 ② 第k次消元 其中 i = k+1,k+2,…,n 2.1 消去与回代计算 其中 i = k+1,k+2,…,n 注：为减少计算量，令 i = k+1,k+2,…,n 2.1 消去与回代计算 当k = n-1时，得 当 时， 称为主元素。 2.1 消去与回代计算 2. 回代 注：若A非奇异。上述过程可行，但无法事先判断。故可能 ① 某一步找不到非零元； ② 到最后一步得 ，无法回代，需在计算中中断。 2.1 消去与回代计算 2.2 gauss消去法的运算量 ① 消元过程:第k步中 除法n–k次, 乘法(n–k)(n–k+1)次，共(n–k)(n–k+2)次，共n – 1步： 以乘除法的次数来衡量算法的运算量 2.1 消去与回代计算 2.2 gauss消去法的运算量 ① 消元过程: ② 回代过程： 求xi中，乘n–i次，除1次，共n–i+1次，i = 1, 2, …, n–1. 共计 2.1 消去与回代计算 2.2 gauss消去法的运算量 ① 消元过程: ② 回代过程： 总次数 MD= 注：当n=20时，约为2670次，比克莱姆法则9.7×1020次大为减少。 §3 Gauss列主元素消去法 3.1 主元素的作用 在消元过程中以主元素 依次作除法计算： 注：① 若消元过程中出现akk(k) = 0，则无法继续 ② 若akk(k)≠0，但较小，则小主元做除数将产生大误差 3.1 主元素的作用 ③ 为避免小主元作除数，在A(k)的第k列主对角线以下元素 中挑选绝对值最大者, 并通过行变换，使之位于主对角线上作为主元素，仍记为 ，然后再进行消元计算。即 ——第k列主对角元以下元素绝对值最大者作主元（该行与第k行对调），称每一步都按列选主元的消去法为gauss列主元素消去法 §3 Gauss列主元素消去法 §3 Gauss列主元素消去法 3.1 主元素的作用 例1 用舍入三位有效数字求解线性方程组 （准确解是x1 = 10.0，x2 = 1.00） 解：1) 不选主元的Gauss消去法计算结果： x1 = -10.0，x2 = 1.01，此解无效； 2) 按列选主元的Gauss消去法计算结果： x1 = 10.0，x2 = 1.00. §3 Gauss列主元素消去法 3.2 算法设计 设Ax=b, A = (aij)n×n 1. 存储方式 ① 用二维数组A(1: n，1: n+1)按行存放增广矩阵(A|b). A(i, j): i=1,2,…,n,j=1,2,…,n+1 ② 计算过程中 §3 Gauss列主元素消去法 3.2 算法设计 1. 存储方式 ② 计算过程中 ③ 解向量存入A的最后一列 ： ——回代。 ④ EXT=1作无解标志； §3 Gauss列主元素消去法 3.2 算法设计 2. 算法 输入A，n，eps For (k = 1 to n - 1) 选主元A(I0，