第3章 解线性方程组及直接法yjs.ppt

解三对角方程组的追赶法 其计算工作量为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三角分解. 定理：（对称正定矩阵的三角分解） 如果A为对称正定矩阵，则存在一个实的非奇异 下三角矩阵L，使A=LLT ，且当限定的对角元素为正时， 这种分解是唯一的。 U = uij = u11 uij / uii 1 1 1 u22 unn 记为 A 对称 即 记 D1/2 = 则 仍是下三角阵 设A为对称正定矩阵, 则有唯一分解A=LU, 且ukk0. 分解A=LLT称为对称正定矩阵的Cholesky分解. 如何求L 对矩阵A进行Cholesky分解，即A=LLT，由矩阵乘法： 对于 i = 1, 2,…, n 计算（一列一列算） 求解下三角形方程组 Ly=b 求解LTx = y Cholesky分解法缺点及优点 2. 对于对称正定阵 A ，从 可知对任意k ? i 有 。即 L 的元素不会增大，误差可控，不需选主元。 缺点：存在开方运算，可能会出现根号下负数。 优点：1.可以减少存储单元。为一般矩阵LU分解计算量的一半。 逐行相乘，并注意到ij，有 由此可得 所以将上述公式改写为 解方程，令LTX = y，先解下三角形方程组LDY = b得 再解上三角形方程组LTX = Y得 算法 3.4：( Cholesky 分解法 ) 步1 输入对称正定矩阵A，右端项b； 步2. Cholesky分解： 对 k=2,3,…,n，计算： 步3.解下三角形方程组LDy = b 步4.解上三角形方程组LTx = y 程序见P53 先看一个例子，设线性方程组 试分析系数矩阵和右端项有微小扰动，解将产生什么样的变化？ 解 方程组的精确解为x=(1,1)T 设系数矩阵有微小扰动 即 若右端项有微小变动 则 若同时对A,b扰动?A ，?b，则 现计算相对误差： x*是方程组(A+△A) x=b+△b (1)的准确解。 在没有舍入误差的解。称方程（1）为方程Ax=b的扰动方程。其中△A， △b为由舍入误差所产生的扰动矩阵和扰动向量。当△A， △b的微小扰动，解得（1）的解与Ax=b的解x的相对误差不大称为良态方程，否则为病态方程。 设线性方程组 Ax=b 1) 系数矩阵是精确的, 常数项有误差?b, 此时记解为x+?x ,则 A(x+?x) =b+?b 于是 A?x =?b 所以 ‖?x ‖= ‖A-1?b‖? ‖A-1‖ ‖?b‖ 又由于 ‖b‖= ‖Ax‖? ‖A‖ ‖x‖ 因此 ‖?x ‖‖b‖? ‖A‖‖A-1‖‖?b‖‖x‖ 即 误差界 2) 再设b是精确的, A有误差?A, 此时记解为x+?x ,则 (A+?A)(x+?x) =b 则有 A?x +?A(x+?x) =0 所以 ?x =-A-1?A(x+?x) 于是 ‖?x‖? ‖A-1‖‖?A‖‖x+?x‖ 也就