计算方法第三章精要.ppt

} 类似可得L的第k列元素计算公式 求得L、U后再求解方程组Ly=b, Ux=y } 用杜利特尔分解法解方程组的计算步骤: (1)三角分解。对k=1,2,…,n，计算 设A=LU= (2)求解。 ①解Ly=b。对k=1,2,…,n 计算 ②解Ux=y。对k=n,n-1,…,1 计算 } 紧凑存贮的杜利特尔分解法 } } 杜利特尔分解法是顺序高斯消去法的一种变形，当计算过程中|ukk|太小、等于０或者近似等于０时，往往不稳定，以致不能进行分解计算，但只要A非奇异，就可以配合按列选主元的技术进行计算 . } 3.6 线性代数方程组的迭代解法 3.6.1 简单迭代法的一般形式 3.6.2 雅可比(Jacobi)迭代法 3.6.3 高斯-赛德尔(Gauss-Seidel)迭代法 线性代数方程组的迭代解法的基本思想： 构造适当的迭代公式，任选一个初始向量x(0)进行迭代计算，使生成的向量序列 收敛于方程组的精确解。 …， ,… 3.6.1 简单迭代法的一般形式 方程组 Ax=b , A非奇异， 化为等价的方程组 } 构造简单迭代公式 其中 任取. M为迭代矩阵. 简单迭代公式的分量形式 } 若存在极限 则称向量序列 收敛于向量 ，并记为 这时，称简单迭代法 是收敛的，否则就是发散的。 收敛时令k→∞，有 ，等价地有Ax*=b . 控制迭代结束的实用标准： 解向量 其中?是指定的精度要求。 } 迭代公式收敛的情形 （图例的准确解是 x1=C1=0.8, x2=C2=0.7, x3=C3=0.9 ) } 迭代公式发散的情形 （图例的准确解是 x1=C1=3, x2=C2=6, x3=C3= -1 ) } 例如 改写为 ，构造迭代公式 取 其矩阵形式为 其中 } } 计算得 设精度要求为ε=0.005 ，计算结果： 根据迭代公式 5.001 1.999 5.005 2.002 4.99 2.01 4.95 1.98 5.1 1.9 5.5 2.2 4 3 0 0 7 6 5 4 3 2 1 0 k 由于 故求得方程组的解为 } 3.6.2 雅可比(Jacobi)迭代法 公式的构造：方程组Ax=b, 设A非奇异，其主对角元aii≠0，(i=1,2,…,n)，并且绝对值相对来说比较大。 从第i个方程 中分离出xi，得到等价的方程组 雅可比迭代公式 } 例 用雅可比迭代法求解方程组 精度要求为?=0.005。 解 等价的方程组为 构造雅可比迭代公式 } 计算结果: 迭代计算 } … 0.99960 0.998 0.992 0.94 0.8 0 x3(k) … 0.99964 0.998 0.99 0.96 0.7 0 x2(k) … 0.99980 0.9992 0.994 0.98 0.9 0 x1(k) … 5 4 3 2 1 0 k } 所以x=x(5)满足精度要求，得 x1=0.999 90 , x2=0.99974 , x3=0.999 81 } (2)回代过程。若ann＝０，则系数矩阵为奇异，停止计算并转结束；否则计算 ③消元计算：对 i=k+1,k+2, …,n, j=k+1,k+2, …,n+1 计算 } 2、全主元高斯消去法 简称全主元消去法，在第k列(k=1,2, …，n-1)消元时，从系数矩阵的右下角(n-k+1)阶子矩阵中，选取绝对值最大的元素 作为主元素。 } 例 用全主元高斯消去法求解 方程组(用三位有效数字计算) 解 选主元 选主元 } 得到等价的上三角形方程组 回代求得 } 3.2.2 对算法的说明 1. 系数矩阵为对称正定或是严格对角占优的情形 n阶矩阵A为对称正定的充要条件是AT=A且其一切顺序主子式均大于0。 n阶矩阵A为严格对角占优矩阵是指其每个主对角元的绝对值大于同一行其他元素绝对值之和，即 一阶严格对角占优矩阵指一个非零数。 } 例如 由于第1行 所以B是严格对角占优矩阵。 C不是严格对角占优矩阵 当方程组的系数矩阵为对称正定或严格对角占优矩阵时，可用顺序高斯消去法求解。 第2行 第3行 } 2.病态方程组 如果一个方程组的系数或者常数项只要有一点微小的误差就会引起解的巨大变化，这样的方程组就称为是病态的。 例 由于右端系数有误差，方程组变成 解就变成 解的最大误差达200％。 对于病态方程组，前面的算法都可能失效，需要使用一些改进精度的算法。 3.3 行列式和逆矩阵的计算 ⒈行列式的计算： } 例