关于线性方程组的解法—高斯消去法的讨论材料学院1012208026张超高.DOC

关于线性方程组的解法—高斯消去法的讨论 材料学院 1012208026 张超 高斯消去法是求解线性代数方程组的直接法，是目前计算机上解低阶稠密矩阵方程组的常用有效方法。它是建立在逐次消去未知元的基础上，把待求解方程组用逐次消去未知元的方法化为与其等价的（或者说具有相同解的）上三角形方程组。消去法完全可以通过增广矩阵的行的初等变换来实现，具体的计算方法常分为顺序高斯消去法、列主元素高斯消去法和全主元素高斯消去法。下面具体做以介绍。 我们考察n阶线性代数方程组， 一、首先讨论，顺序高斯消去法： 第一，消元过程： 记系数矩阵，右端向量 (1)设，记，(i=2,3,….,n)，然后将增广矩阵中第i行减去第1行乘以，(i=2,3,…,n)，完成第一次消元，得： 其中，，(i,j=2,3,….,n)，，(i=2,3,….,n)， 这样，就得到了与原方程组同解的方程组。 (2)设，记，(i=3,4,…,n)，然后将矩阵中第i行减去第2行乘以(i=3,4,…,n)，完成第二次消元，得到了与同解的，(i,j=3,4,…,n)，(i=3,4,…,n)。继续重复上述过程，按此推广至第n-1次消元，这样就得到了与原方程组同解的上三角形方程组，即： 第二，回带过程： 由上述三角方程组自下而上按未知元xi的足标逆序回代得到原方程组的解。称(k=1,2,…,n)为各布消元的主元素，称为各步消元的乘数。 二、再次讨论，列主元素高斯消去法 对于方程组，仍的顺序依次进行高斯消元，只是在每一步消元之前增加一步按列选主元的工作，具体叙述如下： 第一步：在消元前先在系数矩阵的第一列元素中选取绝对值最大的元素为主元素 (不能为0，否则第一列都是0，即系数矩阵为奇异的，与假设是非奇异矛盾)。 若最大元素不在(1,1)位置上，则可通过交换第一行与第行，把换到(1,1)位置上作为主元素，记行交换后的增广矩阵为，然后再进行之前顺序高斯消去法的第一次消元，得矩阵。 第二步：在矩阵的第二列元素中选取绝对值最大的元素为第二步主元素，同理该主元素也不为0，若其不在(2,2)位置，则将增广的第二行与第行交换，把换到(2,2)位置，做为主元素，再进行顺序高斯消去法的第二次消元，得到矩阵。 如此进行下去，在每步消元前都要按列选取主元素，直到完成n-1次消元为止，这时原方程组已经被化成等价的上三角方程组，然后用顺序高斯消去法的回代过程，求得方程组的全部解。上述选主元素的方法，由于每一步是按照列选择主元素，因此称为列主元素高斯消去法。 三、顺序高斯消去法以及列主元素消去法的区别（以及Matlab编程举例说明） 顺序高斯消去法是按照方程组的未知元给定的自然顺序逐个消元的，需要满足的前提条件是各步消元的主元素全不能为0，否则即使的行列式值不等于0，消元过程无法进行。如果主元素绝对值很小没用它做为除数在回代过程求解，则其倒数绝对值很大，由于计算机的精度是有限的，这样会引入很大误差，使方程组的解精度受到严重影响。这样的解决方法就是要选择主元必须要有技巧性，列主元素高斯消去法就是一种带有选主元技巧的消去方法。 在编程实现过程中，列主元素高斯消去法与顺序高斯消去法基本相同，主要分为消元过程以及回代求解过程，前者比后者只是增加了列主元选取的程序代码内容，下面举两个例子体现上述内容。 例1 解线性代数方程组： 首先编写顺序高斯消去法的M文件，名为shunxuxiaoqu.m： clear format rat A=input(输入增广矩阵A=) [m,n]=size(A); for i=1:(m-1) numb=int2str(i); disp([第,numb,次消元后的增广矩阵]) for j=(i+1):m A(j,:)=A(j,:)-A(i,:)*A(j,i)/A(i,i); end A end %回代过程 disp(回代求解) x(m)=A(m,n)/A(m,m); for i=(m-1):-1:1 x(i)=(A(i,n)-A(i,i+1:m)*x(i+1:m))/A(i,i); end x 运行该文件，在命令窗口出现： “输入增广矩阵A=” 按要求输入增广矩阵如下以及回车后运行结果为： 输入增广矩阵A=[3,1,-1,2,6;-5,1,3,-4,-12;2,0,1,-1,1;1,-5,3,-3,3] A = 3 1 -1 2 6 -5 1 3 -4