3.2 LU分解.pptVIP

计算方法 矩阵分解及其在解方程组中的应用 南京信息工程大学数学系 3.2.1 矩阵分解的基本概念 定义1 给定n阶矩阵A,如果存在一个下三角形矩阵L与上三角形U使得A=LU,就说A有三角分解。 这里，L和U的形式如下： 由于(3)、(4)是两个特殊的方程组，(3)是下三角方程组，(4)是上三角方程组，其算法非常简单，对(3)其递推计算式为： 对(4)其递推计算式为: 3.2.3 直接分解的本质或理论来源 我们知道，对n阶线性方程组，高斯消元法要做n-1次消元，观察下面的式子： 分析第1次消元，目的要将第1列除第1行元素之外的所有元素消去，其步骤是： 1、计算各行的消元因子： 2、计算第i行： 从矩阵变换的角度来看，这次消元的实质是对矩阵A左乘矩阵L1,即 实际上，做第k次消元，类似于第1次，即要计算： 1、计算各行的消元因子： 2、计算第i行： 从矩阵变换的角度来看，这次消元相当于： 另一种变化是减少前代过程，减少前代过程是基于考察高斯消元法的结果。事实上，如果在对A做直接分解的同时，对b相量也作相应的计算，于是得到Doolittle分解的一个改进递推公式如下： [二.紧凑LU分解算法] S1:初始化数组:a[n][n+1] /*为方程组的增广矩阵 S2:循环k[分别计算U及L的第k行和第k列]：k=0到n-1 S21:计算U的第k行 S211[循环计算U的第j列]j=k to n a[k][j]= a[k][j]- S22:计算L的第k列 S221[循环计算L的第i列]i=k+1 to n-1 a[i][k]= (a[i][k]- )/a[k][k] a[i][n]= a[i][n]- a[i][k]* a[k][n] S3:输出解： a[n-1][n]=a[n-1][n]/a[n-1][n] 循环：i=n-2 to 0 a[i][n]=(a[i][n]- )/a[i][i] 输出a[i][n](i=0 to n-1) S4:结束 3.2.6 改进的平方根法 一、求解对象是当A为对称正定矩阵时。此时由A的正定性质，则A可以分解为: 进一步有： 这里L为下三角矩阵。 二、递推公式： 此时，由于U=LT，从而计算量减少了一半，也即只要求解L即可。可以推得递推公式如下： 例 1 利用LLT分解法求解方程组 定理3：若对称矩阵A各阶顺序主子式不为零时, 则A可以唯一分解为A= LDLT ，这里 当A有LDLT分解时，利用矩阵运算法则及相等原理易得计算lij,dk,下面给出推算过程： 由于： A=LDLT 也即： NANJING UNIVERSITY OF INFORMATION SICENCE TECHNOLOGY 第3章 解线性方程组的直接方法 数理学院 陈文兵 南京信息工程大学数学系 3.2 矩阵的直接分解法 本节要点 掌握矩阵分解的概念 直接分解解方程组的步骤 直接分解的本质或理论来源 直接分解的递推计算公式 例子与算法 0 0 3.2.2 直接分解解方程组的步骤 假定我们能把矩阵A写成下列两个矩阵相乘的形式：A=LU ，其中L为下三角矩阵，U为上三角矩阵。这样我们可以把线性方程组 Ax=b (1) 等价地写成 Ax=(LU)x=L(Ux ) = b (2) 令 Ux=y,则原线性方程组Ax=b等价于下面两个方程组的求解： Ly=b (3) Ux=y (4) 于是可首先求解向量 y,使Ly=b; 然后求解 Ux=y,从而达到求解线性方程组 Ax=b的目的. 是否将求解问题复杂化？ 0 第K列 第n-1次消元之后，实际