数值分析第5篇.pptx

1 第5章 解线性方程组的直接方法 5.1 引言与预备知识 5.2 高斯消去法 5.3 矩阵三角分解法 5.4 向量和矩阵的范数 5.5 误差分析 2 5.1 引言与预备知识 5.1.1 引言 线性方程组的数值解法一般有两类： 1. 直接法 经过有限步算术运算，可求得方程组精确解的方法(若 计算过程中没有舍入误差).但实际计算中由于舍入误差的存 在和影响，这种方法也只能求得线性方程组的近似解. 2. 迭代法 是用某种极限过程去逐步逼近线性方程组精确解的方法. 3 5.1.2 向量和矩阵 4 也可写成行向量的形式 矩阵可以写成列向量的形式 5 对于方程组 6 矩阵的基本运算： 7 (6) 非奇异矩阵 8 行列式性质： (7) 矩阵的行列式 9 5.1.3 矩阵的特征值与谱半径 定义1 设 ，若存在数 （实数或复数） 和非零向量 ，使 （1.1） （1.2） 10 （1.3） 11 因为 次代数方程 在复数域中有 个根 由(1.3)中的行列式展开可得 故 12 并且 及 13 14 5.1.4 特殊矩阵 15 16 (10) 初等置换阵 (11) 置换阵 由初等置换阵的乘积得到的矩阵. 17 18 其中 19 20 5.2.1 高斯消去法 5.2 高斯消去法 21 例2 用消去法解方程组 22 得到与原方程组等价的三角形方程组 显然，方程组(2.6)是容易求解的，解为 上述过程相当于 23 上述过程就是用行的初等变换将原方程组系数矩阵化 为简单形式(上三角矩阵)，从而将求解原方程组(2.1)的 问题转化为求解简单方程组的问题. 24 下面讨论求解一般线性方程组的高斯消去法. 其中 25 26 即已计算好与(2.1)等价的方程组 27 28 由(2.1)约化为(2.10)的过程称为消元过程. 29 求解三角形方程组(2.10)，得到求解公式 (2.10)的求解过程(2.11)称为回代. 30 继续这过程，高斯消去法照样可进行计算. 31 且计算公式为： 32 (b) 回代计算 33 34 证明 略 35 5.2.2 矩阵的三角分解 下面借助矩阵理论进一步对消去法作些分析，从而建 立高斯消去法与矩阵因式分解的关系. 36 其中 其中 37 重复这过程，最后得到 38 其中 为单位下三角矩阵. 39 证明 略 40 例3 对于例2，系数矩阵 41 例3 求解方程组 5.2.3 列主元消去法 42 用4位浮点数进行计算. 解法1 用高斯消去法 精确解舍入到4位有效数字为 43 计算解为 显然计算解是一个很坏的结果，不能作为方程组的近似解. 其原因是我们在消元计算时用了小主元0.001，使得约化后的方程组元素数量级大大增长，经再舍入使得在计算(3,3)元素时发生了严重的相消情况，因此经消元后得到的三角形方程组就不准确了. 44 解法2 交换行，避免绝对值小的主元作除数. 得计算解为 45 对一般矩阵来说，最好每一步选取系数矩阵(或消元后 的低阶矩阵)中绝对值最大的元素作为主元素，以使高斯消 去法具有较好的数值稳定性. 这就是全主元素消去法. 在选主元时要花费较多机器时间，目前主要使用的是 列主元消去法. 46 设方程组(2.1)的增广矩阵为 47 48 最后将原方程组化为 49 回代求解 算法1 (列主元素消去法) 50 (1) 按列选主元 (4) 消元计算 51 4. 回代求解 52 例3的解法2用的就是列主元素消去法. 列主元素消去法也可用矩阵运算描述： 53 利用(2.15)得到 则有 54 由(2.16)得到 55 定理8 (列主元素的三角分解定理)如果 为非奇异矩阵, 在实际计算中只能在计算过程中做行的交换. 56 5.3 矩阵三角分解法 5.3.1 直接三角分解法 57 1. 不选主元的三角分解法 58 由(3.1)有： 有 故 59 又由(3.1)有 60 61 例5 用直接三角分解法解 62 从而 63 64 65 2. 选主元的三角分解法 采用与列主元消去法类似的方法，将直接三角分解法修改为(部分)选主元的三角分解法. 66 67 于是有 引进