CN107797451B - 未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法 - Google Patents

未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法 Download PDF

Info

Publication number
CN107797451B
CN107797451B CN201710987844.3A CN201710987844A CN107797451B CN 107797451 B CN107797451 B CN 107797451B CN 201710987844 A CN201710987844 A CN 201710987844A CN 107797451 B CN107797451 B CN 107797451B
Authority
CN
China
Prior art keywords
time
discrete
lag
state
tank
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
CN201710987844.3A
Other languages
English (en)
Other versions
CN107797451A (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.)
Zhejiang Juxin Automotive Electronics Co ltd
Original Assignee
Liaoning Shihua 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 Liaoning Shihua University filed Critical Liaoning Shihua University
Priority to CN201710987844.3A priority Critical patent/CN107797451B/zh
Publication of CN107797451A publication Critical patent/CN107797451A/zh
Application granted granted Critical
Publication of CN107797451B publication Critical patent/CN107797451B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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)
  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)

Abstract

本发明涉及未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法，包括以下步骤：步骤一：建立具有不确定时滞系统的状态空间模型，步骤二：将构建模型转化为扩展时滞不确定状态空间模型，步骤三：设计基于上述扩展模型的控制器，步骤四：采用线性矩阵不等式的形式求解控制器增益。本发明的优点为：采用扩展时滞不确定状态空间模型来设计控制器，在时滞存在的条件下，可以保证控制器单独调节系统的状态和输出误差，提供了更多调节的自由度；给定新的、具有更小保守性和较为简单的上述扩展模型基于LMI形式的时滞依赖稳定条件来求解系统的控制器增益，降低了控制器的保守性；引入H‑infinity性能指标到控制器设计中，克服任意有界干扰，降低了控制成本。

Description

