第7章 解方程组.doc

第 七 章 解 方 程 组 的 数 值 方 法 §1 引 言 在自然科学和工程技术中有很多问题的解决归结为求解线性方程组或者非线性方程组的数学问题。例如，电学中网络问题，用最小二乘法求实验数据的曲线拟合问题，三次样条的插值问题，用有限元素法计算结构力学中一些问题或用差分法解椭圆型偏微分方程边值问题等等。 在工程实际问题中产生的线性方程组，其系数矩阵大致有两种，一种是低阶稠密矩阵（这种矩阵的全部元素都存贮在计算机的存贮器中）；另一类是大型稀疏矩阵（此类矩阵阶数高，但零元素较多）。 本章§2一§6介绍求解线性方程组的直接法。目前这种方法是计算机上解低阶稠密矩阵的有效方法。§8将介绍计算机上解大型稀疏矩阵方程组的迭代法。 §2 高 斯 消 去 法 高斯消去法是一个古老的求解线性方法组的方法，但由它改进得到的选主元的高斯消去法则是目前计算机上常用的解低阶稠密矩阵方程组的有效方法。 例1 用高斯消去法解方程组 （2.1） 解 第1步：将（2.1）的第一个方程分别乘上（）和（），并分别加到第二、第三个方程上去，则得到与原方程组等价的方程组 （2.2） 第 2步 ： 将（2.2）的第二个方程乘上2加到第三个方程，消去其中的未知数，又得到与原方程组等价的三角形方程组 （2.3） 最后由上述方程组，用回代的方法，即可求得原方程组的解为： 若用矩阵来描述消去法的约化过程，即为 这种求解过程，称为具有回代的高斯消去法。 从这个例子看出，用高斯消去法解方程组的基本思想是用矩阵行的初等变换将系数矩阵A约化为具有简单形式的矩阵（上三角形矩阵，单位矩阵等），而三角形方程组是很容易解的。 下面讨论求解一般线性方程组的高斯消去法。设有n个未知数的线性方程组 （2.4） 引进记号 ，， （2.4）可用矩阵形式表示 （2.5） 且为了下面讨论方便， 记，。假设A为非奇异矩阵（即设） 。 第1步：设 计算乘数 ， 用乘上（2.4）第一个方程，加到第个方程上去即施行行的初等变换，，消去第2个方程 ～ 第 n个方程的未知数，得到（2.4）的等价方程组 （2.6） 记为 其中（2.6）式中右上角标为（2）的元素为这一步需要计算的元素，计算公式为 第 k步 ：继续上述消去过程，设第1步 ～ 第k— l步计算已经完成，得到与原方程组等价的方程组 （2.7） 记为 现进行第步消元计算，设，计算乘数 ， 用乘（2.7）的第个方程，加到第个方程上去，消去第个方程的未知数，得到与原方程组等价的方程组（略），简记为 其中 元素计算公式为： （2.8） 与前行元素相同，与前个元素相同。 最后，重复上述约化过程，即且设，共完成步消元计算，得到与原方程组（2.4）等价的三角形方程组 （2.9） 用回代的方法，即可求得（2.9）的解，计算公式为 ， （2.10） 元素称为约化的主元素。将（2.4）约化为（2.9）的过程称为消元过程，（2.9）求解过程（2.10）称为回代过程，由消元过程和回代过程求解线性方程组的方法称为高斯消去法。 定理1（高斯消去法）设，其中。如果约化的主元素，则可通过高斯消去法（不进行交换两行的初等交换）将方程组约化为三角形矩阵方程组（2.9）且消元和求解公式为 消元计算 ， 回代计算 ， 如果A为非奇异矩阵时，但可能有某在第列存在有元素，于是可能通过交换的第 k行和第行元素将调到位置，然后再进行消元计算。于是，在A为非奇异矩阵时，只要引进行交换，则高斯消去法可将约化为三角形方程组（2.9），且通过回代即可求得方程组的解。 高斯消去