CN114237187A - 一种工业过程的约束学习先进控制方法 - Google Patents
一种工业过程的约束学习先进控制方法 Download PDFInfo
- Publication number
- CN114237187A CN114237187A CN202111560384.9A CN202111560384A CN114237187A CN 114237187 A CN114237187 A CN 114237187A CN 202111560384 A CN202111560384 A CN 202111560384A CN 114237187 A CN114237187 A CN 114237187A
- Authority
- CN
- China
- Prior art keywords
- linear
- model
- industrial process
- optimization problem
- control law
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000004519 manufacturing process Methods 0.000 title claims abstract description 24
- 238000000034 method Methods 0.000 title claims abstract description 24
- 238000005457 optimization Methods 0.000 claims abstract description 20
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 claims abstract description 19
- 230000006870 function Effects 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims description 10
- 230000033001 locomotion Effects 0.000 claims description 6
- 238000012887 quadratic function Methods 0.000 claims description 3
- 238000013461 design Methods 0.000 abstract description 4
- 238000004886 process control Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 229960000074 biopharmaceutical Drugs 0.000 description 1
- 238000003889 chemical engineering Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012847 fine chemical Substances 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 238000004540 process dynamic Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/418—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
- G05B19/41885—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by modeling, simulation of the manufacturing system
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/32—Operator till task planning
- G05B2219/32339—Object oriented modeling, design, analysis, implementation, simulation language
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- General Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种工业过程的约束学习先进控制方法。本发明首先通过分段模型识别方法导出一系列线性仿射模型,从而得到多个线性仿射模型。其次针对工业过程的二维特性,通过引入状态误差和输出误差,构建新的等价状态空间模型,此基础上进行设计,设计了含有抗外界扰动的二次性能指标函数,采用两步优化的方法进一步提高了系统的控制性能。本发明将非线性模型转化成一系列线性仿射模型,不仅使得所考虑的非线性模型转化为形式较简单的线性模型,还使得所转化这个模型更接近非线性系统,同时系统受到不同形式的干扰也可考虑其中,解决了非线性特性给系统带来的影响。
Description
技术领域
本发明属于自动化过程控制领域,涉及一种工业过程的约束学习先进控制方法。
背景技术
近年来,伴随着精细化工,生物制药、金属加工等领域内大量新产品的开发和研究,工业过程的先进控制已经受到了工业界和学术界的广泛关注。
在实际生产过程中,工业过程普遍具有非线性特性,因此线性模型和实际工业过程之间存在较大的不匹配问题,若利用确定性线性模型来简单地描述非线性过程,将使得在实际应用中很难达到最佳的控制效果。另一方面,工业过程还具有二维特性,即批次维度特性和时间维度特性,寻找与两个维度上相关的控制算法对提高工业过程中的控制性能也极为重要。
作为一种先进的控制算法,模型预测控制方法大多在一维方向上设计控制律,但只考虑时间方向,那么批次方向只是单独的重复,控制性能无法随着批次方向的递增而提高;只考虑批次方向则会产生无法确定初值等问题。对于以上问题,不仅影响系统的控制性能,甚至会影响系统的稳定性。
因此,针对工业过程的非线性和二维特性,提出一种非线性工业过程受限鲁棒模型预测控制最小-最大的综合优化方法,在保证系统稳定的基本条件下,提高系统的控制性能具有重要意义。
发明内容
本发明的目的是为针对非线性工业过程以及工业过程具有的二维特性,提出了一种非线性工业过程受限鲁棒模型预测控制最小-最大的综合优化方法,使得系统稳定运行并实现良好的控制性能。
本发明首先通过分段模型识别方法导出一系列线性仿射模型,从而得到多个线性仿射模型。其次针对工业过程的二维特性,通过引入状态误差和输出误差,构建新的等价状态空间模型,此基础上进行设计,设计了含有抗外界扰动的二次性能指标函数,采用两步优化的方法进一步提高了系统的控制性能。
本发明的具体实现步骤包括:
步骤一:将原始的非线性工业过程分为一系列仿射运算区域,通过对仿射运算区域相对于多个平衡操作点线性化,得到一系列线性仿射模型。并且为了在实际应用中可以描述批次内或批次间的模型不确定性,通常将线性仿射模型转化为多面体不确定性模型。
步骤二:基于工业过程的二维特性,设计迭代学习控制律,并通过引入状态误差和输出跟踪误差,选取包含状态变量和跟踪误差的新的状态空间变量,构建新的等价状态空间模型。
步骤三:针对新的等价状态空间模型,设计含有抗外界扰动的二次性能指标函数。将性能指标函数在整个无限时域范围内分为两部分:前N步和N+1步及以后。对于前N步是一个有限时域约束的最小-最大优化问题,N+1步及以后是一个无限时域约束的最小-最大优化问题。根据线性矩阵不等式约束以及利用性能指标函数的凸优化问题,对优化问题进行求解,得到系统的控制律。
进一步说,步骤三中,在整个无限范围内,本发明没有将输入参数化为单个线性状态反馈控制律,而是在线性反馈控制律之前添加了N个自由控制运动。对N+1步及以后的无限时域约束的最小-最大优化问题,建立了一个参数相关的二次函数来保证系统的鲁棒稳定性,在考虑有干扰的条件下,采用线性矩阵不等式的形式对迭代更新律的增益进行求解,进一步计算得到线性状态反馈控制律。系统的控制律构成为前N步自由控制运动加上线性状态反馈控制律。
本发明与现有技术相比有益效果为:
(1)本发明将非线性模型转化成一系列线性仿射模型,不仅使得所考虑的非线性模型转化为形式较简单的线性模型,还使得所转化的模型更接近非线性系统,同时系统受到不同形式的干扰也可考虑其中,解决了非线性特性给系统带来的影响。
(2)本发明提出了一种非线性工业过程受限鲁棒模型预测控制最小-最大的综合优化方法,它提供了更多的自由度来调节系统控制性能,从而使系统在满足稳定性的基本要求下获得了更好的控制性能。
具体实施方式
下面通过具体实施例对本发明做进一步的说明。
非线性连续搅拌釜反应过程控制是一个典型的工业过程,考虑高度非线性连续搅拌釜反应器(CSTR),目标是让反应器温度遵循给定的曲线,具体实现包括以下步骤:
步骤一:将高度非线性连续搅拌釜反应器分为一系列分段仿射运算区域:
其中,t和k分别表示时间和批次指数,x(t,k),y(t,k),u(t,k)分别表示系统在批次k时间t的状态,输入和输出量,Ω表示实际的输出响应空间,包括由Ωi(s=1,2,...,m)表示的m个分段仿射运算区域,TP是“循环时间”,表示每个循环的初始重置条件,并且为了便于控制设计,每个循环可以相对于每个仿射操作区域重置为零。
为了描述批次内或批次间的模型不确定性,具有多面体不确定性的可观察规范形式的离散时间模型结构可以写成:
s=1,2,...,m;
步骤二:根据二维迭代学习控制律引入状态误差δk(x[s](t,k))=x[s](t,k)-x[s](t,k-1)和跟踪误差e[s](t,k)=yr [s](t)-y[s](t,k),表示为迭代过程的初值,r[s](t,k)表示为迭代学习更新率,表示为设定值。选取包含状态变量和跟踪误差的状态空间变量z[s](t,k),构建新的等价状态空间模型:
具有适当维度的零矩阵,I表示为具有适当维度的单位矩阵。
步骤三:针对新构建的等价状态空间模型,设计含有抗外界扰动的二次性能指标函数,使反应器温度遵循给定的曲线。首先选取相关的性能指标函数:
约束为:
将性能指标函数分为两部分:
对于i=0,1,...,N-1的有限水平约束最小-最大优化问题是在线性反馈控制律之前添加了N个自由控制动作,由于包含跟踪误差和状态误差,可以改进系统的控制性能。
对于i≥N的无限水平约束最小-最大优化问题,引入线性状态反馈控制律。定义二次函数Vi [s](t,k))=z[s](t+i,k)TPi [s](t,k)z[s](t+i,k),i≥N,其中Pi [s](t,k)>0表示为待定的正定矩阵。保证满足以下鲁棒稳定性约束,使优化问题可解:
通过从i=N到∞进行求和可以得到:
应用Schur引理,将约束转换为线性矩阵不等式的形式,以线性矩阵不等式的形式对迭代学习更新律r[s](t,k)=-K[s](t,k)z[s](t,k)的增益K[s]进行求解,进一步计算得到线性状态反馈控制律。对于整个优化问题的控制律构成为前N步自由控制运动加上线性状态反馈控制律。通过这种方式,将产生更多的最优控制输入,提高了系统的控制性能。
Claims (2)
1.一种工业过程的约束学习先进控制方法,其特征在于该方法包括以下步骤:
步骤一:将原始的非线性工业过程分为一系列仿射运算区域,通过对仿射运算区域相对于多个平衡操作点线性化,得到一系列线性仿射模型;
将线性仿射模型转化为多面体不确定性模型;
步骤二:基于工业过程的二维特性,设计迭代学习控制律,并通过引入状态误差和输出跟踪误差,选取包含状态变量和跟踪误差的新的状态空间变量,构建新的等价状态空间模型;
步骤三:针对新的等价状态空间模型,设计含有抗外界扰动的二次性能指标函数;
将性能指标函数在整个无限时域范围内分为两部分:前N步和N+1步及以后;对于前N步是一个有限时域约束的最小-最大优化问题,N+1步及以后是一个无限时域约束的最小-最大优化问题;
根据线性矩阵不等式约束以及利用性能指标函数的凸优化问题,对优化问题进行求解,得到系统的控制律。
2.根据权利要求1所述的一种工业过程的约束学习先进控制方法,其特征在于:步骤三中,在整个无限范围内,在线性反馈控制律之前添加了N个自由控制运动;对N+1步及以后的无限时域约束的最小-最大优化问题,建立了一个参数相关的二次函数来保证系统的鲁棒稳定性,在考虑有干扰的条件下,采用线性矩阵不等式的形式对迭代更新律的增益进行求解,进一步计算得到线性状态反馈控制律;系统的控制律构成为前N步自由控制运动加上线性状态反馈控制律。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111560384.9A CN114237187A (zh) | 2021-12-20 | 2021-12-20 | 一种工业过程的约束学习先进控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111560384.9A CN114237187A (zh) | 2021-12-20 | 2021-12-20 | 一种工业过程的约束学习先进控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114237187A true CN114237187A (zh) | 2022-03-25 |
Family
ID=80759163
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111560384.9A Pending CN114237187A (zh) | 2021-12-20 | 2021-12-20 | 一种工业过程的约束学习先进控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114237187A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114721268A (zh) * | 2022-04-08 | 2022-07-08 | 江南大学 | 注塑成型喷嘴压力鲁棒启发式迭代学习控制方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050107895A1 (en) * | 2001-05-25 | 2005-05-19 | Efstratios Pistikopoulos | Process control |
CN108227494A (zh) * | 2018-01-05 | 2018-06-29 | 海南师范大学 | 非线性批次过程2d最优约束模糊容错控制方法 |
CN110597055A (zh) * | 2019-08-05 | 2019-12-20 | 海南师范大学 | 抗不确定性的2d分段仿射间歇过程最小-最大优化的预测控制方法 |
CN113093536A (zh) * | 2021-03-23 | 2021-07-09 | 南京航空航天大学 | 非线性系统执行器故障的自适应增量优化容错控制方法 |
US20210373513A1 (en) * | 2020-05-29 | 2021-12-02 | Mitsubishi Electric Research Laboratories, Inc. | Nonlinear Optimization Method for Stochastic Predictive Control |
-
2021
- 2021-12-20 CN CN202111560384.9A patent/CN114237187A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050107895A1 (en) * | 2001-05-25 | 2005-05-19 | Efstratios Pistikopoulos | Process control |
CN108227494A (zh) * | 2018-01-05 | 2018-06-29 | 海南师范大学 | 非线性批次过程2d最优约束模糊容错控制方法 |
CN110597055A (zh) * | 2019-08-05 | 2019-12-20 | 海南师范大学 | 抗不确定性的2d分段仿射间歇过程最小-最大优化的预测控制方法 |
US20210373513A1 (en) * | 2020-05-29 | 2021-12-02 | Mitsubishi Electric Research Laboratories, Inc. | Nonlinear Optimization Method for Stochastic Predictive Control |
CN113093536A (zh) * | 2021-03-23 | 2021-07-09 | 南京航空航天大学 | 非线性系统执行器故障的自适应增量优化容错控制方法 |
Non-Patent Citations (2)
Title |
---|
罗文平: "具有执行器故障的间歇过程2D迭代学习预测控制策略研究", 《中国优秀硕士学位论文全文数据库 信息科技辑》, pages 1 - 2 * |
邹志云等: "一种新型线性约束系统预测控制算法", 《化工学报》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114721268A (zh) * | 2022-04-08 | 2022-07-08 | 江南大学 | 注塑成型喷嘴压力鲁棒启发式迭代学习控制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ji et al. | Parameter estimation for block‐oriented nonlinear systems using the key term separation | |
Zhao et al. | Intelligent tracking control for a class of uncertain high-order nonlinear systems | |
Zhang et al. | Variance‐constrained state estimation for networked multi‐rate systems with measurement quantization and probabilistic sensor failures | |
Manzano et al. | Output feedback MPC based on smoothed projected kinky inference | |
Hou et al. | Local learning‐based model‐free adaptive predictive control for adjustment of oxygen concentration in syngas manufacturing industry | |
Li et al. | Maximum likelihood identification of dual‐rate Hammerstein output‐error moving average system | |
Liu et al. | Event‐based model‐free adaptive control for discrete‐time non‐linear processes | |
Zhang et al. | Data‐driven design of two‐degree‐of‐freedom controllers using reinforcement learning techniques | |
Shen et al. | Nonfragile fuzzy control for nonlinear fast sampling singularly perturbed systems subject to Markov jumping parameters | |
Si et al. | Adaptive neural control for stochastic pure‐feedback non‐linear time‐delay systems with output constraint and asymmetric input saturation | |
CN114237187A (zh) | 一种工业过程的约束学习先进控制方法 | |
Chi et al. | Compensation‐based data‐driven ILC with input and output package dropouts | |
CN110597055B (zh) | 抗不确定性的2d分段仿射间歇过程最小-最大优化的预测控制方法 | |
Esterhuizen et al. | Control design with guaranteed transient performance: An approach with polyhedral target tubes | |
Belhaj et al. | MIMO PI controllers for LTI systems with multiple time delays based on ILMIs and sensitivity functions | |
Ikonen et al. | Process control using finite Markov chains with iterative clustering | |
Hast et al. | Optimal proportional–integral–derivative set‐point weighting and tuning rules for proportional set‐point weights | |
Yi et al. | An optimal fuzzy control method for nonlinear time-delayed batch processes | |
Luan et al. | Compensator design based on inverted decoupling for non‐square processes | |
Malik et al. | State and output feedback local control schemes for nonlinear discrete-time 2-D Roesser systems under saturation, quantization and slope restricted input | |
Cheng | A grey CES production function model and its application in Calculating the contribution rate of economic growth factors | |
Luo et al. | 2D fuzzy constrained fault-tolerant predictive control of nonlinear batch processes | |
de Castro et al. | Unrestricted horizon predictive control applied to a nonlinear SISO system | |
CN107194181A (zh) | 基于四元数和最小平均峭度准则的多维时间序列预测方法 | |
Xia et al. | Adaptive DSC of stochastic non‐linear systems with input unmodelled dynamics |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |