则矩阵a可进行三角分解.ppt

列主元消去法 在第k 步消元前，在系数矩阵第k 列的对角线以下的元素中找出绝对值最大的元。 一般计算公式 LU 分解求解线性方程组 Cholesky方法 回顾：对称正定阵A的几个重要性质 （1）A?1 亦对称正定，且 aii 0 （2）A 的顺序主子阵 Ak 亦对称正定 （3）A 的特征值 ?i 0 （4）A 的全部顺序主子式 det ( Ak ) 0 定理 设矩阵A对称正定，则存在唯一的对角元全为正的下三角阵G 使得 A＝GGT 追赶法 * 直接法: 经过有限步运算后可求得方程组精确解的方法(不计 舍入误差!) 迭代法：从解的某个近似值出发，通过构造一个无穷序列去逼 近精确解的方法。 解线性方程组的两类方法 将增广矩阵的第 i 行 + li1 ? 第1行，得到： 消去过程： 第一步：设 ，计算因子 其中 第k步：设 ，计算因子 将增广矩阵的第 i 行 + lik ? 第k行，得到： 其中 定理：若A的所有顺序主子式 均不为0，则高斯消去法能顺序进行消元，得到唯一解。 回代过程： 共进行 n ? 1步，得到 在高斯消去法消去过程中可能出现 的情况，这时 高斯消去法将无法进行；即使主元素 但很小， 其作除数 ，也会导致其它元素数量级的严重增长和舍入 误差的扩散 例：单精度解方程组 /* 精确解为 和 */ 8个 8个 用Gauss消去法计算： 8个 小主元 /* Small pivot element */ 可能导致计算失败。 若p≠k，交换第k个与第p个方程后，再继续消去计算. 这种方法称为列主元Gauss消去法。 列主元Gauss消去法保证了｜lik｜≤1(i=k+1,k+2,…，n). 全主元消去法 在第k步消去前， 在系数矩阵右下角的n-k+1阶主子阵中，选绝对值最大的元素作为主元素。 (1) 如果p ? k，则交换第 k 行与第p行; 如果q ? k，则交换第 k 列与第q列; (2) 消元 注：列交换改变了 xi 的顺序，须记录交换次序， 解完后再换回来。 直接三角分解 每一步消去过程相当于左乘初等变换矩阵Lk A 的 LU 分解 定理 如果Gauss消去法能顺序进行消去，则 矩阵A可进行三角分解，即A＝LU Doolittle分解法 ： 通过比较法直接导出L 和 U 的计算公式。 定理 若A的所有顺序主子式 均不为0，则 A 的 LU 分解唯一（其中 L 为单位下三角阵）。 证明：不难得到，LU 分解存在，下面证明唯一性。 若不唯一，则可设 A = L1U1 = L2U2 ，推出 计算格式为