第三章线性方程组数值解法1.ppt

第3章 线性方程组的数值解法 §1 消元法 §2 矩阵的直接分解法 §3 范数 §4 迭代法 ? §5 迭代法的收敛性 线性方程组的概念 例1（续） Gauss消元法的基本步骤1（4阶） Gauss消元法的基本步骤2（4阶） Gauss消元法的基本步骤3（4阶） Gauss消元法的基本步骤4（4阶） Gauss消元法的消元过程1、2（n阶） Gauss消元法的消元过程3（n阶） Gauss消元法的消元过程3（n阶） Gauss消元法的回代过程（n阶） Gauss消元法的计算量 Gauss法与Cramer法则的计算量比较 高斯主元素消元法是消去法的一种改进。它的基本思想是在逐次消元时总是选绝对值最大的元素(称之为主元)做除数,按消元法的步骤消元。 这里主要介绍求解线性方程组最常用的列主元素消去法 列主元素消去法 所谓列主元素消去法就是在每一步消元过程中取系数子矩阵的第一列元素中绝对值最大者作主元。对线性方程组(3―1)进行n-1次消元后,可得到上三角形方程组 例2（续1） 例2（续2） 例2两种解法的误差分析 总结 1、高斯消去法的计算步骤为: 1）消元过程 对于k=1,2,…,n-1,若按顺序有某一ark≠0,r≥k,则交换k与r行,然后计算 2、列主元素消元法的计算过程为: (1)消元过程。对于k=1,2,…,n-1进行下述运算: ①选主元,确定r,使得 若ark=0,说明系数矩阵为奇异,则停止计算否则进行下一步。 ②交换A的r、k两行 ③消元：对i=k+1,k+2,…,n,j=k+1,k+2,…,n+1计算 算法： Step1:输入系数矩阵A; Step2:for k = 1 to n-1 反复做step3~step6 Step3: max = |akk|; //找最大元 for i = k+1 to n 做 if(max |aik|) 则 max = |aik| ; r = i; //r记录最大元的行号 Step4: if (ark = =0) ,则打印“没有唯一解”，并转step10，否则做step5 Step5:交换矩阵的第r行,第k行：//交换行 if (r≠k)则 for (j= k to n+1) akj←→arj Step6:for i=k+1 to n //消元，i控制行号，j控制列号 for j = k+1 to n+1 Step7:if ann ==0则打印“没有唯一解” ，转step10 Step8: xn = ann+1/ann , for k = n-1 to 1 做 Step9:输出解：x1,x2,……,xn Step 10 : 结束 计算|A|算法： Step1:输入系数矩阵A; det = 1; Step2:for k = 1 to n-1 反复做step3~step6 Step3: max = |akk|; for i = k to n 做 if(max |aik|) 则 max = |aik| ; r = i; Step4: if (ark = =0) ,则打印“没有唯一解”，并转step10，否则做step5 Step5:交换矩阵的第r行,第k行： if (r≠k)则 for (j= k to n) akj←→arj Step6:for i=k+1 to n for j = k+1 to n Step7: for k = 1 to n det = det*akk Step8:输出 det Step 9: 结束 3.3 矩阵分解法 杜利特尔（Doolittle）分解——LU分解 对A进行LU分解 对A进行LU分解的具体步骤 矩阵A的LU分解举例 三角分解的紧凑格式 三角分解的紧凑格式举例 直接三角分解法 三角分解的紧凑格式举例 紧凑格式解线性方程组举例 三角分解法的几点说明 3.3.2 平方根法与改进的平方根法 3.3.2 平方根法 Choleskg分解1 Choleskg分解2 平方根法举例 改进的平方根法 改进的平方根法说明 改进的平方根法举例 3.4 向量范数和矩阵范数 常用的向量范数 对2范数 Rn中范数的等价性 向量的误差 矩阵范数 常用的矩阵范数 最大