未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法
技术领域
本发明属于工业过程的先进控制领域，具体涉及一种未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法。
背景技术
随着现代工业的快速发展，工业产品社会需求量日益增加，自动控制系统规模逐渐扩大、复杂性日益增加，但对产品质量和系统性能的要求日益严格，这势必要求生产设备运行在“卡边”条件下。为此，如何既能保障系统可靠性和安全性又能提高系统的性能是一个开放性的问题。在实际工业过程中，大多数过程尤其是石油化工过程具有不确定性和时变时滞，并且受输入输出约束和外界干扰的影响。这些因素可能会导致系统性能恶化甚至使系统不稳定。为此，研究一种针对工业过程具有不确定性、未知干扰、状态时变时滞和输入输出约束的先进控制方法是非常必要的。以往的技术或方法采用将时滞扩展到模型中，把有时滞的问题转化为无时滞的问题，但这种方法较大的扩大了模型维数，从而增加了计算量。另外也有直接采用具有时滞的模型，但需要引进一些松弛变量，增加了一定的保守性。
目前针对三容水箱系统液位的工业过程具有不确定性、未知干扰、状态时变时滞和输入输出约束等特性，现有的技术和方法所的设计的控制器具有较高的计算量和一定的保守性，这势必会增加三容水箱控制系统的计算负担和降低系统的性能。因此，为了实现提高产品质量、增加产品收率、降低能源消耗和提升经济效益等目标要求，迫切需要提出一种新的先进控制方法来保证系统高效平稳运行。
发明内容
本发明正是针对具有不确定性、未知干扰、状态时变时滞和输入输出约束的三容水箱，提出了具有未知干扰的不确定性系统的时滞依赖鲁棒约束预测控制方法，该方法可以降低三容水箱控制系统的计算量、优化系统的控制性能、降低控制器的保守性。
本发明是通过以下技术方案实现的：
该方法首先将具有不确定性、未知干扰、状态时变时滞和输入输出约束的离散系统表示为三容水箱的状态空间形式，并将输出误差扩展到该模型中，从而形成新的扩展状态空间模型，并给出了一种新的、具有更小保守性和较为简单的上述扩展模型基于LMI形式的时滞依赖稳定条件，从而求解三容水箱系统的控制律，同时为了克服任意未知干扰，将H-infinity性能指标引入到控制器设计中。
未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法，其特征在于，包括以下步骤：
步骤一：建立三容水箱的单输入单输出二阶模型，：
Figure GDA0002457057860000021
式中，h1,h3分别为T1和T3的液位高度；Qin为注水管1的流量，作为操纵变量；Q13为罐1向罐3流入液体的流量；Qout为水箱3至底部水槽的流量，其中：
Figure GDA0002457057860000022
Figure GDA0002457057860000023
其中，az1＝0.48，az2＝0.58；S1＝Sn＝5×10-5m2，为连接管道的截面积；S＝0.154m2，为罐的截面积；Hmax＝0.6m是液位高度上限；sgn(·)为符号函数，h1,h3的初值为0，则状态变量和状态输入分别为：
Figure GDA0002457057860000031
u(k)＝Qin(k) (5)
步骤二：将操作点0.33Hmax的局部线性化，得到三容水箱系统具有不确定状态时的变时滞离散状态空间模型：
Figure GDA0002457057860000032
式中，1≤d(k)≤3,
Figure GDA0002457057860000033
Figure GDA0002457057860000034
C＝[1 0]，
Figure GDA0002457057860000035
w(k)＝(0.0005Δ3 0.0005Δ4)T，其中Δ1,Δ2,Δ3,Δ4是[-1 1]之间的随机数；
步骤三：设计基于上述三容水箱系统扩展模型的控制器：
Figure GDA0002457057860000036
式中，
Figure GDA0002457057860000037
是本发明方法所求的控制器增益，则方程(7)可以转化为：
Figure GDA0002457057860000038
式中，
Figure GDA0002457057860000039
该发明方法将系统(1)转化为扩展的模型(8)，利用该模型来设计三容水箱系统的控制器，基于上述扩展模型(8)，系统优化可以描述为如下min-max优化问题：
Figure GDA00024570578600000310
Figure GDA00024570578600000311
式中，
Figure GDA00024570578600000312
和
Figure GDA00024570578600000313
分别是系统跟踪误差和控制输入相应维数的加权矩阵；Δu(k+i|k)和
Figure GDA0002457057860000041
分别是离散k+i时刻的增量控制输入和状态变量；ΔuM和ΔyM分别是指系统增量的输入和输出的边界值；
步骤四：采用线性矩阵不等式的形式求解控制器增益
Figure GDA0002457057860000042
给定时滞依赖的稳定判据条件，根据线性矩阵不等式约束，应用MATLAB软件LMI工具箱，求解控制器增益
Figure GDA0002457057860000043
该时滞依赖的稳定判据条件如下：
Figure GDA0002457057860000044
Figure GDA0002457057860000045
Figure GDA0002457057860000046
Figure GDA0002457057860000047
其中，
Figure GDA0002457057860000048
均为正定对称矩阵，矩阵
Figure GDA0002457057860000049
以及标量
Figure GDA00024570578600000410
γ＞0,θ＞0,0≤dm≤dM，
Figure GDA00024570578600000411
Figure GDA00024570578600000412
Figure GDA00024570578600000413
Figure GDA00024570578600000414
Figure GDA0002457057860000051
所述步骤一中的三容水箱具体型号为TTS20，其中三个透明水箱分别称之为T1、T2和T3，且T1与T2之间、T2与T3之间均通过阀门依次连接，T3后侧安装有排水的阀门，且上述三个阀门的横截面积相同。
三个所述水箱外壁设有泄水阀。
三个所述水箱的内腔均安装有一压力液位传感器，且液位传感器用来作为系统的测量元件。
基于LMI形式的时滞依赖稳定条件(9)-(12)，可以通过MATLAB软件求解三容水箱系统的控制增益，保证三容水箱系统的鲁棒性是渐进稳定的，且鲁棒H∞性能指标小于γ。不同以往的技术和方法，本发明方法采用一种差分的方法来构建李亚普诺夫函数，该函数没有引入一些松弛变量，而是充分利用时滞的上下界的信息，在推导时滞依赖稳定条件时避免了采用差分不等式对于交叉项的边界和模型转化，从而降低了控制器的保守性。同时，为了克服任意未知干扰，本发明方法将H-infinity性能指标γ引入到控制器设计中，γ值越小，表明系统的抗干扰能力越强。
与现有技术相比，本发明的有益效果为：本发明针对三容水箱系统具有不确定性、未知干扰、状态时变时滞和输入输出约束等特性，发明了一种时滞依赖的鲁棒约束模型预测控制方法，使被控过程能更加稳定和高效的运行，提高了系统的性能。建立了具有时滞的扩展状态空间模型，在时滞存在的条件下，保证增量状态和输出误差同时收敛，并为控制器的调节提供了更多的自由度，改善了系统性能；采用一种新的、具有更小保守性和较为简单的上述扩展模型基于LMI形式的时滞依赖稳定条件来求解系统的控制器增益，降低了控制器的保守性；引入H-infinity性能指标到控制器设计中，可以克服任意有界干扰，减少了外界扰动对系统影响，降低了控制成本。
附图说明
图1为提出的方法与传统的方法的(a)输出响应和(b)控制输入。
图2为本发明步骤流程图。
图3为水箱TTS20实物图。
图4为TTS20结构和过程流程图。
图5为定义TTS20变量和参数的原理图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步的说明。
本发明提供的未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法，具体步骤如下：
步骤一：建立具有不确定时滞系统的状态空间模型：
Figure GDA0002457057860000061
式中，
Figure GDA0002457057860000062
w(k)是表示在离散k时刻的系统状态、输入、输出和未知外界干扰，d(k)是依赖于离散k时刻的时变时滞，满足：
dm≤d(k)≤dM (16)
式中，dM和dm分别是时滞的上界和下界，
A(k)＝A+Δa(k),Ad(k)＝Ad+Δd(k),B(k)＝B+Δb(k),A,Ad,B和C是相应维数的常数矩阵，并且Δa(k),Δd(k)和Δb(k)是在离散k时刻的不确定摄动，可以表示为：
[Δa(k) Δd(k) Δb(k)]＝NΔ(k)[H Hd Hb] (17)
且ΔT(k)Δ(k)≤I；
式中，N,H,Hd和Hb是相应维数的已知常数矩阵，Δ(k)是依赖于离散时间k的不确定摄动；
步骤二：将构建的时滞状态空间模型转化为扩展时滞不确定状态空间模型，在方程(1)左右两边分别乘以后移算子Δ，为：
Figure GDA0002457057860000071
式中，Δ＝1-q-1，
Figure GDA0002457057860000072
Figure GDA0002457057860000073
定义设定值为c(k)，则跟踪误差可以表示为：
e(k)＝y(k)-c(k) (19)
综合方程(4)和(5)，得
Figure GDA0002457057860000074
通过结合跟踪误差和增量的状态变量，新扩展的具有不确定性、未知干扰、状态时变时滞的状态空间模型可以表示为：
Figure GDA0002457057860000075
式中，
Figure GDA0002457057860000076
步骤三：设计基于上述扩展模型的控制器：
Figure GDA0002457057860000077
式中，
Figure GDA0002457057860000078
是所发明方法的控制器增益，可以通过如下步骤四计算，则方程(7)可以转化为：
Figure GDA0002457057860000081
式中，
Figure GDA0002457057860000082
该发明方法将系统(1)转化为扩展的模型(9)，利用该模型来设计系统的控制器，此模型的优点是当时滞存在的条件下，可以保证控制器单独调节系统的状态和输出误差，从而使得其分别收敛，并为控制器提供了更多调节的自由度；
为此，基于上述扩展模型(8)，系统优化可以描述为如下min-max优化问题：
Figure GDA0002457057860000083
Figure GDA0002457057860000084
式中，
Figure GDA0002457057860000085
和
Figure GDA0002457057860000086
分别是系统跟踪误差和控制输入相应维数的加权矩阵；Δu(k+i|k)和
Figure GDA0002457057860000087
分别是离散k+i时刻的增量控制输入和状态变量；ΔuM和ΔyM分别是指系统增量的输入和输出的边界值；
步骤四：计算控制器增益
Figure GDA0002457057860000088
给定时滞依赖的稳定判据条件，根据线性矩阵不等式约束，应用MATLAB软件LMI工具箱，求解控制器增益
Figure GDA0002457057860000089
该时滞依赖的稳定判据条件如下：
Figure GDA00024570578600000810
Figure GDA0002457057860000091
Figure GDA0002457057860000092
Figure GDA0002457057860000093
其中，
Figure GDA0002457057860000094
均为正定对称矩阵，矩阵
Figure GDA0002457057860000095
以及标量
Figure GDA0002457057860000096
γ＞0,θ＞0,0≤dm≤dM，
Figure GDA0002457057860000097
Figure GDA0002457057860000098
Figure GDA0002457057860000099
Figure GDA00024570578600000910
Figure GDA00024570578600000911
实施例
如图1和2所示，本发明为针对TTS20型三容水箱系统的液位进行仿真对比的研究，从图中可以看出，利用扩展的具有时滞的状态空间模型，提出的发明方法能够更好的跟踪阶跃设定值，具有更好的跟踪性能，而采用传统的方法输出具有一定超调。同时通过控制输入曲线可以看出，提出的方法输入更加平滑，可以更好的抵抗未知干扰和不确定性，体现良好的抗干扰能力。因此，提出的发明方法通过高性能的控制可以更好的满足控制需求，从而保证系统在“卡边”控制的稳定性，提高系统的可靠性和安全性。
作为工业生产过程中诸多被控对象的抽象模型，三容水箱因其在非线性、时滞性等方面具有典型的代表性，因此，实际工业系统中的诸多被控对象，其整体或者局部都可以看作是一种三容水箱模型，本发明方法以三容水箱作为仿真对象对于实际的工业系统而言具有较为典型的代表意义。
如图3所示，为实验室常见的TTS20三容水箱，TTS20型水箱的结构和整体流程如图4所示。其中，中间的三个透明水箱分别称之为Tank1(T1)，Tank2(T2)，Tank3(T3)，其中，T1、T2之间、T2、T3之间均通过其底部的阀门依次连接，在T3后装有一个专门用来排水的阀门，可以将T3中的水排放到贮水槽中，且上述三个阀门横截面相同。此外，在三个水箱又各自拥有一个独立的泄水阀，可以将水箱的水直接泄露到下面的贮水槽之中，通过泄水阀及泄露到贮水槽的水流量来对水箱的故障信息进行描述。
左侧上下两个泵分别称之为Pump1(P1)、Pump2(P2)，他们可以分别将贮水槽中的水抽出并分别通过一个线性比例电磁阀将水注入到T1、T3中。而上述从底部阀门及T3排水阀门排除到贮水槽中的水供P1、P2使用，因此便构成了一个回路。其中，T1、T2、T3中均安装有一个测量压力液位传感器作为系统的测量元件，系统中的每一个阀门都可以通过给定的信号完成阀门的开关动作。
通过开启和关闭连接阀门及泄露阀，水箱可以很方便的改造为单输入单输出，多输入多输出，三阶，二阶，一阶模型。定义TTS20水箱变量和参数的原理图如图5所示。这针对水箱的单入单出二阶模型进行仿真，水箱的改造方法为：仅打开T1与T3之间连接阀门以及T3与水箱底部水槽之间的阀门，把泵1的流量作为控制输入，把T1中液体的液位高度作为系统输出。
未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法，其特征在于，包括以下步骤：
步骤一：建立三容水箱的单输入单输出二阶模型：
Figure GDA0002457057860000111
式中，h1,h3分别为T1和T3的液位高度；Qin为注水管1的流量，作为操纵变量；Q13为罐1向罐3流入液体的流量；Qout为水箱3至底部水槽的流量，其中：
Figure GDA0002457057860000112
Figure GDA0002457057860000113
其中，az1＝0.48，az2＝0.58；S1＝Sn＝5×10-5m2，为连接管道的截面积；S＝0.154m2，为罐的截面积；Hmax＝0.6m是液位高度上限；sgn(·)为符号函数，h1,h3的初值为0，则状态变量和状态输入分别为：
Figure GDA0002457057860000114
u(k)＝Qin(k) (5)
步骤二：将操作点0.33Hmax的局部线性化，得到三容水箱系统具有不确定状态时的变时滞离散状态空间模型：
Figure GDA0002457057860000115
式中，1≤d(k)≤3,
Figure GDA0002457057860000116
Figure GDA0002457057860000117
C＝[1 0]，
Figure GDA0002457057860000118
w(k)＝(0.0005Δ3 0.0005Δ4)T，其中Δ1,Δ2,Δ3,Δ4是[-1 1]之间的随机数；
步骤三：设计基于上述三容水箱系统扩展模型的控制器：
Figure GDA0002457057860000119
式中，
Figure GDA00024570578600001110
是本发明方法所求的控制器增益，则方程(7)可以转化为：
Figure GDA0002457057860000121
式中，
Figure GDA0002457057860000122
该发明方法将系统(1)转化为扩展的模型(8)，利用该模型来设计三容水箱系统的控制器，基于上述扩展模型(8)，系统优化可以描述为如下min-max优化问题：
Figure GDA0002457057860000123
Figure GDA0002457057860000124
式中，
Figure GDA0002457057860000125
和
Figure GDA0002457057860000126
分别是系统跟踪误差和控制输入相应维数的加权矩阵；Δu(k+i|k)和
Figure GDA0002457057860000127
分别是离散k+i时刻的增量控制输入和状态变量；ΔuM和ΔyM分别是指系统增量的输入和输出的边界值；
步骤四：采用线性矩阵不等式的形式求解控制器增益
Figure GDA0002457057860000128
给定时滞依赖的稳定判据条件，根据线性矩阵不等式约束，应用MATLAB软件LMI工具箱，求解控制器增益
Figure GDA0002457057860000129
该时滞依赖的稳定判据条件如下：
Figure GDA00024570578600001210
Figure GDA00024570578600001211
Figure GDA00024570578600001212
Figure GDA0002457057860000131
其中，
Figure GDA0002457057860000132
均为正定对称矩阵，矩阵
Figure GDA0002457057860000133
以及标量
Figure GDA0002457057860000134
γ＞0,θ＞0,0≤dm≤dM，
Figure GDA0002457057860000135
Figure GDA0002457057860000136
Figure GDA0002457057860000137
Figure GDA0002457057860000138
Figure GDA0002457057860000139
采用传统的鲁棒约束预测控制和本文所提的发明方法进行对比，控制器的模型选为(14)，控制参数统一选为：Q＝1，R＝0.1，
Figure GDA00024570578600001310
Figure GDA00024570578600001311
输入和输出约束为：
Figure GDA00024570578600001312
设定值取为：
Figure GDA00024570578600001313
综上，本发明以TTS20型水箱的液位控制设计为例，来验证本发明所提出的控制方法的有效性和可行性。仿真结果表明系统在具有不确定性、未知干扰、状态时变时滞和输入输出约束情况下，可以更好的跟踪液位设定值和抵抗未知随机干扰，具有较好的跟踪性能和抗干扰的能力，可以使闭环系统在最优和稳定的条件下运行且具有良好的控制性能。因此，这种发明方法的提出，从长远来看，可以保证系统高效、安全和平稳运行，从而可以提高产品质量、增加产品收率、降低能源消耗和提升经济效益等。

