第六章 解线性方程组的消去法.pptVIP

§6.1 约当(Jordan)消去法 二、若当消去法： §6.2 高斯(Gauss)消去法 §6.3 追赶法 由（6.18）的第一个方程我们得 ——6.19 由（6.18）的第二个方程又得 ——6.20 于是（6.20）成为 以此类推。一般地，令 为了能兼容（6.19）式，规定 ——6.21 为此先由（6.18）的最后一个方程和（6.21）中k=n-1的方程联立可解得 * * 引 言 在实际问题中，存在大量的解线性方程组的问题：如求一块方形平板的二维稳定状态，导热过程的温度分布，给了25个节点，就需解一个9阶方程组。另外很多数值方法到最后也会涉及到线性方程组的求解问题：如 三次样条插值最后的三对角线性方程组，曲线拟合的法方程，方程组的Newton迭代等问题。因此掌握线性方程组的求解方法是很有必要的。 ※ 线性代数中常用克莱姆法则求解线性方程组，但其 不适合于求解大型方程组。 ※ 工程计算中的线性方程组，其系数阵大致有两种， 从而解线性方程组的方法可以分为2类： ①低阶稠密阵：常用消去法求解，其特点是结果准， 可靠，理论上得到的解是精确的。 ②大型稀疏阵：由于此类矩阵阶数高，但0元素多，故 常用迭代法求解，特点是求解速度快，但结果有误差。 ※消去法的基本思想： ※ 消去法有多种计算格式，这里我们只介绍较为实用的约当消去法和高斯消去法。 将一个方程乘或除以某个常数，或者将两个方程相加减，通过这两种手续逐步减少方程当中变元的数目，从而得到所求的解。 一、［引例］解方程组： 从而得到方程组（6.3），其中只有第一个方程含有 这就是所求的解。 上述方法就是所谓若当消去法，其消元过程可归纳如下： … … … 以此类推，直到每个方程仅有一个变元为止。 对一般形式的线性方程组 ——6.5 用约当消元法第k－1步消元后得到等价形式如下： 符号说明： 消元过程的第k步要做两项工作： ——6.7 注意到(6.6)的第i个方程中 为此需要进行的计算是 ——6.9 算式(6.7)、(6.9)是关于下标k的递推公式，对第一步，即k=1时，所要得到的方程组(6.6)的原始形式是所给方程组(6.5)，因此令系数 ——6.10 于是得到所求解。 三、算法与流程图： 1.算法的总计算量： ——6.11 次乘除法。 2.算法的压缩存储： 3.流程图：（略） ※ 高斯方法是约当方法的一种改进，它在约当方法的基础上减少了计算量。由其改进得到的选主元的高斯消元法是目前计算机上常用的解低阶稠密矩阵方程的有效方法。 一、［引例］解方程组： [解] 最后对上述方程组，用回代法，求得原方程组的解为： ♀上述求解过程的矩阵描述： 二、基本思想： 用矩阵的初等变换将系数矩阵A约化为具有简单形式的上三角形矩阵，而三角形矩阵是容易解的。 三、解法： ——6.5 设有n个未知数的线性方程组 ——6.13 进一步加工将与约当法不同。 高斯法的第二步，不动方程（6.12），将它留到回代时再处理，第二步将针对方程组（6.13）继续消元。 如此继续下去，第k步所要加工的方程组是： ——6.14 由于高斯法的第k步计算并不涉及到前k-1个方程，反映在算式上，i应从k+1算起，于是高斯法的算式为： ——6.15 ——6.16 按照上述手续进行n步后，不像约当法那样立即得到方程组之解，而将原方程组加工成如下形状： ——6.12 这种形式的方程组是三对角形方程组，其系数矩阵为 主对角线以下元素全为零。 消元过程归结出的（6.12）通过逆“回代过程”很容易求解，回代公式是 ——6.17 即 ♀ 小结： 高斯方法可分为消元过程和回代过程两个环节。 消元过程将所给方程（6.5）加工成三角方程组（6.12），其系数按（6.15）、（6.16）计算。 对三角方程组（6.12）通过回代可求解，按求解公式（6.17）计算。 四、算法与流程图： 1.算法的总计算量：（加减运算可忽略） （1）消元计算： （2）回代计算： （3）总计算量： ♀：高斯方法由于添加了回代过程，算法结构较约当法稍许复杂了一点，但因此节省了约六分之一的计算量。 2.流程图：（略） 五、选主元素的高斯消元法： 单精度解方程组 /* 精确解为 和 */ 8个 8个 用Gaussian 消元法计算： 8个 小主元可能导致计算失败。 ［例］ ? ［改］ 进行行交换，避免小主元。 1.列主元素高斯消元法： 用高斯消去法解方程组时，小主元可能导致计算失败。故在计算时应该选主元，即在消去的每一步应该在系数矩阵或消元后的低阶矩阵中选