第四节 矩阵分解新.pptVIP

消元法的实质是对增广矩阵作初等行交换,而初等变换是可以用矩阵运算来描述的.消元过程相当于下述矩阵乘法运算。;令;消元过程相当于分解 A为单位下三角阵L与上三角阵U的乘积，解方程组Ly=b，回代过程就是解方程组Ux=y。其中 ;定理1 若矩阵A非奇异，则A能分解为LU 的充分必要条件是A的顺序主子式不为0，即 ;二 LU 分解的计算公式;比较式 A=LU 两端的元素, 按下图所示顺序逐框进行,先求 ukj 后求 lik . 由第一框可得;得;有了矩阵 A 的LU分解计算公式 , 解线性方程组 Ax=b 就转化为依次解下三角方程组 Ly=b 与上三角方程组 Ux=y 其计算公式如下： ;解 用紧凑格式;所以;例 3 将矩阵A 进行LU 分解 ;; 当矩阵 A 对称且 其顺序主子式均不为0时由定理2可知其LU分解必存在。且;于是 A=LU=LDM. 又据A的对称性有;同 A 的 LU 分解过程, 可得LDLT分解的计算公式为;此时, 解方程组 Ax=b 等价于???LDLTx=b，而后者可转化为依次解两个三角形方程组 Ly=b , LTx=D-1y，其计算公式如下 ;四 平方根法（LLT 分解法）;定理4 设矩阵 A 对称正定,则必存在下三角阵 L, 使 A=