高斯消元法解线性方程组

.z.高斯消元法解线性方程组在工程技术和工程管理中有许多问题经常可以归结为线性方程组类型的数学模型，这些模型中方程和未知量个数常常有多个，而且方程个数与未知量个数也不一定一样。则这样的线性方程组是否有解呢？如果有解，解是否唯一？假设解不唯一，解的构造如何呢？这就是下面要讨论的问题。一、线性方程组设含有n个未知量、有m个方程式组成的方程组〔3.1〕其中系数，常数都是数，是未知量〔也称为未知数〕。当右端常数项,,…,不全为0时，称方程组〔3.1〕为非齐次线性方程组；当==…==0时，即〔3.2〕称为齐次线性方程组。由n个数,,…,组成的一个有序数组〔,,…,〕，如果将它们依次代入方程组〔3.1〕中的,,…,后，〔3.1〕中的每个方程都变成恒等式，则称这个有序数组〔,,…,〕为方程组〔3.1〕的一个解。显然由=0,=0,…,=0组成的有序数组〔0,0,…,0〕是齐次线性方程组〔3.2〕的一个解，称之为齐次线性方程组〔3.2〕的零解，而当齐次线性方程组的未知量取值不全为零时，称之为非零解。〔利用矩阵来讨论线性方程组的解的情况或求线性方程组的解是很方便的。因此，我们先给出线性方程组的矩阵表示形式。〕非齐次线性方程组〔3.1〕的矩阵表示形式为：A*=B其中A=，*=，B=称A为方程组〔3.1〕的系数矩阵，*为未知矩阵，B为常数矩阵。将系数矩阵A和常数矩阵B放在一起构成的矩阵=称为方程组〔3.1〕的增广矩阵。齐次线性方程组〔3.2〕的矩阵表示形式为：A*=O二、高斯消元法〔下面介绍利用矩阵求解方程组的方法，则矩阵初等行变换会不会改变方程组的解呢？我们先看一个定理。〕定理3.1假设用初等行变换将增广矩阵化为，则A*=B与C*=D是同解方程组。证由定理3.1可知，存在初等矩阵,,…,，使…=记…=P，则P可逆，即存在。设为方程组A*=B的解，即A=B在上式两边左乘P，得PA=PB即C=D说明也是方程组C*=D的解。反之，设为方程组C*=D的解，即C=D在上式两边左乘，得C=D即A=B说明也是方程组A*=B的解。因此，方程组A*=B与C*=D的解一样，即它们是同解方程组。〔证毕〕(由定理3.1可知，求方程组〔3.1〕的解，可以利用初等行变换将其增广矩阵化简。又有第二章定理2.10可知，通过初等行变换可以将化成阶梯形矩阵。因此，我们得到了求解线性方程组〔3.1〕的一般方法：)用初等行变换将方程组〔3.1〕的增广矩阵化成阶梯形矩阵，再写出该阶梯形矩阵所对应的方程组，逐步回代，求出方程组的解。因为它们为同解方程组，所以也就得到了原方程组〔3.1〕的解。这种方法被称为高斯消元法，〔下面举例说明用消元法求一般线性方程组解的方法和步骤。〕例1解线性方程组〔3.3〕解先写出增广矩阵，再用初等行变换将其逐步化成阶梯形矩阵，即=上述四个增广矩阵所表示的四个线性方程组是同解方程组，最后一个增广矩阵表示的线性方程组为将最后一个方程乘，再将项移至等号的右端，得将其代入第二个方程，解得再将代入第一个方程组，解得因此，方程组〔3.3〕的解为〔3.4〕其中可以任意取值。由于未知量的取值是任意实数，故方程组〔3.3〕的解有无穷多个。由此可知，表示式〔3.4〕表示了方程组〔3.3〕的所有解。表示式〔3.4〕中等号右端的未知量称为自由未知量，用自由未知量表示其它未知量的表示式〔3.4〕称为方程组〔3.3〕的一般解，当表示式〔3.4〕中的未知量取定一个值〔如=1〕，得到方程组〔3.3〕的一个解〔如，，，〕，称之为方程组〔3.3〕的特解。注意，自由未知量的选取不是唯一的，如例1也可以将取作自由未知量。如果将表示式〔3.4〕中的自由未知量取一任意常数k，即令=k，则方程组〔3.3〕的一般解为，其中k为任意常数。用矩阵形式表示为=〔3.5〕其中k为任意常数。称表示式〔3.5〕为方程组〔3.3〕的全部解。〔用消元法解线性方程组的过程中，当增广矩阵经过初等行变换化成阶梯形矩阵后，要写出相应的方程组，然后再用回代的方法求出解。如果用矩阵将回代的过程表示出来，我们可以发现，这个过程实际上就是对阶梯形矩阵进一步简化，使其最终化成一个特殊的矩阵，从这个特殊矩阵中，就可以直接解出或“读出〞方程组的解。例如，〕对例1中的阶梯形矩阵进一步化简，上述矩阵对应的方程组为将此方程组中含的项移到等号的右端，就得到原方程组〔3.3〕的一般解，〔3.4〕其中可以任意取值。例2解线性方程组解利用初等行变换，将方程组的增广矩阵化成阶梯阵，再求解。即=一般解为例3解线性方程组解利用初等行变换，将方程组的增广矩阵化成阶梯阵，再求解。即=阶梯形矩阵的第三行“0,0,0,-2〞所表示的方程为：，由该方程可知，无论，，取何值，都不能满足这个方程。因此，原方程组无解。三、线性方程组的解的判定前面介绍了用高斯消元法解线性方程组的方法，通过例题可知，线性方程组的解的情况有三种：无穷多解、唯一解和无解。从求解过程可以看出，方程组〔3.1〕是否有解，关键在于增广矩阵[AB]化成阶梯非零行的行数与系数矩阵A化成阶梯形矩阵后非零行的行数是否相等。因此，线性方程组是否有解，就可以用其系数矩阵和增广矩阵的秩来描述了。定理3.9线性方程组〔3.1〕有解的充分必要是=。证设系数矩阵A的秩为r，即=r。利用初等行变换将增广矩阵[AB]化成阶梯阵：[AB]=[CD]故A*=B与C*=D是同解方程组，因此A*=B有解=0==r即==r。〔证毕〕推论1线性方程组有唯一解的充分必要条件是＝=。推论2线性方程组有无穷多解的充分必要条件是＝。〔将上述结论应用到齐次线性方程组〔3.2〕上，则总有＝。因此齐次线性方程组一定有解。并且有〕例4判别以下方程组是否有解？假设有解，是有唯一解还是有无穷多解？(1)(2)(3)解(1)用初等行变换将增广矩阵化成阶梯阵，即[AB]=因为=4，=3，两者不等，所以方程组无解。(2)用初等行变换将增广矩阵化成阶梯阵，