解线性方程组的直接法4.3.ppt

EAST CHINA INSTITUTE OF TECHNOLOGY 理学院 4.3 Gauss列主元消去法 例1. 用Gauss消去法解线性方程组(用3位十进制浮 点数计算） 解: 本方程组的精度较高的解为 用Gauss消去法求解(用3位十进制浮点数计算) 一、Gauss列主元消去法的引入 T x ). 1 ,. 0 ( * = -9999 回代后得到 与精确解相比,该结果相当糟糕。 究其原因,在求行乘数时用了很小的数0.0001作除数 -9998 如果在求解时将1,2行交换,即 0.9999 回代后得到 这是一个相当不错的结果。 0.9998 例2. 解线性方程组(用8位十进制尾数的浮点数计算) 解: 这个方程组和例1一样,若用Gauss消去法计算会有 小数作除数的现象,若采用换行的技巧,则可避免 绝对值最大 不需换行 经过回代后可得 事实上,方程组的准确解为 例2所用的方法是在Gauss消去法的基础上,利用换行 避免小主元作除数,该方法称为Gauss列主元消去法 二、Gauss消元过程与系数矩阵的分解 1.Gauss消去法消元过程的矩阵描述 行变换相 当于左乘 初等矩阵 由于 令 则 显然若令 则有 因此 从而 故 即 且 定义1. 不带行交换的Gauss 消去法的消元过程,产生 一个单位下三角矩阵L和一个上三角矩阵U,即 该过程称之为 由上述分析不难得到 Gauss消去法 可以执行 证明略 定理1 2.Gauss列主元消去法消元过程的矩阵描述 由于Gauss列主元消去法每一步都要选取列主元,因 此不可避免要进行行交换 即 初等矩阵 因此,Gauss列主元消去法的消元过程为 : 显然 上三 角阵 仍然为单位 下三角矩阵 消元过程为 时 当 例如 , 4 , = n 初等矩阵的乘积,称为排列阵 则 推广到一般情形 令 仍然为单位 下三角矩阵 则 单位下三角阵与上三角阵的乘积 综合以上讨论,有 定理2 开始 输出无解信息 消元 换行 停机 回代求解 Gauss列主元消去法的算法设计 (一) 流程图 (二) 自然语言 选主元 换行 消元