数值方法PPT课件.ppt

第三章 线性代数方程组的解法; 大量的科学与工程实际问题常常可以归结为求解含有多个未知量x1，x2，…，xn的线性代数方程组求解。;其矩阵表示形式为：; 常用的线性代数方程组解法可分为直接法和迭代法两种:;一．三角形方程组的解法 对于（3.2）所示的下三角方程组; 对于（3.3）所示的下三角方程组;二、 Gauss消去法;Gauss消去法具体过程：;Step1：;按上述做法，做完n-1步消元，原方程可化为同解的上三角方程组：;最后，若 ，逐步回代可得到原方程组的解：;Remark2：在消去过程中，也可以采用Gauss-Jordan消去法，将方程组化为对角形方程组，进而化为单位阵，则右端列向量就化为方程组的解向量。该方法不需回代过程，但总的计算量比Gauss顺序消去法略有增加。;Remark3：在消去过程中，消去过程能够进行的前提条件是 。当detA= 时方程组存在唯一解，但未必能满足 的???件。要使Gauss顺序消去法能够求得方程组的解，应满足如下的定理： ; 此时高斯顺序消去法能进行下去，但不能求出唯一解。; 在计算过程中舍入误差增大迅速，造成计算解与真解相差甚远，这一方法就是不稳定的方法，反之，在计算过程中的舍入误差增大能得到控制，该方法就是稳定的。小主元是不稳定的根源，这就需要采用“选主元素”技术，即选取绝对值最大的元素作为主元。;B= ;再考虑 右下角矩阵，选取绝对值最大的元素作为主元素,经过行的对换把主元素移到　 ， ;例题：用列主元消去法解线性方程组 ;彻摩危唇缆月铁擞州欺酣釜屯锨惰汁蹄禹携脆蜘祸弧街完眶椒息比蓟贯递数值方法PPT课件数值方法PPT课件;2)全主元消去法;都 ;Remark;一.Gauss顺序消去法的矩阵形式;锅甘镇显蛀锤狭耙雍告坡借伪依茹须腑皱瓶弧射史蔚焊惺卤花倔鞘疏八展数值方法PPT课件数值方法PPT课件;则n-1步消元过程为;即 ;令L= ; 即 ;二.矩阵的三角分解及条件 ;Remark：实际中对A进行三角分解，不是利用初等变换矩阵，而是直接使用矩阵乘法得到。（若不加说明，后面我们讲到的三角分解一律指Doolittle型分解。）;三.直接三角分解法 ;Step2： ;可以算出： ;的求解过程为:;对比计算 和 公式， ;康摇凹绰斟捡镇卞近靶蛤缄腺喜留霉唯宪桔巫珊药锡获嘶牲辖室萤溉镍驭数值方法PPT课件数值方法PPT课件;攻答睬夸干御胺办槛陌盒檬拦凸懒测倦彦由虾憎温铃霓睬池梨湘妻理双捡数值方法PPT课件数值方法PPT课件;四.平方根法（Cholesky分解） ;平方根法的计算公式 （为简单起见设 ，L为下三角矩阵） ;Remark1:由于在上式分解过程中有n次开方运算，故Choledsky分解法也称为平方根法。 ;Remark3:可以证明，若用Gauss顺序消去法求解对称正 定的方程组Ax＝b，则有;五、解三对角方程组的追赶法;2.追赶法解三对角方程组; 当A为上述三对角阵时，A有更特殊的三角分解形式;即 ;由 得 ;1．向量的范数;（1）;2．矩阵的范数;常用的矩阵范数：;迭代法的一般格式为;将该方程组等价变形为 构造简单迭代格 式 ， 。若 收敛于确定的向量 ，则 就是方程组的解。此时称简单迭代法 ， 关于初始向量 收敛。;③? 变形为 的方式不唯一。;2.简单迭代法的收敛性; (2)收敛的充分条件 ;相减有 ;因为 ，所以;3.Jacobi迭代法;建立迭代格式;雅可比迭代式成为：; (2) Jacobi迭代法的收敛条件;二、Seidel迭代法;与简单迭代法相应的Seidel??代可化为;的分量形式;有 上式对i=1,2,……，n都成立，;所以 。当 k 则;3.与Jacobi方法对应的Seidel迭代格式;因 ,故 存在，于是Seidel迭代式为; （2）收敛条件; 三、逐次超松弛迭代法 （SOR法-Successive Over Relaxation）;这就是逐次超松驰迭代法（SOR方法)， 称为松驰 因子。;2.SOR方法的收敛性;Remark：能使SOR方法收敛最快的松弛因子叫做最佳松弛因子，记为?opt。 对某些特殊类型