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

第三 章 线性方程组 直接解法 第三章目录 线性方程组的数值解法 §1 Gauss消元法 例1（续） Gauss消元法的基本步骤3（4阶） Gauss消元法的消元过程1、2（n阶） Gauss消元法的回代过程（n阶） Gauss消元法的计算量 Gauss法与Cramer法则的计算量比较 §2 主元素法 例2（续1） 例2（续2） 例2两种解法的误差分析 2.2 列主元素法 列主元素法 2.3 全主元素法 2.4 解三对角方程组的追赶法 追赶法举例 §3 矩阵分解法 Gauss消元法的矩阵形式 Gauss消元法的矩阵形式（续） 杜利特尔（Doolittle）分解——LU分解 3.2 矩阵的三角分解 对A进行LU分解 对A进行LU分解 （r 行r 列） 对A进行LU分解的具体步骤 矩阵A的LU分解举例 3.3 直接三角分解法 直接三角分解法（续） 解线性方程组举例 三角分解法的几点说明 三角分解法的几点说明（续） 解线性方程组举例 §4 平方根法与改进的平方根法 4.1 平方根法 定理 2.2（续） Choleskg分解1 Choleskg分解2 平方根法举例 4.2 改进的平方根法 改进的平方根法说明 改进的平方根法举例 经过n?1步消元后得到： 因为Lk (k = 1,2,…, n?1) 均为非奇异阵，故它们 的逆矩阵存在。容易求出： 这说明：在 的条件下，消元过程实际上是把系数矩阵A分解成单位下三角阵与上三角矩阵的乘积的过程。 事实上，只要A满足一定条件，由上述结论，它一定可以分解成两个三角形矩阵的乘积，即：A=LU 。 上述分解称为杜利特尔（Doolittle）分解，也称为LU分解，当系数矩阵完成三角分解后，对于求解方程组：Ax = b 。 消元过程相当于分解A = LU及求解三角形方程组Ly = b，回代过程则是求解另一个三角形方程组Ux = y，因此，解线性方程组问题可转化为矩阵的三角分解问题。 其中：L为单位下三角矩阵， U为上三角矩阵： 正如Gauss消元法要在一定条件下才能进行到底一样，矩阵A也必须满足一定条件才能进行三角分解。 设A为n阶方阵，若A的顺序主子式Ai (i = 1,2,…,n )均不为零,则矩阵A存在唯一的Doolittle分解。 定理2.1 下面讨论如何对A进行LU分解： （1行1列） 由于两个矩阵相等就是它们的对应元素都相等，因此通过比较A与LU的对应元素，即可得到直接计算L、U的元素的公式：A=LU 即： 紧接下屏： A L U 由矩阵乘法规则及比较（2-11）两端的元素，得： 即可由（2-12）求出U的第一行，L的第一列。 下面讨论一般情况，即：U的第i 行，L的第j列： 一般情况下，可由： 计算过程应按U第1行、L第1列（第1框），U第2行、L第2列（第2框），……的顺序 具体分解步骤见下屏： 得到计算uij 和 lij 的公式： 1. 计算U的第1行，L的第1列，亦称为计算 第1框； 2. 计算U的第r 行，L的第r 列（r =2,…,n）， 即第r 框 ： 解：按分解公式（2-13）， 一框一框分解，每框计算时先行后列 ： 所以： 例5 若线性方程组Ax = b的系数矩阵 A完成三角分解，A = LU，那么解方程组Ax = b等价于求解两个三角形方程组Ly = b，Ux = y，即由： 再由： 可解得： 容易看出，式（2-14）与式（2-13）的运算规律相同， 或者由： 例6 解 ：按表2-3计算： 1、用三角分解法求线性方程的乘除法运算量也是n3/3 数量级。由于在求出uij, lij和yi后，aij和bi无需保留， 故上机计算时，可把L，U和y存在A，b所占的单元， 回代时x取代y，整个计算过程中不需要增加新的存 贮单元。 3、完成A=LU分解后可以较容易地求出行列式|A|的值： 2、从三角分解法的推导及例中可以看出，系数矩阵的 三 角分解与右端项无关。因 而在计算多个系数矩阵 为A而右端不同的线性方程组系时，用三角分解法更 为 简便（如可用于求逆矩阵）。 6、分解法的优点除上述2、3外，还有： a. 可求解A2z = b,因为算A2计算量大，可用 b. 可根据A的形状设计算法，当A为大型稀疏，且非零 元素有规律如带状，三对角等，作分解时能充分利用A的 特点，L，U能