第三章线性方程组直接解法.pptVIP

§3矩阵分解法如果用矩阵形式表示，Gauss消元法的消元过程是对方程组（2-1）的增广矩阵（A、b）进行一系列的初等行变换，将系数矩阵A化成上三角矩阵的过程，也等价于用一串初等变换阵去左乘增广矩阵，因此，消元过程可以通过矩阵运算来实现。紧接下屏：3.1Gauss消元法的矩阵形式事实上，Gauss消元法的第一次消元相当于用初等矩阵：第29页，共53页，星期日，2025年，2月5日第二次消元相当于用初等矩阵：第k次消元相当于用初等矩阵：Gauss消元法的矩阵形式第30页，共53页，星期日，2025年，2月5日经过n?1步消元后得到：因为Lk(k=1,2,…,n?1)均为非奇异阵，故它们的逆矩阵存在。容易求出：这说明：在的条件下，消元过程实际上是把系数矩阵A分解成单位下三角阵与上三角矩阵的乘积的过程。Gauss消元法的矩阵形式（续）第31页，共53页，星期日，2025年，2月5日杜利特尔（Doolittle）分解——LU分解事实上，只要A满足一定条件，由上述结论，它一定可以分解成两个三角形矩阵的乘积，即：A=LU。上述分解称为杜利特尔（Doolittle）分解，也称为LU分解，当系数矩阵完成三角分解后，对于求解方程组：Ax=b。消元过程相当于分解A=LU及求解三角形方程组Ly=b，回代过程则是求解另一个三角形方程组Ux=y，因此，解线性方程组问题可转化为矩阵的三角分解问题。其中：L为单位下三角矩阵，U为上三角矩阵：第32页，共53页，星期日，2025年，2月5日3.2矩阵的三角分解正如Gauss消元法要在一定条件下才能进行到底一样，矩阵A也必须满足一定条件才能进行三角分解。设A为n阶方阵，若A的顺序主子式Ai(i=1,2,…,n)均不为零,则矩阵A存在唯一的Doolittle分解。定理2.1下面讨论如何对A进行LU分解：（1行1列）由于两个矩阵相等就是它们的对应元素都相等，因此通过比较A与LU的对应元素，即可得到直接计算L、U的元素的公式：A=LU即：紧接下屏：第33页，共53页，星期日，2025年，2月5日对A进行LU分解ALU由矩阵乘法规则及比较（2-11）两端的元素，得：即可由（2-12）求出U的第一行，L的第一列。下面讨论一般情况，即：U的第i行，L的第j列：第34页，共53页，星期日，2025年，2月5日一般情况下，可由：对A进行LU分解（r行r列）计算过程应按U第1行、L第1列（第1框），U第2行、L第2列（第2框），

……的顺序具体分解步骤见下屏：得到计算uij和lij的公式：第35页，共53页，星期日，2025年，2月5日对A进行LU分解的具体步骤1.计算U的第1行，L的第1列，亦称为计算

第1框；2.计算U的第r行，L的第r列（r=2,…,n），

即第r框：第36页，共53页，星期日，2025年，2月5日矩阵A的LU分解举例解：按分解公式（2-13），一框一框分解，每框计算时先行后列：所以：例5第37页，共53页，星期日，2025年，2月5日3.3直接三角分解法若线性方程组Ax=b的系数矩阵A完成三角分解，A=LU，那么解方程组Ax=b等价于求解两个三角形方程组Ly=b，Ux=y，即由：再由：可解得：第38页，共53页，星期日，2025年，2月5日第三章线性方程组直接解法返回前进第三章线性方程组直接解法第1页，共53页，星期日，2025年，2月5日第三章目录§1.Gauus消元法§2.主元素法2.1引入主元素法的必要性2.2列主元素法2.3全主元素法2.4解三对角方程组的追赶法§3.矩阵分解法3.1Gauss消去法的矩阵形式3.2矩阵的三角分解3.3直接三角分解法§4.平方根法与改进的平方根法§5.矩阵求逆§6.方程组的性态和条件数第2页，共53页，星期日，2025年，2月5日设n阶线性方程组：其矩阵形式为：Ax=b(2-2)其中：在科学研究和工程技术中所提出的计算问题中，线性方程组的求解问题是基本的，常见的，很多问题如插值函数，最小