Claims (4)

1.未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法，其特征在于，包括以下步骤：
步骤一：建立三容水箱的单输入单输出二阶模型：
Figure FDA0002752331950000011
式中，h1,h3分别为水槽T1和水槽T3的液位高度；Qin为注水管1的流量，作为操纵变量；Q13为罐1向罐3流入液体的流量；Qout为水箱3至底部水槽的流量，其中：
Figure FDA0002752331950000012
Figure FDA0002752331950000013
其中，g为当地重力加速度，h2为水槽T2的液位高度，az1＝0.48，az2＝0.58；S1＝Sn＝5×10-5m2，为连接管道的截面积；S＝0.154m2，为罐的截面积；Hmax＝0.6m是液位高度上限；sgn(·)为符号函数，h1,h3的初值为0，则状态变量和状态输入分别为：
Figure FDA0002752331950000014
u(k)＝Qin(k) (5)
其中，x1(k)为离散k时刻水槽T1的液位状态变量，x2(k)为离散k时刻水槽T3的液位状态变量，h1(k)为离散k时刻水槽T1的液位高度，h3(k)为离散k时刻水槽T3的液位高度，Qin(k)为离散k时刻注水管1的流量；
步骤二：将操作点0.33Hmax的局部线性化，得到三容水箱系统具有不确定状态时的变时滞离散状态空间模型：
Figure FDA0002752331950000021
式中，x(k),u(k),y(k),w(k)是表示在离散k时刻的系统状态、输入、输出和未知外界干扰，x(k+1)表示在离散k+1时刻的系统状态，A(k)＝A+Δa(k)为离散k时刻系统的状态矩阵，Ad(k)＝Ad+Δd(k)为离散k时刻系统的状态时滞矩阵，B(k)＝B+Δb(k)为离散k时刻系统的输入矩阵，Δa(k),Δd(k)和Δb(k)是在离散k时刻的不确定摄动，可以表示为[Δa(k) Δd(k) Δb(k)]＝NΔ(k)[H Hd Hb]，d(k)是依赖于离散k时刻的时变时滞，满足
Figure FDA0002752331950000022
Figure FDA0002752331950000023
Figure FDA0002752331950000024
C＝[1 0]，
Figure FDA0002752331950000025
w(k)＝(0.0005Δ3 0.0005Δ4)T，其中Δ1,Δ2,Δ3,Δ4是[-1 1]之间的随机数；
步骤三：设计基于上述三容水箱系统扩展模型的控制器：
Figure FDA0002752331950000026
式中，Δu(k)为离散k时刻系统的输入增量，Δx(k)为离散k时刻系统的状态增量，e(k)为离散k时刻系统的输出跟踪误差，
Figure FDA0002752331950000027
为离散k时刻系统的扩展状态变量，
Figure FDA0002752331950000028
是本发明方法所求的控制器增益，则方程(7)可以转化为：
Figure FDA0002752331950000029
式中，
Figure FDA00027523319500000210
Δy(k)为离散k时刻系统的输出增量，z(k)为离散k时刻系统的输出跟踪误差，
Figure FDA00027523319500000211
Figure FDA0002752331950000031
Figure FDA0002752331950000032
Figure FDA0002752331950000033
Figure FDA0002752331950000034
为离散k-1-d(k-1)时刻系统的状态增量，1≤d(k-1)≤3,
Figure FDA0002752331950000035
Δ(k-1)为-1到1之间的随机数，
Figure FDA0002752331950000036
为输出矩阵C和不确定性矩阵N组合成的不确定性扩展矩阵， x(k-1)为离散k-1时刻系统的状态，u(k-1)为离散k-1时刻系统的控制输入，Δw(k)为离散k时刻系统的干扰增量；
该发明方法将系统(1)转化为扩展的模型(8)，利用该模型来设计三容水箱系统的控制器，基于上述扩展模型(8)，系统优化可以描述为如下min-max优化问题：
Figure FDA0002752331950000037
Figure FDA0002752331950000038
式中，
Figure FDA0002752331950000039
为最优性能指标，
Figure FDA00027523319500000310
为离散k时刻最优性能指标；i指在离散k时刻基础上累加的时间，取为0到∞；
Figure FDA00027523319500000311
和
Figure FDA00027523319500000312
分别是系统跟踪误差和控制输入相应维数的加权矩阵；Δu(k+i|k)和
Figure FDA00027523319500000313
分别是离散k+i时刻的增量控制输入和状态变量；Δy(k+i)是离散k+i时刻的控制输出；ΔuM和ΔyM分别是指系统增量的输入和输出的边界值；
步骤四：采用线性矩阵不等式的形式求解控制器增益
Figure FDA00027523319500000314
给定时滞依赖的稳定判据条件，根据线性矩阵不等式约束，应用MATLAB软件LMI工具箱，求解控制器增益
Figure FDA00027523319500000315
该时滞依赖的稳定判据条件如下：
Figure FDA0002752331950000041
Figure FDA0002752331950000042
Figure FDA0002752331950000043
Figure FDA0002752331950000044
其中，
Figure FDA00027523319500000417
均为正定对称矩阵，矩阵
Figure FDA0002752331950000047
以及标量
Figure FDA0002752331950000048
γ＞0,θ＞0,0≤dm≤dM，
Figure FDA0002752331950000049
Figure FDA00027523319500000410
为误差变量矩阵
Figure FDA00027523319500000418
的转置，
Figure FDA00027523319500000411
为系统控制输入加权矩阵
Figure FDA00027523319500000412
的
Figure FDA00027523319500000413
次幂，
Figure FDA00027523319500000414
Figure FDA00027523319500000415
Figure FDA00027523319500000416
Figure FDA0002752331950000051
Figure FDA0002752331950000052
Figure FDA0002752331950000053
r∈(k-3,k),
Figure FDA0002752331950000054
为离散时间在k-3到k之间变化的系统状态变量，
Figure FDA0002752331950000055
Figure FDA0002752331950000056
为系统输入增量上界的平方值，
Figure FDA0002752331950000057
为输出增量上界的平方值。
2.根据权利要求1所述的未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法，其特征在于，所述步骤一中的三容水箱具体型号为TTS20，其中三个透明水箱分别称之为T1、T2和T3，且T1与T2之间、T2与T3之间均通过阀门依次连接，T3后侧安装有排水的阀门，且上述三个阀门的横截面积相同。
3.根据权利要求2所述的未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法，其特征在于，三个所述水箱外壁设有泄水阀。
4.根据权利要求2所述的未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法，其特征在于，三个所述水箱的内腔均安装有一压力液位传感器，且压力液位传感器用来作为系统的测量元件。
CN201710987844.3A 2017-10-21 2017-10-21 未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法 Active CN107797451B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710987844.3A CN107797451B (zh) 2017-10-21 2017-10-21 未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710987844.3A CN107797451B (zh) 2017-10-21 2017-10-21 未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法

