计算方法第六章解线性方程组的直接法.pptVIP

考虑在整个矩阵范围选主元，这就是所谓的全主元 消去法，此时要注意的是，在做列的变换时，要同 时记录当前变量的次序，以免自变量的含义不清。 范数等价的定义: Step n ? 1: Lk = 其中 记为 L 记 U = 由上述讨论可知，高斯消去法实质上产生了一个将系数 矩阵A分解为上三角阵与下三角阵相乘的因式分解。 若A的所有顺序主子式 均不为0，则 A 的 LU 分解唯一（其中 L 为单位下三角阵）。 设有方程组AX=b,并设A=LU,于是 AX=LUX=b 令UX=Y, 则 LY=b. 于是求解AX=b的问题等价于求解两个方程组UX=Y和LY=b （1）利用顺推过程解LY=b, 其计算公式为: （2）利用回代过程解UX=Y , 其计算公式为: 定理1：（矩阵的三角分解）设A为n ? n实矩阵，如果 解AX = b用高斯消去法能够完成（限制不进行行的交 换，即 ），则矩阵A可分解 为单位下三角矩阵L与上三角知阵U的乘积。 A = LU 且这种分解是唯一的。 定理2：约化主元素（ , i = 1, 2, …, k） 充要条件是矩阵A的顺序主子式 矩阵的三角分解 通过比较法直接导出L 和 U 的计算公式。 思路 （1）对i=1,2,…,n （2）计算 U 的第 r 行, L 的第 r 列元素 对r =2,3…,n 直接三角分解法解AX = b的计算公式 对于r = 2, 3, …, n计算 （2）计算U的第r行元素 （3）计算L的第r 列元素 (r ? n) (1) (4) (5) §4 平方根法 1．矩阵的LDR分解 定理3：如果n阶矩阵A的所有顺序主子式均不等于零， 则矩阵A存在唯一的分解式A = LDR其中L和R分别是 n阶单位下三角阵和单位上三角阵，D是n阶对角元素 的不为零的对角阵，上述分解也称为A的LDR分解。 2．平方根法 如果A为对称正定矩阵，则存在一个实的非奇异 下三角矩阵，使A=LLT ，且当限定的对角元素为正时， 这种分解是唯一的。 定理4：（对称正定矩阵的三角分解） 将对称 正定阵 A 做 LU 分解 U = uij = u11 uij / uii 1 1 1 u22 unn 记为 A 对称 即 记 D1/2 = 则 仍是下三角阵 定理 设矩阵A对称正定，则存在非奇异下三角阵 使得 。若限定 L 对角元为正，则分解唯一。 注： 对于对称正定阵 A ，从 可知对任意k ? i 有 。即 L 的元素不会增大，误差可控，不需选主元。 用平方根法解线性代数方程组的算法 （1）对矩阵A进行Cholesky分解，即A=LLT，由矩阵乘法： 对于 i = 1, 2,…, n 计算 （2）求解下三角形方程组 （3）求解LTX = y 3．改进平方根法 其中 改进平方根法解对称正定方程组的算法 令LTX = y，先解下三角形方程组LDY = b得 解上三角形方程组LTX = Y得 §5 向量和矩阵的范数 1．向量的范数 定义1：设X ? R n，??Ｘ?? 表示定义在Rn上的一个实值函数， 称之为X的范数，它具有下列性质： （3）三角不等式：即对任意两个向量X、Y? R n，恒有 (1) 非负性：即对一切X ? R n，X ? 0, ??Ｘ?? 0 (2) 齐次性：即对任何实数a ? R，X ? R n， --------(1) --------(2) --------(3) --------(4) 常用的向量x的范数有 下述范数的几何意义是: 例3 求下列向量的各种常用范数 例3 求下列向量的各种常用范数 解: 1*4≤9≤9/4*4=9 即 本例中 显然 , 2 1 1 ￥ ￥ ￡ ￡ x c x x c , 设‖·‖A 和‖·‖B 是R上任意两种范数，若存在 常数 C1、C2 0 使得 则称 和 等价。 Rn 上一切范数都等价。 定理1 * AX = b （3.1） 第三章 解线性方程组的直接法 §1 高斯消去法 1．三角形方程组的解法 （3.2） （3.3） 首先将A化为上三角阵 ，再回代求解 。 = (一) 高斯消去法的求解过程,可分为两个阶段: 首先,把原方程组化为上三角形方程组,称之为 “消元”过程; 然后,用逆次序逐一求出三角方程