高斯消去法一.PPTVIP

§1 引言与预备知识 §2 高斯消去法 §3 矩阵的三角分解法 三、解三对角方程组的追赶法 解三对角方程组的追赶法 §4 向量和矩阵的范数 §5 误差分析 为了研究线性方程组的近似解的误差估计和迭代法的 收敛性, 我们需要对Rn中的向量(或Rnⅹn中的矩阵)的 大小引进某种度量——向量(或矩阵)的范数. 先考虑Rn中向量的长度, 然后可定义向量(或矩阵)的范数. 定义2 在Rn中，对 ?=(a1,a2, ?,an)T， ? =(b1,b2, ?,bn)T， 数量积: (?, ?)= ?T?= a1b1+a2b2+ ? +anbn. 欧氏范数: ||?||2 = (?, ?)1/2 . 在Cn中，(?, ?)= ?H?. (1) 正定性： 等号当且仅当 时成立； (2) 齐次性： (3) 三角不等式： 则称 为向量 的范数或模. 由(3)得 (4) 几种常用范数 (无穷范数) (1-范数) (2-范数) (p-范数) 可以验证它们都是范数. 易见前三种范数是p-范数的特殊情况 例6 计算向量 的几种常用范数 定义5(矩阵的范数) (1) 正定性： 等号当且仅当 时成立； (2) 齐次性： (3) 三角不等式： 则称 为 矩阵 的范数或模。 常用矩阵算子范数有： 例7 计算矩阵 的几种常用范数 一、矩阵的条件数 考虑线性方程组 AX=b 系数矩阵A和右端b的小扰动对解X的影响. 例8 方程组 准确解为 常数项微小变化后 准确解 定义7 如果矩阵A或常数项b的微小变化,引起线性方程组AX=b的解的巨大变化,则称此方程组为病态方程组，矩阵A称为病态矩阵,否则称方程组为良态方程组,矩阵A为良态矩阵. 例9 求Hilbert矩阵H3的条件数. * 第5章 解线性代数方程组的直接法 一、引言 线性方程组的来源 线性方程组的两类解法： 1、直接法 2、迭代法 二、向量和矩阵(略) 三、矩阵特征值与谱半径(略) 四、特殊矩阵 对角矩阵 三对角矩阵 上三角矩阵 上海森伯格(Hessenberg)阵 对称矩阵 埃尔米特矩阵 对称正定矩阵 正交矩阵 酉矩阵 初等置换阵 置换阵 定理1 设A∈Rnⅹn, A非奇异?…? 定理2 若A∈Rnⅹn对称正定矩阵，则?…? 定理3 若A∈Rnⅹn对称矩阵，则对称正定矩阵=…? 定理4(若当标准型)… 其中 对角化的条件：1）…；2） …. 一、顺序高斯消去法 设有线性方程组：AX=b 一般地,顺序高斯消去法: （1）消元过程 其中 第一步：若 用 乘第一行 加到第i行中，得到 第二步：若 用… …. … … 第k步：若 用 乘第k行 加到第i行中，得到 其中 第n-1步： … … （2）回代过程 若 则 例（略） 说明: 若线性方程组的系数矩阵非奇异，则它总可以通过带行交换的高斯消去法进行求解。 定理5 (1) 可通过顺序高斯消去法求解. (2)系数矩阵非奇异，总可以通过带行交换的高斯消去 法进行求解。 算法. 二、矩阵的三角分解 下面建立高斯消去法与矩阵的因式分解的关系. 例3’采用3位十进制，用顺序消元法求解 解法1： 三、 列主元素消去法 解法2： 设有线性方程组：AX=b 第一步:先在A的第一列选取绝对值最大的元素作主元素, 然后交换第1行和第i1行(当i1≠1时),再进行第1次消元. … … 第k步选主元素, 然后交换第k行和第ik行(当ik≠k时),再进行第k次消元. … … 第n-1步 回代求解 算法(列主元消去法). … … (见P150) 设有线性方程组：AX=b 一、直接三角分解法 当A的各阶顺序主子式非零时，由顺序高斯消去法知 则有如下的三角分解算法： 于是，可以通过求解两个三角形方程组 得到原方程组的解, 求解方程组计算公式: 练习 利用LU(Doolittle)分解法求解方程组 二、平方根法 应用有限元法解结构力学问题时，最后归结为求解线性代数方程组，系数矩阵往往对称正定。平方根法是一种对称正定矩阵的三角分解法，广泛用于求解系数矩阵为对称正定的线性代数方程组。 设A为对称矩阵,且顺序主子式不为零，则 若A为对称正定矩阵,则 解AX=b的平方根法:… 解AX=b的改进平方根法(略去) 例： 在数值求解常微分方程边值问题、热传导方程和建立三次样条函数时，都会要解三对角方程组：AX=b 并且满足 下面讨论三角分解 比较两边得到 练习 用追赶法求解三对角线性方程组AX=b，其中 *