一种多变量半自适应预测控制方法与流程

本发明涉及优化控制领域，具体地说是一种多变量半自适应预测控制方法。

背景技术：

模型预测控制(modelpredictivecontrol，mpc)是一种能够处理多变量、约束、纯滞后等问题的控制算法，目前在工业的各个领域中应用广泛。其通过预测模型预测未来时刻的输出，并通过二次性能指标得到最优的控制作用。预测模型的精确度是预测控制实施成功与否的关键因素之一，不准确的预测模型会导致大的预测偏差，进而影响控制效果。另外，目前工业应用中广泛采用双层结构的mpc，其稳态目标计算(steady-statetargetcalculation,sstc)的稳态模型来自动态优化(dynamicoptimization,do)层动态模型的稳态增益，稳态增益直接影响最优的经济工作点。因此，模型的精确度不仅仅影响控制效果，同样影响生产效益。

在mpc的实施过程中，主要包括如下阶段：

1)效益分析与初步设计；

2)预测试；

3)动态响应测试与模型辨识；

4)控制器离线仿真与参数整定；

5)控制器投运与培训；

6)控制器维护。

通常预测试采用开环阶跃测试辨识过程模型，过程测试和模型辨识占整个工程实施50％以上的时间.由于生产装置的生命周期特征，设备磨损、结垢、产品特性以及操作条件改变等因素，导致过程动态特性具有慢时变性。而模型与装置之间的误差，是影响预测控制器性能的主要因素之一。为了确保mpc效果和经济效益，有必要在模型预测出现较大偏差时，通过再辨识更新预测模型。再辨识过程中，断开控制器，产品质量不能达到要求，同时危及生产安全。

技术实现要素：

针对现有技术的不足，本发明提供一种多变量半自适应预测控制方法，实现在保证生产安全的同时，解决由于模型失配导致的控制器性能下降问题。

本发明为实现上述目的所采用的技术方案是：

一种多变量半自适应预测控制方法，包括以下步骤：

步骤1：初始化多变量过程的测试信号，建立辨识所需要的种子模型和控制器性能评判标准；

步骤2：判断当前控制器是否满足控制器性能评判标准，如果满足，则继续执行步骤2，实时监测控制器性能；否则执行步骤3；

步骤3：判断当前的预测控制器结构，如果当前的预测控制器结构为单层结构或，当前的预测控制器结构为双层结构且满足经济效益放松条件，则执行步骤4；

步骤4：修改设定点预测控制为区间预测控制，以多变量测试信号作为前馈输入，建立具有测试和控制功能的区间预测控制器，并进行测试，采集测试数据；

步骤5：对测试数据进行辨识，得到辨识后的模型，用该模型更新原有预测控制器模型，并返回步骤2。

所述种子模型为运行的预测控制器模型或由生产过程中被控变量和操作变量的历史数据分析获取。

所述控制器性能评判标准采用基于低扰动实验的mpc模型诊断方法建立。

所述多变量过程的测试信号为正交四值信号。

所述经济效益放松条件为：

步骤1：计算最优效益目标值jopt＝cuss(n)，其中，c＝[c1c2…cm]分别表示代价系数；uss(k)表示当前时刻操作变量的稳态值；下角标ss(steady-state)表示稳态；

步骤2：计算当前效益值j(k)＝cu(k)，其中，u(k)表示当前时刻操作变量值；

步骤3：计算当前效益差δj(k)＝j(k)-rejopt，其中，re为效益放松系数，通过初始化设置；

步骤4：如果δj(k)≥0，当前的系统状态满足经济效益放松条件，否则不满足，返回步骤1。

所述区间预测控制包括以下优化形式：

s.t.

ull≤um(k)+δum(k)≤uhl

δull≤δum(k)≤δuhl

其中，y∈rⁿ为被控变量；u∈r^m为操作变量；v∈r^q为可测扰动变量；m和p分别表示控制时域和预测时域；和分别表示第i个被控变量初始预测值和预测值；δum(k)＝[δu1,m(k)…δum,m(k)]^t，δuj,m(k)＝[δuj(k),…,δuj(k)]^t；δv(k)＝[δv1(k)…δvq(k)]^t，因为扰动无法预测，因此通常假设vl(k)＝,…,＝vl(k+p)，δvl(k)＝vl(k)-vl(k-1)；ysp(k)表示设定值；ll和hl分别表示上下限，yll表示被控变量的下限，q和r分别表示被控变量和控制权矩阵；g和h分别为输入输出动态矩阵和可测扰动输出动态矩阵；“*”表示hadamard积；目标函数中第一项和第二项分别表示对被控变量开环预测值低于区间下限和高于区间上限的惩罚，ql和qh为对应的权矩阵；l¹和l²分别违反区间约束的预测值的系数向量，即，当被控变量的预测值满足区间约束时，对应的系数(和)为0；否则，设置为1。当ysp,hl＝ysp,ll，区间预测控制等价于设定值预测控制。

