3.1线性代数方程组的解法——列主元高斯消去算法.ppt

线性方程组的解法 Gauss消去法 顺序Gauss消去算法的实现 3 计算实例 4 进一步讨论 5 列主元Gauss消去法 Mathematical Methods its Applications J. G. Liu School of Math. Phys. North China Elec. P.U. 考虑如下线性方程组 写成矩阵形式 其中 , (1) 由克莱姆法则可知 (1)有唯一的解，而且解为： 其中 若行列式用按行（列）展开的方法计算 ， 当方程组阶数较高时，计算量很大，因此克莱姆法则通常仅有理论上的应用价值，计算线性方程组的解还要考虑数值解法。 求解线性方程组的数值方法分为直接解法和迭代解法： 直接解法 若计算过程没有舍入误差，经过有限次算术运算就能求出方程组（1）精确解的数值方法。 用克莱姆法则求解（1）的计算量为: 迭代解法 若计算过程没有舍入误差，也不能经过有限次算术运算求得方程组（1）的精确解，而只能是逐步逼近的数值方法。 我们先考虑直接方法。 1 基本思想 我们知道，下面有3种方程的解我们可以直接求出： ① n次运算 ② (n＋1)n/2次运算 ③ (n＋1)n/2次运算 众所周知，对方程组作如下变换，解不变， ①交换两个方程的次序； ②一个方程的两边同时乘以一个非零的数； ③一个方程的两边同时乘以一个数，加到另一个方程上。 因此，对应的，对增广矩阵(A，b)，作如下的行变换，解不变， ①交换矩阵的两行； ②某一行乘以一个非零的数； ③某一行乘以一个数，加到另一行。 消去法 就是对增广矩阵作上述行变换，化为可以直接求解的3种方程之一，而后求解。 考虑到： 方程组（1） 增广矩阵 (A，b) 思路: Gauss消去法就是先将(1)的系数矩阵A化为上三角阵，再回代求解。 其中 ——主元素 2 计算步骤 例1 求解方程组 解: 下面看求解n元线性方程组的一般过程， 第一步： 运算量： (n-1)*(n+1) 其中 即 若 ， 第二步： 运算量： (n-2)*n 即 其中 若 ， 类似的做下去， 第k步： 运算量： (n－k)*(n－k＋2)。 则 若 ， 则 n－1步以后，我们可以得到变换后的增广矩阵为： 回代可得原方程的解！ 因此，总的计算量为： 加上 解上三角形方程组的计算量(n+1)n/2，总共为： 注： 这是计算行列式的一种有效数值方法！ (1) (3) 该方法又称顺序高斯消去法。 (2) 求逆阵； 消元过程 回代过程 输入A和b，定义并初始化相关变量； 一个变量即可 例2 用（顺序）高斯消去法求解线性方程组 解： 注意到，计算过程中 处在被除的位置， 所以，Gauss消去法的可行条件为： A的所有顺序主子式均不为0， 因此整个计算过程要保证它不为0！ 8个 8个 如果某个 很小的话，会引入大的误差！ 下面看一个例子： 因此，有些线性方程组，不能用Gauss消去求解! (为什么？) 例3 单精度解方程组 /* 精确解为： 8个 用Gauss 消去法计算： 小主元导致计算失败。 小主元会导致： （1）大数吃掉小数；或者 （2）舍入误差被放大！ 在Gauss消去法第k步之前，做如下的事情（选列主元）： 若 交换k行和j行 行的交换，不改变方程组的解，同时又有效地克服了顺序Gauss消去法的缺陷，具有较好的数值稳定性。 例3 ? —— 列主元Gauss消去法