数值分析——线性方程组直接解法.ppt

第五章 线性方程组的直接解法 第五章 线性方程组的直接解法 §1 Gauss消去法 1.1 顺序Gauss消去法 1.2 列主元Gauss消去法 §2 直接三角分解方法? 2.1 Gauss消去法的矩阵运算 2.2 Doolittle分解法 2.3 平方根法 2.4 追赶法 §1 Gauss消去法 一、顺序Gauss消去法 二 列主元Gauss消去法 §2 直接三角分解方法 一 、 Gauss消去法的矩阵运算 二、Doolittle分解法 例7-3 利用Doolittle三角分解法分解矩阵 例7-4 利用Doolittle三角分解方法解线性方程 例7-5 利用Doolittle三角分解方法解线性方程组 三、平方根法 四、追赶法 五、三角分解方法的优点 第七章 总 结 关于方程组 Ax=b , 如果对系数矩阵进行了平方根分解 A＝LLT，则将方程组化为： Ly＝b , LTx＝y 解得 于是，关于系数矩阵是对称正定矩阵的线性方程组Ax=b的求解，分两步进行： 第一步：系数矩阵的平方根分解 第二步：解等价方程组 例7-6 用平方根法求解对称正定方程组 解: 首先进行A 的Cholesky 分解 A＝LLT 2 -0.5 0.5 2 1.5 1 2 -0.5 0.5 2 1.5 1 得 y1＝2，y2＝3.5 ，y3＝1 得 x1＝1，x2＝1，x3＝1 求解Ly＝b: 再求解 LTx＝y: 关于对称正定方程组 也可以用一般的三角分解法求解，这时由 4 -1 1 4 -0.25 0.25 4 3 7 0.75 1 1 求得 x1＝1，x2＝1，x3＝1 追赶法是专门用于求解三对角方程组的。这类方程组经常出现于用差分方法或有限元方法求解二阶常微分方程边值问题、热传导问题及三次样条函数插值等问题，三对角方程组Ax＝b的系数矩阵具有如下形式： 设A为一个三对角矩阵，那么它的顺序主子式均不为零的一个充分条件是： 在此条件下，可对A进行三角分解，设 比较矩阵的对应元素，根据矩阵乘法规则，可得到 i-1列 i 行 第3列 若 ，令 ，i＝2，3，…，n，则有 施行第二步消元，我们得到 如此下去，施行第n－1步消元，得到 由此可见，在顺序Gauss消去法的过程中，系数矩阵A＝A（1）经过一系列单位下三角矩阵的左乘运算约化为上三角矩阵A（n），即 这时 由 得 令 容易验证 则从顺序Gauss消去法的矩阵运算表示式可知，系数矩阵A可分解为一个单位下三角矩阵L和一个上三角矩阵U的乘积，即 其中 第一个方程组的系数矩阵为下三角矩阵，第二个方程组的系数矩阵为上三角矩阵，两个方程组都非常容易求解，具体求解结果如下： 我们将A=LU 称为矩阵A的三角分解,这时线性方程组为： 令 则有 对于 由 解得 对于 由 求得 可以看出对于方程组： 只要对系数矩阵作了三角分解： 由这个简单的计算过程可知，系数矩阵的三角分解很关键，如何进行三角分解更容易？下面介绍几种方法。 通过如下两组公式很容易求解： 前已述及，若在顺序Gauss消去法的过程中，每步消元的主元素 akk(k)≠0，则矩阵A可分解为A＝LU，L为单位下三角矩阵，U为上三角矩阵，此分解称为A的 Doolittle（杜利特尔）分解。可以证明akk(k) ≠0的充要条件是A的各阶顺序主子式不为零，于是有如下定理。 定理7.1 设n阶方阵A的各阶顺序主子式不为零，则存在惟一单位下三角矩阵L和上三角矩阵U使A＝LU。 下面介绍矩阵三角分解的Doolittle分解方法。 根据 A=LU 有等式成立： 比较等式两端对应元素，有 n n i i 可以解得： 当 i=1 时 当 j=1 时 当 i 1 时 当 j 1 时 于是，对于矩阵的三角分解： 可按照以下公式进行： 对于 i=2,3, … ，n, 计算 （7.2） （7.3） （7.4） 用计算公式(7.3)、(7.4)对矩阵A作的分解(7.2),称作Doolittle分解。 下面，我们对具体矩阵进行Doolittle 三角分解。 为了表示和存储方便，可以将分解后的两个矩阵用一 个矩阵表示 解：分解时用到如下公式 1 2 3 4 1 1 1 2 6 12 3 7 6 24 6 24 1 2 3 4 1 1 1 2 6 12 3 7 6 24 6 24 可以写成：