第热能系统分析与最优综合四章6,7.ppt

第六节联立方程法中非线性方程的求解 式中 ——状态变量（相关变量）向量； ——决策变量（独立变量）向量； ——需同时求解的所有方程 F v 联立方程法是把描述流程的所有方程式收集起来，形 成一大组需要联立求解的线形或非线性代数方程组。 在数学上可表示为 一、联立方程的建立及其稀疏性 基本方程——所有单元模型均可由有限个种类的方 程构成，这些方程称为基本方程。 等效方程——一类特殊的线性方程，用它们可以直 接求出非线性方程不能直接解出的未 知量的近似值，作为求解时的初值。 基本方程 种类与用途 混合与分流的物料衡算方程 脱氧器物料衡算方程 泵、透平的能量衡算方程 绝热等熵膨胀方程 加热／冷却器的能量衡算方程 部分基本方程一览表 在建立联立方程时，要了解需要的变量、方程式以及自由度的数目，从而正确地选择决策变量和给定必要的参数值，使自由度为零。 　　　决策变量的选择，对系统的模拟解算过程也有很大影响。若决策变量选择是恰当的，可以尽量消除再循环迭代。 　　　模型方程通常是一组大型稀疏方程组。所谓稀疏方程组就是每个方程中只含有少数几个非零系数。 　　　对于线性方程组，可以用矩阵形式表示: 　　　 　AX=b 　系数矩阵A中的大部分元素为零。 　　　对于大型线性方程组的求解已有几种成熟方法，最常用的是高斯消元法。然而，非线性广泛存在于热能工程领域，联立方程法的核心问题是求解大型稀疏非线性方程。 二、大型稀疏非线性方程组的分隔降价解法 　 利用稀疏的特点，可以将复杂的大型稀疏非线性方程组分解若干个较为简单的小型稀疏方程组，然后依次求解——分隔降阶解法 其基本手段是采用分隔和断裂技术： 1. 方程组的分隔 当用一组方程式来描述系统时，也可以用布尔矩阵来表达，称为事件矩阵。事件矩阵是描述每个方程式对各变量的依赖关系。 事件矩阵的特点： （1）事件矩阵的每一行对应一个方程式，每一列 对应一个系统变量； （2）矩阵中的元素 按下列规定只取布尔数1或0； 事件矩阵是表达信息流结构的一种简洁方法，可以用它来进行方程组的分隔。 1 当变量j出现于方程式i中时 0 变量j在方程式i中不出现时 如下方程组 重新 排列 和 中只包含变量 和 ， 和 中只包含 和 ，可以分别单独联立求解 对于n阶必须单独联立求解的稀疏方程组，还可以进一步分隔和排序。n阶方程中可以找到一个包含　个变量的　阶子方程，这个　阶子方程可以单独求解，得到　个变量。余下　　　阶方程的按上述方法继续求解。 例如 原来的五元联立方程可以分隔成对角线上3个可按顺序求解的子方程 联立求解 2.不可分隔稀疏方程组的断裂降阶解法 无法分隔的联立方程组，可以通过选择断裂变量，把高阶方程组降阶成为低阶方程组进行迭代求解。 如 相应事件方程 选择 中 的为断裂变量，赋予初值，由 解出 这样可以分别联立求解 和 ，将解出的变量代入 中，检验是否满足。若不满足，修改断裂变量，重复上述计算，直到满足为止。 计算结果 Yes No 解算过程图 三、非线性方程组求解的牛顿-拉夫森迭代法 由于热能系统模型中的大部分方程是线性的，只有少数是非线性方程，因此，可以考虑把非线性方程用某种方式线性化，然后再联立求解，称为拟线性解法。 1.牛顿-拉夫森迭代法 这种方法是先用函数的偏导数来获得非线性方程的线性近似式，然后将其并入线性方程组，联立求解包含少量近似式的线性方程组。 对一个n维非线性方程组 … … … 给出一组初始近似解 … … … 泰勒级数展开 只取一阶线性项 通过解线性方程组，得到： 重复上述过程，直到满足精度（ ）为止，得到原方程组的近似解。 方程组也可以用向量式表示 … … … = 其中 雅可比矩阵 例 求解下列非线性方程组 解 求出雅可比矩阵 取初值 k 0 2 10 5 18.4142 48.0000 0.9505 1 0.7084 8.3416 4.0334 1.4871 5.3527 -0.1837 2 0.9336 8.0185 3.9932 0.0014 0.1711 -0.0247 3 0.9995 8.00