计算方法第6章矩阵及线性方程组.ppt

定理　设有三对角线方程组Ax=f，且A满足上面条件，则A为非奇异矩阵，且A的所有顺序主子式都不为零。 即 将A进行三角分解，A=LU 通过矩阵的乘法可得计算待定系数{ }{ }{ }的计算公式，即 (2) (1) (3) 得到解方程Ax=f的追赶法公式 （1）分解计算公式A＝LU （2）求解Ly=f的递推公式 求解Ax =f?求解 （a）Ly =f求y (b) Ux = y求x 上一页 下一页 返回首页 6.1 线性方程组的直接解法 第6章 矩阵及线性方程组 6.2矩阵的三角分解 6.3解三角方程组的追赶法 6.4矩阵求逆 6.1 线性方程组的直接解法 6.1.1 直接法概述 设有线线方程组Ax = b 其中 当det A 0 时，方程组的解存在且唯一，对增广矩阵（A，b）施行行初等变换，化A为上三角形矩阵 ，同时b化为 即: 方法1 对应的线性方程组 消元过程 从最后一个方程解得 回代过程 方法2 是直接将方程组Ax = b的系数矩阵A分解为两个形式简单的三角形矩阵L和U的乘积，即A = LU。 即: 从而求解A x = b的问题转化为解三角形方程组Ly = b和U x = y 。 6.1.2三角形线性方程组的解法 则Lx =b为下三角形方程组，它的第i个方程为 假定 ，按 的顺序解得： 上三角形方程Ux =b的第i个方程为 假定 ，按 的顺序求解得： 6.1.3 Gauss消去法 1.消元与回代计算 设det A 0,记Ax = b的增广阵 假定 0，第一步消元是从方程组的第2至n个方程中消去未知数 0 0 为此，作如下变换 第i行元素减去第一行对应元素乘以 ，使 第一列主对角元以下元素为零，同时有 记 整个过程可描述如下： Step 1：设 ，计算因子 将增广矩阵 第 i 行 ? mi1 ? 第1行，得到 其中 Step k：设 ，计算因子 且计算 共进行 ？ 步 n ? 1 回代过程 如果矩阵A非奇异，总可以通过带有行交换或不带行交换的消元过程，将A化为非奇异上三角形矩阵 ，因此，回代求解过程也可以进行到底。 2.Gauss消去法的运算量 第k步消元计算中，做乘法 次，除法 次，完成 步消元需做乘除法的总数为 次。 回代过程共做乘除法的次数为 于是Gauss消去法中乘除法总数为 6.1.4 选主元素的Gauss消去法 用Gauss消去法解Ax = b时，其中设A为非奇异矩阵，可能出现 = 0情况下，这时必须进行行带交换的Gauss消去法。但在实际计算中即使 但其绝对值很小时，用 作除数，会导致中间矩阵 元素数量级严重增长和舍入误差的扩散，使得最后的计算结果不可靠。 例 设有方程组 解 精确解为 [方法1] 用Gauss消去法求解（用具有舍入的3位浮点数进行计算）。 回代得到计算解 比较 ╳ [方法2]　用具有行交换的高斯消去法（避免小主元）。 方法1计算失败的原因，是用了一个绝对值很小的数作除数，乘数很大，引起约化中间结果数量级严重增长，再舍入就使得计算结果不可靠了。 回代求得 比较 √ 这个例子告诉我们，在采用高斯消去法解方程组时，小主元可能导致计算失败，故在消去法中应避免采用绝对值很小的主元素。 通常在元过程中进行选主元，即在高斯消去法的每一步应该在系数矩阵或消元后的低阶矩阵中选取绝对值最大的元素作为主元素，保持乘数 ≤1，以便减少计算过程中舍入误差对计算解的影响。 1．完全主元素消去法 线性方程Ax=b，其中A为非奇异矩阵,增广矩阵为: 第k步选主元区域 第k步计算可描述如下: Step k: ① 选取 ② If ik ? k then 交换第 k 行与第 ik 行; If jk ? k then 交换第 k 列与第 jk 列; ③ 消元 经过上面的过程，即从第1步到 步完成选主元，交换两行，交换两列，消元计算，原方程组约化为 注：列交换改变了 xi 的顺序，须记录交换次序，解完后再换回来。 回代求解 这里的 为未知数 调换后的次序。 2．列元素消去法 完全主元素方法在选主元时要花费一定的计