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

图 3.5 例4 用克劳特分解方法求解下列方程组 解： 令 利用矩阵乘法可得到 这样原方程组就化为依次求下列两个三角形方程组 代入第二个方程组可求得原方程组的解为 2. 杜利特尔分解方法 杜利特尔分解方法是令A=LU 这里 例5 用杜利特尔分解方法求解下列方程组的解: 解： 利用式(3―50)、(3―51)可求得 求得原方程组的解为 §5 行列式和逆矩阵的计算 5.1 行列式的计算 1. 高斯消去法 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 设 求A-1。 解:将方程组写成同一个增广矩阵,即 先对第一列消元,将a11化为1,a21、a31化为0,有 再将第二列a32化为0,得 将a33化为1,并将a23、a13化为0,得 最后将a13化为0,得 所以 第3章 线性方程组的数值解法 第3章 线性方程组的数值解法 《计算方法与实习》 第3章 线性代数计算方法 §1 高斯消去法 §2 高斯―约当消去法 §3 解实三对角线性方程组的追赶法 §4 矩阵的三角分解 §5 行列式和逆矩阵的计算 §6 迭代法 ? §7 迭代法的收敛性 解线性方程组的两类方法: 直接法: 经过有限次运算后可求得方程组精确解的方法(不计舍入误差!) 迭代法：从解的某个近似值出发，通过构造一个无穷序列去逼近精确解的方法(一般有限步内得不到精确解) §1 高斯消去法 高斯消去法是一个古老的直接法，由它改进的方法是目前计算机上常用的求低阶稠密矩阵方程组的有效方法。 思路 首先将方程组Ax=b 化为上三角方程组，此过程称为消去过程，再求解上三角方程组，此过程称为回代过程. 一、顺序消去法 解方程组 ① ② ③ (3―6) §1 高斯消去法 作②-①消去②中的x1,作③-①×4消去③中的x1,则方程组(3―6)化为 ① ② ③ (3―6′) 对方程组(3―6′)作③-②× ,得到三 角形方程组 ① ② ③ (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―11) 2）回代过程 对于k=n,n-1,…,2,1,计算 (3―12) 二、主元素消去法 1.列主元素消去法 所谓列主元素消去法就是在每一步消元过程中取系数子矩阵的第一列元素中绝对值最大者作主元。 取四位有效数字计算。 解：第一列消元时，②中-18为主元,交换②和①得 ① ② ③ ① ② ③ 例1 用列主元素消去法解方程组 ②+①×12/18,③+①×1/18得 ① ② ③ 第二列消元时,主元为1.167,交换方程②和③得 ① ② ③ ③+②×1/1.167得 ① ② ③ 回代求得 x1=1.000, x2=2.000, x3=3.001 方程组的实际解 x1=1, x2=2, x3=3 列主元素消去法的计算过程: (1)消元过程。对于k=1,2,…,n-1进行下述运算: ①选主元,确定r,使得 若ark=0,说明系数矩阵为奇异,则停止计算；否则进行下一步。 ②交换A的r、k两行 ③对i=