数值分析5.2.pptVIP

湘潭大学数学与计算科学学院 * 即得 又可以解出 一般地，对 有 由 的正定性可知平方根中值 为正的 湘潭大学数学与计算科学学院 * 由矩阵乘法解得 例 湘潭大学数学与计算科学学院 * 设线性方程组 的系数矩阵 为三对角矩阵 当 的所有顺序主子矩阵非奇异时可作如下分解 2 解三对角方程组的追赶法 湘潭大学数学与计算科学学院 * 追赶法： （1） 分解（Doolittle 分解） 对 计算 湘潭大学数学与计算科学学院 * （2）解 （“追”过程） 对 计算 （3）解 （“赶”过程） 对于 计算 湘潭大学数学与计算科学学院 * 同样可以作Crout分解： 对 计算 湘潭大学数学与计算科学学院 * 湘潭大学数学与计算科学学院 * 上一页 下一页 湘潭大学数学与计算科学学院 * §5.2 Gauss消去法、矩阵分解 湘潭大学数学与计算科学学院 * 2.1 Gauss消去法 下面通过简单例子导出一般算法。 设给定方程组 (1) 湘潭大学数学与计算科学学院 * 乘以第一个方程，这样方程组（1） (2) 化为 其中： 显然方程组（2）和原方程组（1）等价 若 ，则以第 个方程减去 (1) 湘潭大学数学与计算科学学院 * 得到 (3) 其中 依此方法继续下去，得到 以（2）的第 个方程 减去 (2) 湘潭大学数学与计算科学学院 * （4） 其中 从（4）的最后一个方程组得到 湘潭大学数学与计算科学学院 * 再将 代入（4）倒数第二个方程，可得： 类似地，得到： 我们称将方程组（1）按以上步骤化为等价方程组 （4）的过程为解线性方程组的消元过程 从（4）中得出解的过程称为高斯消去法的回代过程 （4） 湘潭大学数学与计算科学学院 * 一般情形 对于一般的 阶线性代数方程组 即 1. 消元过程 首先消去第一列除 之外的所有元素， 湘潭大学数学与计算科学学院 * 设 总可由消元过程得到系数矩阵为上三角阵的线性代数 方程组，其第 步的结果为 湘潭大学数学与计算科学学院 * 其中 这里取 2. 回代过程 若通过消元过程原方程组已化为等价的三角形 方程组 湘潭大学数学与计算科学学院 * 且 , 则逐步回代可得原方程组的解 湘潭大学数学与计算科学学院 * Gauss逐步消去法有如下的缺点： 任一主元 ，就无法做下去 任一 绝对值很小时,也不行（舍入误差的影响大） 2.2 Gauss主元素消去法 下面我们讨论列主元消去法. 假设Gauss消去法的消元过程进行到 第 步， 设 湘潭大学数学与计算科学学院 * 并令 为达到最大值 的最小行标 ， 则交换 和 中的第 行和第 行， 再进行消元过程的第 步. 这时每个乘子 都满足 可以防止有效数字大量丢失而产生误差. 湘潭大学数学与计算科学学院 * 例 用列主元消去法解如下方程组 解 对增广矩阵按列选主元再进行高斯消元 湘潭大学数学与计算科学学院 * 湘潭大学数学与计算科学学院 * 回代求解得 湘潭大学数学与计算科学学院 * 全主元消去法 定义 此时交换 和 的行及A的列，使主元位置的元素 的绝对值具有给出的最大值 ， 然后进行第 步消元过程 注意：因为有列的交换，因此未知量的次序有改变，待消元过程结束时必须还原.多使用列主元消去法. 湘潭大学数学与计算科学学院 * Gauss消去法的实质是将矩阵 分解为 其中 --单位下三角矩阵， --上三角矩阵. 事实上，线性方程组 经过 步消元过程后，有等价方程组 其中： ，而 和 的形式为： 2.3 矩阵的三角分解与Gauss消去法的变形 湘潭大学数学与计算科学学院 * （1） 可以直接验证 ， 湘潭大学数学与计算科学学院 * 湘潭大学数学与计算科学学院 * 其中 湘潭大学数学与计算科学学院 * 乘积 是下三角阵，且对角元全部等