计算方法31.pptVIP

线性方程组的数值解法大致可分为两种： 直接法（精确解法）：在没有舍入误差的条件下，经过有限次四则运算求得方程组的精确解的方法. 优点:可以事先估计计算量. 适用于：低阶稠密矩阵方程组以及某些大型稀疏方程组. 迭代法：用某种极限过程去逼近方程组的精确解的方法. 优点：简单，所需计算机的存储单元少,易编制程序. 适用于：大型稀疏系数矩阵的方程组. 顺序高斯消去法 列主元高斯消去法 高斯－约当消去法 小结 了解直接法解方程组的基本思想； 掌握顺序Gauss消去法的使用条件与一般步骤； 会用顺序Gauss消去法解方程组. 第三章 线性代数方程组的数值解法 高斯消去法 矩阵的三角分解法 向量和矩阵的范数 方程组的性态和误差分析 迭代法 迭代的收敛性 序 本章主要讨论n阶线性代数方程组 的数值解法. 未知向量 常数向量 系数矩阵 其矩阵形式为 Ax=b ,其中 §3.1 高斯消去法 一、顺序高斯消元法 自然顺序依次消元,把原方程组转化为同解的上三角形方程组. 基本思想： 按照方程自上而下、未知元从左到右的 当 时 ,方程组有唯一解. 自下而上回代，即可求出方程组的解. 系数矩阵A 是上三角形矩阵. 1.消元过程：用矩阵的初等行变换把原方程组的 2.回代过程：对上三角形方程组按自下而上的顺序回代求解. 基本步骤： 的顺序把主对角元下方的元素化为0. 依自上而下,从左到右 增广矩阵化为上三角形矩阵. 注意：1. 顺序高斯消去法不作行交换； 2. 化零时，用行变换 ri－krj ； 例1 用高斯消去法解方程组 解 （消元过程用增广矩阵的初等行变换来实现） 第一次消元 第二次消元 然后回代，解得： ， 过程分析： 第1次消元时， 第一个主行的主元素 然后计算消去行的乘数： 第2次消元时， 系数行列式的计算： 消元过程 det A=1×2×(-7)= -14 把上述过程一般化，对线性方程组 Ax=b ,记为 主元为1，2，-7， 第二个主行的主元素 然后计算消去行的乘数： 对线性方程组 Ax=b的增广矩阵施行行初等变换: ① 当 时， 计算乘数: 作行变换 ，得： 第一次消元 于是： ② 当 时， 计算乘数: 作行变换 ，得： 于是： 第二次消元 重复以上过程...... , 当 时， 经过n－1次消元后，就把原方程组化为同解的三角形方程组 ② 系数行列式的计算： 说明： 当 时，自下向上回代，即得原方程组的解. ① 顺序高斯消去法的使用条件是主元 原因：A的k阶顺序主子式： ③ 系数矩阵A的n个顺序主子式全不为零，则顺序高斯消去法能实现方程组的求解.(定理3-1) ④ A是严格对角占优矩阵,则可实现顺序高斯消去法.(定理3-2) 严格行对角占优矩阵是指矩阵的每个主对角元的绝对值大于同一行其他元素绝对值之和，即: (定义3-1) 顺序高斯消去法的运算量 包含加减乘除运算的算法，一般只关注乘除的运算量. 除法次数(n-k) 乘法次数(n-k) 第 k 步 消 元 第k步消元的乘除法运算次数为：(n-k)(n-k+2)次； 乘法次数(n-k)(n-k) 完成全部n-1步消元的乘除法运算次数为： 全部回代过程需作乘除法运算的次数为： 于是顺序高斯消去法的乘除法运算总的次数为： 乘法次数（n-i） 除法次数 1 回代 与n2同阶大 二、列主元高斯消去法 例如： 用顺序消去法解方程组 第一次消元 ① 若第k步中主元素 ② 就可能造成误差的严重扩散，从而造成结果严重失真。 若其绝对值相对较小(此时称 为小主元),它作除数， 顺序消去法中可能会遇到下面两个问题： 则消元过程就无法进行； 如果选用2.000作为主元素（即先交换两个方程的位置）， 然后再进行消元，则有 交换两行 第一次消元 则得 主元素，这种消元法称为主元素消元法。 由此，在每次消元前，应选取绝对值尽可能大的元素作为消元 优点在于没有增加求解过程的运算量，但大大减小误差. 完全主元消去法： 每次消元前，在系数矩阵的某一行中， 选择绝对值最大的元素作为主元的消去法。 每次消元前，在系数矩阵的某一列中， 选择绝对值最大的元素作为主元的消去法。 行主元消去法： 列主元消去法： 每次消元前，在系数矩阵的所有元素中， 选择绝对值最大的元素作为主元的消去法。 按选取主元素的方法不同，主元素消去法可分为以下几种： 列主元消去法最为常用，其一般步骤为： 进行第k次消元前有增广矩阵