高斯列主元消去法解线性方程组教案.pdf

线性方程组的数值解法 ——高斯列主元消去法解线性方程组的MATLAB 实现 班级：MATH 20XX 学号：xxxxxxxxxx 姓名：LI Wen 高斯列主元消去法解线性方程组的MATLAB 实现 高斯列主元消去法解线性方程组的MATLAB 实现 摘 要 在自然科学和工程技术中许多问题的解决常常归结为解线性代数方程组，并且在实 际计算中很多线性方程组是无法求其精确解的，所以只能求他的数值解。而解线性方程 组的数值解法一般有两类：直接法和迭代法。其中直接法就是经过有限的算术运算，求 得线性方程解的方法。但在实际计算中由于舍入误差的存在和影响，这种方法也只能求 的线性方程组的近似解。本论文将阐述这类方法中的最基本的高斯消去法及其变形而来 的列主元消去法，并运用MATLAB 程序来实现它。由于高斯消去法在计算的过程中，会 出现某个主元等于零的情况，从而无法计算；或主元绝对值很小的情况，也使得结果与 真实结果相差甚远。为了避免高斯消去法消元过程中出现的上述两个问题，一般在实际 计算中采用所谓的选择主元法。目前计算机上常用的是按列选主元的方法。因此本文中 做的也是列选主元高斯消去法。下面本文将详细介绍高斯消去法和列主元消去法的计算 过程，并给出列主元消去法的程序算法。 关键词：线性方程组，消去法，列主元，MATLAB，数值计算 -1- By Lyvnee 高斯列主元消去法解线性方程组的MATLAB 实现 引言 在自然科学和工程技术中许多问题的解决常常归结为解线性代数方程组，例如电学 中的网络问题，船体数学放样中建立三次样条函数问题，用最小二乘法求实验数据的曲 线拟合问题，解非线性方程组问题，用差分法或者有限元法解常微分方程，偏微分方程 边值问题等都导致求解线性方程组，而且后面几种情况常常归结为求解大型线性方程组。 关于解线性方程组的数值解法一般有两类：直接法和迭代法。直接法就是经过有限 的算术运算，可求的线性方程组精确解的方法（若计算过程中没有舍入误差）。但在实 际计算中由于舍入误差的存在和影响，这种方法也只能求的线性方程组的近似解。本论 文将阐述这类方法中的最基本的高斯消去法及其变形而来的列主元消去法。 一、问题的提出 我们都知道，高斯列主元素消去法是计算机上常用来求解线性方程组的一种直接的 方法。就是在不考虑舍入误差的情况下，经过有限步的四则运算可以得到线性方程组的 准确解的一类方法。实际运算的时候因为只能有限小数去计算，因此只能得到近似值。 在实际运算的时候，我们很多时候也常用高斯消去法。但是高斯消去法在计算机中运算 的时候常会碰到两个问题。 1）一旦遇到某个主元等于0 ，消元过程便无法进行下去。 2 ）在长期使用中还发现，即使消元过程能进行下去，但是当某个主元的绝对值很 小时，求解出的结果与真实结果相差甚远。 为了避免高斯消去法消元过程中出现的上述两个问题，一般采用所谓的选择主元法。 其中又可以分为列选主元和全面选主元两种方法。目前计算机上常用的按列选主元的方 法。因此我在这里做的也是按列选主元高斯消去法。 二、高斯消去法的基本思想 大家知道，如果一个线性方程组的系数矩阵是上三角矩阵时，即这种方程组我们称 -2- By Lyvnee 高斯列主元消去法解线性方程组的MATLAB 实现 之为上三角方程组，它是很容易求解的。我们只要把方程组的最下面的一个方程求解出 来，在把求