总结第3章 节线性方程组的数值解法.ppt

第3章 线性代数计算方法 ;§1 高斯消去法 ; 一、顺序消去法 解方程组 ; 作②-①消去②中的x1,作③-①×4消去③中的x1,则方程组(3―6)化为 ; 从方程组(3―6“)的方程③解出x3,将所得的结果代入方程②求出x2,再把x3、x2同时代入方程①解出x1。这样可求出方程组的解为 ; 2、顺序高斯消去法的计算步骤： 在计算机上实现时,我们常把方程组右端的常数项排于系数矩阵的第n+1列， 1）消元过程 对于k=1,2,…,n-1列,若按顺序有某一ark≠0,r≥k,则交换k与r行,然后计算 ;3、计算量;二、主元素消去法 1.列主元素消去法 所谓列主元素消去法就是在每一步消元过程中取系数子矩阵的第一列元素中绝对值最大者作主元。; 取四位有效数字计算。 解：第一列消元时，②中-18为主元,交换②和①得 ; ②+①×12/18,③+①×1/18得 ; ③+②×1/1.167得 ;列主元素消去法的计算过程:;图 3.1 ;图 3.1 ; 2. 全主元素消去法 所谓全主元素消去法,就是每步消元时选取系数子矩阵中绝对值最大的元素作主元。 说明：若矩阵第i列与第j列对调,则未知数xi与xj也相应地对调了,xi的结果实质上为xj的结果。 ; 例2 用全主元素消去法求解方程组 ; 再全选主元,主元为2.333,交换x2和x3所在的两列,同时改变两未知数的排列号得 ;全主元素消去法的计算过程: (1)消元过程。对k=1,2,…,n-1进行下列运算: ; 图 3.2 ; 图 3.2 ;§2 高斯―约当消去法 ;§3 解实三对角线性方程组的追赶法 ; ;§4 矩阵的三角分解 ?;1.克劳特分解方法; 这样,L、U中的元素都已求出。计算L的各列与U的各行的次序如下图所示 。 ;克劳特分解求解线性方程组的计算过程: ① LU分解过程:对于k=1,2,…,n依次计算 ; 图 3.5 ; 图 3.5 ; 例4 用克劳特分解方法求解下列方程组 ; 利用矩阵乘法可得到 ; 这样原方程组就化为依次求下列两个三角形方程组 ; 2. 杜利特尔分解方法 杜利特尔分解方法是令A=LU 这里 ; 例5 用杜利特尔分解方法求解下列方程组的解: ;§5 行列式和逆矩阵的计算 ;2. LU分解法 ①克劳特分解法: detA=detL·detU=detL=l11l22…lnn ②杜利特尔分解法: detA=detL·detU=detU=u11u22…unn ③乔累斯基分解法: 这里A为正定矩阵。 detA=detL·detLT=l211l222…l2nn或者 detA=detL·detD·detLT=detD=d1d2…dn;5.2逆矩阵的计算 ; 例7 设 ; 先对第一列消元,将a11化为1,a21、a31化为0,有 ; 将a33化为1,并将a23、a13化为0,得 ;§6 迭代法 ;例如将A分解为两个矩阵之差 A=B-C 其中矩阵B可逆, 于是方程组成为 Ax=Bx-Cx=b Bx=Cx+b x=B-1 Cx+B-1b 令M=B-1 C, f=B-1 b 则x（k+1）=Mx（k）+f ;雅可比（Jacobi）迭代法 取 特别当M的对角元素均不为零且按绝对值来说较大时, ; 例8 用简单迭代法解下列方程组 ; 取初始值x(0)= 0,按迭代公式 ; 图 3.8 ; 6.2 赛德尔(Seidel)迭代法 设方程组的等价形式为 x=Mx+f 将矩阵M分解为M=B+C， 其中 ;于是 x=Bx+Cx+f x(k+1)=Bx(k+1)+Cx(k)+f, k=0,1,2,… 即： 故;