南京工业大学 计算机在材料中应用 3.线性方程组数值解法.pptVIP

3 线性方程组数值解法; ?a11 a12 ……a1n ? ?X1? ?b1? A=?a21 a22 ……a2n? X= ?X2? B= ?b2? ? …… ? ?... ? ?…? ?an1 an2 ……ann? ?Xn? ?bn? A称为方程组的系数矩阵，当是n阶非奇异矩阵时，既?A ??0，此时方程组有唯一解。 X阵是解向量，B阵是常向量。 在线性代数中学过用克莱姆法则求解，它是一种直接法（属于解析法），但随着n??计算工作量?; ?直接法：经有限次计算得准确解（在无舍 ?入误差下），实际上舍入误差客观存 ?在，得到的依然还是近似解。方法有主元 |消去法、三角??解法。计算量小、精度较 |高，但程序设计复杂。 解线性方程组? 的数值解法 | ?迭代法：将问题构成一个无穷序列，逼近 |准确解。方法有雅可比法、高斯-赛得尔法 | 算法简单、易于编程，但存在收敛性及收敛 |速度的 |问题，只对具有某些性质的系数矩阵的方程 ? 组才适用。;3.2 消去法; ?2 4 -2? 设A（1）=?1 -1 5? =（aij(1)） =(aij)=A (i=1,2,3, j=1,2,3) ?4 1 -2? ?6? B(1)=?0?=(bi) (i=1,2,3) ?2? 令li1=ai1(1)/a11(1) (i=2,3)，称li1为消元因子，用li1消去除第一个方程外其余各方程的x1，即将第一个方程乘以- li1后加到第i个方程上去（i=2，3），则有; ?2 4 -2 ? ?x1? ? 6 ? ?0 -3 6 ? ?x2?= ? - 3 ? ?0 -7 2 ? ?x3? ?-10? 此时记为A（2）X=B(2) 其中 aij(2)=aij(1)-li1a1j(1) (i=j=2,3) bi(2)=bi(1)-li1b1(1) 令li2=ai2(2)/a22(2) (i=3)，将第二个方程乘以-li2后加到第i个方程上去 （i=3） 则有 ?2 4 -2? ?x1? ? 6 ? ?0 -3 6? ?x2?= ? - 3 ? ?0 0 -12? ?x3? ? 3 ? ;记为A（3）X=B(3) 其中 aij(3)=aij(2)-li2a2j(2) (i=j=3) bi(3)=bi(2)-li2b2(2) (i=3) 现在A（3）已成上三角阵，上述过程称之为消元过程，下面进行回代过程，所以从第三个方程可解出x3后，回代入第二个方程得x2，再将x2、 x3回代到第一个方程解出x1 ，解分别为0.25 、1.5 、 0.25。 在上述过程中要求aii?0，若aii=0，则可通过行、列变换A来满足aii?0， aii称为主元素。 ?aii ?越小，作除数引起的误差就越大，所以尽可能选用max ? aii ?，这样就有了全主元、列主元消去法。;现讨论一般n阶方程组 AX=B 记A（1）X=B(1) ，其中A(1)=(aij(1))=（aij） =A B(1)