所述具有测试和控制功能的区间预测控制器为：

s.t.

δum(k)＝δucon,m(k)+λδuident,m(k)

ull≤um(k)+δum(k)≤uhl

δull≤δum(k)≤δuhl

其中，δucon,m(k)表示未来m步预测输出增量值；δuident,m(k)表示未来m步激励信号增量值；λ为幅值强度；和分别表示第i个被控变量初始预测值和预测值；δum(k)＝[δu1,m(k)…δum,m(k)]^t，ysp(k)表示设定值；ll表示下限；；hl表示上限；yll表示被控变量的下限；q和r分别表示被控变量和控制权矩阵；g为输入输出动态矩阵；“*”表示hadamard积；ql和qh为对被控变量开环预测值低于区间下限和高于区间上限惩罚的权矩阵；l¹和l²分别违反区间约束的预测值的系数向量。

当幅值强度λ＝0时，控制器处于控制模式，当设定值下限ysp,ll与设定值上限ysp,hl相等时，控制目标为对设定点的实时跟踪；

当幅值强度λ≠0时，控制器处于综合测试模式，δuident,m(k)作为激励信号，可视为系统的可测扰动，实现辨识需要的持续激励。δucon,m(k)保证系统在控制区间的约束满足。

所述控制模式为通过区间预测控制或设定值预测控制保障生产过程中产品质量和生产的安全。

所述综合测试模式为在模型测试过程中，在实现辨识过程充分激励的同时，通过区间预测控制或设定值预测控制保障生产过程中产品质量和生产的安全。

本发明具有以下有益效果及优点：

1.本发明解决辨识所需的充分激励信号与系统平稳运行的矛盾；

2.本发明解决辨识输入信号与不可测噪声的相关性问题；

3.本发明不需要增加新的约束条件，且不转变优化形式，无需额外计算量；

4.本发明减少测试时间；

5.本发明保证了测试过程中，生产过程的安全持续运行；

6.本发明在线解决因模型因素造成的控制器性能下降的问题。

附图说明

图1是本发明的方法流程图；

图2是本发明的初始化流程图；

图3是本发明的性能检测过程流程图；

图4是本发明的测试过程流程图；

图5是本发明的辨识过程流程图；

图6是本发明的系统架构图；

图7是本发明的面向约束mpc的测试模块结构图；

图8是本发明的面向双层mpc的测试模块结构图。

具体实施方式

下面结合附图及实施例对本发明做进一步的详细说明。

如图1为本发明的方法流程图。

包括如下步骤：

步骤s1初始化；

步骤s2为控制性能监测过程；

步骤s3为测试过程；

步骤s4，辨识过程。

如图2所示是本发明的初始化流程图。

步骤s1进一步包括：

步骤s11，建立辨识所需要的种子模型。其可为投运的预测控制器模型；也可由生产过程中被控变量和操作变量的历史数据分析获取。通常其具有合理的延迟时间和稳态时间。

步骤s12，建立多变量过程的测试信号。针对多输入输出系统，为了同时激励所有输入节省辨识时间。本发明采用一种改进的正交四值信号，其由二个周期信号乘积构成。

ui(k)＝hi(k)p(k)

其中，p(k)表示周期为np，幅值为[a,-a×p]的prbs信号，信号hi(k)的周期为nh，选自阶次为nh＝2^m-1的hadamard矩阵i行元素。因此，单个输入信号的周期长度为nh×np。

步骤s13，建立控制器性能评判标准

步骤s14，设置效益放松系数re。

如图3所示是本发明的性能检测过程流程图。

步骤s2为控制性能监测过程。

具体的，步骤s21根据步骤s13建立的控制器性能评判标准实时判断当前控制器是否满足评判标准，如满足，则继续进行控制器的性能判断；否则，进入下一步骤(步骤s3，测试)。

如图4所示是本发明的测试过程流程图。

步骤s3为测试过程，进一步包括：

步骤s31判断当前的mpc控制器结构是否为双层结构。如果是，进入步骤s32，否则进入步骤s36。

步骤s32计算最优效益目标值jopt＝cuss(n)，其中，c＝[c1c2…cm]分别表示代价系数；uss(k)表示当前时刻操作变量的稳态值；下角标ss表示稳态。

步骤s33计算当前效益值j(k)＝cu(k)，其中，u(k)表示当前时刻操作变量值。

步骤s34计算当前效益差δj(k)＝j(k)-rejopt。

步骤s35判读当前时刻的效益差δj(k)是否大于0。如果δj(k)≥0，说明当前的系统状态满足步骤s14设置的经济效益放松条件，进入步骤s36，进行测试；否则进入s32。

