线性方程组的直接解法.ppt

Gauss-Jordan消元法求矩阵的逆Gauss-Jordan消元法求矩阵的逆（续1）Gauss-Jordan消元法求矩阵的逆（续2）Gauss-Jordan消元法求矩阵的逆（续3）Gauss-Jordan消元法求逆阵(续4)Gauss-Jordan消元法求逆阵(续5)Gauss-Jordan消元法求逆阵(续6)Gauss-Jordan消元法求逆阵(续7)Gauss-Jordan消元法求逆阵(续8)Gauss-Jordan法求逆阵的具体步骤Gauss-Jordan法求逆阵的具体步骤(续)Gauss-Jordan法求逆阵举例1……设经过k–1步后得到：其中：完成n步消元后，A?1放在原A的位置。按上述紧缩存贮原则，可节省存贮单元，同时还使得整个计算更简单了。可总结求逆步骤如下：上述1，2是求第k列元素，构成Mk（即求主列）（计算其他元素，但少k列，k行）用上述Gauss-Jordan法求逆阵，计算量约为n3，是Gauss消元法的3倍，为保证方法稳定性，还可选列主元，若仍按上述紧缩存贮原则，则最后需按行交换的相反次序作列交换才能得到A?1。例解：按紧缩存贮方式，逐次计算结果与存贮如下：第一步：k=1，

在第一列中选

主元，交换1，

2行，得：第三步：k=3以a33=1/6为主元，消元后得：交换2、3列最后：按行交换的相反次序进行列交换：先交换2，3列，再

交换1，2列得A?1。交换1，2列第二步：k=2在第二列对角元下选主元，交换2，3行由1，2

先计算第2列，由3计算其他元素（除2列2行外）而由4计算剩

下的第2行的元素（这里k=2的第2列第行称为主列，主行）**线性方程组的直接解法Gauss消去法直接三角分解方法方程组的性态与误差估计*在自然科学与工程技术中，很多问题的解决常常归结为解线性方程组的问题：很多数值计算方法到最后也涉及到线性方程组的求解问题：如电学中的网络问题，机械和建筑结构的设计和计算等。如求样条插值的M和m的关系式，解曲线拟合的法方程，求矩阵特征值的反幂法等问题。电子计算机线性方程组{稠密和稀疏（按系数矩阵含零元多少分）高阶和低阶（按阶数的高低分）对称正定、对角占优等（按系数矩阵的形状性质分）基本解法{直接法（通过有限步计算得到精确解，适用于低阶、大型带型阵）迭代法（通过逐次迭代逼近得到近似解，适用于大型稀疏、非带型阵）线性方程组及方法分类*求解线性方程组：对此方程组进行求解有两种方法：采用Cramer法则、消元法。*Cramer(克莱姆)法则对于20阶的线性方程组,若用Cramer法则求解,其乘、除运算次数为9.7*1020,用一亿次/秒的计算机,要30.8万年！若用高斯消去法进行数值求解，乘、除运算只需约3060次。定理：如果线性方程组则方程组有唯一解：其中Ak是将A的第k列元素依次换成常数项b1,…,bn得到的行列式。的系数行列式非零，即计算量大*Gauss消去法一、高斯顺序消去法思路首先将A化为上三角阵/*upper-triangularmatrix*/，再回代求解/*backwardsubstitution*/。=是一种古老的求解线性方程组的方法,按自然顺序进行消元的方法.*例1解方程组解step1消元*Step2对上三角形方程进行回代求解,得下面我们来一般性地讨论求解n阶线性方程组的高斯顺序消去法.*消元Step1：设，计算因子将增广矩阵/*augmentedmatrix*/第i行?li1?第1行，得到与(1)式等价的方程组*Step2：一般第k次消元(1≤k≤n-1)**Step3：继续上述过程,且设aii（i-1）≠0(i=1,2,…,n-1),直到完成第n-1次消元,最后得到与A(0)x=b(0)等价的三角形方程组A(n-1)x=b(n-1).将(1)式化为(2)式的过程称为消元过程.forforforGauss消去法的消元过程算法Gauss消去法工作量为*回代定理若A的所有顺序主子式/*determinantofleadingprincipalsubmatrices*/均不为0，则高斯消元无需换行即可进行到底，得到唯一解。注：事实上，