数值计算方法课件-CH2 解线性方程组的直接法—2.1～2.3 Gauss消去法.pptVIP

华长生制作 华长生制作 回代 停机 . 8 输出无解信息 . 7 . )), ( , ), 2 ( ), 1 ( ( . 6 并转8 输出解 n x x x x T L = ) ( ) , ( / k x k k A S T S S n k A T - + ) 1 , ( S j x j k A S T + ) ( * ) , ( n k j , , 1 L + = 对于 注意：为提高算法效率，下列矩阵和向量共用了相同的存储空间： * * 第二章 解线性方程组的 直接法 线性方程组 { 稠密和稀疏（按系数矩阵含零元多少分） 高阶和低阶（按阶数的高低分） 对称正定、三对角占优等（按系数矩阵的形 状性质分） 基 本 解 法 { 直接法（通过有限步计算得到精确解，适用于低阶线性 方程组，第二章内容） 迭代法（通过逐次迭代逼近得到近似解，适用于大型线 性方程组和非线性方程，第六章内容） §2.1 直接法与三角形方程组求解 线性方程组及方法分类 （2.1） 求解线性方程组： AX = b 根据Cramer(克莱姆)法则, 若 则方程组　　 有唯一解。 三角形线性方程组解法 （2.3） 上三角形 §2.2 高斯消去法 Gaussian Elimination 高斯消去法 消元 回代 ＋ ＝ 化上三角方程组的过程 思路 首先将A化为上三角阵 ，再回代求解 = 自下而上解上三角形方程组的过程 经过n-1次 以上求解线性方程组的方法称为高斯(Gauss )消去法 利用 同 解 假设 ，对增广矩阵作行初等变换: 执行过程 上三角阵 通过解 得原方程组 的解。 1. 不改变行列式的值 2. 不改变未知量次序 设 ，计算因子 将增广矩阵 第i 行 ? mi1?第1行，得到： 其中 记 消元过程 行乘数 Step 1： 设 ，计算因子 Step k： 其中 行乘数 将增广矩阵 第i 行 ? mik?第k行，得到： 共进行 ？ 步 n ? 1 ? 以上消元过程均假定 注: 如果 怎么办 由于 ，则 A 的第一列中至少有一元素不为零。 如果 ，则将 的第1行与第 行交换后再 消元，得 假设 例： 当 时，采取类似的处理措施。 回代过程 由于 ，所以 即得线性方程组 Ax = b 解： 因此，上三角形方程组 有唯一解。 No unique solution exists. Then we must find the integer k ? i with , and interchange the k-th row with the i-th row. No unique solution exists. 注: 事实上, 只要 A 非奇异, 即 A?1 存在, 则可通过逐次消元及行交换, 将方程组化为三角形方程组, 求出唯一解。 定理 若A的所有顺序主子式 均不为0，则高斯消元无需换行即可进行到底，得到唯一解。 What if ? ? What if ? ? What if we can’t find such k ? ? 由于计算机中乘除运算的时间远远超过加减运算的时间，故估计某种算法的运算量时，往往只估计 乘除次数，而且通常以乘除次数的 最高次幂 为运算量的 数量级。 (n ? k) 次 Step k：设 ，计算因子 且计算 共进行 n ? 1 步 (n ? k)2 次 (n ? k) (n ? k + 2) 次 消元乘除次数： 1 次 (n ? i +1) 次 回代乘除次数： Gauss消去法的运算量 (n ? k) 次 n-1 步回代过程需作乘除运算总次数为： Gauss消去法的乘除运算总次数为： 完成 n-1 步消元需作乘除运算总次数为： 例如 而如果用 Cramer 法则的乘除法运算次数约为： 当 n 很大时， 当 n = 20 时，Gauss