数值分析-线性方程组直接解法().pptVIP

第三 章 线性方程组 直接解法 （上） 第三章目录 线性方程组的概念 线性方程组的概念（续） 线性方程组的数值解法 §1 Gauss消元法 例1（续） Gauss消元法的基本步骤1（4阶） Gauss消元法的基本步骤2（4阶） Gauss消元法的基本步骤3（4阶） Gauss消元法的基本步骤4（4阶） Gauss消元法的消元过程1、2（n阶） Gauss消元法的消元过程3（n阶） Gauss消元法的消元过程3（n阶） Gauss消元法的回代过程（n阶） Gauss消元法的计算量 Gauss法与Cramer法则的计算量比较 §2 主元素法 例2（续1） 例2（续2） 例2两种解法的误差分析 2.2 列主元素法 列主元素法 2.3 全主元素法 主元素法举例 2.4 解三对角方程组的追赶法 追赶法的解题步骤 追赶法举例 §3 矩阵分解法 Gauss消元法的矩阵形式 Gauss消元法的矩阵形式（续） 杜利特尔（Doolittle）分解——LU分解 3.2 矩阵的三角分解 对A进行LU分解 对A进行LU分解 （r 行r 列） 对A进行LU分解的具体步骤 矩阵A的LU分解举例 三角分解的紧凑格式 三角分解的紧凑格式 三角分解的紧凑格式举例 3.3 直接三角分解法 直接三角分解法（续） 紧凑格式解线性方程组举例 紧凑格式解线性方程组举例（续） 三角分解法的几点说明 三角分解法的几点说明（续） §4 平方根法与改进的平方根法 4.1 平方根法 定理 2.2（续） Choleskg分解1 Choleskg分解2 平方根法举例 4.2 改进的平方根法 改进的平方根法说明 改进的平方根法举例 事实上，只要A非奇异，由上述结论，它一定可以分解成两个三角形矩阵的乘积，即：A=LU 。 消元过程相当于分解A = LU及求解三角形方程组Ly = b，回代过程则是求解另一个三角形方程组Ux = y，因此，解线性方程组问题可转化为矩阵的三角分解问题。 其中：L为单位下三角矩阵， U为上三角矩阵： 上述分解称为杜利特尔（Doolittle）分解，也称为LU分解，当系数矩阵完成三角分解后，对于求解方程组：Ax = b 。 正如Gauss消元法要在一定条件下才能进行到底一样，矩阵A也必须满足一定条件才能进行三角分解。 设A为n阶方阵，若A的顺序主子式Ai (i = 1,2,…,n ? 1)均不为零,则矩阵A存在唯一的Doolittle分解。 定理2.1 存在性证明 唯一性证明 下面讨论如何对A进行LU分解： （1行1列） 由于两个矩阵相等就是它们的对应元素都相等，因此通过比较A与LU的对应元素，即可得到直接计算L、U的元素的公式：A=LU 即： 紧接下屏： A L U 由矩阵乘法规则及比较（2-11）两端的元素，得： 即可由（2-12）求出U的第一行，L的第一列。 下面讨论一般情况，即：U的第i 行，L的第j 列： 一般情况下，可由： 计算过程应按U第1行、L第1列（第1框），U第2行、L第2列（第2框），……的顺序 具体分解步骤见下屏： 得到计算uij 和 lij 的公式： 1. 计算U的第1行，L的第1列，亦称为计算 第1框； 2. 计算U的第r 行，L的第r 列（r =2,…,n）， 即第r 框 ： 解：按分解公式（2-13）， 一框一框分解，每框计算时先行后列 ： 所以： 例5 根据式（2-13）的特点，矩阵的三角分解可按以下格式及顺序进行。这种格式既便于记忆，又便于计算，称为紧凑格式。见下表2-1： …… …… …… ? ? 表2-1 (a11)u11 (a12)u12 (a13)u13 (a1n)u1n …… (a21)l21 (a22)u22 (a23)u23 (a2n)u2n …… (an1)ln1 (an2)ln2 (an3)ln3 (ann)unn …… (a31)l31 (a32)l32 (a33)u33 (a3n)u3n …… … … … … … （1）计算顺序：将aij ，uij ，lij 按表2-1列好，计算 时按框从外到内进行， 每一框中先算行。从 左向右依次计算 uij ；再算列，自上而下求 lij ； （2）计算方法：按行计算时，需将所求元 uij 的 对应元aij 逐次减去 uij 所在行左面各框的元 素 lij 乘以 uij 所在列上面各框相应的元 uij 。 按列计算 lij 时，在作上述运算后还需除以 lij