化工系统第2章讲解.ppt

* * 4. 求雅可比矩阵的逆矩阵时工作量大 在求解过程中，雅可比矩阵有可能是奇异矩阵。两种原因： 建模型时出错，此时需对模型进行检查修正； 迭代过程中出现的（相当于单变量情况下 时的情形），此时可重设初始值。 * * 5. 解析求导困难时，可用下式进行数值求导： Δxi 为差分步长，xi 一般可取为0.001 * * 例：对串联的油换热器组进行最优设计时，得到如下方程组 求油换热器进出口的温度T1 和 T2，已知初始值是（180，292）。要求精度为0.01 解：（1）写出 即： * * （2）计算 （3）求出 J * * （4）计算 （5）求逆 （6）迭代 （7）重复(2)至(6)步，直到满足计算精度要求 计算 若：J 的行列式为非奇异（满秩）方阵，即 则 J 的逆矩阵为： * * J 的伴随矩阵 矩阵求逆 求 J*：先求 中元素的代数余子式： 则 * * 三、布洛伊顿（Broyden）拟牛顿法 牛顿-拉夫森法优点：有一定理论基础，收敛速度较快，效果较好。 缺点：迭代公式中含有一个待解方程组中函数向量 的雅可比矩阵的逆矩阵 雅可比矩阵的求取，是同方程组数学形式直接有关的运算，当方程组的规模较大时，相当麻烦。 用解析求导，工作量很大，有时还有困难；用数值求导也并不简单。 * * 比照着牛顿-拉夫森法的迭代公式把迭代公式写成如下形式： （2-31） 如果 ，（2-31）式就是牛顿-拉夫森迭代公式。如果不按照牛顿-拉夫森方法取值，则这样的求解方法就叫做拟牛顿法（Quasi-Newton Method）。 对于迭代矩阵的具体构成，可以提出各种不同的方案，这样就可形成不同的拟牛顿法。 * * 各种拟牛顿法中的迭代矩阵，一般并不与待解方程组的数学形式直接挂钩。因此在应用拟牛顿法时，只需逐轮进行各方程的函数值的计算，就可使迭代进行下去，方便，适用范围更广。 布洛伊顿法是应用得最广泛的拟牛顿法。在某些文献中，在称拟牛顿法时，指的就是布洛伊顿法。 * * 自学布洛伊顿法，例2-4 作业：P 250，习题 4 迭代法总结： 1.确定迭代方法 2.转换方程形式 3.写出迭代通式 4.列表计算 5.终止迭代 建议-多借助软件 Matlab Mathematic Maple 四、方程组的分块和切割1.稀疏方程组 一个方程组共含有n个变量和n个方程，但这并不意味着，n个方程中的每一个，都必含有n个变量的全体。很可能是某些方程中只含有n个变量中的某几个，另一些方程只含有n个变量中的另外几个。这叫做方程组具有“稀疏性(Sparseness)”。这样的方程组也被称为稀疏方程组。 化工系统的模型方程，通常都是稀疏方程组。而且稀疏程度一般很高。 * * 方程组的分块和变量切割，就是针对求解稀疏方程组的一种方法。 方程组的分块和变量的切割称为方程组的分解。它与第六章要介绍的化工系统的分解，本质上是一致的，做法也类似。 仅介绍方程组分块和变量切割的基本概念。 * * 2.方程组的分块（Partitioning） 例：稀疏方程组 * * 将方程组分成单独求解的维数较低的子方程组，并确定子方程组求解顺序的过程 x1、 x4 f1 、f4 x2 f3 x3、 x5 f2 、f5 1 2 3 3. 变量的切割（Tearing） 稀疏方程组分块后，其子方程组还可能是稀疏的，但又不能再通过分解降维。可以设想以下求解办法： 首先选择几个变量并给以估计值，然后利用稀疏方程算出这些变量的计算值，利用前面介绍的迭代计算方法，计算出这些变量的解。 这一方法叫做切割法，被选择的少数几个变量叫做切割变量（Tearing Variables）。 * * 稀疏方程的切割解法 * * 选出x3和x5为切割变量 x2 x4 x1 x5′ x3′ 注 意 切割法实际上只是在求解方程组时，可供选用的一种求解策略，它并不是和本章以上所讲的各种解法地位相当的一种方法，它不能取代上面的解法。 采取了切割法以后，只是把需要迭代求解的变量数减少了，相当于只须去求解一个规模缩小了的方程组。规模尽管缩小了，要去求解，仍然还是要选用上面讲过的某种解法。 * * 注意：这个要去求解的较小方程组，不再具有通常方程组所具有的那种数学形式。其求解过程中的每一轮，是从这少数几个变量已取得了估计值（初值）的情况开始，通过对原方程组中含有的某些方程进行一系列运算，直至又得出这些变量的计算值，再据此生成下一轮这些变量的估计值。 进行这样的求解最宜选用的解法，是只涉及函数值，不涉及函数