chap-1直接法.pptVIP

方程组可写成 若令 则上式变成 则求解原方程组Ax=b便转化为 上述矩阵分解法称为Doolittle分解,计算公式为 其中 得 其余全为0。 例2 ：对于例1，由增广矩阵表示消元过程，有 故有 高斯主元素消去法 例1 在只取小数点后四位数字的条件下，用Gauss顺序消去法解方程组 解： 由第二个方程可求出 代入第一个方程可求出 而方程组的准确解为 设 由第一个方程知 因此 而 于是，有 小主元素是算法不稳定的主要原因。 在LU分解法中，若 时，计算无法 进行，或者 绝对值很小时，按分解公式 计算时可能会引起舍入误差的增大，因此， 与Gauss消去法一样，为了保证运算能顺利 进行以及方法的稳定性，三角分解一般采用 选列主元的技术，并称之为Crout分解 。 Gauss-Jordan消去法 Crammer法则 Gauss顺序消去法 Gauss-Jordan LU分解法 运算量比较 线性代数方程组的数值解 Gauss 消去法 矩阵三角分解法 对称矩阵的平方根法 三对角方程组的追赶法 向量与矩阵范数及方程组的性态 解线性方程组的迭代法 分快迭代法 §1.引言 n阶线性方程组： 可以表示成矩阵形式： 解线性方程组的两类方法： 直接法: 经过有限次运算后可求得方程组精确解的方法(不计舍入误差! ） 迭代法：从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法。（一般有限步内得不到精确解) §2. Gauss 消去法 简单消去法 Gauss顺序消去法的可行性及计算量 矩阵的三角分解法 主元素消去法 Gauss-Jordan列主元消去法 一、简单消去法 将n阶线性方程组转化为等价（或同解）的 三角形方程组 的过程称为消元过程，逐次求出 的步骤称为回代过程。 Gauss 消去法计算过程: 统一记号 原方程为 相当于第i个方程-第一个方程×数→新的第i方程—同解！第一方程不动！ Step1 假设 上述消元过程除第一个方程不变以外,第2—第n个方程全消去了变量 ?1，而系数和常数项全得到新值： 得到同解方程组为 其中 Step2.若 ，保留方程组中第一及第二个方程，消去其余方程中变量x3 ,得同解方程组 记作 ，其中 系数矩阵与常数项： 计算出 的过程为消元过程。 消去过程算法 回代过程算法 对于 例1. 用消去法解方程组 解：用增广矩阵表示求解过程 二、Gauss顺序消去法的可行性及计算量 Gauss消去和回代过程都要求元素 称之为主元素。若 ，则计算无法 进行。但在系数矩阵保持一定条件时，则 可保证主元素非零，使Gauss顺序消去法在 计算机上顺利执行。 定理1 若矩阵A的各阶顺序主子式均不为零，即 ，则 定理2 若矩阵A对称正定，则 定理3 若矩阵A严格对角占优，则 矩阵 严格对角占优，若 消去第一列的 n-1 个系数要计算(n-1)+n*(n-1） 次乘法。 Gauss消去法乘法计算量 消去第二列的 n-2 个系数要计算(n-2)+(n-1）*(n-2） 次乘法。 …… 乘法计算量 回代过程总计算量 Gauss消去法乘除法总运算量为 三、 矩阵三角分解法 其中 每一步消去过程相当于左乘初等变换矩阵Lk ，记 (1) 交换矩阵的两行 (2) 以一个非零的数 乘矩阵的某一行(第 行乘数); (3) 把矩阵的某一行的 倍加到另一行(第 行乘 加到 行). 依此类推，得矩阵的LU分解形式 其中