CN110825051B - 一种基于gap metric的不确定系统的多模型控制方法 - Google Patents

一种基于gap metric的不确定系统的多模型控制方法 Download PDF

Info

Publication number
CN110825051B
CN110825051B CN201911113348.0A CN201911113348A CN110825051B CN 110825051 B CN110825051 B CN 110825051B CN 201911113348 A CN201911113348 A CN 201911113348A CN 110825051 B CN110825051 B CN 110825051B
Authority
CN
China
Prior art keywords
model
max
linear
gap
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.)
Active
Application number
CN201911113348.0A
Other languages
English (en)
Other versions
CN110825051A (zh
Inventor
杜静静
陈俊风
贾静
李建
姜学平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changzhou Campus of Hohai University
Original Assignee
Changzhou Campus of Hohai University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changzhou Campus of Hohai University filed Critical Changzhou Campus of Hohai University
Priority to CN201911113348.0A priority Critical patent/CN110825051B/zh
Publication of CN110825051A publication Critical patent/CN110825051A/zh
Application granted granted Critical
Publication of CN110825051B publication Critical patent/CN110825051B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/418Total 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/41885Total 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
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/32Operator till task planning
    • G05B2219/32339Object oriented modeling, design, analysis, implementation, simulation language
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total 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

本发明公开一种基于gap metric的不确定系统的多模型控制方法。首先,对系统的不确定参数网格化,在每个网格点建立线性模型。然后,利用gap metric对这些线性模型进行归类,得到不确定系统的子模型集。最后基于子模型集设计多模型控制器,对原不确定系统进行优化控制。本发明的方法可以有效减少子模型集的建立对先验知识的依赖,减少子模型集的冗余,简化多模型控制器的结构。这对于提高多模型控制系统的性能,改善不确定性系统的闭环性能大有裨益。

Description

一种基于gap metric的不确定系统的多模型控制方法
技术领域
本发明公开了一种基于gap metric的不确定系统的多模型控制方法,属于不确定系统的多模型控制领域。
背景技术
实际工业生产中往往需要考虑很多不确定因素,比如系统故障,参数变化,外部干扰等等。由于这些不确定因素的存在,使得对系统的控制变得十分复杂。由于传统控制方法是基于被控对象的精确模型对系统进行分析和设计的,因此难以对具有不确定性的系统进行控制和优化。
发明内容
为了克服上述不足,本发明提供了一种基于gap metric的不确定系统的多模型控制方法,多模型方法的特点是分而治之,即将复杂的系统分解成一系列简单的问题,通过对简单问题的求解,得到实现对原始复杂问题的求解。多模型方法有效的简化了控制问题的复杂性,降低了控制难度。因此多模型控制方法在控制领域得到了广泛的关注,非常适用于解决具有不确定性的系统的控制问题。
多模型控制方法主要包含:多模型分解(子模型集建立)、子控制器设计、以及多模型合成三个关键步骤。传统的多模型分解方法有:按照系统的物理组成单元进行分解;按照系统内在的物理、化学规律分解;按照系统数学模型的展开式进行分解;按照系统的控制目标进行分解;按照系统的工况进行分解。此类分解方法往往特别依赖系统的定性知识、操作经验等先验知识,或者需要从系统的输入-输出数据中获取信息,缺少系统的指导原则。近年来,学者们将gap metric引入多模型控制,用于对非线性系统进行分解并设计多模型控制器,可有效减少对先验知识的过分依赖,减少模型冗余,简化模型控制器的结构。
本发明的技术方案如下:
一种基于gap metric的不确定系统的多模型控制方法,包括以下步骤:
(1)、对系统的不确定参数网格化,在每个网格点建立线性模型;
(2)、然后,利用gap metric对这些线性模型进行归类,得到近似原不确定系统的子模型集;
(3)、最后基于子模型集设计多模型控制器,对原不确定系统进行优化控制。
上述步骤(2)的具体步骤如下:
(2-1)、假如网格点为n,在每个网格点建立线性模型,根据公式(1)计算线性模型之间的gap matrix,[δij]n×n
Figure BDA0002273379330000021
其中Pi表示第i个线性模型,Pj表示第j个线性模型,M1、N1是Pi的规范化互质分解因子,M2、N2是Pj的规范化互质分解因子,Q是H空间的有理数,inf表示下确界;
(2-2)、选择分解的阈值γ(范围0.4~0.6),令i=1;
(2-3)、如果i≤n,令j=i,并跳转到步骤(2-4),否则跳转到步骤(2-11)
(2-4)、根据min-max原则计算子模型P*,min-max原则如公式(2)所示:
P*:={Ph:mini≤h≤j(maxi≤l≤j(δ(Ph,Pl)))} (2)
h和l是i到j之间的一个自然数,Ph表示第h个线性模型,Pl表示第l个线性模型;maxi≤l≤j(·)表示在h不变的情况下,当l从i变化到j时,求解和线性模型Ph的gap metric距离最大的那个线性型及其距离;mini≤h≤j(·)表示在l不变的情况下,当h从i变化到j时,求解和线性模型Pl的gap metric距离最小的那个线性模型及其距离;
(2-5)、根据公式(3)计算最大gap
δmax:=maxi≤h≤j(δ(P*,Ph)) (3)
maxi≤h≤j(δ(P*,Ph))表示i到j之间的线性模型与P*距离最大的gap metric距离的值;
(2-6)、如果δmax≤γ且j<n,则令j=j+1并返回步骤(2-4)反之,若δmax>γ,令Setj=j–1,进入步骤(2-7);
(2-7)、更新子模型,即根据min-max原则更新第i到第j个线性模型中的子模型;
(2-8)、子模型P*记为Gmc,其对应的网格点记为OPmc.将Gmc放入队列pQc,并将OPmc放入队列oQc
(2-9)、令i=j+1,并返回步骤(2-3);
(2-10)、分解结束;
(2-11)、得到近似原不确定系统的子模型集。
本发明的有益效果如下:
基于gap metric的归类方法可以有效减少模型个数,同时减少对先验知识的依赖,并避免模型冗余,简化多模型控制器的结构。然后,基于子模型集设计多模型控制器对系统进行优化控制。这对于改善不确定系统的多模型控制效果大有裨益。
本发明的方法可以有效减少子模型集的建立对先验知识的依赖,减少子模型集的冗余,简化多模型控制器的结构。这对于提高多模型控制系统的性能,改善不确定性系统的闭环性能大有裨益。
附图说明
图1是实施例中30个线性化模型两两之间的gap距离;
图2是实施例中在A=20时在本发明的多模型控制器下的控制效果;
图3是实施例中在A=10时在本发明的多模型控制器下的控制效果。
具体实施方式
为使得本发明的发明目的、特征、优点能够更加的明显和易懂,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,下面所描述的实施例仅仅是本发明一部分实施例,而非全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
下面结合附图并通过具体实施方式来进一步说明本发明的技术方案。
实施例:
对于系统
Figure BDA0002273379330000031
A的变化范围是[1,30],设计两个线性PID控制器,构成多模型控制器对其进行控制。
S1:选择分解阈值γ=0.4。
S2:对不确定参数网格化,A被平均分成30等分。在各个网格点建立系统的线性模型。并计算线性模型之间的gap矩阵。
S3:令i=1。
S4:如果i≤30,令j=i,并跳转到S5。否则跳转到S11。
S5:根据min-max原则计算子模型。
S6:计算最大gap。
S7:如果δmax≤γ且j<30,则令j=j+1并返回S4.反之,若δmax>γ,令Set j=j–1,进入S8。
S8:更新子模型,即根据min-max原则更新第i到第j个线性模型中的子模型。
S9:子模型P*记为Gmc,其对应的网格点记为OPmc.将Gmc放入队列pQc,并将OPmc放入队列oQc
S10:令i=j+1,并返回S4。
S11:分解结束。
S12:基于得到的子模型集设计多模型PID控制器对被控系统进行控制。
最后,系统被分成了2个子系统。对应的标称模型分别是第6个和第21个模型。分解结果详见表1。
表1.案例分解的结果
Figure BDA0002273379330000041
30个线性模型之间的gap距离如图1所示。
基于此分解结果设计多模型PID控制器,闭环控制效果如图2和图3所示。两种情况下输出y都可以在20s之内整定完毕。图2系统输出有超调,但超调量不超过5%,而图3系统超调非常小。因此闭环系统的控制性能非常好。

Claims (1)

1.一种基于gap metric的不确定系统的多模型控制方法,其特征在于,包括以下步骤:
(1)、对系统的不确定参数网格化,在每个网格点建立线性模型;
(2)、然后,利用gap metric对这些线性模型进行归类,得到近似原不确定系统的子模型集;
(3)、最后基于子模型集设计多模型控制器,对原不确定系统进行优化控制;
所述步骤(2)的具体步骤如下:
(2-1)、假如网格点为n,在每个网格点建立线性模型,根据公式(1)计算线性模型之间的gap matrix,[δij]n×n
Figure FDA0003931495810000011
其中Pi表示第i个线性模型,Pj表示第j个线性模型,M1、N1是Pi的规范化互质分解因子,M2、N2是Pj的规范化互质分解因子,Q是H空间的有理数,inf表示下确界;
(2-2)、选择分解的阈值γ,γ的取值范围:0.4~0.6,令i=1;
(2-3)、如果i≤n,令j=i,并跳转到步骤(2-4),否则跳转到步骤(2-11);
(2-4)、根据min-max原则计算子模型P*,min-max原则如公式(2)所示:
P*:={Ph:mini≤h≤j(maxi≤l≤j(δ(Ph,Pl)))} (2)
h和l是i到j之间的一个自然数,Ph表示第h个线性模型,Pl表示第l个线性模型;maxi≤l≤j(·)表示在h不变的情况下,当l从i变化到j时,求解和线性模型Ph的gap metric距离最大的那个线性型及其距离;mini≤h≤j(·)表示在l不变的情况下,当h从i变化到j时,求解和线性模型Pl的gap metric距离最小的那个线性模型及其距离;
(2-5)、根据公式(3)计算最大gap
δmax:=maxi≤h≤j(δ(P*,Ph)) (3)
maxi≤h≤j(δ(P*,Ph))表示i到j之间的线性模型与P*距离最大的gap metric距离的值;
(2-6)、如果δmax≤γ且j<n,则令j=j+1并返回步骤(2-4)反之,若δmax>γ,令Set j=j–1,进入步骤(2-7);
(2-7)、更新子模型,即根据min-max原则更新第i到第j个线性模型中的子模型;
(2-8)、子模型P*记为Gmc,其对应的网格点记为OPmc,将Gmc放入队列pQc,并将OPmc放入队列oQc
(2-9)、令i=j+1,并返回步骤(2-3);
(2-10)、分解结束;
(2-11)、得到近似原不确定系统的子模型集。
CN201911113348.0A 2019-11-14 2019-11-14 一种基于gap metric的不确定系统的多模型控制方法 Active CN110825051B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911113348.0A CN110825051B (zh) 2019-11-14 2019-11-14 一种基于gap metric的不确定系统的多模型控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911113348.0A CN110825051B (zh) 2019-11-14 2019-11-14 一种基于gap metric的不确定系统的多模型控制方法

Publications (2)

Publication Number Publication Date
CN110825051A CN110825051A (zh) 2020-02-21
CN110825051B true CN110825051B (zh) 2023-02-14

Family

ID=69555340

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911113348.0A Active CN110825051B (zh) 2019-11-14 2019-11-14 一种基于gap metric的不确定系统的多模型控制方法

Country Status (1)

Country Link
CN (1) CN110825051B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111694279B (zh) * 2020-06-30 2022-09-23 河海大学常州校区 多变量非线性系统自适应均衡多模型分解控制方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104981699A (zh) * 2012-12-10 2015-10-14 安帕奇蒙股份有限公司 用于测量垂直风分量的方法和系统
CN105867119A (zh) * 2016-01-15 2016-08-17 南京航空航天大学 一种采用保护映射理论的空天飞行器大包线切换控制方法
CN108614431A (zh) * 2018-06-08 2018-10-02 河海大学常州校区 一种基于夹角的Hammerstein-Wiener系统多模型分解及控制方法
CN109063277A (zh) * 2018-07-12 2018-12-21 佛山科学技术学院 一种基于间隙度量的动态模式识别方法及装置

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8494988B2 (en) * 2009-05-15 2013-07-23 The Aerospace Corporation Systems and methods for generating feasible solutions from two parents for an evolutionary process

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104981699A (zh) * 2012-12-10 2015-10-14 安帕奇蒙股份有限公司 用于测量垂直风分量的方法和系统
CN105867119A (zh) * 2016-01-15 2016-08-17 南京航空航天大学 一种采用保护映射理论的空天飞行器大包线切换控制方法
CN108614431A (zh) * 2018-06-08 2018-10-02 河海大学常州校区 一种基于夹角的Hammerstein-Wiener系统多模型分解及控制方法
CN109063277A (zh) * 2018-07-12 2018-12-21 佛山科学技术学院 一种基于间隙度量的动态模式识别方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
网络信息模式下大范围工况运行系统的多模型协同预测控制;亓晓雯等;《信息与控制》;20181231;第7页 *

Also Published As

Publication number Publication date
CN110825051A (zh) 2020-02-21

Similar Documents

Publication Publication Date Title
CN102662321B (zh) 一种主元分析监控模型的在线更新方法
CN108227494A (zh) 非线性批次过程2d最优约束模糊容错控制方法
EP2196872A1 (en) A system for gain scheduling control
Chen et al. H/sub/spl infin//output feedback control of discrete-time fuzzy systems with application to chaos control
Karg et al. Learning-based approximation of robust nonlinear predictive control with state estimation applied to a towing kite
CN107976942A (zh) 无穷时域优化的间歇过程2d约束容错控制方法
CN108107723A (zh) 非线性间歇过程的2d最优模糊控制器设计方法
CN110503251A (zh) 一种基于Stacking算法的非节假日负荷预测方法
Liu et al. Prescribed-time containment control with prescribed performance for uncertain nonlinear multi-agent systems
CN111427266A (zh) 一种针对带扰动的非线性系统辨识方法
DE102015203210A1 (de) Verfahren zum Regeln einer Regelstrecke, Vorrichtung zur Erzeugung von Reglerparametern und Steuergerät
CN110825051B (zh) 一种基于gap metric的不确定系统的多模型控制方法
CN106094524A (zh) 基于输入趋势补偿的快速模型预测控制方法
Lin et al. PID control strategy for UAV flight control system based on improved genetic algorithm optimization
CN110658722B (zh) 一种基于gap的自均衡多模型分解方法及系统
CN110888323A (zh) 一种用于切换系统智能优化的控制方法
CN103439964B (zh) 一种ocsvm监控模型的在线更新系统及方法
CN104155876A (zh) 一种pid控制器的分离实现方法
CN114399078A (zh) 一种负荷日前预测与超短期预测的协同预测方法及系统
CN113705929A (zh) 一种基于负荷特征曲线及典型特征值融合的春节假期负荷预测方法
CN110134095A (zh) 火电厂热工模拟量控制系统优化的方法及终端设备
da_Fonseca Neto et al. On the convergence of DLQR control system design and recurrences of riccati and lyapunov in dynamic programming strategies
CN112415886B (zh) 一种基于pid的一体化多模型控制方法
Kiknadze et al. INNOVATIVE MODEL DESIGN FOR THE MANAGEMENT OF REGIONAL SUSTAINABLE DEVELOPMENT
CN103064286A (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
GR01 Patent grant
GR01 Patent grant