步骤s36改变当前的mpc设定点控制模式为区间控制模式。

设定点控制模式的mpc目标函数和约束条件如下：

s.t.

ull≤um(k)+δum(k)≤uhl

δull≤δum(k)≤δuhl

其中，y∈rⁿ为被控变量；u∈r^m为操作变量；v∈r^q为可测扰动变量；m和p分别表示控制时域和预测时域；和分别表示第i个被控变量初始预测值和预测值；δum(k)＝[δu1,m(k)…δum,m(k)]^t，δuj,m(k)＝[δuj(k),…,δuj(k)]^t；δv(k)＝[δv1(k)…δvq(k)]^t，因为扰动无法预测，因此通常假设vl(k)＝,…,＝vl(k+p)，δvl(k)＝vl(k)-vl(k-1)；ysp(k)表示设定值；ll和hl分别表示上下限，yll表示被控变量的下限，q和r分别表示被控变量和控制权矩阵；g和h分别为输入输出动态矩阵和可测扰动输出动态矩阵。

区间mpc的优化形式如下：

s.t.

ull≤um(k)+δum(k)≤uhl

δull≤δum(k)≤δuhl

其中，“*”表示hadamard积。目标函数中第一项和第二项分别表示对被控变量开环预测值低于区间下限和高于区间上限的惩罚，ql和qh为对应的权矩阵；l¹和l²分别违反区间约束的预测值的系数向量，即，当被控变量的预测值满足区间约束时，对应的系数({和)为0；否则，设置为1。当ysp,hl＝ysp,ll，区间mpc等价于设定值mpc。

步骤s37将s12建立的测试信号，δuident(k)作为区间控制器的前馈输入。则式(2)的优化命题变为：

s.t.

ull≤um(k)+δum(k)≤uhl

δull≤δum(k)≤δuhl

其中，δucon,m(k)表示未来m步预测输出增量值；δuident,m(k)表示未来m步激励信号增量值；λ为幅值强度。

当幅值强度λ＝0时，控制器处于控制模式，当设定值下限ysp,ll与设定值上限ysp,hl相等时，控制目标为对设定点的实时跟踪。

当幅值强度λ≠0时，控制器处于综合测试模式，δuident,m(k)作为激励信号，可视为系统的可测扰动，实现辨识需要的持续激励。δucon,m(k)保证系统在控制区间的约束满足。

注：随着幅值强度λ的增加，相应的操纵变量变化幅度及被控变量cv变化幅度也随之增加，测试效果也将越好，但同时被控变量违反区间约束的可能性也将越大。工程师可在目标跟踪控制与测试数据质量之间作出优化权衡。

步骤s38采集过程数据，待后期辨识所用。

如图5所示是本发明的辨识过程流程图。

步骤s4为辨识过程，进一步包括：

步骤s41通过步骤s38得到的数据(可以为实时或离线数据)进行辨识，得到新的模型。

步骤s42将步骤s41得到的新的模型作为mpc的控制器模型。更新mpc控制器模型后，此时，如果性能检测模块检测控制器性能满足要求，则停止测试和辨识；否则，继续测试和辨识，直至更新模型后的控制器性能满足要求。

图6为本发明涉及的多变量半自适应预测控制系统架构简图。

主要包括性能监测模块、辨识模块和测试模块。性能监测模块实时检测当前控制器是否满足控制性能指标。当性能监测模块检测到控制器性能指标不满足要求时，测试模块中测试信号实现对过程的充分激励，获得辨识模块所需的数据。辨识模块根据mpc的结构分为面向约束mpc的测试模块和面向双层结构mpc的测试模块。辨识模块通过测试模块获得的数据，进行模型辨识，并将辨识后的模型作为新的控制器模块。更新后的控制器如果满足控制器性能指标，结束测试和辨识过程，否则，继续测试和辨识过程，直至控制器性能指标得到满足。

图7为本发明中面向约束mpc的测试模块结构图。

主要包括区间预测控制器、多变量测试信号。其中多变量测试信号实现测试过程的充分激励，区间预测控制器实现在测试过程中，保证过程的被控变量满足区间约束。

图8为本发明中面向双层mpc的测试模块结构图。

主要包括效益放松条件判断、区间预测控制器、多变量测试信号。其中，效益放松条件判断表示对当期状态下的经济效益是否满足设定要求。在该结构下，如果当期状态下的经济效益满足设定要求，则进行面向约束mpc的测试，即在充分激励测试的同时，保证被控变量满足区间约束。否则，继续进行判断，直至满足设定要求。

以上公开的仅为本申请的具体实施例，但并不能因此理解为对本发明专利范围的限制。应当指出，在不脱离本发明构思的前提下，任何本领域的技术人员能思之的变化，都应落在本申请的保护范围内。