Publications (2)

Publication Number Publication Date
CN107797451A CN107797451A (zh) 2018-03-13
CN107797451B true CN107797451B (zh) 2021-01-05

Family

ID=61533760

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710987844.3A Active CN107797451B (zh) 2017-10-21 2017-10-21 未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法

Country Status (1)

Country Link
CN (1) CN107797451B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108803347B (zh) * 2018-08-02 2021-05-11 杭州电子科技大学 一种反向过程鲁棒控制方法
CN109212974B (zh) * 2018-11-12 2021-07-13 辽宁石油化工大学 区间时变时滞系统的鲁棒模糊预测容错控制方法
CN110703607B (zh) * 2019-11-07 2022-05-06 辽宁石油化工大学 具有执行器故障的区间时变时滞系统的随机鲁棒预测容错控制方法
CN111221252B (zh) * 2020-02-25 2023-07-14 西安建筑科技大学 一种针对带分数阶滞后过程工业系统的预测控制器参数解析整定方法
CN111506033B (zh) * 2020-05-08 2023-03-28 辽宁石油化工大学 基于喷嘴压力的注塑机保压容错切换控制方法
CN113223721B (zh) * 2021-03-23 2022-07-12 杭州电子科技大学 一种新型冠状病毒肺炎的预测控制模型

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102944994A (zh) * 2012-12-09 2013-02-27 冶金自动化研究设计院 基于不确定性离散模型的液压活套鲁棒模糊控制方法
CN104407515A (zh) * 2014-11-13 2015-03-11 华南理工大学 一种基于不确定模型的LMIs状态反馈系统控制方法
CN105388758A (zh) * 2015-11-10 2016-03-09 南京航空航天大学 一种基于Delta算子的液位控制系统的自适应滑模控制方法
CN107077106A (zh) * 2014-07-29 2017-08-18 莱恩斯特里姆技术有限公司 自抗扰控制的优化参数化
CN107102553A (zh) * 2017-06-16 2017-08-29 国网重庆市电力公司电力科学研究院 基于时滞鲁棒模型算法的并网逆变器控制方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102944994A (zh) * 2012-12-09 2013-02-27 冶金自动化研究设计院 基于不确定性离散模型的液压活套鲁棒模糊控制方法
CN107077106A (zh) * 2014-07-29 2017-08-18 莱恩斯特里姆技术有限公司 自抗扰控制的优化参数化
CN104407515A (zh) * 2014-11-13 2015-03-11 华南理工大学 一种基于不确定模型的LMIs状态反馈系统控制方法
CN105388758A (zh) * 2015-11-10 2016-03-09 南京航空航天大学 一种基于Delta算子的液位控制系统的自适应滑模控制方法
CN107102553A (zh) * 2017-06-16 2017-08-29 国网重庆市电力公司电力科学研究院 基于时滞鲁棒模型算法的并网逆变器控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《Incremental multivariable predictive functional control and its application in a gas fractionation unit》;SHI Hui-yuan.ect;《Journal of Central South University of Technology》;20151231;第4653-4668页 *
《基于T-S模型的多变量非线性自适应预测函数控制》;杨宇等;《Automation & Instrumentation》;20161231(第7期);第6-11页 *
《增量式多变量预测函数控制及其在气体分馏装置中的应用》;苏成利等;《南京理工大学学报》;20151031(第5期);第625-631页 *

