计算方法（潘林强）NA06b.pptVIP

* §2 三角分解法 /* Matrix Factorization */ ? 高斯消元法的矩阵形式 /* Matrix Form of G.E. */： Step 1: 记 L1 = ，则 Step n ? 1: 其中 Lk = §2 Matrix Factorization – Matrix Form of G.E. 记为 L 单位下三角阵 /* unitary lower-triangular matrix */ 记 U = A 的 LU 分解 /* LU factorization */ Hey hasn’t GE given me enough headache? Why do I have to know its matrix form??! When you have to solve the system for different with a fixed A. Could you be more specific, please? Factorize A first, then for every you only have to solve two simple triangular systems and . §2 Matrix Factorization – Matrix Form of G.E. 定理 若A的所有顺序主子式 /* determinant of leading principal submatrices */ 均不为0，则 A 的 LU 分解唯一（其中 L 为单位下三角阵）。 证明：由§1中定理可知，LU 分解存在。下面证明唯一性。 若不唯一，则可设 A = L1U1 = L2U2 ，推出 Upper-triangular Lower-triangular With diagonal entries 1 ? 注： L 为一般下三角阵而 U 为单位上三角阵的分解称为Crout 分解。 实际上只要考虑 A* 的 LU 分解，即 ，则 即是 A 的 Crout 分解。 §2 Matrix Factorization – Doolittle ? 道立特分解法 /* Doolittle Factorization */： —— LU 分解的紧凑格式 /* compact form */ 反复计算, 很浪费哦 …… 通过比较法直接导出L 和 U 的计算公式。 思路 §2 Matrix Factorization – Doolittle 固定 i : 对 j = i, i+1, …, n 有 lii = 1 a 将 i ，j 对换，对 j = i, i+1, …, n 有 b Algorithm: Doolittle Factorization Step 1: u1j = a1j; lj1 = aj1 / u11; ( j = 1, …, n ) Step 2: compute and for i = 2, …, n?1; Step 3: a b 一般采用列主元 法增强稳定性。但注意 也必须做相应的 行交换。 §2 Matrix Factorization – Choleski ? 平方根法 /* Choleski’s Method */： ——对称 /* symmetric */ 正定 /* positive definite */ 矩阵的分解法 定义 一个矩阵 A = ( aij )n?n 称为对称阵，如果 aij = aji 。 定义 一个矩阵 A 称为正定阵，如果 对任意非零向量 都成立。 ?回顾：对称正定阵的几个重要性质 ? A?1 亦对称正定，且 aii 0 ? 若不然，则 存在非零解，即 存在非零解。 ? ? 对任意 , 存在 , 使得 , 即 。 ? 其中 第 i 位 ? A 的顺序主子阵 /* leading principal submatrices */ Ak 亦对 称正定 对称性显然。对任意 有 , 其中 。 ? A 的特征值 /*