三角分解法.ppt

(2)假设对n2阶方阵成立，则对n+1阶方阵充分性还成立吗？ 为此在n+1阶方阵A的第一列中选主元素，如果主元素在第r1行则交换A的第一行与第r行，这相当于用置换阵P P1去左乘A即： 由于A的可逆性故PA的第一行第一列元素不为0，从而可进行一次消元即将beta消成0，为此构造Gauss矩阵： 它还是可逆方阵？ 显然上面分解中的n-1阶方阵B一定是可逆的，由归纳假设，存在n-1阶置换阵PB、单位下三角阵LB、上三角阵UB使 并且LB的元素模不大于1。 其中 易证四项相乘恰好是一个单位下三角阵! 记为P它正好是一个置换阵 记为L 对k=1:n-1 (1)确定主元素： 带主元的LU分解的算法 (2)换行：第k行与第r行交换，并保存p(k)=r; vr=A(k,k:n);A(k,k:n)=A(r,k:n);A(r,k:n)=vr; (3)消元： 带主元的LU分解的算法 查看ＭＡＴＬＡＢ代码 PLU分解的例子 线性方程组的列主元消去法 其中P是置换阵： 定理(Cholesky分解定理)对称正定方阵A必有正对角元的下三角阵L使A=LLT。其中L称为Cholesky因子。 证明：A对称正定则它的全体顺序主子式大于0，故必有LU分解： 对称正定方阵的Cholesky分解法 其中U一波是单位上三角形阵，但A是对称正定的故： 由此得A的分解可表示为： 从而得： 这是单位下三角阵! 这是单位上三角阵! 由A的对称 正定性知D是对角正定阵 修正的Cholesky分解： 构造矩阵： 则得Cholesky分解： 其中L是单位下三角阵，D是对角阵。 Cholesky分解法计算方法 对比两边的元素可建立公式! 首先： 由(k,k)元素得： Cholesky分解法计算公式 再由(i,k)得： 以上公式表达L的计算顺序是： Cholesky分解法计算公式 这是将L的元素按列计算的方法! (5) (4) (3) (2) (1) 修正Cholesky分解法计算公式 其中L是单位下三角阵，D是对角阵。 可以避免开平方! * 如果一个方程组系数矩阵具有三角形形状比如： 三角形方程组求解 则可以用前代法求解下三角形方程组Ly=b： 三角形方程组求解 则可以用回代法求解上三角形方程组Ｕx=y： 如果能将一个系数方阵分解成两个三角形方阵的乘积的形式即： 方阵的三角分解 这种分解的好处是：可将方程组 Ax=b化成两个三角方程组去解： 这节我们将涉及到一种较特殊的下三角阵： 单位下三角阵 它的主对角线元素为1 现在设Gauss顺序消去法已经进行完第k-1次消元得到矩阵k=1,2,……,n-1 Gauss消去法的再分析 将n-k个消元因子： 放在单位阵的第k列构造一个矩阵,称为Gauss矩阵 它的主对角线元素为1，第k列第k+1~n元素为负消元因子，其它元素为0。 比如前面哪个例子中产生的消元因子可以生成如下Gauss矩阵： Gauss矩阵有很好的性质，比如求逆： 它的逆阵很好求!而且逆阵也是Gauss矩阵! 第k步消元实际上是用Lk－1左乘A(k-1)得到 这正好是上三角形阵! 经过一系列的消元得： 比如第一步消元实际上是用L1-1左乘(A,b)得到: 合并上面的一系列消元过程中的矩阵得： 可以证明：它们的乘积正好是一个单位下三角阵! 这时U是上三角阵!称这种分解为系数方阵A的三角分解! 取A(0)，A(n-1)的前面n列得 Doolittle分解 我们课本上用的就是这种分解!当然还有其它分解形式，如Crout分解。 Doolittle分解分解的求法 第一种方法：保存Gauss顺序消去法产生的消元因子． 第二种方法：紧凑格式法 通过比较如下分解两边的位于(i,j)的元素，可得： 由这两个情况可以建立直接求三角分解的紧凑格式法公式! 当列数j不小于行数i时即：j=i,i+1,……,n 当列数j小于行数i时即：j=1,……,i-1 紧凑格式法的计算公式 紧凑格式法的计算图示：5*5的方阵 (8) (7) (6) (5) (4) (3) (2) (1) 这种图示表明(1)上三角阵U按行从上到下求出； (2)单位下三角阵按列从左到右求出。 (3)求U的行与L的列交叉进行。 分析上面的计算公式可得可得这样的计算顺序! 紧凑格式法的计算图示：5*5的方阵 如果仔细分析上面公式的计算顺序，还可得另一种计算顺序!它的计算公式自己建立! (9) (8) (7) (6) (5) (4) (3) (2) (1) 这种图示表明(1)上三角阵U按行从上到下求出； (2)单位下三角阵也按行从左到右求出。 三角分解的例子：要求会求 你可以使用Gauss顺序消去法，也可以用紧凑格式法。 这是用Gauss顺序消去法啊!比较容易掌握!