数学线性代数方程组.pptxVIP

第3章 线性方程组求解的数值方法;1 n阶线性代数方程组的一般形式;在没有舍入误差的情况下，经过有限次 运算可以得到方程组的精确解的方法。;§3.1 Gauss消去与矩阵LU分解;2 上三角方程组与回代过程;3 下三角方程组与前推过程;4 Gauss消去过程;则Gauss消去过程如下：;;顺序高斯消去法必须保证第k 步的akk≠0，因为它在消去过程和回代过程中起关键作用，所以称它为主元素。;例：用基本Gauss消元法求解下列方程组;;?基本Gauss消元法的工作量;基本Gauss消元法的实现条件;小主元 可能导致计算失败;主元素要作除数，其绝对值相对于被除数过小会影响到计算的精度，所以，我们可以采用;?列主元消去法/* Column Pivoting Strategies */; 算法: Gauss列主元消去算法 求方程组Ax=b 的解. 输入:增广矩阵An?(n+1)=（A|b）. 输出: 近似解 xk=ak,n+1(k=1,2,…,n) 或失败信息. 消元过程 for k = 1,2,…,n-1 do Step 1 - Step 4 Step 1 寻找行号 ik , 使得 Step 2 如果 ，则交换第k行和ik行； 否则转Step 7 ; 算法: Gauss列主元消去算法（续） Step 3 for i=k+1,…,n 计算 Step 4 for j=k+1,…,n+1 计算 回代过程 Step 5 Step 6 for i=n-1,…,1 计算 Step 7 Output (系数矩阵奇异); /*不成功 */ STOP. ;例：用Gauss列主元消去法求解下列方程组;Step 1;;;二 LU分解;解方程组的直接三角分解法有3种形式： （1）A=LU;2 解方程组的直接三角分解法;Doolittle分解法; 为单位下三角;设;;?方程组求解的计算公式;例8：用Doolittle分 解法求解下列方程组;§3.2 乔列斯基（Cholesky）分解法(平方根法);（Cholesky分解）;单位上三角;思想;Cholesky分解公式;因对称性无需存储;例10：用Cholesky分 解法求解下列方程组;求解方程组;? Cholesky分解法求解方程组中需说明的几个问题;;;例6：用改进的平方 根法求解下列方程组;求解方程组;;先求解方程组;§3.3 向量范数与矩阵范数;例：设x=(2,-4,3)T,求||x||1,||x||2和||x||∞;称为矩阵ATA的谱半径; §3.4线性代数方程组的迭代解法 /*Iterative Method for Solving Linear Algebraic Systems*/;二、Jacobi迭代法;2 J法的迭代矩阵;二 Gauss-Seidel迭代;2、GS法迭代矩阵;J法适用于并行计算，GS法适用于串行计算，并且所需要的存储空间较小。;§3.5 迭代法的分析;由该定理可知,只要迭代矩阵的谱半径小于1,则J法和GS法收敛.;定理3.2：如果A 是严格对角占优矩阵，则它一定非奇异。;定理3.3：如果方程组的系数矩阵A 是严格对角占优的，则对任意初始近似值，Jacobi迭代法和G-S迭代法均收敛。;§3.6 超松弛(SOR)迭代法/*Overrelaxation Iteration*/;?超松弛(SOR)迭代法的分量形式：;超松弛(SOR)迭代法的迭代公式;解：;;求解方程组 的SOR法收敛的必要 条件是：;§3.7 线性方程组的条件;如果对系数矩阵作”微小”改变;如何评价方程组的“好坏”？;分别考虑右端项及系数矩阵变化对解的影响;定义3.1 设A为非奇异矩阵，称 cond(A)=||A-1||||A|| 为矩阵A的条件数。;定理3.6 对任意的非奇异矩阵 A，cond(A)是由任意的矩阵范数定义的条件数，则;Thank you