2.3直接三角分解法.pptVIP

华长生制作 矩阵的三角分解 一、LU分解(Doolittle分解) 2.存储问题和全部流程 用LU分解来解线性方程组的全部流程如下: 例1:利用LU分解求解线性方程组 矩阵的三种形式的分解: Doolittle分解: A = LU (单位下三角与上三角) Crout分解: (下三角与单位上三角) LDU分解: A = LDU (单位下三角, 对角及单位上三角) 二、追赶法 定理3:设上述三对角矩阵A满足 4.追赶法的计算流程 总结 事实上，追赶法的求解过程就是将系数矩阵分解两个简单的二对角线矩阵，从而归结为求解两个简单三角形方程组的过程。 例2:用追赶法求解方程组 三、解正定矩阵方程组的平方根法 * * 第二章 解线性方程组的直接法 2.4 直接三角分解法 1. LU分解 用n=3来举例说明。三元方程组Ax=b的增广矩阵为 ,用初等行变换变换为 这相当于将A矩阵左乘初等行变换矩阵 和 　　　　　　　　　　　　　　 使M2M1[A b]=[A(1) b(1)],其中l21=a21/a11, l31=a31/a11 再将 进行初等行变换,得到 这相当于左乘初等行变换矩阵 ,其中 ,则 。记 ,则U是一个上三角矩阵,有 。则 。其中 , , 。 记 L是下三角矩阵。这就是矩阵A的LU分解:A=LU。 3. 分解的充要条件 定理1:A能分解为　　的充分必要条件是A的各阶顺序主子式均不为零。即 , , 。 4、LU分解的紧凑格式 　　我们不必在高斯消去法过程中产生L和U,而是直接用矩阵A来进行LU分解。 下面导出计算公式 因A=LU则 4、LU分解的紧凑格式 　　我们不必在高斯消去法过程中产生L和U,而是直接用矩阵A来进行LU分解。 下面导出计算公式 因A=LU则 利用矩阵乘法,得 综合以上分析,有 因此可以推导出 U的第一行 L的第一列 ------(1) ------(2) U的第r行 L的第r列 ------(3) ------(4) 称上述(1) ~ (4)式所表示的分解过程为Doolittle分解 由此可以得到 和 的计算公式: 。在U的第 行和L的第　　　　 列计算出来后再 先计算U的第一行 计算L的第j列: 在U的第 行和L的 再计算U的第i行: 计算的过程是: U第1行,L第1列,U第2行,L第2列,……顺序计算。 第 列均计算出来后 1.计算流程 U可以存储在A的上三角部分,L可以存储在A的下三角部分,对角元存U的对角元,A的对角元为1而不必存储,如下图: 解:第一步,A的LU分解 第二步,求解 解得: , , , 。 第二步,求解 解得: , , , 。 1. 三对角矩阵:在实际问题中,经常会遇到如下形式的方程组, 其系数矩阵为 没有写的部分均为零(下同),对角元为 后次 对角线上的元素为 ,前次对角线上的元素为 。 我们称矩阵A为三对角矩阵,相应的方程组称三对角方程组。如果对A进行LU分解,则将有如下形式: 并称它们为二对角矩阵。 且 则对矩阵A的LU分解能进行,且分解是唯一的。 3.追赶法的计算公式 利用矩阵乘法可得: 从而可以得到: 解 得: 解 得: , , 第一个循环称之为追的过程,相当于消元过程; 第二个循环称之为赶的过程,相当于回代过程。 追赶法的原理和高斯消去法相同，但考虑到 方程组的特点，计算时会把大量零元素撇开， 从而大大节省计算量。 也称Thomas法 解:追的过程: , ; , , ; , , ; , , 。 赶的过程: , , , 。 定理1. (Cholesky分解) 且该分解式唯一。 这种关于对称正定矩阵的分解称为Cholesky分解 如果方程组的系数矩阵A的对称正定矩阵,可以证明: A可以唯一分解为 ,其中L是下三角矩阵, 是L的转置, 即 A= =L 。 由矩阵乘法可知,在对角元上 ; 在第 列均计算完后得 ;在 时, ; * * * *