数值分析线性方程组直接法.pptVIP

数 值 分 析 Numerical Analysis 机械与汽车工程学院 主讲人：孔胜利 kongsl@spu.edu.cn 2012-09-01 第5章 线性方程组的直接法 5.1 线性方程组的简介 5.2 高斯消去法 5.3 高斯主元素消去法 5.4 矩阵三角分解法 5.1 线性方程组的简介 线性方程组的两类解法 直接法: 经过有限步运算后可求得方程组精确解的方法 (不计舍入误差) 迭代法：从解的某个近似值出发，通过构造一个无穷序 列去逼近精确解的方法。分为两类： 逐次逼近法(一般有限步内得不到精确解) 共轭斜量法（不考虑计算过程的舍入误差，只用有限步就收敛于方程组的精确解） 5.2高斯（GAUSS ）消去法 高斯（GAUSS ）消去法 设有线性方程组 或写成矩阵的形式 简记为 例题：用高斯消去法求解线性方程组 其增广矩阵为 回代即可得 5.3 主元素消去法 Gauss消去法能求解 的情况。 主元素消去法 对于一般矩阵来说，最好每一步选取系数矩阵（或消 元后的低阶矩阵）中绝对值最大的元素作为主元素，以 使高斯消去法具有较好的数值稳定性，这就是全主元素 消去法。 目前主要使用的是列主元素消去法和完全主元素消去 法。 高斯--若当消去法：高斯消去法是消去对角线下方的 元素，而高斯—若当消去法是消去对角线上方和下方的 元素。 例题 按列主元消去法求解方程组 解 例题 按完全主元消去法求解方程组 解 首先我们指出，完全主元消去法是有列交换的，如果在增广矩阵的第一行上方再加一行 ，系数阵作列交换时未知数也跟着变换。求解过程如下： 5.4 矩阵三角分解法 直接三角分解法 平方根法 追赶法 习题 习题 习题 习题 习题 习题 习题 习题 习题 习题 习题 习题 习题 *数值分析 Numerical Analysis 5.* r1*(-1)+r2 r1*(-1)+r3 r1*(-4)+r4 r1*(-4)+r5 r2*(-3)+r3 r2*(-0.5)+r4 r2*(-5)+r5 r3*(-0.5)+r4 r3*(-2)+r5 r4*(2)+r5 但是，对于可能出现的 的情况，消去 法将无法进行； 另外，即使主元素 但很小时，用其作 除数，也会导致其它元素数量级的严重增长和舍入误差 的扩散，导致解的不稳定和不可靠。 因此，Gauss消去法适用范围受限。 *数值分析 Numerical Analysis 5.*