高斯消元法(完整).docVIP

高斯消元法解线性方程组 在工程技术和工程管理中有许多问题经常可以归结为线性方程组类型的数学模型，这些模型中方程和未知量个数常常有多个，而且方程个数与未知量个数也不一定相同。那么这样的线性方程组是否有解呢？如果有解，解是否唯一？若解不唯一，解的结构如何呢？这就是下面要讨论的问题。 一、线性方程组 设含有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）的矩阵表示形式为： AX B 其中 A ，X ，B 称A为方程组（3.1）的系数矩阵，X为未知矩阵，B为常数矩阵。将系数矩阵A和常数矩阵B放在一起构成的矩阵 称为方程组（3.1）的增广矩阵。 齐次线性方程组（3.2）的矩阵表示形式为：AX O 二、高斯消元法 （下面介绍利用矩阵求解方程组的方法，那么矩阵初等行变换会不会改变方程组的解呢？我们先看一个定理。） 定理3.1 若用初等行变换将增广矩阵化为，则AX B与CX D是同解方程组。 证 由定理3.1可知，存在初等矩阵, , …, ，使 … 记… P，则P可逆，即存在。 设为方程组A X B的解，即 A B 在上式两边左乘P，得 P A PB 即 C D 说明也是方程组C X D的解。反之，设为方程组C X D的解，即 C D 在上式两边左乘，得 C D 即 A B 说明也是方程组AX B的解。 因此，方程组A X B与C X 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 解线性方程组 解 利用初等行变换，将方程组的增广矩阵化成阶梯阵，再求解。即