Also Published As

Publication number Publication date
CN107797451A (zh) 2018-03-13

Similar Documents

Publication Publication Date Title
CN107797451B (zh) 未知干扰下不确定性系统时滞依赖鲁棒约束预测控制方法
CN108764539B (zh) 一种梯级电站的上下游水位预测方法
Zhou et al. RBF-ARX model-based MPC strategies with application to a water tank system
Segovia et al. Model predictive control and moving horizon estimation for water level regulation in inland waterways
Andrade et al. A practical NMPC with robustness of stability applied to distributed solar power plants
CN110703607B (zh) 具有执行器故障的区间时变时滞系统的随机鲁棒预测容错控制方法
CN104656448A (zh) 一种基于解耦和扰动观测的超临界机组预测控制方法
Vega et al. Integrated design and control of chemical processes–Part II: An illustrative example
CN105974953A (zh) 一种反应釜负压精馏模糊控制方法
CN114692993B (zh) 融合季节ARIMA和BiLSTM的水利设施形变预测方法
Gouta et al. Experimental sensorless control for a coupled two-tank system using high gain adaptive observer and nonlinear generalized predictive strategy
Chi et al. A model predictive control approach with relevant identification in dynamic PLS framework
Yang et al. Nonlinear programming properties for stable and robust NMPC
CN105676645A (zh) 一种基于函数型权rbf-arx模型的双回路水箱液位预测控制方法
CN105068422A (zh) 一种基于三角区间软约束的模型预测控制方法
Cui et al. Digital twin for power system steady-state modelling, simulation, and analysis
Ben et al. A feedback simulation procedure for real-time control of urban drainage systems
CN102662324A (zh) 槽式反应器基于在线支持向量机的非线性模型预测控制方法
Ye et al. An active approach to space-reduced NCO tracking and output feedback optimizing control for batch processes with parametric uncertainty
Wang et al. Robust model predictive control with bi-level optimization for boiler-turbine system
Gerkšič et al. Feedback action in predictive control: an experimental case study
Grancharova et al. Dual-Mode distributed model predictive control of a quadruple-tank system
CN103345161B (zh) 废塑料裂解炉余热烘干装置压力控制方法
Hafeez et al. Design and testing of decentralized PI (PID) controller for MIMO coupled tank system
Zhang et al. Backstepping-enhanced decentralised PID control for MIMO processes with an experimental study

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
CB03 Change of inventor or designer information
CB03 Change of inventor or designer information

Inventor after: Shi Huiyuan

Inventor after: Li Ping

Inventor after: Su Chengli

Inventor after: Jiang Xueying

Inventor after: Peng Bo

Inventor after: Guo Ying

Inventor after: Cao Jiangtao

Inventor before: Shi Huiyuan

Inventor before: Li Ping

Inventor before: Su Chengli

Inventor before: Guo Ying

Inventor before: Cao Jiangtao

GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240104

Address after: 311100 No.158 Wuchang Avenue, Wuchang Street, Yuhang District, Hangzhou City, Zhejiang Province

Patentee after: Zhejiang Juxin Automotive Electronics Co.,Ltd.

Address before: 113000 No.1, west section of Dandong Road, Wanghua District, Fushun City, Liaoning Province

Patentee before: Liaoming Petrochemical University