基于连续搅拌反应釜的多目标分层预测控制方法
专利汇可以提供基于连续搅拌反应釜的多目标分层预测控制方法专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种基于连续搅拌反应釜的多目标分层预测控制方法。该发明通过引入双层预测控制 算法 ，改变传统连续搅拌釜反应釜(CSTR)控制过程中，只考虑动态控制目标，忽略经济指标的现状。在控制过程中，采用两级分层预测控制结构，动态控制层引入 扩展卡 尔曼观测器对状态进行无偏估计，稳态优化层采用多目标分层优化，优先满足重要的优化目标。同时，整个控制系统采用有限元 正交 配置法将非线性优化控制命题转化为非线性规划问题，模型方程作为约束条件，配置点上的值作为预测值，不必通过线性化或阶跃响应获取 预测模型 ，从而避免线性化误差，存在积分对象需单独处理等问题。，下面是基于连续搅拌反应釜的多目标分层预测控制方法专利的具体信息内容。
1.一种基于连续搅拌反应釜的多目标分层预测控制方法，其特征在于，包括以下步骤：
(1)采集状态变量x(t)和输入变量u(t)的初始值，所述状态变量x(t)包括出口浓度、反应器体积、反应器温度，所述输入变量u(t)包括反应器之间的传热量、反应器出口阀位；
(2)将步骤1中的状态变量x(t)和输入变量u(t)的初始代入到CSTR动态
模型 中。根据动态模型得到CSTR稳态优化模型
为x(t)的导数；y(t)表示输出变量，所述输出
变量包括反应器混合物出口浓度、反应器温度，f、g表示模型方程；
(3)对CSTR稳态优化模型、CSTR动态模型采用联立法中的有限元正交配置法进行离散化处理，确定输入变量u(t)的预测值和输出变量y(t)的预测值，并确定使工厂经济效益最大化的输出最优值 输入最优值 为动态控制层提供最优控制目标，所述 通过以下方法得到：
(3.1)可稳性判断，确定系统稳定时输入变量允许波动的范围u′min,u′max；
yij＝g(xij)
xiK＝x(i+1)K＝xs,i＝1...NE
uiK＝u(i+1)K＝us,i＝1...NE
umin≤us≤umax
其中，J1表示可稳性优化目标，y(kt+l|kt)、y(kt+l-1|kt)分别表示kt时刻对kt+l、kt+l-1时刻的输出预测值，wy表示输出变量的权重矩阵，τij表示第i个有限元第j个插值点，φk(τij)表示状态变量的拉格朗日插值函数，xik表示第i个有限元第k个插值点处状态变量的值，NE表示有限元个数，K表示插值点的个数，hi表示第i个有限元的长度，f(xij,uij)表示模型状态方程，yij表示通过公式yij＝g(xij)计算得到的第i个有限元第j个插值点处对应的输出变量的值，xiK、x(i+1)K、xs分别表示第i个有限元第K个插值点的状态值、第i+1个有限元第K个插值点的状态值、稳态时的状态值，uiK、u(i+1)K、us分别表示第i个有限元第K个插值点的输入值、第i+1个有限元第K个插值点的输入值、稳态时的输入变量值，umin、umax分别表示初始要求允许的输入变量的最小、最大值，若在初始要求允许的范围内，J1≠0，系统的可稳性条件不满足，控制器停止运行，将umin、umax分别调整为u′ min,u′max，使得J1＝
0。
(3.2)可行性分析，确定系统稳定时输出变量允许波动的范围y′min,y′max；
yij＝g(xij)
xiK＝x(i+1)K＝xs,i＝1...NE
uiK＝u(i+1)K＝us,i＝1...NE
y(kt+l|kt)-y(kt+l-1|kt)＝0
u′min≤us≤umax
ymin-δy≤ys≤ymax+δy
δy≥0
其中，J2表示可行性优化目标，δy用来衡量实际输出值是否在控制要求允许波动范围内，us表示输入变量稳态值，wy是输出变量y对应的优先权重矩阵；ymin、ymax分别为初始要求允许的输出变量的最小值和最大值；若在初始要求允许的范围内，J2≠0，不满足可行性条件，将ymin、ymax按照以下方式进行调整，使得J2＝0。
y′min←inf(g(xNEK),ymin) y′max＝sup(g(xNEK),ymax)
其中xNEK表示第NE个有限元，第K个配置点处的状态变量值。
(3.3)经济优化，确定全厂经济效益最大化时的输入变量值 和输出变量值
yij＝g(xij)
xiK＝x(i+1)K＝xs,i＝1...NE
uiK＝u(i+1)K＝us,i＝1...NE
y(kt+l|kt)-y(kt+l-1|kt)＝0
u′min≤us≤u′ max
y′min≤ys≤y′ max
其中，J3表示经济优化目标，IRVu、IRVy由全厂优化层给出，a1、a2分别为us，ys的二次权重；在u′min≤us≤u′max和y′min≤ys≤y′max条件下，求出使得J3值最小的并将u′min，u′max，y′min，y′max做如下修改：
(3.4)工作点优化，确定扰动条件下的最优输入变量值
yij＝g(xij)
xiK＝x(i+1)K＝xs,i＝1...NE
uiK＝u(i+1)K＝us,i＝1...NE
其中，J4表示工作点优化目标， 是kt-1时刻操纵变量的最佳稳态值，wu是输入变量u对应的优先权重矩阵。在 和 条件下，求出使得J4值最小的
(4)根据稳态优化层提供的控制目标 和 确定动态优化命题：
其中，JA为动态优化命题，jj表示采样时刻，y(kt+jj|kt)表示根据当前kt时刻的输出值y(kt)和动态离散化模型对未来kt+jj时刻的输出预测值，y(kt)＝g(x(kt))，x(kt)通过扩展卡尔曼观测器得到；u(kt+jj)表示未来kt+jj时刻的操作变量值，Δu(kt+jj)表示输入变量改变量，
Δu(kt+jj)＝u(kt+jj)-u(kt+jj-1)，Qjj、Rjj、Tjj分别表示输出变量权重矩阵、输入变量权重矩阵、输入变量改变量权重矩阵，P、M分别表示预测时域、控制时域。
(5)对输出预测值y(kt+jj|kt)通过以下方式进行校正：
其中,JB为校正后的动态优化命题， 表示kt+jj时刻的真实值y(kt+jj)与预测值y(kt+jj|kt)之间的差，采用内点法求解器IPOPT对改进后的动态优化命题进行求解，确定当前kt时刻操作变量 的值。
(6)不断重复上述步骤，使得y(kt)不断接近于 直到当前kt时刻操作变量 的值与最优稳态输入值 的值相同，JB＝0，y(kt)等于 完成多目标分层预测的控制。
2.根据权利要求1所述的方法，其特征在于，所述步骤2中的有限元正交配置法具体为：
把计算域剖分为NE个互不重叠且相互连接的单元，在有限元上选取K个点作为配置点，选取拉格朗日多项式为插值函数，对步骤1中的动态模型和稳态模型分别进行离散化，得到如下动态离散模型：
i＝1...NE,j＝1...K
yij＝g(xij)
x10-x0＝0
x(i+1)0＝xiK
得到如下稳态离散模型：
hif(xij,uij)＝0,i＝1...NE,j＝1...K
yij＝g(xij)
x10-x0＝0
x(i+1)0＝xiK
hi表示每个有限元的长度，φk表示状态变量对应的拉格朗日插值函数，τij代表配置点，xik表示第i个有限元上第j个配置点的状态变量值，uij表示第i个有限元上第j个配置点的操作变量值，x10为第一个有限元的第0个配置点的值，即状态的初值，x(i+1)0表示第i+1个有限元的第0个配置点的状态值，xiK表示第i个有限元的第K个配置点的状态变量值。
3.根据权利要求1所述的方法，其特征在于，所述步骤4中的根据当前kt时刻的输出值y(kt)和动态离散化模型预测未来kt+jj时刻的输出预测值y(kt+jj|kt)的预测方法具体为：
通过扩展卡尔曼观测器得到kt时刻状态估计值x(kt)，作为第一个有限元第0个配置点的初值x10，根据动态离散模型 i＝1...NE,j＝1...K得到
x11...x1K的值，再结合yij＝g(xij)得到输出预测值y(kt+jj|kt)。
基于连续搅拌反应釜的多目标分层预测控制方法
技术领域
背景技术
发明内容
出变量包括反应器混合物出口浓度、反应器温度，f、g表示模型方程；
具体实施方式
A反应釜体积V1 200
A反应釜出料产品浓度q1 0.0357
A反应釜温度T1 446.471
B反应釜体积V2 100
B反应釜出料产品浓度q2 0.0018
B反应釜温度T2 453.2585
B反应釜的出口阀位u1 2
传递到A反应釜的热量u2 0
B反应釜的出口阀位u1 0-1
传递到A反应釜的热量u2 0-0.6
A反应釜的出料产品浓度Ca1 0.004-0.006
B反应釜温度T2 480-495
A反应釜的出料产品浓度Ca1 0.0042
B反应釜温度T2 489.3997
B反应釜的出口阀位u1 0.9013
传递到A反应釜的热量u2 0.5602
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种基于长短期记忆循环神经网络的空调能耗预测方法 | 2020-10-26 | 0 |
| 一种儿童近视风险预测方法、装置及系统 | 2021-12-29 | 0 |
| 一种高炉炼铁过程一氧化碳利用率的建模方法 | 2022-07-07 | 0 |
| 一种基于iPLS-PA算法的硝酸盐浓度预测方法 | 2020-05-20 | 2 |
| 一种语音评价方法及装置 | 2020-12-10 | 0 |
| 一种基于选择性隐朴素贝叶斯分类器的网络故障诊断方法 | 2021-06-19 | 1 |
| 心脑血管疾病风险预测方法及系统 | 2022-04-22 | 1 |
| 容量资源的预测方法及装置 | 2020-08-03 | 0 |
| 一种基于历史数据的汽车行驶轨迹及里程预测方法 | 2021-05-22 | 1 |
| 基于结构磁共振的带状疱疹性神经痛疗效预测方法及系统 | 2020-10-13 | 1 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
