线性方程组直接解法在化学中的应用.docVIP

线性方程组直接解法在化学中的应用 一、前言 线性方程组直接法在化学化工中应用很普遍，线性方程组直接法可分为Gauss消元法和直接三角分解法（又可称LU分解法），详细介绍并对比这两种方法。结合其在化学数值求解方面的应用，求解部分化工实际问题。 在化学工程中需要求解线性方程组的场合很多，如多组分体系的物料衡算，计算各种化合物的物理化学性质以及稳态动力学计算等。在化工中涉及到的非线性方程及非线性方程组的问题，如各种形式的真实气体状态方程，多组分混合溶液的沸点，饱和蒸汽压计算，流体在管道中阻力的计算，多组分多平衡级分离操作模拟计算，多效蒸发过程设计，化工流程最优化计算等。事实上这些问题很多数值法最终也都可以归结为解线性方程组问题，所以说线性代数计算方法是数值计算的重要基础。主要讨论高斯消去法和直接三角分解法在化工计算中的优劣。 二、算法简述 Gauss消去法主要是消元和回代两个过程组成，可以分为顺序Gauss消去法和列主元素Gauss消去法。顺序Gauss消去法，其算法如下： 记aij(1)= aij (i,j=1,2,…,n), bi(1)= bi (i=1,2,…,n) 消元过程 对于k=1,2,…,n-1执行 （1）如果akk(k)=0,则算法失效，停止计算；否则转（2）。 （2）对于i=k+1,k+2,…,n计算 mik= aik(k)/ akk(k) aij(k+1)= aij(k)- mik akj(k), j=k+1,k+2,…,n bi(k+1)= bi(k)- mik bk(k) 回代过程 ，k=n-1,n-2,…1 列主元素Gauss消去法，其算法如下： 记aij(1)= aij (i,j=1,2,…,n), bi(1)= bi (i=1,2,…,n) 消元过程 对于k=1,2,…,n-1执行 （1）选行号ik,使 （2）交换与（j=k,k+1,…,n）与所含的数值。 （3）对于i=k+1,k+2,…,n计算 mik= aik(k)/ akk(k) aij(k+1)= aij(k)- mik akj(k), j=k+1,k+2,…,n bi(k+1)= bi(k)- mik bk(k) 回代过程 ，k=n-1,n-2,…1 这两种方法计算量相同但是后者精度显著要高于前者。 直接三角分解法又可称LU分解法。如果方程组系数矩阵A能分解成 A=LU 根据分解式中L是单位下三角矩阵，U是上三角矩阵称为Doolittle分解；如果L是下三角矩阵，U是单位上三角矩阵称为Crout分解。可以证明LU分解法与高斯消去法是等价的，计算工作量也基本相同。但LU分解法有其特点，首先节省计算机的存储空间，LU分解所得元素可存入A矩阵元素各存贮单元内；如果用LU分解法求解具有相同系数矩阵的多个方程组，则可节省大量机时和内存空间。 三、上述算法的化学运用 1、Gauss消去法的运用 假设物系服从Beer定律，试确定下列混合物四种组分的浓度。设光程长度为1cm，观测数据如下表。 摩尔吸收率/（L/mol·㎝） 波长 对二甲苯 间二甲苯 邻二甲苯 乙苯 总吸收率 12.5 1.502 0.0514 0 0.0408 0.1013 13.0 0.0261 1.1516 0 0.0820 0.09943 13.4 0.0343 0.0355 2.532 0.2933 0.2194 14.3 0.0340 0.0684 0 0.03470 0.0339 解：因为服从Beer定律，则有Atoti= 其中 Atoti—波长为 时观测到的总吸收率； —波长为 时第j个组分的摩尔吸收率; —混合物中第j组分的摩尔浓度。 根据题意可建立线性方程组，写成矩阵形式如下 = 经第一次消元后，得到 = 第二次消元后，得 = 对于此题仅经两次消元即变为三角形方程组，然后回代可解出 =0.026/0.3413=0.0762(mol/L) =(0.2142-0.2899×0.0762)/2.532=0.0759(mol/L) =(0.0977-0.0813×0.0762)/1.1507=0.0795(mol/L) =(0.1013-0.0514×0.0795-0.0408×0.0762)/1.502=0.0627(mol/L) 2、LU分解法的运用。 乙炔的摩尔热容与温度的经验关联式，用最小二乘法拟合实测数据得到如下正规方程组 105.21=8a+28b+140c 402.29=28a+140b+784c 2070.29=140a+784b+4676c 试用LU分解法求解方程组，确定参数a,b,c 解：首先将方程组的系数矩