数值分析_第5章haha.ppt

 在自然科学和工程技术中,很多问题归结为解线性方程组.有的问题的数学模型中虽不直接表现为含线性方程组,但它的数值解法中将问题“离散化”或“线性化”为线性方程组.因此线性方程组的求解是数值分析课程中最基本的内容之一. §2 高斯消去法 §3 高斯主元素消去法 §4 矩阵的三角分解法 三、解三对角方程组的追赶法 解三对角方程组的追赶法 §5 向量和矩阵的范数 §6 误差分析 如何发现判断矩阵是病态的? 如何解决和处理? 预处理方法. 例10 设 则 化为 则 二、迭代改善法(略去) 结束 * 结束 * 结束 * 结束 * 结束 * 结束 * 结束 * 结束 * 练习 利用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，其中 为了研究线性方程组的近似解的误差估计和迭代法的 收敛性, 我们需要对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 计算矩阵 的几种常用范数 一、矩阵的条件数 考虑线性方程组 AX=b 系数矩阵A和右端b的小扰动所产生的相对误差. 例8 方程组 准确解为 常数项微小变化后 准确解 定义7 如果矩阵A或常数项b的微小变化,引起线性方程组AX=b的解的巨大变化,则称此方程组为病态方程组矩阵A称为病态矩阵,否则称方程组为良态方程组,矩阵A为良态矩阵. 例9 求Hilbert矩阵H3的条件数. * 线性方程组: (5.1) 结束 * 第5章 解线性方程组的直接方法 1 引言与预备知识 常记为矩阵形式 Ax=b (5.2) 此时A是一个n×n方阵,x 和 b是n维列向量. 根据线性代数知识若｜A | ≠0,(5.2)的解存在且唯一. 结束 * 关于线性方程组的解法一般分为两大类,一类是直接法,即经过有限次的算术运算,可以求得(5.1)的精确解(假定计算过程没有舍入误差).如线性代数课程中提到的克莱姆算法就是一种直接法.但该法对高阶方程组计算量太大,不是一种实用的算法.实用的直接法中具有代表性的算法是高斯消元法,其它算法都是它的变形和应用. 另一类是迭代法,它将(5.1)变形为某种迭代公式,给出初 始解 x0 ,用迭代公式得到近似解的序列｛xk｝,k=0,1,2, ?,在一定的条件下 xk→x* (精确解).迭代法显然有一个收敛条件和收敛速度问题. 这两种解法都有广泛的应用,我们将分别讨论,本章介绍 直接法. 结束 * 一、高斯消去法 设有线性方程组：AX=b 例1 解方程组 (5.3) (5.4)