数值线性代数课程设计高斯消去法.docxVIP

数值线性代数课程设计线性方程组的直接解法数理学院 0940501120 沈骁摘要：如何利用电子计算机来快速、有效的求解线性方程组的问题是数值线性代数的核心问题。本文将主要介绍解线性方程组的基本的直接法——高斯消去法，平方根法，并用实例来验证此方法的有效性。关键字：高斯消去法，顺序消去法，选主元消去法，平方根法，消元过程，回代过程,主元数和乘数引言：因为各种各样的科学与工程问题往往最终都要归结为一个线性方程组的求解问题。本文在比较着几个方法的基础上，通过一道实例来得到最方便最有效的方法。基本原理：工程计算和科学研究中的许多问题，最终归结为线性代数方程组的求解。求解的方法也有很多，如高斯消去法（顺序消去法，选主元消去法），平方根法。高斯消去法是目前求解中小规模线性方程组最常用的方法；平方根法是求解对称正定线性方程组最常用的方法之一。为了更快速、更方便的求解线性方程组，下面我们比较一下这几种方法哪种更好。高斯（Causs）消去法就是逐步消去变元的系数，将原方程组化为系数矩阵为三角形的等价方程组，然后求解系数矩阵为三角形的方程组而得出原方程组解的方法。把逐步消元去变元的系数，将方程组化为以系数矩阵为三角形的等价方程组的过程称为小院过程；把求系数矩阵为三角形的方程组解的过程称为回代过程。最初求解方程组的高斯消去法也称为顺序消去法，它由消元过程和回代过程组成。顺序消去法消元过程考虑一般方程组，为了推导过程方便，记系数矩阵的元素为，右端向量的元素记为，于是方程组 （1.1）成为假设，将第1个方程乘以加到第个方程,得到第1个导出方程组其中：，，。由于因子不止一次地用到，常记为。再假设，由第1个导出方程组的第2个方程乘以加到第个方程，得到第2个导出方程组类似地记：，则第2个导出方程组的元素，，。重复上述过程次，得到个导出方程组（1.2）其中第个导出方程组的元素的递推关系是， （1.3） 。由于上述消元过程只是将原方程组的系数矩阵和右端进行初等变换，因此，第个导出方程组与原方程组等价，即通过消元过程将方程组（1.1）化成了等价的上三角方程组（1.2）。由第个导出方程组的计算公式（1.3），容易得到消元过程所需要的乘法和加法次数为除法次数为2.回代过程 回代过程就是求等价三角形方程组（1.2）的解。只要，就可以最后一个方程得到的值，再从第个方程得的值。在一般情形，可求得的回代递推公式（1.4） 由公式（1.4）可知，回代过程需要次除法，其乘法和加法的次数同为。所以回代过程需次加法，次乘、除法。3.顺序消去法的运算次数与计算步骤由消元过程和回代过程的运算次数可知，顺序消去法的加法次数为，乘除法次数为。消元过程在编程上机运算时，需采用三重循环，即对于计算；对于计算。回代过程只需要二重循环，即计算，对于对于计算4.主元数和乘数由顺序消去法的推导过程可知，无论是消元过程还是回代过程的都不需要对未知元作真正的运算，而仅需要对方程组 的系数矩阵的元素和右端作运算。因此，在实际运算中，总是将方程组的系数矩阵和右端合在一起，记成增广矩阵。由消元过程（1.3）可以看到，元素为“主元素“，另外，在消元过程中不止一次用到数，这个数称为消元过程的乘数。二．平方根法平方根法又叫Cholesky分解法，是求解对称线性方程组最常用的方法之一。对于一般方阵，为了消除分解的局限性和误差的过分积累，而采用了选主元的方法。但对于对称正定矩阵而言，选主元却是完全不必要的。设是对称正定的，即满足而且0对一切的非零向量成立.此时,由定理容易推出Cholesky分解定理：若对称正定，则存在一个对角元均为正数的下三角阵,求得.上式中的称作的Cholesky因子。若线性方程组的系数矩阵是对称正定的，则我们自然可按如下的步骤求其解：（1）求的Cholesky分解：;(2)求解得；（3）求解得。实验：1.用高斯（Causs）消去法求解下面的84阶方程组。解：2.用你编写的程序求解对称求解对称正定方程组，其中随机的选取，系数矩阵为100阶。解:改进的平方根方法：结论：通过用两种方法对题目的求解，我们得到高斯消去法适合用于中小规模线性方程的求解，它一般用于系数矩阵稠密而又没有任何特殊结构的线性方程组，由它改进后得到的选主元消去法是目前计算机上常用的有效方法；而对于一般方阵，为了消除分解的局限性和误差的过分积累，而采用了选主元的方法。但对于对称正定矩阵而言，选主元却是完全不必要的。通过上面两道题目的解答让我们对求解线性方程组的求法也有了具体的深刻的了解，更好的认识到了高斯（Causs）消元法和平方根法所适合的范围，也知道了这两种方法各有各的优点。因此在我们做题时应该综合考虑各种因素来确定我们所要用的方法。这样才能使题目得到最方便、最有效的解答。参考文献：徐树方 高