用高斯消元法解线性方程组.ppt

用高斯消元法 解线性方程组 北京景山学校 何江舟 GPA排名系统（CTSC2001） 高等院校往往采用GPA来评价学生的学术表现。传统的排名方式是求每一个学生的平均成绩，以平均成绩作为依据进行排名。对于不同的课程，选课学生的平均成绩会受到课程的难易程度等因素的影响，因此这种排名方式不够合理。 为此，我们需要对排名系统进行这样的改进：对第i门课的每一个学生的成绩加上一个特定的修正值di（调整后的成绩不按照百分制），使得经过调整后，该课的平均分等于选该课的所有学生的所有课的平均分。对每一门课都这样调整，使得上述条件对所有课程都满足。 你的任务是根据一个年级学生某学年的成绩，通过上述调整，得出他们的排名。 简要分析 Ai：选修第i门课的学生的集合 Bj：第j个学生选修课程的集合 Gi,j：第j个学生第I门课的成绩 di：第i门课的修正值 对于第p门课，可列出如下关系式： 这是关于di（i=1,2,…,n）的线性方程，我们可以整理出n个这样的方程。 线性方程组的一般形式 a1,1x1+a1,2x2+……+a1,nxn=b1 a2,1x1+a2,2x2+……+a2,nxn=b2 …… an,1x1+an,2x2+……+an,nxn=bn 下面是n元线性方程组的一般形式： 我们可以把它表示为增广矩阵的形式： a1,1 a1,2 …… a1,n b1 a2,1 a2,2 …… a2,n b2 …… an,1 an,2 …… an,n bn 先看一个例子 2 -1 3 1 4 2 5 4 1 2 0 7 2 -1 3 1 4 -1 2 2.5 -1.5 6.5 2 -1 3 1 4 -1 2 -0.875 5.25 ×2 ×0.5 ×2.5 得出： x3=5.25/(-0.875)=-6 x2=(2-(-1)x3)/4=-1 x1=(1-(-1)x2-3x3)/2=9 消元过程 a1,1(1) a1,2 (1) …… a1,n (1) b1 (1) a2,1(1) a2,2 (1) …… a2,n (1) b2 (1) …… an,1(1) an,2 (1) …… an,n (1) bn (1) 注：用上标(k)表示第k次消元前的状态 第1次消元，第1行的乘数： （i=2,3,…,n） a1,1(1) a1,2 (1) …… a1,n (1) b1 (1) a2,2 (2) …… a2,n (2) b2 (2) …… an,2 (2) …… an,n (2) bn (2) 得到新的增广矩阵： ai,j(2)=ai,j(1)-mi,1a1,j(1) bi(2)=bi(1)-mi,1b1(1) （i,j=2,3,…,n） 第k次消元，第k行的乘数： （i=k+1,k+2,…,n） 消元过程 a1,1(1) a1,2 (1) ………… a1,n (1) b1 (1) a2,2 (2) ………… a2,n (2) b2 (2) ………… ak,k(k) …… ak,n (k) bk (k) …… an,k(k) …… an,n (k) bn (k) 第k次消元前的增广矩阵： ai,j(k+1)=ai,j(k) -mi,kak,j(k) bi(k+1)=bi(k) -mi,kbk(k) 增广矩阵的变化： （i,j=k+1,k+2,…,n） 第k步消元的主行 第k步消元的主元素 回代过程 a1,1(1) a1,2 (1) …… a1,n (1) b1 (1) a2,2 (2) …… a2,n (2) b2 (2) …… an,n (n) bn (n) 最后得到的增广矩阵： 最终结果的计算： 为什么要选主元素 前面介绍的消元法都是按照自然顺序，即x1、x2、……、xn的顺序消元的。有： 所以每一次消元的主元素都不能为0。如果按照自然顺序消元的过程中出现的ak,k(k)=0，那么消元无法继续进行下去。或者| ak,k(k) |很小，也会严重影响计算精度。 为什么要选主元素 例如（假设运算过程中使用单精度实数）： 10-10 1 1 1 1 2 10-10 1 1 -1010 -1010 解得：x1=0，x2=1 这个解与第二个方程差异很大。究其原因，因为消元过程中第一个方程所乘的系数过大，使得上式“吃掉”了下式，所以在结果中根本无法体现下式。 但如果调整一下顺序： 1 1 2 10-10 1 1 1 1 2 1 1 解得：x1=1，x2=1，这个解基本符合原方程 所以每次消元的主元素的绝对值应该