化工数值计算-5.pptVIP

? 线性代数方程组的数值解;化学化工中,用线性代数方法描述复杂的化学反应、定义多原子和多分子体系的线性空间、计算各种物质的物理化学性质的加和方法、物料衡算和热量衡算、确定反应体系的独立反应数、化工量纲分析中确定独立变量数以及某些分析化学计算等,都常常要运用线性代数方程组的计算。;线性代数方程组的解法大体上可分为两大类,即迭代法与直接法。所谓迭代法是将方程组写成迭代形式,由给定的一组初值,代入迭代方程组,通过多次迭代来得到一个可逼近方程组真实解的近似解。它的优点是计算简单,编制程序较容易;其缺点是要求方程组的系数矩阵具有某种特殊性质(如具有主对角线占优),否则不收敛或收敛;很慢,使计算量很大,失去使用意义。 直接法又称精确法,由于运算过程中有舍入误差,实际上所求得的解是近似解。它的优点是计算量较小,在机器的内存容量等许可的情况下,求解线性方程组以采用直接法为多。 5.1????解线性代数方程组的迭代法;? 　　设有线性代数方程组 ?;其矩阵形式为: ? 式中,A = (aij)n×n为系数矩阵;X = (x1,x2,…,xn)T为n维向量;B = (b1,b2,…,bn)T为等式右端常数向量。现讨论其迭代法求解方法。;5.1.1　简单迭代法 对于线性方程组AX = B,即 ? 若A为非奇异矩阵,且aii ≠ 0(i=1,2,…,n),则可将其转化为下列等价方程组:;? 若对式(5-3)实施迭代法,则其迭代格式为 ? 令cij=?,di=?,则上式变为;? 这种迭代法又称为雅可比(Jacobi)迭代法。如果 ? 存在,则称迭代式是收敛的,这时极限值?(i=1,2,…,n)就是线 性方程组的解。下面来讨论其迭代收敛的条件。;方程组的精确解?与第k+1次迭代近似值?之差为 ? k+1次迭代近似值的误差为;? 于是由式(5-6)有 ? 式中;? 由式(5-7)可得 ? 由此可见,如果L1,则k→∞,必有ek+1→0,亦即有 ?;这说明对任意给定初值?(i=1,2,…,n),迭代式(5-5)收敛。 由以上分析可知, ? 式(5-11)是迭代式(5-5)收敛的充分条件。此即要求各行应有;? 或 ? 这样的方程组系数矩阵就是严格对角占优矩阵。;对于迭代收敛的情形,可根据精度要求ε,用条件 ? 来控制迭代过程的结束。 5.1.2　赛德尔迭代法 5.1.2.1　方法概述;在迭代式(5-5)的迭代计算中,第k次迭代用的是第k-1次迭代的近似值?。然而,实际计算?时,它前面的其他未知数的本次迭代 近似值已经计算出来了。;一般来说,新值总比旧值更接近真实值,因此应该优先使用它们,这可能加速迭代过程的收敛。这就是赛德尔(Seidel)迭代法的基本思想。可见,赛德尔迭代法是简单迭代法的一种改进。 对于线性方程组AX = B,由式(5-5)可得赛德尔迭代法的迭代格式为 ? 式中;? 可以证明,迭代式(5-14)收敛的充分条件仍然是式(5-12)。一般情况下,赛德尔迭代法比简单迭代法的收敛速度快,但是也有相反的情况。 5.1.2.2　程序框图 图5-1为赛德尔迭代法的通用计算程序框图,其中主要变量如下:;N　　方程组中变量的个数,也即系数矩阵的维数 E　　给定精度 A　　二维数组,存放线性方程组的系数矩阵(N×N) B　　一维数组,存放线性方程组等号右端的常数 X　　一维数组,存放线性方程组解的初始值或迭代值 R　　迭代最大误差;5.1.2.3　计算实例;　　解　由方程组的系数可以看出,该方程组为严格对角占优矩阵,不必进行转化就可采用赛德尔迭代法按图5-1编写源程序。 源程序: ’******************************************** ’　　Example 5-1------Eg5-1.frm ’******************************************** DefDbl A-H,O-Z Private Sub Command1_Click() Dim X(50),A(50,50),B(50) Dim XA As Variant,AB As Variant Cls N = 3: E = 0.00001: ID = 100 XA = Array(0,0,0) For J = 1 To N: X(J) = XA(J - 1): Next J AB = Array(10,-2,-1,3,-2,10,-1,15,-1,-2,5,10) K = 0 　　For I = 1 To N 　　　　For J = 1 To N 　　　　　　A(I,J) = AB(K): K = K + 1 　　　　 Next J 　　　　　　B(I) = AB(K): K = K + 1