化工过程分析与合成 化工过程系统的优化培训教案.ppt

化工过程分析与合成Analysis and Synthesis of Chemical Process  掌握线性规划问题的求解方法； 掌握拉格朗日乘子法及罚函数法； 了解过程系统优化的种类及优化的基本概念； 了解过程系统参数模拟优化方法； 了解变分法。 4.1 概述 在多种可行方案中找出费用最小的流程结构，保证该方案满足安全、环保、易操作等方面的要求 确定冷、热物流的匹配方式，以便充分利用系统内部热量，降低公用工程消耗 动态系统参数的最优化又称连续系统最优化，这是由于优化问题的解是时间t的连续函数 动态系统参数优化问题的一般模型 可以看出，目标函数随状态变量和决策变量的不同而不同，就是说明目标函数是函数的函数。 在数学上，这种函数的称为泛函，求泛值的问题称为变分问题。因此，连续系统的最优化问题就是一个变分问题。 一 无约束连续系统的最优化 当函数 y (x) 经微小改变后变为y1 (x)，则称 为函数的变分，表示了y(x) 的微小改变。也写为 式中 是一个连续可微的任意函数，ε是一个很小的正数，当 时， 若y(x)的微小改变要求在区间[x1，x2]两端固定，即保持y(x1)=y1，y(x2)=y2，则 应满足 或记为： 相应于函数y(x)的微小改变，泛函I [ y (x)]的改变量为： 将大括号内的函数的用泰勒级数展开，并略去ε的高次项得: 式中 称为泛函I [ y (x)]的第一变分。泛函 I [ y (x)]取极值的必要条件为： 将前式等号右侧第二项分部积分得到： 在端点固定（即y(x1)=y1及y(x2)=y2）的条件下， 所以上式等号右侧第一项为零。因此 由于 为任意函数，不恒等于零，所以要式4-32成立，必然有 上式就是使泛函I [ y (x)]取极值的必要条件，称为欧拉方程。有了这个方程，求泛函极值的问题就转化为求解微分方程的问题了 对于（4-27）式给出的优化问题，其极值为下列偏微分方程组的解： (2) 泛函极值的充分条件 要判断满足欧拉方程的函数是使泛函极大还是极小，需计算第二变分 对 作泰勒级数展开，忽略 以上的高次项，得到： 式中等号右侧第二项记为 ，称为第二变分： 对于满足欧拉方程的函数y(x)， ，因此，若函数y(x)使第二变分 ,则有 即函数y(x)使泛函I [ y (x)]取极小； 若函数使第二变分 ，则 即函数y(x)使泛函I [ y (x)]取极大。 我们把第二变分写成矩阵的形式： 其中矩阵： 就是Hesse矩阵，这时大于零（或小于零）与 Hesse的正定（或负定）是一致的，二者都 可作为判定泛函极值的充分条件 二 有约束连续系统的最优化 (1) 目标函数不含终态 该问题是求函数向量x(t)和w(t)，使之满足常微分形式的约束条件，且使泛函 极小的问题 先考虑一维情况。引进拉格朗日泛函： 式中 称做伴随函数或拉格朗日乘子函数 利用泛函极值必要条件欧拉方程可得： 即： 边界条件： 由于 时x是自由的，所以由欧拉方程的推导过程（4-33)式可知， 。此时只有 才能使欧拉方程成立。因此有： 由（4-43）式有 所以， 应满足： 这样，原变分问题便转化成下面微分方程组的边值问题： 边界条件： 若给定终端条件 的话，只需用其替代式中的边界条件即可。 和一般函数一样，上述方法称为拉格朗日方法。也可以推广到多维状态函数和决策函数的情形 边界条件： 同样可定义拉格朗日泛函： 由此可导出相应的两点边值问题： 边界条件： 这里共有2m+r个方程，可解出2m+r个未知函数，作为原问题的最优解 例 4-8 有回流的连续反应系统的最优控制 浓度为xi的物料以常流量q流入混合器，与分离器回流的物料混合后，以浓度x0及流量q+r进入反应器。反应后浓度为x?，在分离器中分离出浓度为xf的产品以流量q流出。另一部分浓度为xr的物料以流量r回流至混合器 反应器内进行 不可逆吸热反应，反应动力学方程： (1) 式中 T 为反应温度。分离器的分割比s定义为 (2) 令C1表示纯产品的单价，C2表示反应器升高单位温度所消耗的能量费用，C3表示原料的单位成本。问题是应如何控制反应温度T，才能在[0,? ]时间内获得最大的利润 分离器物料平衡为： （4） 将2式代入4得： （5） 混合器物料平衡为： （6） 将2式代入6得： （7