CN108681317A - 一种化工工业过程鲁棒学习控制方法 - Google Patents
一种化工工业过程鲁棒学习控制方法 Download PDFInfo
- Publication number
- CN108681317A CN108681317A CN201810760548.4A CN201810760548A CN108681317A CN 108681317 A CN108681317 A CN 108681317A CN 201810760548 A CN201810760548 A CN 201810760548A CN 108681317 A CN108681317 A CN 108681317A
- Authority
- CN
- China
- Prior art keywords
- moment
- kth
- period
- error
- batch
- 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
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], 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], 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]
Abstract
本发明公开了一种化工批次过程鲁棒学习控制方法,包括如下步骤:步骤1、建立以预测值形式表示的闭环系统模型;步骤2、设计被控对象的批次处理控制器。该方法首先根据给定的系统模型设计鲁棒学习控制量,然后通过引入状态误差和输出误差,将系统的原始模型转化成以预测值形式表示的闭环系统模型,进一步根据设计的最优化性能指标和李雅普诺夫稳定性理论,给出了线性矩阵不等式表示的闭环系统鲁棒渐进稳定的充分条件和最优控制律表达形式。不同于传统的控制方法,本方法设计了无限时域的性能指标,使系统能够在不超过给定值的情况下克服由最小状态能量和输入增量引起的最大扰动。
Description
技术领域
本发明属于自动化工业过程控制领域,涉及到一种化工批次过程鲁棒学习控制方法。
背景技术
在现代制造业中,批次处理作为重要的生产模式已广泛应用于小规模高附加值产品制造,与此同时,对于批次处理的理论研究也取得了巨大的突破。在批次处理的一些常规程序中,会表现出明显的重复性和不确定性,而导致系统不能稳定的运行。为了优化实际生产过程中的重复性和不确定性,特别是系统内部参数扰动和外部干扰,结合复合的鲁棒学习控制策略,提高系统的鲁棒性,使系统运行更加稳定。
发明内容
本发明的目的是优化批次处理过程中产生的重复性和不确定性。该方法首先根据给定的系统模型设计鲁棒学习控制量,然后通过引入状态误差和输出误差,将系统的原始模型转化成以预测值形式表示的闭环系统模型,进一步根据设计的最优化性能指标和李雅普诺夫稳定性理论,给出了线性矩阵不等式表示的闭环系统鲁棒渐进稳定的充分条件和最优控制律表达形式。不同于传统的控制方法,本方法设计了无限时域的性能指标,使系统能够在不超过给定值的情况下克服由最小状态能量和输入增量引起的最大扰动。
本发明的技术方案是通过模型建立、控制器设计、算法设计、优化等手段,提出了一种化工批次过程鲁棒学习控制方法,有效的弥补了传统批次处理中稳定性的不足,满足了实际工业过程的需求。
本发明的方法步骤包括:
步骤1、建立以预测值形式表示的闭环系统模型,具体步骤是:
1-1.建立具有不确定参数扰动的离散时间模型,形式如下:
其中,t和k分别是运行时刻和运行周期,x0,k是第k周期批次处理运行的初始条件,
u(t,k)和y(t,k)分别是t时刻第k周期系统的输入和输出,x(t,k)、x(t+1,k)分别是t时刻、 t+1时刻第k周期的系统状态,w(t,k)是t时刻第k周期未知的外部干扰,B、
C是具有适当维数的常数矩阵。
1-2.结合步骤1-1,引入鲁棒学习控制量:
u(t,k)=u(t,k-1)+r(t,k)u(t,0)=0,t=0,1,2,…
其中,u(t,0)是t时刻系统的初始输入,u(t,k)、u(t,k-1)分别是t时刻第k周期、第k-1周期的系统输入,r(t,k)是t时刻第k周期的更新律。
1-3.定义输出跟踪误差和误差的批次方向函数:
e(t,k)=y(t,k)-yr(t)
δf(t,k)=f(t,k)-f(t,k-1)
其中,yr(t)是输出的参考值,e(t,k)是t时刻第k周期的输出跟踪误差,f(t,k)是t时刻第k周期的系统状态、输出变量或未知的外部干扰,f(t,k-1)是t时刻第k-1周期的系统状态、输出变量或未知的外部干扰,δf(t,k)是t时刻第k周期误差的批次方向函数。
1-4.结合步骤1-1至步骤1-3,可以得到:
其中,δx(t,k)、δx(t+1,k)是t时刻、t+1时刻第k周期批次方向的系统状态误差,δw(t,k)是t时刻第k周期批次方向的外部干扰误差且 是不确定的系统矩阵,x(t,k-1)是t时刻第k-1周期的系统状态,δw(t,k)是t时刻第k周期的外部干扰,e(t+1,k)是t+1时刻第k周期的输出跟踪误差,e(t+1,k-1)是 t+1时刻第k-1周期的输出跟踪误差。
1-5.得到鲁棒渐进稳定的系统模型,形式如下:
其中, z(t+1,k)、z(t,k)是 t+1时刻、t时刻第k周期综合批次方向的系统状态误差和输出跟踪误差的扩展状态, z(t+1,k-1)是t+1时刻第k-1周期综合批次方向的系统状态误差和输出跟踪误差的扩展状态,δy(t,k)是t时刻第k周期批次方向的输出误差,A、I、E、F是适当维度的常量矩阵,G(t,k)是t时刻第k周期的常量矩阵。
1-6.以预测值形式表示的闭环渐进稳定系统模型,形式如下:
其中,j是预测值,H1、H2是对应的更新律增益,z(t+j+1|t,k)、 z(t+j+1|t,k-1)分别是t时刻到t+j+1时刻第k周期、第k-1周期的预测状态, z(t+j|t,k)是t时刻到t+j时刻第k周期的预测状态,是t时刻到t+j时刻第k周期批次方向的外部干扰误差,δx(t+j|t,k)、δy(t+j|t,k)分别是t时刻到t+j时刻第k周期批次方向的系统状态误差、输出误差
步骤2、设计被控对象的批次处理控制器,具体是:
2-1.基于步骤1,鲁棒性能保证控制下的更新律形式如下:
2-2.在重复性和非重复性扰动下,均可得到增益矩阵控制律的形式如下所示:
H1=γ-1Y1P,H2=γ-1Y2P
其中,P是满足系统渐进稳定的对称正定矩阵,Y1,Y2是满足系统条件的矩阵,γ-1>0是满足系统条件的系数。结合步骤1-2至步骤2-1即可得到最优系统输入u(t,k)。
2-3.在下一时刻,重复步骤1-6到2-2得到新的最优系统输入u(t,k),再将其作用于控制对象,并依次循环。
具体实施方式
以注塑成型工艺为例:
这里以注塑成型过程中喷嘴压力控制为例加以描述,调节手段是控制比例阀的阀门开度。
步骤1、建立以预测值形式表示注塑成型的闭环系统模型,具体步骤是:
1-1.建立注塑成型不确定参数扰动的离散时间模型,形式如下:
其中,t和k分别是注塑成型运行时刻和运行周期,x0,k是第k周期注塑成型批次处理运行的初始条件,u(t,k)和y(t,k)分别是t时刻第k周期注塑成型的阀门开度和喷嘴压力,x(t,k)、x(t+1,k)分别是t时刻、t+1时刻第k周期注塑成型的系统状态,w(t,k)是 t时刻第k周期注塑成型外部干扰,B、C是具有适当维数的常数矩阵。
1-2.结合步骤1-1,首先设计注塑成型过程比例阀的阀门开度:
u(t,k)=u(t,k-1)+r(t,k)u(t,0)=0,t=0,1,2,…)
其中,u(t,0)是t时刻比例阀的阀门开度的初始值,u(t,k)、u(t,k-1)分别是t时刻第k周期、第k-1周期比例阀的阀门开度,r(t,k)是t时刻第k周期的更新律。
1-3.定义注塑成型过程的输出跟踪误差和误差的批次方向函数:
e(t,k)=y(t,k)-yr(t)
δf(t,k)=f(t,k)-f(t,k-1)
其中,yr(t)是输出喷嘴压力的参考值,e(t,k)是t时刻第k周期注塑成型过的输出跟踪误差,f(t,k)是t时刻第k周期注塑成型的系统状态、输出变量或未知的外部干扰, f(t,k-1)是t时刻第k-1周期注塑成型的系统状态、输出变量或未知的外部干扰,δf(t,k)是t时刻第k周期注塑成型误差的批次方向函数。
1-4.结合步骤1-1至步骤1-3,可以得到:
其中,δx(t,k)、δx(t+1,k)是t时刻、t+1时刻第k周期注塑成型的批次方向系统状态误差,是t时刻第k周期注塑成型的批次方向外部干扰误差且ΔA(t,k)是不确定的系统矩阵,x(t,k-1)是t时刻第 k-1周期注塑成型的系统状态,δw(t,k)是t时刻第k周期注塑成型的外部干扰, e(t+1,k)是t+1时刻第k周期注塑成型的输出跟踪误差,e(t+1,k-1)是t+1时刻第 k-1周期注塑成型的输出跟踪误差。
1-5.得到注塑成型鲁棒渐进稳定的系统模型,形式如下:
其中,C1=[C 0], 是恒等于,z(t+1,k)、 z(t,k)是t+1时刻、t时刻第k周期注塑成型综合批次方向的系统状态误差和输出跟踪误差的扩展状态,z(t+1,k-1)是t+1时刻第k-1周期注塑成型综合批次方向的系统状态误差和输出跟踪误差的扩展状态,δy(t,k)是t时刻第k周期注塑成型批次方向的输出误差, A、I、E、F是适当维度的常量矩阵,G(t,k)是t时刻第k周期的常量矩阵。
1-6.以预测值形式表示的注塑成型闭环渐进稳定系统模型,形式如下:
其中,j是预测值,H1、H2是对应的更新律增益,z(t+j+1|t,k)、 z(t+j+1|t,k-1)分别是t时刻到t+j+1时刻第k周期、第k-1周期注塑成型的预测状态, z(t+j|t,k)是t时刻到t+j时刻第k周期注塑成型的预测状态,是t时刻到 t+j时刻第k周期注塑成型的批次方向外部干扰误差,δx(t+j|t,k)、δy(t+j|t,k)分别是t时刻到t+j时刻第k周期注塑成型的批次方向系统状态误差、输出误差。
步骤2、设计喷嘴压力的注塑成型过程控制器,具体是:
2-1.基于步骤1,鲁棒性能保证控制下的注塑成型过程的更新律形式如下:
2-2.在重复性和非重复性扰动下,均可得到注塑成型过程的增益矩阵控制律的形式如下所示:
H1=γ-1Y1P,H2=γ-1Y2P
其中,P是满足系统渐进稳定的对称正定矩阵,Y1,Y2是满足系统条件的矩阵,γ-1>0是满足系统条件的系数。结合步骤1-2和步骤2-1即可得到最优系统输入u(t,k)。
2-3.在下一时刻,重复步骤1.6到2.2得到比例阀的阀门开度u(t,k),再将其作用于喷嘴压力,并依次循环。
Claims (3)
1.一种化工批次过程鲁棒学习控制方法,包括如下步骤:
步骤1、建立以预测值形式表示的闭环系统模型;
步骤2、设计被控对象的批次处理控制器。
2.如权利要求1所述的化工批次过程鲁棒学习控制方法,其特征在于:
步骤1具体如下:
1-1.建立具有不确定参数扰动的离散时间模型,形式如下:
其中,t和k分别是运行时刻和运行周期,x0,k是第k周期批次处理运行的初始条件,u(t,k)和y(t,k)分别是t时刻第k周期系统的输入和输出,x(t,k)、x(t+1,k)分别是t时刻、t+1时刻第k周期的系统状态,w(t,k)是t时刻第k周期未知的外部干扰,B、C是具有适当维数的常数矩阵;
1-2.结合步骤1-1,引入鲁棒学习控制量:
u(t,k)=u(t,k-1)+r(t,k) u(t,0)=0,t=0,1,2,…
其中,u(t,0)是t时刻系统的初始输入,u(t,k)、u(t,k-1)分别是t时刻第k周期、第k-1周期的系统输入,r(t,k)是t时刻第k周期的更新律;
1-3.定义输出跟踪误差和误差的批次方向函数:
e(t,k)=y(t,k)-yr(t)
δf(t,k)=f(t,k)-f(t,k-1)
其中,yr(t)是输出的参考值,e(t,k)是t时刻第k周期的输出跟踪误差,f(t,k)是t时刻第k周期的系统状态、输出变量或未知的外部干扰,f(t,k-1)是t时刻第k-1周期的系统状态、输出变量或未知的外部干扰,δf(t,k)是t时刻第k周期误差的批次方向函数;
1-4.结合步骤1-1至步骤1-3,得到:
其中,δx(t,k)、δx(t+1,k)是t时刻、t+1时刻第k周期批次方向的系统状态误差,是t时刻第k周期批次方向的外部干扰误差且ΔA(t,k)是不确定的系统矩阵,x(t,k-1)是t时刻第k-1周期的系统状态,δw(t,k)是t时刻第k周期的外部干扰,e(t+1,k)是t+1时刻第k周期的输出跟踪误差,e(t+1,k-1)是t+1时刻第k-1周期的输出跟踪误差;
1-5.得到鲁棒渐进稳定的系统模型,形式如下:
其中,C1=[C 0],z(t+1,k)、z(t,k)是t+1时刻、t时刻第k周期综合批次方向的系统状态误差和输出跟踪误差的扩展状态,z(t+1,k-1)是t+1时刻第k-1周期综合批次方向的系统状态误差和输出跟踪误差的扩展状态,δy(t,k)是t时刻第k周期批次方向的输出误差,A、I、E、F是适当维度的常量矩阵,G(t,k)是t时刻第k周期的常量矩阵;
1-6.以预测值形式表示的闭环渐进稳定系统模型,形式如下:
其中,j是预测值,H1、H2是对应的更新律增益,z(t+j+1|t,k)、z(t+j+1|t,k-1)分别是t时刻到t+j+1时刻第k周期、第k-1周期的预测状态,z(t+j|t,k)是t时刻到t+j时刻第k周期的预测状态,是t时刻到t+j时刻第k周期批次方向的外部干扰误差,δx(t+j|t,k)、δy(t+j|t,k)分别是t时刻到t+j时刻第k周期批次方向的系统状态误差、输出误差。
3.如权利要求2所述的化工批次过程鲁棒学习控制方法,其特征在于:
步骤2具体如下:
2-1.基于步骤1,鲁棒性能保证控制下的更新律形式如下:
2-2.在重复性和非重复性扰动下,均可得到增益矩阵控制律的形式如下所示:
H1=γ-1Y1P,H2=γ-1Y2P
其中,P是满足系统渐进稳定的对称正定矩阵,Y1,Y2是满足系统条件的矩阵,γ-1>0是满足系统条件的系数。结合步骤1-2至步骤2-1即可得到最优系统输入u(t,k);
2-3.在下一时刻,重复步骤1-6到2-2得到新的最优系统输入u(t,k),再将其作用于控制对象,并依次循环。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810760548.4A CN108681317A (zh) | 2018-07-11 | 2018-07-11 | 一种化工工业过程鲁棒学习控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810760548.4A CN108681317A (zh) | 2018-07-11 | 2018-07-11 | 一种化工工业过程鲁棒学习控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108681317A true CN108681317A (zh) | 2018-10-19 |
Family
ID=63814039
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810760548.4A Pending CN108681317A (zh) | 2018-07-11 | 2018-07-11 | 一种化工工业过程鲁棒学习控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108681317A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110826021A (zh) * | 2019-10-31 | 2020-02-21 | 哈尔滨工业大学 | 一种非线性工业过程鲁棒辨识和输出估计方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107966902A (zh) * | 2017-11-27 | 2018-04-27 | 辽宁石油化工大学 | 一种不确定性间歇过程的约束2d跟踪控制方法 |
-
2018
- 2018-07-11 CN CN201810760548.4A patent/CN108681317A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107966902A (zh) * | 2017-11-27 | 2018-04-27 | 辽宁石油化工大学 | 一种不确定性间歇过程的约束2d跟踪控制方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110826021A (zh) * | 2019-10-31 | 2020-02-21 | 哈尔滨工业大学 | 一种非线性工业过程鲁棒辨识和输出估计方法 |
CN110826021B (zh) * | 2019-10-31 | 2021-03-12 | 哈尔滨工业大学 | 一种非线性工业过程鲁棒辨识和输出估计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107966902B (zh) | 一种不确定性间歇过程的约束2d跟踪控制方法 | |
CN107942667B (zh) | 基于时变时滞和干扰的注塑过程混杂2d跟踪控制方法 | |
CN105911868B (zh) | 多批次间歇反应器二维迭代学习反馈控制方法 | |
CN104765350B (zh) | 基于组合模型预测控制技术的水泥分解炉控制方法及系统 | |
CN109557810B (zh) | 一种基于新型二自由度内模pid的加热炉温度控制方法 | |
CN101334676B (zh) | 小型聚合釜dcs生产聚氯乙烯反应温度的控制方法 | |
CN109407512B (zh) | 依赖时滞的间歇过程2d输入输出约束控制方法 | |
CN110764414A (zh) | 针对多种干扰的多阶段批次异步切换过程的鲁棒预测控制方法 | |
CN108681317A (zh) | 一种化工工业过程鲁棒学习控制方法 | |
CN103412486A (zh) | 一种聚氯乙烯汽提过程智能控制方法 | |
CN105353618B (zh) | 一种批次注塑过程的约束跟踪控制方法 | |
CN104527009A (zh) | 基于高斯预测的注塑成型工艺参数优化调控方法 | |
CN109212972B (zh) | 间歇过程的受限滚动时域混杂2d跟踪控制方法 | |
CN110412873A (zh) | 基于终端约束的时滞间歇过程2d迭代学习预测控制方法 | |
CN113641101B (zh) | 一种基于数值仿真的多渠池控制参数寻优算法 | |
CN110794672A (zh) | 一种水泥生产过程分解炉炉温显式控制方法 | |
CN106610588B (zh) | 一种串级预测控制系统及方法 | |
CN108803315B (zh) | 一种化工间歇过程的智能跟踪控制方法 | |
CN108803314B (zh) | 一种化工间歇过程的新型复合跟踪控制方法 | |
CN109491245A (zh) | 一种cstr系统的扰动补偿控制方法 | |
CN112506050A (zh) | 一种基于潜变量过程迁移模型的间歇过程集成优化方法 | |
CN108897219B (zh) | 一种化工不确定工业过程约束预测控制方法 | |
CN117193407B (zh) | 甲醇氧化制甲醛的控制方法、电子设备及存储介质 | |
CN108008653A (zh) | 一种面向重复运动的双向控制系统及加工设备 | |
CN114311574B (zh) | 一种注塑机注射速度的优化控制方法、系统及装置 |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20181019 |