计算方法Gauss消去法ch05ar解答.pptVIP

第五章解线性方程组的直接方法 本章内容 本讲内容 线性方程组直接解法 预备知识 预备知识 特殊矩阵 性质 Gauss 消去法 Gauss 消去法 Gauss 消去法 Gauss 消去法 几点注记 运算量 LU 分解 LU 分解 LU 分解存在唯一性 列主元 Gauss 消去法 列主元 Gauss 消去法 PLU 分解 全主元Gauss消去法 作业 * * 计算方法 —— Gauss 消去法 预备知识 Gauss 消去法 矩阵三角分解法 误差分析 向量与矩阵 特征值与谱半径 一些特殊矩阵 预备知识 一般过程 对应的矩阵三角分解 列主元Gauss消去法 Gauss 消去法 自然科学和工程计算中, 很多问题最终都需要求解一个线性代数方程组 (1) 直接法: 适合低阶方程组或某些特殊大型稀疏方程组 目前使用的数值解法: (2) 迭代法: 解大型稀疏方程组的主流算法 在本章中，我们总是假定 A 是 n 阶方阵 向量与矩阵 预备知识 特征值与特征向量 矩阵的谱： ? (A) = { A 的所有特征值 } 矩阵的迹： 矩阵的谱半径： 相关性质 A 与 AT 有相同的特征值 一些特殊矩阵 对角矩阵、三角矩阵、三对角矩阵 对称矩阵、Hermite对称矩阵、对称正定矩阵 正交矩阵、酉矩阵 初等置换阵、置换阵（排列阵） 上 Hessenberg 矩阵 定理 1 （解的存在唯一性，教材 141 页） 定理 2 （对称正定矩阵的性质，教材 141 页） 定理 3 （对称正定矩阵的充分条件，教材 141 页） 定理 4 （Jordan 标准型，教材 141 页） 例：直接法解线性方程组 解： 高斯消去法的主要思路： 将系数矩阵 A 化为上三角矩阵，然后回代求解。 考虑 n 阶线性方程组： 矩阵形式 = 第一步：消去第一列 依次将增广矩阵的 第 i 行 - mi1 ? 第 1 行，得 设 ，计算 其中 第二步：消去第二列 依次将上述矩阵的 第 i 行 - mi2 ? 第 2 行，得 其中 设 ，计算 记 ，即 。 高斯消去法 第 k 步：消去第 k 列 依此类推，直到第 n-1 步，原方程化为 设 ，计算 回代求解： 计算 ( i = k+1, …, n ) ( i = n-1, …, 1 ) 主元： Gauss 消去法能进行到底的条件：主元全不为 0 定理： （i=1, 2, ..., n）的充要条件是 A 的顺序主子式不为零，即 推论： 第 k 步：消第 k 列 计算 计算 ( i, j = k+1, …, n ) 回代求解： ( i = n-1, …, 2,1 ) n – k 次 (n – k)2 次 n – k 次 n (n+1)/2 次 Gauss 消去法的乘除运算量为： 计算机中做乘除运算的时间远远超过做加减运算时间，故我们只估计 乘除运算 的次数 Gauss消去过程其实就是一个矩阵的三角分解过程 则 A(k) 与 A(k+1) 之间的关系式可以表示为： 其中： ( i = k + 1, …, n ) 将 Gauss 消去过程中第 k-1 步消元后的系数矩阵记为： ( k = 1, …, n-1) LU 分解 记： ，则 其中：L --- 单位下三角矩阵，U --- 上三角矩阵 LU 分解 于是有： 容易验证： ( k = 1, …, n-1) (杜利脱尔Doolittle分解) LU 分解存在 高斯消去法不被中断 定理：若 A 的所有顺序主子式不为零，则 A 存在唯一的 LU分解 （证：自学） 所有顺序主子式不为零 Gauss 消去法有效的条件是： 主元全不为零 例：解线性方程组 ① 先选取列主元： ② if ik ? k then 交换第 k 行和第 ik 行 ③ 消元 列主元 Gauss 消去法 在第 k 步消元时，在第 k 列的剩余部分选取主元 算法 (列主元Gauss消去法 ) for k=1 to n-1 if then stop if ik ? k then swap k-th and ik-th row (including b) for i=k+1 to n end end 列主元 Gauss 消去法对应的矩阵分解为 PLU 定理：若 A 非奇异，则存在排列矩阵 P，使得 PA = LU 其中 L 为单位下三角矩阵，U 为上三角矩阵 列主元 Ga