CN109100940A - 一种基于gap metric加权函数的多模型预测控制方法 - Google Patents
一种基于gap metric加权函数的多模型预测控制方法 Download PDFInfo
- Publication number
- CN109100940A CN109100940A CN201811137302.8A CN201811137302A CN109100940A CN 109100940 A CN109100940 A CN 109100940A CN 201811137302 A CN201811137302 A CN 201811137302A CN 109100940 A CN109100940 A CN 109100940A
- Authority
- CN
- China
- Prior art keywords
- model
- weighting function
- model predictive
- gap metric
- sub
- 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
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 238000000034 method Methods 0.000 claims description 22
- 238000005457 optimization Methods 0.000 claims description 4
- 238000010586 diagram Methods 0.000 description 6
- 238000003786 synthesis reaction Methods 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 239000002826 coolant Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 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
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于gap metric加权函数的多模型预测控制方法,利用gap metric提出一种系统、有效的多模型加权算法,对各个子模型预测控制器进行合成,得到多模型预测控制器,对非线性系统进行优化控制。和传统加权函数,比如梯形加权函数相比,该加权函数只有一个整定参数,大大减少了参数整定的工作量,而且权值可离线计算,并存储在查询表里备用,大大减低了在线计算量;本发明提供的算法保留了传统加权算法的优点,而克服了其缺点,是一种有效的多模型加权算法,对提高多模型预测控制器的性能大有裨益。
Description
技术领域
本发明涉及一种基于gap metric加权函数的多模型预测控制方法,属于非线性系统多模型控制技术领域。
背景技术
近年来,多模型预测控制方法越来越多的用于具有宽操作空间的非线性系统的控制。基于分解-合成原则的多模型预测控制方法融合了多模型方法和预测控制的优点。一方面,多模型预测控制方法将复杂的非线性控制问题转化为一系列简单的线性控制问题;另一方面,可以将控制系统中的约束条件直接包含在目标函数中。此外,由于采用了滚动优化机制,多模型预测控制的闭环控制性能得到了进一步改善。
多模型预测控制方法中,通常采用硬切换和加权两种方式进行子控制器的合成。虽然硬切换方式可以选择最合适的模型控制器切入到闭环系统中,然而硬切换很容易导致系统输出抖动,以至于不满足闭环控制要求。相比之下,加权合成方式可以使得系统输出平滑,避免输出抖动。因此加权函数合成方式更受欢迎。
传统的加权函数,比如梯形加权函数、高斯加权函数等虽然形式简单,且可离线计算,然而整定参数的个数随着子模型个数的增加而增加,给参数整定带来巨大的挑战。Gapmetric近年来越来越多的用于多模型控制方法中,比如用于选择模型集及子模型个数精简等,取得了非常好的效果。因此,充分发挥gap metric的优点,构造简单方便的加权函数对子模型预测控制进行加权合成,将有利于多模型预测控制方法性能的进一步提高。
发明内容
目的:为了克服现有多模型预测控制方法中加权函数存在的缺点和不足,本发明提供一种基于gap metric的加权函数,对各个子模型预测控制器进行加权,合成多模型预测控制器,对非线性系统进行优化控制。
技术方案:为解决上述技术问题,本发明采用的技术方案为:
一种基于gap metric加权函数的多模型预测控制方法,包括如下步骤,
S1.根据非线性系统的特性,选择合适的能够反映系统操作条件的调度变量θ;
S2.选择m个线性子模型P1,P2,…,Pm构成模型集对非线性系统进行近似;
S3.基于每个子模型P1,P2,…,Pm设计对应的子模型预测控制器K1,K2,…,Km;
S4.t时刻,非线性系统的模型记为nPt,其线性化模型记为P(θt);则非线性系统nPt与线性子模型Pi之间的gap metric距离定义为:
γi(θt)=δ(Pi,P(θt)),i=1,...,m 公式一;
函数δ(P1,P2)表示线性系统P1和P2之间的gap metric距离;
第i个子控制器Ki在t时刻时的加权函数根据公式二计算,
其中,kg是加权函数唯一的整定参数;
S5.为了减少在线计算量,对步骤S4定义的加权函数,采用离线计算的方式,并将计算得到的权值存储在查询表中备用;
S5.1对调度变量θ进行网格化,假设共得到Ng个网格点,即得到θ1,θ2,...,θj,...,θNg;计算非线性系统在每个网格点的线性化模型LP1,LP2,...,LPj,...,LPNg;
S5.2计算网格点θj处非线性系统与第i个子模型Pi的gap metric距离γi(θj)=δ(Pi,LPj),进而计算第i个线性模型预测控制器Ki在系统运行到θj处时对应的权值其中j=1,2,...,Ng,且i=1,...,m;
S5.3遍历θj,j=1,2,...,Ng,计算每个控制器在整个调度变量范围内的权值,并以查询表的形式储存,得到查询表Tab_w1,Tab_w2,...,Tab_wm。
S6.子模型预测控制器K1,K2,...,Km在t时刻的输出记为u1(t),u2(t),...,um(t);
S7.多模型预测控制器在t时刻的输出根据公式三计算得出,并对非线性系统进行优化控制;
其中,θt为t时刻时调度变量对应的值,可直接在查询表中查询θ=θt时m个预测控制器对应的权值w1(θt),w2(θt),...,wm(θt);若查询表中没有值和θt对应,则可以采用线性差值的方法求取对应的权值。
有益效果:与现有技术相比,本发明提供的技术方案具有以下优点和积极效果:
1)加权函数只有一个整定参数,在线整定方便快捷;
2)权值可离线计算并储存在查询表中备用,大大减低了在线计算量;
3)提高了闭环系统的控制性能。
附图说明
图1为本发明提供的基于gap metric的多模型加权函数的示意图(假设3个子模型);
图2为基于gap metric加权函数的多模型预测控制结构图;
图3为CSTR系统在基于梯形加权函数的多模型预测控制方法下的闭环仿真图;
图4为CSTR系统在本发明提供的多模型预测控制方法下的闭环仿真图。
具体实施方式
下面结合附图和实施例对本发明作更进一步的说明。
实施例1,针对一个连续搅拌反应釜系统(CSTR)进行仿真并加以分析。
CSTR系统:
其中,x1和x2是无量纲的反应率和反应釜内温度,输入变量u是无量纲的冷却剂温度。方程中的常数Da=0.072,γ=20,B=8及β=0.3。此系统非线性程度很强,单个线性控制器无法满足要求。
利用本发明提供的基于gap metric加权函数的多模型预测控制方法对CSTR系统进行控制,步骤如下:
S1.CSTR系统的输出y的变化直接反映了系统的非线性程度,因此将y选为调度变量;该系统的操作空间为{y|y∈[0,6]}。
S2.选择3个线性模型对CSTR系统进行近似,分别为
S3.基于子模型P1,P2,P3分别设计线性模型预测控制器K1,K2,K3。
S4步骤省略,直接进行步骤S5的离线计算并建立CSTR系统的权值查询表Tab_w1,Tab_w2,Tab_w3。
S5.1对CSTR系统的调度变量y,利用42个静态点对其进行网格化,并线性化得到每个网格点处的线性化模型LPj,j=1,2,...,42。
S5.2计算网格点θj处非线性系统与第i个子模型Pi的gap metric距离γi(θj)=δ(Pi,LPj),进而计算第i个线性模型预测控制器Ki在系统运行到θj处时对应的权值其中j=1,2,...,42,i=1,2,3,kg=6。
S5.3遍历θj,j=1,2,...,42,建立CSTR系统的权值的查询表Tab_w1,Tab_w2,Tab_w3。
S6.子模型预测控制器K1,K2,K3在t时刻的输出记为u1(t),u2(t),u3(t)。
S7.多模型预测控制器在t时刻的输出根据公式计算得出,并对非线性系统进行优化控制。
图1为本发明提供的一种基于gap metric的多模型加权函数的示意图(假设3个子模型),图2为基于gap metric加权函数的多模型预测控制结构图。
图3给出了CSTR系统在基于梯形加权函数的多模型预测控制器下的输入和输出值,显然系统的输出不能够很好的在整个操作范围内跟踪参考信号的变化,存在较明显的静态误差,且存在输出抖动的情况。图4是CSTR系统在本发明提供的基于gap metric加权函数的多模型预测控制器下的输入和输出,显然输出跟踪参考信号的效果比前者好很多,且跟踪精度高、速度快,输出既平滑又准确。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (6)
1.一种基于gap metric加权函数的多模型预测控制方法,其特征在于:包括如下步骤,
S1.根据非线性系统的特性,选择合适的能够反映系统操作条件的调度变量θ;
S2.选择m个线性子模型P1,P2,…,Pm构成模型集对非线性系统进行近似;
S3.基于每个子模型P1,P2,…,Pm设计对应的子模型预测控制器K1,K2,…,Km;
S4.t时刻,非线性系统的模型记为nPt,其线性化模型记为P(θt);则非线性系统nPt与线性子模型Pi之间的gap metric距离定义为:
γi(θt)=δ(Pi,P(θt)),i=1,…,m 公式一;
第i个子控制器Ki在t时刻时的加权函数根据公式二计算,
S5.为了减少在线计算量,对步骤S4定义的加权函数,采用离线计算的方式,并将计算得到的权值存储在查询表中备用;
S6.子模型预测控制器K1,K2,…,Km在t时刻的输出记为u1(t),u2(t),…,um(t);
S7.多模型预测控制器在t时刻的输出根据公式三计算得出,并对非线性系统进行优化控制;
其中,θt为t时刻时调度变量对应的值,可直接在查询表中查询θ=θt时m个预测控制器对应的权值w1(θt),w2(θt),…,wm(θt)。
2.根据权利要求1所述的基于gap metric加权函数的多模型预测控制方法,其特征在于:步骤S5所述的离线计算方式的步骤如下,
S5.1对调度变量θ进行网格化,假设共得到Ng个网格点,即得到θ1,θ2,…,θj,…,θNg;计算非线性系统在每个网格点的线性化模型LP1,LP2,…,LPj,…,LPNg;
S5.2计算网格点θj处非线性系统与第i个子模型Pi的gap metric距离γi(θj)=δ(Pi,LPj),进而计算第i个线性模型预测控制器Ki在系统运行到θj处时对应的权值
S5.3遍历θj,j=1,2,…,Ng,计算每个控制器在整个调度变量范围内的权值,并以查询表的形式储存,得到查询表Tab_w1,Tab_w2,…,Tab_wm。
3.根据权利要求1所述的基于gap metric加权函数的多模型预测控制方法,其特征在于:步骤S4中函数δ(P1,P2)表示线性系统P1和P2之间的gap metric距离。
4.根据权利要求1所述的基于gap metric加权函数的多模型预测控制方法,其特征在于:步骤S4的公式二中,kg是加权函数唯一的整定参数。
5.根据权利要求2所述的基于gap metric加权函数的多模型预测控制方法,其特征在于:步骤S5.2中,且i=1,…,m。
6.根据权利要求1所述的基于gap metric加权函数的多模型预测控制方法,其特征在于:步骤S7中,若查询表中没有值和θt对应,则可以采用线性差值的方法求取对应的权值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811137302.8A CN109100940A (zh) | 2018-09-28 | 2018-09-28 | 一种基于gap metric加权函数的多模型预测控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811137302.8A CN109100940A (zh) | 2018-09-28 | 2018-09-28 | 一种基于gap metric加权函数的多模型预测控制方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109100940A true CN109100940A (zh) | 2018-12-28 |
Family
ID=64867563
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811137302.8A Pending CN109100940A (zh) | 2018-09-28 | 2018-09-28 | 一种基于gap metric加权函数的多模型预测控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109100940A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110442027A (zh) * | 2019-08-16 | 2019-11-12 | 河海大学常州校区 | 一种gap多模型加权函数参数自整定方法 |
CN110658722A (zh) * | 2019-10-18 | 2020-01-07 | 河海大学常州校区 | 一种基于gap的自均衡多模型分解方法及系统 |
CN111694279A (zh) * | 2020-06-30 | 2020-09-22 | 河海大学常州校区 | 多变量非线性系统自适应均衡多模型分解及控制方法 |
CN112415886A (zh) * | 2019-08-21 | 2021-02-26 | 河海大学常州校区 | 一种基于pid的一体化多模型控制方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106054618A (zh) * | 2016-08-20 | 2016-10-26 | 南京理工大学 | 含输入间隙的电机伺服系统自适应鲁棒控制方法 |
US20170351226A1 (en) * | 2016-06-01 | 2017-12-07 | Rockwell Automation Technologies, Inc. | Industrial machine diagnosis and maintenance using a cloud platform |
CN107450325A (zh) * | 2017-09-06 | 2017-12-08 | 东南大学 | 一种燃烧后co2捕集系统的多模型预测控制方法 |
-
2018
- 2018-09-28 CN CN201811137302.8A patent/CN109100940A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170351226A1 (en) * | 2016-06-01 | 2017-12-07 | Rockwell Automation Technologies, Inc. | Industrial machine diagnosis and maintenance using a cloud platform |
CN106054618A (zh) * | 2016-08-20 | 2016-10-26 | 南京理工大学 | 含输入间隙的电机伺服系统自适应鲁棒控制方法 |
CN107450325A (zh) * | 2017-09-06 | 2017-12-08 | 东南大学 | 一种燃烧后co2捕集系统的多模型预测控制方法 |
Non-Patent Citations (2)
Title |
---|
杜静静: "Comparison of two gap-metric-based weighting methods in multilinear", 《2017 36TH CHINESE CONTROL CONFERENCE(CCC)》 * |
陈垣君: "基于LPV模型非线性预测控制的精细化建模及精确求解", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110442027A (zh) * | 2019-08-16 | 2019-11-12 | 河海大学常州校区 | 一种gap多模型加权函数参数自整定方法 |
CN112415886A (zh) * | 2019-08-21 | 2021-02-26 | 河海大学常州校区 | 一种基于pid的一体化多模型控制方法 |
CN112415886B (zh) * | 2019-08-21 | 2023-08-29 | 河海大学常州校区 | 一种基于pid的一体化多模型控制方法 |
CN110658722A (zh) * | 2019-10-18 | 2020-01-07 | 河海大学常州校区 | 一种基于gap的自均衡多模型分解方法及系统 |
CN110658722B (zh) * | 2019-10-18 | 2022-04-26 | 河海大学常州校区 | 一种基于gap的自均衡多模型分解方法及系统 |
CN111694279A (zh) * | 2020-06-30 | 2020-09-22 | 河海大学常州校区 | 多变量非线性系统自适应均衡多模型分解及控制方法 |
CN111694279B (zh) * | 2020-06-30 | 2022-09-23 | 河海大学常州校区 | 多变量非线性系统自适应均衡多模型分解控制方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109100940A (zh) | 一种基于gap metric加权函数的多模型预测控制方法 | |
CN102023574B (zh) | 一阶反应连续搅拌釜式反应器的混杂模型优化控制方法 | |
Odgaard et al. | On using Pareto optimality to tune a linear model predictive controller for wind turbines | |
US10281897B2 (en) | Model predictive control with uncertainties | |
CN104375475B (zh) | 间歇反应器中间歇反应过程的优化控制方法 | |
CN110376900B (zh) | 参数优化方法、终端设备以及计算机存储介质 | |
CN111123871B (zh) | 针对化工过程遗传算法优化的预测函数控制方法 | |
CN104698842B (zh) | 一种基于内点法的lpv模型非线性预测控制方法 | |
CN111781835B (zh) | 一种镇定二阶惯性加纯滞后系统的线性自抗扰控制器设计方法 | |
WO2019176370A1 (ja) | 将来状態推定装置および将来状態推定方法 | |
CN105388764A (zh) | 基于动态矩阵前馈预测的电液伺服pid控制方法及系统 | |
CN101452261B (zh) | 一种聚丙烯装置牌号切换控制方法 | |
CN106527125A (zh) | 智能控制中的无模型控制方法 | |
CN105404144A (zh) | 连续搅拌反应釜的多模型自适应控制方法及系统 | |
Tian et al. | Hybrid modeling of molten steel temperature prediction in LF | |
CN107678276B (zh) | 一种基于转台控制的自适应复合控制方法 | |
CN110094838B (zh) | 一种基于空调系统的可变参数无模型自适应控制方法 | |
Peng et al. | Robust switched predictive control for multi-phase batch processes with uncertainties and unknown disturbances | |
Xiao et al. | A new discrete variable structure control algorithm based on sliding mode prediction | |
CN104111605B (zh) | 单输入单输出非自衡生产过程的控制器及控制方法 | |
CN113110519A (zh) | 一种舰船用非增量型无模型自适应艏向控制方法 | |
CN106647250B (zh) | 基于离线优化/在线查表方式的双层结构预测控制方法 | |
CN110442027A (zh) | 一种gap多模型加权函数参数自整定方法 | |
CN111694279B (zh) | 多变量非线性系统自适应均衡多模型分解控制方法 | |
CN112415886B (zh) | 一种基于pid的一体化多模型控制方法 |
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: 20181228 |