数值剖析ch05三角分解法.pptVIP

第五章 解线性方程组的直接方法 矩阵三角分解法  直接三角分解法 Matrix factorization method 口直接三角分解的基本思想 如果将线性方程组Ax=b的系数矩阵A分解成两个 三角矩阵L和U,即A=LU(也称为LU分解),则: Ax=b LUx =b Ly=b ●问题: 对于一般非奇异矩阵是否存在LU分解?  LU分解 口Gaus消去过程其实就是一个矩阵的三角分解过程 将 Gauss消去过程中第k-1步消元后的系数矩阵记为 (k=1,,n-1 则A(与A1之间的关系式可以表示为:A+1=LA) 其中:「1 i a/a k+1,k (i=k+1,…,n) 3  于是有:A=Ln-1…L2L1A →A=A=(L2LL)A 容易验证:「1 (k=1,n-1) 记:L=L1L2…Ln,U=A,则 A=LU→LU分解(Do0ittl分解) 其中:L-单位下三角矩阵,U-上三角矩阵  分解存在唯一性 LD分解存在(》高斯消去法不被中断 所有顺序主子式不为零〈a≠0 定理:若A的所有顺序主子式不为零,则A存在 唯一的L分解  PLU分解 列主元 Gauss消去法对应的矩阵分解为PLU分解 定理:若A非奇异,则存在排列矩阵P,使得 PA= LU 其中L为单位下三角矩阵,U为上三角矩阵 6  1