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

*第五章线性方程组的直接解法*比较对应元素：第62页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*当i=j时当ji时固定j按列求第63页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*于是，根据计算公式可以对对称正定矩阵进行平方根分解l11…l21l22l31l32l33………ln1ln2ln3…lnn第64页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*按列计算元素平方根法不必选主元，精度高，稳定性好，是最有效算法之一。第65页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*由于这两种方法的精度差不多，且全主元Gauss消去法程序设计复杂占用机器时间较多，实际应用中一般采用列主元Gauss消去法，它既简单又能保证计算精度。有时候在消元过程中可以在系数矩阵所有元素中选择绝对值最大的元素作为主元素（换行换列），这样的Gauss消去法和叫做全主元Gauss消去法。无回代的列主元消去法：高斯—约当消去法（初等变换法求逆矩阵的规范算法），计算量大约是n3/2次乘除法。第30页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*§2直接三角分解方法一、Gauss消去法的矩阵运算从§1中讨论可知，Gauss消去法的消元过程是将增广矩阵[A,b]＝[A（1）,b（1）]逐步约化为矩阵[A（n）,b（n）]。现在说明，在消元过程中，系数矩阵A＝A（1）是如何经矩阵运算约化为上三角矩阵A（n）,即用矩阵运算的观点来看，消元的每一步计算等价于用一个单位下三角矩阵左乘前一步约化得到的矩阵。第31页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*若，令，i＝2，3，…，n，得到下三角矩阵施行第一步消元，我们得到第32页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*若，令，i＝2，3，…，n，则有施行第二步消元，我们得到初等行变换第33页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*如此下去，直至施行第n－1步消元，得到初等行变换第34页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*由此可见，在顺序Gauss消去法的过程中，系数矩阵A＝A（1）经过一系列单位下三角矩阵的左乘运算约化为上三角矩阵A（n），即这时由得令第35页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*容易验证则从顺序Gauss消去法的矩阵运算表示式可知，系数矩阵A可分解为一个单位下三角矩阵L和一个上三角矩阵U的乘积，即其中第36页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*第一个方程组的系数矩阵为下三角矩阵，第二个方程组的系数矩阵为上三角矩阵，两个方程组都非常容易求解，具体求解结果如下：我们将A=LU称为矩阵A的三角分解,这时线性方程组为：令则有高斯消元法的实质是把A分解为两个三角矩阵乘积第37页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*对于由解得第38页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*对于由求得第39页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*可以看出对于方程组：只要对系数矩阵作了三角分解：由这个简单的计算过程可知，系数矩阵的三角分解很关键，如何进行三角分解更容易？下面介绍几种方法。通过如下两组公式很容易求解：第40页，共92页，星期日，2025年，2月5日*第五章线性方程组的直接解法*二、Doolittle分解法前已述及，若在顺序Gauss消去法的过程中，每步消元的主元素akk(k)≠0，则矩阵A可分解为A＝LU，L为单位下三角矩阵，U为上三角矩阵，此分解称为A的Doolittle（杜利特尔）分解。已经证明主元akk(k)≠0的充要条件是A的各阶顺序主子式不为零，于是有如下定