计算方法2_线性方程组直接法总结.ppt

高斯消去法的乘除总运算分析如下： 消元次数k 消元乘法次数 消元除法次数 回代乘除法总次数 1 n(n-1) n-1 2 (n-1)(n-2) n-2 . . . k (n-k+1)(n-k) n-k . . n-1 2*1 1 n(n+1)/2 故高斯消去法的计算量为 N=n(n2-1)/3+n(n-1)/2+n(n+1)/2=n3/3+n2-n/3 当 N 充分大时为 n3/3 2.4 高斯消去法的计算量 方法的特点 全主元素法的精度优于主元素法,这是由于全主元素是在全体系数中选主元,故它对控制舍入误差十分有效.但全主元素法在计算过程中,需同时作行与列的互换,因而程序比较复杂,计算时间较长。列主元素法的精度虽稍低于全主元素法,但其计算简单工作量大为减少,且计算经验与理论分析均表明,它与全主元素法同样具有良好的数值稳定性,列主元素法是求解中小型稠密性方程组的最好方法之一。 选主元消去法(包括解线性方程组的所有直接的方法)比较适用于中小型方程组.对高阶方程组,即使系数矩阵是稀疏的,但在计算中很难保持稀疏性,因而有存储量大，程序复杂等不足,所幸的是这一缺点可用迭代法解决。 另外,高斯选主元消去法还可技巧性的解决一些特殊线性方程组。 在计算过程中，由于计算机字长的有限性，不可避免地产生舍入误差。同时，由于所求问题的初始数据（例如线形方程组的系数矩阵和右端项系数）往往是带有一定误差的。因此计算结果总是不可避免地带有误差，或者说，如果初始数据有扰动，势必将带来具有一定误差的计算结果。就拿Ax = b来说，由于观测或计算等原因，线性方程组两端的系数A和b都带有误差?A和?b，这样实际建立的方程组是近似方程组(A+?A)(x+?x)=b+?b。对近似方程组求出的解是原问题的真解x加上误差?x,即x+?x。而?x是由?A及?b引起的，它的大小将直接影响所求解的可靠性。 这种解依赖于方程组系数的误差?A及?b的问题，称为线性方程组解对系数的敏感性。 2.5 线性方程组解对系数的敏感性 方程组 此方程组的准确解为x1=0, x2=-1。现将其右端加以微小的扰动使之变为： 经计算可得准确解为x1=2, x2=-3. 这两个方程组的解相差很大，说明方程组的解对常数项b的扰动很敏感。 病态方程组： 如果方程组AX=b由于A或b的小扰动而导致解严重失真，则此方程组称为病态方程组，否则称为良态方程组。 判定一个病态方程组的简单方法； 病态方程组一般不能用解方程组的常用方法求解，而采用“迭代求解法”来计算 2.6 LU分解 当A的所有顺序主子式均不为零时，矩阵A可唯一地分解为两个三角矩阵的乘积 A＝LU 其中，L是单位下三角矩阵，U是上三角矩阵 (2.13) 定义 叫 的三角（因子）分解，其中 是 是上三角。 下三角, 为单位下三角阵（对角元全为1）， 为上三角阵，则称 为Doolittle分解; 若 是下三角， 是单位上三角，则称 定理 n阶阵 有唯一Doolittle分解(Crout） 的前n-1个顺序主子式不为0.（证略） 三角分解不唯一,为此引入 定义 若 为Crout分解。 为什么要讨论三角分解？ 若在消元法进行前能实现三角分解A＝LU，则 容易回代求解 (2.14) L是单位下三角矩阵，因此 (2.15) 其中，yi是向量Y的分量，Y＝(y1,y2,…,yn)T,再从UX＝Y中解出X (2.16) 1.直接三角分解法（以Doolittle分解为例） 设 　　 　 ＝ 由矩阵乘法 ……………………… (2.18) (2.17) [解]先LU分解系数矩阵，由2.17式得 再由2.18式得 对第一个方程组，由2.15式得 再由2.16式得 对第二个方程组，由2.15式得 再由2.16式得 2．平方根法 定理 设A对称正定，则有非奇异下三角阵L，使 ---- 理论基础 (证略） 分解方法：设 ( choleskey分解) 3 追赶法 追赶法仍然保持LU分解特性，它是一种特殊的LU分解。充分利用了系数矩阵的特点，而且使之分解更简单，得到对三对角线性方程组的快速解法。 因三对角矩阵的非零元素呈“带状”，我们也因此将它叫做带状矩阵。