一种数据驱动的系统关键性能指标自适应调节方法与流程

本发明涉及一种系统关键性能指标自适应调节方法，属于跟踪控制技术领域。

背景技术：

近年来随着计算机技术、信息化水平日益提高，工业生产过程如化工、冶金、机械等对关键性能指标(如产品质量)的要求越来越高，高精度的跟踪控制得到了工业界和学术界的关注。现有预测控制技术一般要求被控对象的模型一致，但实际应用中由于测量噪声和复杂的系统机理等问题，精确的系统模型难以获取。另一方面随着传感器和计算机储存技术的提高，系统存在大量未被充分利用的测量数据。如何利用这些数据，得到在系统约束下使关键性能指标维持在期望值的最优控制序列，得到了学术界的广泛关注。

另一方面，现有数据驱动的系统关键性能指标自适应调节技术通常需要对大量的参数进行调节，这些参数不仅很大程度影响了性能(控制收敛性、跟踪误差)，而且也为其实际应用带来了很多不便和困难。如利用pid(比例proportion、积分integral、微分differential)控制进行系统关键性能指标自适应调节的技术，不同的参数下的跟踪误差不同，而低跟踪误差对应的参数较难确定，因此综合来说跟踪误差较高。提出一种需要极少先验知识的系统关键性能指标自适应调节技术，提高其泛用性、简化在应用中的实际操作，也是亟待解决的关键问题。

技术实现要素：

本发明为解决现有数据驱动的系统关键性能指标自适应调节技术通常需要对大量的参数进行调节，导致跟踪控制性能差，跟踪误差较高的问题，提供了一种数据驱动的系统关键性能指标自适应调节方法。

本发明所述一种数据驱动的系统关键性能指标自适应调节方法，通过以下技术方案实现：

步骤一、将系统输入设为随机方波信号采集数据，进行所述方法的初始化：

步骤a1、设定时间窗长度，采集被控对象在随机方波激励信号下，过去时间窗长度的输入输出数据，构造截止上一时刻输入矩阵和截止上一时刻输出矩阵，并计算截止上一时刻输入矩阵的广义逆；

步骤a2、利用改进偏最小二乘法计算回归模型初始值进而得到预测模型；

步骤a3、利用被控对象的期望输出计算上一时刻跟踪误差，进而得到待求解二次型问题的相关参数；然后结合约束矩阵得到上一时刻控制序列的增量；

步骤a4、计算并储存当前时刻被控对象的输入；

步骤二、利用对被控对象的当前时刻测量数据进行所述方法的在线更新：

步骤b1、利用步骤a1中获得的截止上一时刻输入矩阵的广义逆和偏最小二乘算法计算当前时刻输入矩阵的广义逆；

步骤b2、采集被控对象的当前时刻输出的测量数据，并更新回归模型和预测模型；

步骤b3、利用当前时刻跟踪误差计算待求解二次型问题的相关参数；并结合约束矩阵得到当前时刻控制序列的增量；

步骤b4、计算并储存下一时刻被控对象的输入；进而由被控对象的输入确定系统关键性能指标；

步骤三、重复执行步骤二，直至系统运行结束。

本发明最为突出的特点和显著的有益效果是：

本发明所涉及的一种数据驱动的系统关键性能指标自适应调节方法，利用被控对象的测量数据，实现对线性系统的关键性能指标自适应调节。本发明方法过程简单，提出的系统关键性能指标自适应调节方法仅需要很少的先验知识，需要预先假设的参数仅有预测步长np和时间窗长度为n，有很高的泛用性；并且本发明基于系统过程数据，不依赖于系统模型，对线性被控对象有很好的跟踪控制性能，跟踪控制预测收敛性好，仿真实验表明，本发明跟踪误差趋近于0。

附图说明

图1为本发明方法流程图；

图2为实施例中系统的跟踪控制效果曲线图；

图3为实施例中预测收敛性曲线图。

具体实施方式

具体实施方式一：结合图1对本实施方式进行说明，本实施方式给出的一种数据驱动的系统关键性能指标自适应调节方法，具体包括以下步骤：

步骤一、将系统输入设为随机方波信号采集数据，进行所述方法的初始化：

步骤a1、设定时间窗长度n，采集被控对象在随机方波激励信号下，过去时间窗长度的输入输出数据，构造截止上一时刻输入矩阵un(k-1)和截止上一时刻输出矩阵yn(k-1)，并计算截止上一时刻输入矩阵的广义逆

步骤a2、利用改进偏最小二乘法计算回归模型初始值m(k-1)进而得到预测模型

步骤a3、利用被控对象的期望输出计算上一时刻跟踪误差e(k-1)，进而得到待求解二次型问题的相关参数hk-1和fk-1；然后结合约束矩阵得到上一时刻控制序列的增量δu(k-1)；

步骤a4、计算并储存当前时刻被控对象的输入(控制器输出)uk；

步骤二、利用对被控对象的当前时刻测量数据进行所述方法的在线更新：

步骤b1、利用步骤a1中获得的截止上一时刻输入矩阵的广义逆和偏最小二乘算法计算当前时刻输入矩阵的广义逆

步骤b2、采集被控对象的当前时刻输出的测量数据yk，并更新回归模型m(k)和预测模型

步骤b3、利用当前时刻跟踪误差计算待求解二次型问题的相关参数hk和fk；并结合约束矩阵得到当前时刻控制序列的增量δu(k)；

步骤b4、计算并储存下一时刻被控对象的输入uk+1；进而由被控对象的输入确定系统关键性能指标；

