[数学]数值分析-第五章.pptVIP

例9 求Hilbert矩阵H3的条件数. 如何发现判断矩阵是病态的? 如何解决和处理? 预处理方法. 例10 设 则 化为 则 二、迭代改善法(略去) * * * * * * * * * * * * * * * * * * * * 就有,不选主元的三角分解算法： 于是，可以通过求解两个三角形方程组 得到原方程组的解, 求解方程组计算公式: 说明: L和U的存放;计算∑aibi; 运算量n3/3; Doolittle分解法;求解方程组优点. 练习 利用LU(Doolittle)分解法求解方程组 2、选主元直接三角分解法 当需选主元时，PA=LU. 设第r-1步已完成,就有 选主元三角分解算法： 求解Ly=Pb及Ux=y的算法： 二、平方根法 应用有限元法解结构力学问题时，最后归结为求解线性代数方程组，系数矩阵往往对称正定。平方根法是一种对称正定矩阵的三角分解法，广泛用于求解系数矩阵为对称正定的线性代数方程组。 设A为对称矩阵,且顺序主子式不为零，则 若A为对称正定矩阵,则 解AX=b的平方根法:… 解AX=b的改进平方根法(略去) 三、解三对角方程组的追赶法 在数值求解常微分方程边值问题、热传导方程和建立三次样条函数时，都会要解三对角方程组：AX=b 并且满足 条件(i)保证方程组不能降阶，条件(ii)保证三角分解可做到底。 下面讨论三角分解 比较两边得到 解三对角方程组的追赶法 说明: 稳定性；运算量5n-4次乘除法；存贮. 练习 用追赶法求解三对角线性方程组AX=b，其中 §4 向量和矩阵的范数 为了研究线性方程组的近似解的误差估计和迭代法的 收敛性, 我们需要对Rn中的向量(或Rnⅹn中的矩阵)的 大小引进某种度量——向量(或矩阵)的范数. 先考虑Rn中向量的长度, 然后可定义向量(或矩阵)的范数. 定义1 在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 计算向量 的几种常用范数 定义4(矩阵的范数) (1) 正定性： 等号当且仅当 时成立； (2) 齐次性： (3) 三角不等式： 则称 为 矩阵 的范数或模。 诱导出的常用范数有： 它们满足如下相容关系： 例7 计算矩阵 的几种常用范数 §5 误差分析 一、矩阵的条件数 考虑线性方程组 AX=b 系数矩阵A和右端b的小扰动所产生的相对误差. 例8 方程组 准确解为 常数项微小变化后 准确解 定义7 如果矩阵A或常数项b的微小变化,引起线性方程组AX=b的解的巨大变化,则称此方程组为病态方程组矩阵A称为病态矩阵,否则称方程组为良态方程组,矩阵A为良态矩阵. §1 引言与预备知识 第5章 解线性代数方程组的直接法 一、引言 线性方程组的来源 线性方程组的分类 线性方程组的两类解法： 1、直接法 经过有限步算术运算，可求得方程组精确解的方法(若计算过程中没有舍入误差). 2、迭代法 是用某种极限过程去逐步逼近线性方程组精确解的方法. 二、 向量和矩阵 用 表示全部 实矩阵的向量空间， 表 示全部 复矩阵的向量空间. 这种实数排成的矩形表，称为 行 列矩阵. 称为 维列向量. 其中 为 的第 列. 其中 为 的第 行. 也可写成行向量的形式 写成列向量的形式 (5) 单位矩阵 矩阵的基本运算： (1) 矩阵加法 (2) 矩阵与标量的乘法 (3) 矩阵与矩阵乘法 (4) 转置矩阵 (6) 非奇异矩阵 设 如果 存在， 则称 为非奇异矩阵. 如果 均为非奇异矩阵， 其中 如果 则称 是 的逆矩阵， 记为 且 则 (7) 矩阵的行列式 设 则 的行列式可按任一行(或列)展开， 其中 为 的代数余子式， 行列式性质： 即 的余子式. 为元素 三、矩阵的特征值与谱半径 四、特殊矩阵 对角矩阵 三对角矩阵 上三角矩阵