第5篇 线性方程组的直接解法0.ppt

第五章 线性方程组的直接解法; 在科学计算中，经常需要求解含有n个未知量 的n个方程构成的线性方程组; 根据 Gramer（克莱姆）法则，求解方程组（5.1）时，要计算大量的行列式，所需乘法次数大约为; 直接方法的特点是: 如果不考虑计算过程中的舍入误差，运用此类方法经过有限次算术运算就能求出线性方程组的精确解。;假设系数行列式不为零，则解存在惟一，可用克莱姆法则求解;§1 Gauss消去法;一、顺序Gauss消去法;这样，对于方程组;顺序Gauss消去法的消元过程可表述如下： ;其中 ;其中 ;增广矩阵[A（n），b（n）]对应如下上三角形方程组 ;对于等价方程组;首先写出增广矩阵;然后进行消元，采用公式;在编程计算时，最后的增广矩阵存放的元素是：;消元过程： ;消元条件：;用Gauss消去法解方程组，应注意：;二 列主元Gauss消去法;此方程组具有四位有效数字的精确解为;经回代求解得 x3＝5.546，x2＝100.0，x1＝－104.0;写出原方程组的增广矩阵：;求得方程的解为：x3＝5.546，x2＝－45.76，x1＝17.46;下面将列主元Gauss消去法的计算步骤叙述如下：; 3. 再进行第二步消元得到增广矩阵[A（3），b（3）]。按此方法继续进行下去，经过n－1步选主元和消元运算，得到增广矩阵[A（n），b（n）]，它对应的方程组;列主元Gauss消去算法 ;5：消元计算，对于i＝k＋1，…，n，计算 ; 由于这两种方法的精度差不多，且全主元Gauss消去法程序设计复杂占用机器时间较多，实际应用中一般采用列主元Gauss消去法，它既简单又能保证计算精度。;§2 直接三角分解方法;若 ，令 ， i＝2，3，…，n，;若 ，令 ，i＝2，3，…，n，则有 ;如此下去，直至施行第n－1步消元，得到 ; 由此可见，在顺序Gauss消去法的过程中，系数矩阵A＝A（1）经过一系列单位下三角矩阵的左乘运算约化为上三角矩阵A（n），即 ;容易验证 ; 第一个方程组的系数矩阵为下三角矩阵，第二个方程组的系数矩阵为上三角矩阵，两个方程组都非常容易求解，具体求解结果如下：;对于;对于;可以看出对于方程组： ;二、Doolittle分解法;根据 A=LU 有等式成立：;先固定i按行求;于是，对于矩阵的三角分解：; 按框先横后竖计算元素;下面，我们对具体矩阵进行Doolittle 三角分解。;例5-3 利用Doolittle三角分解法分解矩阵;1;如果我们要求解方程组;由;（5.5）;例5-4 利用Doolittle三角分解方法解线性方程;得到;例5-5 利用Doolittle三角分解方法解线性方程组;1 2 3 -4 -2;三、平方根法; 记 Ak（1≤k≤n）为A的 k 阶顺序主子阵，则det（Ak）为A的k阶顺序主子式。由上式，利用矩阵分块运算规则，容易验证 ;即： A=LDM ，其中 DM=U， M=D-1U。;当A＝AT为对称矩阵时，根据 A＝LDM, 得到; 如果对称正定矩阵A具有如下分解 A＝GGT，其中G为下三角矩阵，则称其为对称正定矩阵的 Cholesky（乔列斯基）分解。 ;Ly＝b LTx＝ y;比较对应元素：;当 i = j 时;于是，根据计算公式; 关于方程组 Ax=b , 如果对系数矩阵进行了平方根分解 A＝LLT，则将方程组化为： Ly＝b , LTx＝y; 于是，关于系数矩阵是对称正定矩阵的线性方程组Ax=b的求解，分两步进行：;例5-6 用平方根法求解对称正定方程组 ;得 y1＝2，y2＝3.5 ，y3＝1 ;关于对称正定方程组 ;四、追赶法;在此条件下，可对A进行三角分解，设;i-1列;i-1列;i-1列;于是，由以上结果：;对于三对角方程组Ax＝b，设A的三???分解为A＝LU，则原方程组等价于;由Ly＝b，即;于是，对于三对角矩阵方程组 Ax=b，如下的两组公式便构成了构成了解三对角方程组的追赶法：;例5-7 用追赶法求解三对角方程组 ;Date;求解方程组Ly＝y，即 ; 当三对角矩阵A满足对角占优条件时，追赶法是数值稳定的。 ;1． ;五、三角分解方法的优点;2. 可以用以求可逆矩阵 A 的逆矩阵 A－1;3. 可以用以求矩阵 A 的行列式