chap2_线性方程组直接法1 .pptVIP

第二章 求解线性方程组的 数值解法 (直接法) 解线性方程组的两类方法: 直接法: 经过有限次运算后可求得方程组精确解的方法(不计舍入误差!) 迭代法：从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法(一般有限步内得不到精确解 列主元消去法 在第k步消元前，在系数矩阵第k列的对角线以下的元素中找出绝对值最大的元。 * 问题：考虑线性方程组 ------------(1) 其中 系数矩阵 未知向量 常数项 根据Cramer(克莱姆)法则,若 当n 比较大时，行列式计算比较困难，计算量大。（O(n!)） 基本想法： 利用方程组的同解变形，逐步将原方程组转化为简单易于求解的特殊形式的线性方程组。 determinantal 行列式的记号 上三角形方程组和下三角形组的解法 前推方向 其解为 回 代 方 向 其解为: 如何将方程组化为三角形方程组 1.高斯消去法 2. 非奇异矩阵三角分解 §2.1 解线性方程组的直接法 一、高斯消去法 思路 首先将方程组Ax=b 化为上三角方程组，此过程称为消去过程，再求解上三角方程组，此过程称为回代过程. §2.1.1 高斯消去法和选主元高斯消去法 解 消去(2.4)中的未知数 得到 第2步. 第1步. 例 用消去法解方程组 将方程(2.2)乘上-2加到方程(2.4)上去， 将方程(2.3)加到方程(2.5)上去，消去方程 (2.5)中的未知数x_2 得到与原方程组等价的三角形方程组 显然，方程组(2.6)是容易求解的，解为 上述过程相当于 由此看出，用消去法解方程组的基本思想是用逐次消去未知数的方法把原方程组 Ax=b 化为与其等价的三角形方程组，而求解三角形方程组可用回代的方法. 上述过程就是用行的初等变换将原方程组系数矩阵化为简单形式(上三角矩阵)，从而将求解原方程组Ax=b 的问题转化为求解简单方程组的问题. 若用初等变换法求解,则对其增广矩阵作行初等变换 经过n-1次 即 同解 将增广矩阵的第 i 行 + li1 ? 第1行，得到： 消元过程： 第一步：设 ，计算因子 其中 第k步：设 ，计算因子 且计算 共进行 n ? 1步，得到 定理：若A的所有顺序主子式 均不为0，则高斯消去法能顺序进行消元，得到唯一解。 回代过程： Gauss消去法的运算量 计算机作乘除运算所耗时间要远远多于加减运算 且在一个算法中，加减运算和乘除运算次数大体相当 故在衡量一个算法的运算量时只需统计乘除的运算次数 乘法次数 除法次数： 全部回代过程需作乘除法的总次数为 于是Gauss消去法的乘除法运算总的次数为 二、 选主元消去法 为避免这种情况的发生, 可通过交换方程的次序，选取 绝对值大的元素作主元. 基于这种思想导出了主元素法 在高斯消去法消去过程中可能出现 的情况，这时 高斯消去法将无法进行；即使主因素 但很小， 其作除数 ，也会导致其它元素数量级的严重增长和舍 误差的扩散 若p≠k，交换第k个与第p个方程后，再继续消去计算. 这种方法称为列主元Gauss消去法。 列主元Gauss消去法保证了｜lik｜≤1 (i=k+1,k+2,…，n). 例. 用Gauss消去法解线性方程组(用3位十进制浮 点数计算） 解: 本方程组的精度较高的解为 用Gauss消去法求解(用3位十进制浮点数计算) 回代后得到 与精确解相比,该结果相当糟糕 究其原因,在求乘数时用了很小的数0.0001作除数 主元 如果在求解时将1,2行交换,即 回代后得到 这是一个相当不错的结果