(第章解线性方程组的直接法.pptVIP

解三对角方程组的追赶法 其计算工作量为6n-5次乘除法。工作量小，有以下定理可得证三对角矩阵求解的充分性条件。 LU分解方式 LU分解算法（续） LU分解算法（续） 3.5 对称正定矩阵的Cholesky分解法 A对称：AT=A A正定：A的各阶顺序主子式均大于零。即 改进Cholesky分解法 改进的cholesky分解A=LDLT 改进的cholesky分解 改进的cholesky分解 3.6 线性方程组的性态和舍入误差对解的影响 病态方程组 由于计算机字长限制，在解Ax=b时，舍入误差是不可避免的。因此我们只能得出方程的近似解x*。 病态方程组的处理 病态方程组的计算 1) 用双精度或更高精度计算 2) 采用数值稳定性较好的算法，如全主元法 3) 用迭代改善法 先解 ,得 再解 ,得 解线性方程组Ax=b的Doolittle三角分解法的计算量约为n3/3,与Gauss消去法基本相同.其优点在于求一系列同系数的线性方程组Ax=bk ,(k=1,2,…,m)时,可大大节省运算量. 此外，还有另一种LU分解，称为Crout三角分解. Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 定理：（对称正定矩阵的三角分解） 如果A为对称正定矩阵，则存在一个实的非奇异 下三角矩阵L，使A=LLT ，且当限定的对角元素为正时， 这种分解是唯一的。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. U = uij = u11 uij / uii 1 1 1 u22 unn 记为 A 对称 即 记 D1/2 = 则 仍是下三角阵 设A为对称正定矩阵, 则有唯一分解A=LU, 且ukk0. 分解A=LLT称为对称正定矩阵的Cholesky分解. 如何求L Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 对矩阵A进行Cholesky分解，即A=LLT，由矩阵乘法： 对于 i = 1, 2,…, n 计算（一列一列算） 求解下三角形方程组 Ly=b 求解LTx = y Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. Cholesky分解法缺点及优点 2. 对于对称正定阵 A ，从 可知对任意k ? i 有 。即 L 的元素不会增大，误差可控，不需选主元。 缺点：存在开方运算，可能会出现根号下负数。 优点：1.可以减少存储单元。为一般矩阵LU分解计算量的一半。 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 逐行相乘，并注意到ij，有 由此可得 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profile 5.2.0.0. Copyright 2004-2011 Aspose Pty Ltd. 所以将上述公式改写为 解方程，令LTX = y，先解下三角形方程组LDY = b得 再解上三角形方程组LTX = Y得 Evaluation only. Created with Aspose.Slides for .NET 3.5 Client Profi