第五章直接法-数值积分.pptVIP

第五章. 解线性代数方程组的直接方法 §1.引言 §2. Gauss (高斯)消去法 Gauss 消去法计算过程 §3. 高斯主元素消去法 选主元素的矩阵表示 例如： L相当是由各L-1k (k=1,2,…,n-1)的所有左下元素拼凑 后加上对角元1而得. 当A进行LU分解后,Ax=b就容易解了. 即Ax=b等价于: 这样,Ax=b分解为解两个三角形方程组,三角形方程组是极易求解的.(*)即是高斯消元法的矩阵表述. 定理7 矩阵 An×n,只要A的各阶顺序主子式非零,则A可以分解为一个单位下三角阵L和一个上三角阵U的乘积,即A=LU,且这种分解是唯一的. (P173) L是下三角矩阵,且所有的对角元为1,故称为单位下三角阵.这样,A=LU称为A的LU分解,其中L是单位下三角阵,U是上三角阵. 例1 直接LU分解法 A的LU分解可以用高斯消元法完成,但也可以用矩阵乘法原理推出另一种方法,结果是完全一致的． 由矩阵乘法公式可推出直接LU分解算法,也称杜利特(Doolittle) 紧凑算法,现总结如下: 1.矩阵分解:A=LU，记 2.解 L y = b 3.解 U x = y 杜利特算法可用文字描述如下: 将A的元素划分为形如“┏”的n框, 如A的第一行和第一列元素为第一框,紧靠第一框内部的第二行和第二列元素为第二框,依此类推,ann一个元素为第n框.