第三章线性代数方程组的直接解法讲述.ppt

定理 3.12 (b的扰动对解的影响) 设A非奇异， Ax=b≠0，且 则有 证:设A精确且非奇异,b有扰动δb,使解x有扰动δx, 则 消去Ax=b，有 又 相比较可得 定理 3.13 (A的扰动对解的影响) 设A非奇异， Ax=b≠0，且 若， 则 证 ：见p66 我们还可证明更为一般的结论： 当方程组的系数矩阵A非奇异和常数项b为非零向量时，且同时有扰动δA，δb，满足 ，若x和x+δx分别是方程组 Ax=b 及 的解 则 例3.13 线性方程组 的系数矩阵带误差，成为如下方程组 求方程组系数矩阵的条件数, 并说明方程组的性态 解 因为 所以 因此方程组是良态的 3.7.2 精度分析 求得方程组Ax=b的一个近似解以后,希望判断其精度，检验精度的一个简单办法是将近似解再回代到原方程组去求出余量r. r = b-A 如果r很小，就认为解是相当精确的。 定理3.14 设 是方程组Ax=b的一个近似解,其精确解记为 ,r为 的余量。则有 证明见P68 例3.14 设A为正交矩阵，证明：cond2(A)=1 分析：由正交矩阵和条件数的定义便可推得 解：因为A是正交矩阵， 故ATA= AAT=I， A-1= AT，从而 例3.15 设A，B为n阶矩阵，证明： cond(AB)≤ cond(A) · cond(B) 分析: 由矩阵范数性质和条件数定义 便可证明 证： cond (AB) = || AB || · || (AB)-1 || ≤ ||A|| · ||B|| ||A-1|| · ||B-1|| = ||A|| · ||A-1|| ||B|| · ||B-1|| = cond (A) · cond (B) 例3.16 设A，B为n阶非奇异矩阵，||?||表示 矩阵的任一种范数，证明： || A-1-B-1 || ≤|| A-1 || || B-1 || || A-B || cond(AB)≤ cond(A) · cond(B) 分析: 由矩阵范数的基本性质即可推证 证： A-1-B-1 = A-1(B-A)B-1 ,从而 || A-1-B-1 ||≤|| A-1(B-A)B-1 || ≤ ||A-1|| · ||B-A|| · ||B-1|| ∴ || A-1-B-1 || ≤ ||A-1|| · ||B-A|| · ||B-1|| 本章小结 本章介绍了解线性方程组的直接法。直接法是一种计算量小而精度高的方法。直接法中具有代表性的算法是高斯（Gauss）消去法（在第一章提到的克莱姆算法也是一种直接法，但该算法用于高阶方程组时计算量太大而不实用），其它算法大都是它的变型，这类方法是解具有稠密矩阵或非结构矩阵（零元分布无规律）方程组的有效方法。 选主元的算法有很好的数值稳定性。从计算简单出发实际中多选用列主元法。 解三对角矩阵方程组（A的对角元占优）的追赶法，解对称正定矩阵方程组的平方根法都是三角分解法，且都是数值稳定的方法，这些方法不选主元素，也具有较高的精度。 向量、矩阵的范数、矩阵的条件数和病态方程组的概念，是数值计算中一些基本概念。线性方程组的病态程度是其本身的固有特性，因此即使用数值稳定的方法求解，也难以克服严重病态导致的解的失真。在病态不十分严重时，用双精度求解可减轻病态的影响 在实际应用中如何选择算法是一个重要问题，往往从三个方面考虑： ?????? ① 解的精度高； ????? ② 计算量小； ?????? ③ 所需计算机内存小。 但这些条件相互间是矛盾而不能兼顾的，因此实际计算时应根据问题的特点和要求及所用计算机的性能来选择算法。一般说，系数矩阵为中、小型满矩阵，用直接法较好；当系数矩阵为大型、稀疏矩阵时，有效的解法是第四章要讨论的迭代法。 作业 习题一3.1~3.13 * * This is a placeholder for the demo. It reminds you when to switch over to the de