资源描述：
《计算方法解线性方程组的直接法.ppt》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、第8次线性方程组的直接解法计算方法(NumericalAnalysis)1）高斯消去法2）高斯主元素消去法3）方程组的性态4)高斯消去法算法构造（编程）本讲内容高斯消去法§5.1引言在工程技术、自然科学和社会科学中，许多问题最终都可归结为求解线性方程组的数学问题。线性方程组的求解对于实际问题是极其重要的。解线性方程组的直接法解线性方程组的直接法可简记为Ax=b，其中(6.1)常见的nxn线性方程组，一般形式为线性方程组的数值解法一般有两类：直接法：就是经过有限步算术运算，可求得方程组精确解的方法（若计算过程中没有舍入误差），如克莱姆法则就是一
2、种直接法，直接法中具有代表性的算法是Gauss消去法。迭代法:就是用某种极限过程去逐步逼近线性方程组的精确解的方法。也就是从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法。(一般有限步内得不到精确解)例子：求解如下的上三角线性方程组：解：由（4），得将（4）带入（3），得将结果代入（2），得将结果代入（1），得§5.2高斯消去法§5.2高斯消去法5.2.1高斯消去法的基本思想①②③解:高斯消去法包括如下的消元和迭代的两个过程。先用一个简单实例来说明Gauss法的基本思想例5.1解线性方程组（1）消元过程④⑤第1步:将方程①乘上(-
3、2)加到方程②上去,将方程①乘上加到方程③上去,这样就消去了第2、3个方程的项,于是就得到等价方程组第2步：将方程④乘上加到方程⑤上去，这样就消去了第3个方程的项，于是就得到等价方程组⑥这样，消元过程就是把原方程组化为上三角形方程组，其系数矩阵是上三角矩阵。（2）回代过程将上述三角形方程组自下而上求解得：从而求得原方程组的解：前述的消元过程相当于对原方程组的增广矩阵进行下列行变换同样可得到与原方程组等价的方程组⑥高斯消去法的基本思想：这种求解上三角方程组的方法称为回代,通过一个方程乘或除以某个常数,利用矩阵行的初等变换将原方程组Ax=b系数矩
4、阵化为上三角形矩阵,然后从最后一个方程开始,依次向前代入求出未知变量：将两个方程相加减,逐步减少方程中的变元数,最终将方程组化成上三角方程组,一般将这一过程称为消元,然后再回代求解。…5.2.3高斯消去法的适用条件注2：设系数矩阵A为非奇异矩阵，则若a11=0，则可以通过调换行的方法，使得在第一行的第一个元素非0。其它在消元过程中，kk位置的情形类似处理。则高斯消元法可以进行。注1：设系数矩阵A为非奇异矩阵，直接使用高斯消元法（不进行行的交换）对于某些简单的矩阵可能失败，例如：证明：上三角形方程组是从原方程组出发，通过逐次进行“一行乘一数加到
5、另一行”而得出的，该变换不改变系数矩阵顺序主子式的值。因此，需要对上述的高斯算法进行修改，首先应该研究原来的矩阵A在何条件下能够保证…定理1若方程组系数矩阵的顺序主子式全不为0，则高斯消去法能实现方程组的求解，即：…设方程组系数矩阵，其顺序主子式(m=1,2,…，n)经变换得到的上三角形方程组的顺序主子式所以能实现高斯消去法求解（m=1,2,…，n）…………………...定义5.1设矩阵每一行对角元素的绝对值都大于同行其他元素绝对值之和则称A为严格对角占优矩阵。……上述条件展开以后为：定理1.1若方程组的系数矩阵A为严格对角占优，则用高斯消去法
6、求解时，全不为0。因此，可以使用高斯消去法求解。练习：用高斯消去法求解如下的线性方程组解：增广矩阵为Home高斯主元素消去法使用高斯消去法求解时，在消元过程中可能会出现的情况，这时消去法将无法进行；§5.3高斯主元素消去法即使，但它的绝对值很小时，用其作除数，会导致其他元素数量级的严重增长和舍入误差的扩散，将严重影响计算结果的精度。实际计算时必须避免这类情况的发生。主元素消去法就可弥补这一缺陷。例4求解如下的方程组解法1使用Gauss消去法求解其精确解为（舍入到4位有效数字）：计算解为：解法2。变换行，避免绝对值小的主元做除数解为：这个解比解
7、法1更加接近真实的解。交换原则：通过方程或变量次序的交换，使在对角线位置上获得绝对值尽可能大的系数作为akk(k)，称这样的akk(k)为主元素，并称使用主元素的消元法为主元素法根据主元素选取范围分为：列主元素法行主元素法(不讲)全主元素法(不讲)§5.3高斯主元素消去法（续）5.3.2列主元素法列主元素法就是在待消元的所在列中选取主元（选取一列中绝对值最大的元素当主元），经方程的行交换，置主元于对角线位置后进行消元的方法。例5.4用列主元素法解下列线性方程组解：选择-20作为该列的主元素，交换方程(1)和(2)得m21=10/-20=-0.
8、5m31=1/-20=-0.05(5)-m21(4),(6)-m31(4)得选6为主元素，交换方程（7）与（8），得：(10)-m32(9)得-2.34168x3=
