第二章线性联立方程组技术方案.pptx

第二章 线性联立方程组 2 第二章 线性联立方程组 1、数值分析求解方法 内容： 2、应用实例 一、引言 二、线性方程组的求解方法 1）克莱姆规则 2）迭代法 3）消去法 3 第二章 线性联立方程组——引言 1、引言 一些化学问题要求解线性联立方程组，这包括： 混合物分析 多组分平衡分析 化学反应配平 非线性回归分析 曲线修匀 量子化学计算 参数的估算 等 4 第二章 线性联立方程组——引例 2、引例 配平化学反应： NaOH+H2SO4→Na2SO4+H2O 该反应的化学计量系数为n1, n2, n3, n4 即 n1NaOH+n2H2SO4→n3Na2SO4+n4H2O 配平反应的依据是各原子数反应前后相等。 元素 Na n1-2n3=0 O n1+4n2-4n3-n4=0 H n1+2n2-2n4=0 S n2-n3=0 因此问题转化成了求解上述线性方程组得问题。 5 3、线性方程组的一般形式为： a11x1+a12x2+…+a1nxn=b1 a21x1+a22x2+…+a2nxn=b2 …… am1x1+am2x2+…+amnxn=bm 第二章 线性联立方程组—一般问题 6 4、克莱姆规则 该规则为： 其中 Di为A得第i列被B取代后得系数矩阵对应的行列式。 第二章 克莱姆规则 当n比较小时（如n=2或n=3), 很容易用克莱姆规则求解，但当n较大时，由于求解行列式需要计算n!项，因此通常不采用。 7 5、 Gauss-Seidel 迭代法 1）方法原理 对于AX=B, 给定一组X的初值X(0),计算的X(1), 继续下去，直至第m次，若X(m) 满足方程条件，则X(m)可视为方程的解。 第二章 Gauss-Seidel 迭代法原理 例如：三元联立方程组 a11x1+a12x2+a13x3=b1 a21x1+a22x2+a23x3=b2 a31x1+a32x2+a33x3=b3 若a11, a22, a33≠0, 则可给定初值x1(0)=x2(0)=x3(0)=0进行迭代： 8 第一次迭代： 第二章 Gauss-Seidel 迭代法 第二次迭代： 9 迭代下去直至第m次： 第二章 Gauss-Seidel 迭代法 若已满足规定的收敛条件，则x1(m), x2(m), x3(m)为方程组的解. 10 第二章 11 实行Gauss-Seidel 方法的条件是： 系数矩阵A要非奇异。 A中得对角元不能等于 0。 矩阵A得对角元是严格占优势的。 即 第二章 Gauss-Seidel 迭代法的条件 12 i. ii. iii. m= max (最大迭代次数) Gauss-Seidel 方法的计算步骤： 给定迭代计算的初值xi(0), (I=1, 2, … , n); 计算新的xi值： 第二章 Gauss-Seidel 迭代法计算步骤 判断收敛条件： 13 Gauss-Seidel 方法的计算步骤： 给定迭代计算的初值xi(0), (I=1, 2, … , n); 计算新的xi值： 第二章 14 i. ii. iii. m= max (最大迭代次数) 第二章 判断收敛条件： 15 过程开始 c=eps c1=a[i,j]*x[j]+c1 c=0 c1=0; c2=0 END 打印结果 kitmax i=1 to n j=1 to i-1 j=i+1 to n c2=a[i,j]*x[j]+c2 x1[i]=(b[i]-c1-c2)/a[i,i] i=1 to n c=abs((x1[i]-x0[i])/x1[i])+c cE x0[i]=x1[i] i=1 to n Gauss-Seidel 迭代法流程图 y y N N 程序 16 主程序开始 END 定义变量 a[i,j];b[i] x0[i] 调用过程 Gauss-Seidel 迭代法流程图 应用程序实例 17 6、消去法： 原理 消去法涉及对增广矩阵得以系列运算。增广矩阵（A|B)为n× (n+1)矩阵： a11 a12 … a1n b1 a21 a22 … a2n b2 … am1 am2 … amn bn 第二章 消去法原理 18 第二章 消去法原理 消去法计算过程 1、消元过程 19 第二章 消去法原理 消去法计算过程 20 第二章 消去法原理 消去法计算过程 21 第二章 消去法原理 消去法计算过程 2、回代过程 a11 a12 …