步骤三、重复执行步骤二，直至系统运行结束。

具体实施方式二：本实施方式与具体实施方式一不同的是，步骤a1中所述截止上一时刻输入矩阵的广义逆具体为：

其中，n为时间窗长度，k表示当前时刻，表示截止上一时刻输入矩阵；表示k时刻被控对象的输入，表示维的实数向量空间；表示截止上一时刻输出矩阵，表示k时刻被控对象的输出，表示n维的实数向量空间；上标“t”表示转置，上标“+”表示广义逆。

其他步骤及参数与具体实施方式一相同。

具体实施方式三：本实施方式与具体实施方式一不同的是，所述预测模型的具体计算过程包括：

利用改进偏最小二乘计算回归模型：

其中，m(k)表示k时刻的回归模型，表示截止k时刻输入矩阵的广义逆；表示截止k时刻输出矩阵；表示k时刻被控对象的输出，表示n维的实数向量空间；

给定预测步长为np，计算预测模型：

其中，n表示被控对象的输入数据的维数，m表示被控对象的输出数据的维数。

其他步骤及参数与具体实施方式一或二相同。

具体实施方式四：本实施方式与具体实施方式一不同的是，所述跟踪误差的具体计算过程为：

其中，e(k)表示k时刻的跟踪误差，ysp(k)表示k时刻被控对象的期望输出，np为预测步长，yk表示k时刻被控对象的输出。

其他步骤及参数与具体实施方式一、二或三相同。

具体实施方式五：本实施方式与具体实施方式三不同的是，所述控制序列的增量的具体计算过程包括：

计算待求解二次型问题的相关参数hk和fk：

其中，q和p均为权重矩阵，e(k)表示k时刻的跟踪误差；

给定符合实际情况(机械结构等)的约束矩阵c1和c2，计算控制序列的增量：

minδu(k)δu^t(k)hkδu(k)+fk^tδu(k)(7)

c1δu(k)≤c2(8)

其中，δu(k)表示k时刻控制序列的增量，δu(k)中的元素δuk表示k时刻被控对象的输入的增量。

其他步骤及参数与具体实施方式一、二、三或四相同。

具体实施方式六：本实施方式与具体实施方式五不同的是，所述控制器输出的具体计算过程为：

uk+1＝uk+δuk(9)

其中，uk+1为k+1时刻的控制器输出，uk为k时刻的控制器输出。

其他步骤及参数与具体实施方式一至五相同。

具体实施方式七：本实施方式与具体实施方式二至六不同的是，步骤b1中所述利用步骤a1中获得的截止上一时刻输入矩阵的广义逆和偏最小二乘算法计算当前时刻输入矩阵的广义逆的具体过程包括：

计算偏最小二乘算法的减量形式的中间矩阵s：

利用和偏最小二乘算法的减量形式计算的广义逆：

计算偏最小二乘算法的增量形式的中间矩阵b、c、d：

c＝uk-dun-1(k-1)(13)

利用和偏最小二乘算法的增量形式计算un(k)的广义逆：

其他步骤及参数与具体实施方式一至六相同。

实施例

采用以下实施例验证本发明的有益效果：

本实施例为连续搅拌加热罐标准工业仿真过程，该被控对象是一个五输入三出输出系统。出于方便，本实施例中以对两个罐的热输入(单位为％)作为输入变量，以罐1的温度作为关键性能指标(输出变量)。被控对象的其他输入变量设置为标准工况的输入(分别为60％，55％和50％)

步骤一、初始化。将系统输入设为随即方波信号采集数据，进行控制算法初始化：

步骤a1、给定时间窗长度n＝50，在系统约束范围内将系统的输入置为随机方波激励信号下的50组输入输出数据，构造输入矩阵u50(k-1)和输出矩阵y50(k-1)；

采用公式(1)计算输入矩阵的广义逆；

步骤a2、利用改进偏最小二乘算法计算回归模型初始值m(k-1)；如公式(2)

给定预测步长np＝50，如公式(3)计算预测模型

步骤a3、利用被控对象的期望输出ysp(k)计算上一时刻跟踪误差e(k-1)；

其中设定点设为：

利用公式(5)、(6)得到待求解二次型问题的相关参数hk-1和fk-1；考虑到连续搅拌加热罐的热输入增量在一个采样间隔变动不超过0.2(δuk≤0.2)，热输入上限100下限为0(0k＜uk＜100)，构造对输入增量的约束如公式(8)，然后合约束矩阵c1、c2，得到上一时刻控制序列的增量δu(k-1)；

步骤a4、计算并储存当前时刻被控对象的输入uk(控制器输出)。

步骤二、在线更新。利用在线对被控对象的测量数据进行在线更新：

步骤b1、如公式(15)，计算当前时刻输入矩阵un(k)的广义逆；

步骤b2、采集被控对象的在线输出测量数据yk，并更新回归模型和预测模型；

步骤b3、利用当前时刻跟踪误差e(k)计算待求解二次型问题的相关参数hk和fk；并结合约束矩阵得到当前时刻控制序列的增量δu(k)；

步骤b4、计算并储存下一时刻被控对象的输入(控制器输出)uk+1，进而由被控对象的输入确定系统关键性能指标；

步骤三、若被控对象仍在运行中，则重复步骤二。

如图2所示为系统的跟踪控制效果；预测模型收敛性如图3所示，其中，测量值拟合(模型拟合准确度)为1时表示期望值与测量值完全拟合，跟踪误差等于0，从图中可以看出，本发明方法的跟踪误差非常小，接近0。

本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。