CN109358511B - 一种数据驱动的系统关键性能指标自适应调节方法 - Google Patents
一种数据驱动的系统关键性能指标自适应调节方法 Download PDFInfo
- Publication number
- CN109358511B CN109358511B CN201811519171.XA CN201811519171A CN109358511B CN 109358511 B CN109358511 B CN 109358511B CN 201811519171 A CN201811519171 A CN 201811519171A CN 109358511 B CN109358511 B CN 109358511B
- Authority
- CN
- China
- Prior art keywords
- input
- moment
- data
- calculating
- control system
- 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
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)
Abstract
本发明提供一种数据驱动的系统关键性能指标自适应调节方法,属于跟踪控制技术领域。本发明首先将系统输入设为随机方波信号采集数据,根据采集的数据、被控对象的期望输出,利用改进偏最小二乘法进行初始化,获得并储存当前时刻被控对象的输入;然后利用对被控对象的当前时刻测量数据和偏最小二乘算法进行本发明方法的在线更新,获得下一时刻被控对象的输入;重复执行上一步骤,直至系统运行结束。本发明解决了现有数据驱动的系统关键性能指标自适应调节技术通常需要对大量的参数进行调节,导致控制性能差,跟踪误差较高的问题。本发明可用于系统关键性能指标的自适应调节。
Description
技术领域
本发明涉及一种系统关键性能指标自适应调节方法,属于跟踪控制技术领域。
背景技术
近年来随着计算机技术、信息化水平日益提高,工业生产过程如化工、冶金、机械等对关键性能指标(如产品质量)的要求越来越高,高精度的跟踪控制得到了工业界和学术界的关注。现有预测控制技术一般要求被控对象的模型一致,但实际应用中由于测量噪声和复杂的系统机理等问题,精确的系统模型难以获取。另一方面随着传感器和计算机储存技术的提高,系统存在大量未被充分利用的测量数据。如何利用这些数据,得到在系统约束下使关键性能指标维持在期望值的最优控制序列,得到了学术界的广泛关注。
另一方面,现有数据驱动的系统关键性能指标自适应调节技术通常需要对大量的参数进行调节,这些参数不仅很大程度影响了性能(控制收敛性、跟踪误差),而且也为其实际应用带来了很多不便和困难。如利用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),并计算截止上一时刻输入矩阵的广义逆
步骤A3、利用被控对象的期望输出计算上一时刻跟踪误差E(k-1),进而得到待求解二次型问题的相关参数Hk-1和fk-1;然后结合约束矩阵得到上一时刻控制序列的增量ΔU(k-1);
步骤A4、计算并储存当前时刻被控对象的输入(控制器输出)uk;
步骤二、利用对被控对象的当前时刻测量数据进行所述方法的在线更新:
步骤B3、利用当前时刻跟踪误差计算待求解二次型问题的相关参数Hk和fk;并结合约束矩阵得到当前时刻控制序列的增量ΔU(k);
步骤B4、计算并储存下一时刻被控对象的输入uk+1;进而由被控对象的输入确定系统关键性能指标;
步骤三、重复执行步骤二,直至系统运行结束。
其中,N为时间窗长度,k表示当前时刻,表示截止上一时刻输入矩阵;表示k时刻被控对象的输入,表示m维的实数向量空间;表示截止上一时刻输出矩阵,表示k时刻被控对象的输出,表示n维的实数向量空间;上标“T”表示转置,上标“+”表示广义逆。
其他步骤及参数与具体实施方式一相同。
具体实施方式三:本实施方式与具体实施方式一不同的是,所述预测模型的具体计算过程包括:
利用改进偏最小二乘计算回归模型:
给定预测步长为np,计算预测模型:
其中,n表示被控对象的输入数据的维数,m表示被控对象的输出数据的维数。
其他步骤及参数与具体实施方式一或二相同。
具体实施方式四:本实施方式与具体实施方式一不同的是,所述跟踪误差的具体计算过程为:
其中,E(k)表示k时刻的跟踪误差,ysp(k)表示k时刻被控对象的期望输出,np为预测步长,yk表示k时刻被控对象的输出。
其他步骤及参数与具体实施方式一、二或三相同。
具体实施方式五:本实施方式与具体实施方式三不同的是,所述控制序列的增量的具体计算过程包括:
计算待求解二次型问题的相关参数Hk和fk:
其中,Q和P均为权重矩阵,E(k)表示k时刻的跟踪误差;
给定符合实际情况(机械结构等)的约束矩阵C1和C2,计算控制序列的增量:
C1ΔU(k)≤C2 (8)
其他步骤及参数与具体实施方式一、二、三或四相同。
具体实施方式六:本实施方式与具体实施方式五不同的是,所述控制器输出的具体计算过程为:
uk+1=uk+Δuk (9)
其中,uk+1为k+1时刻的控制器输出,uk为k时刻的控制器输出。
其他步骤及参数与具体实施方式一至五相同。
具体实施方式七:本实施方式与具体实施方式二至六不同的是,步骤B1中所述利用步骤A1中获得的截止上一时刻输入矩阵的广义逆和偏最小二乘算法计算当前时刻输入矩阵的广义逆的具体过程包括:
计算偏最小二乘算法的减量形式的中间矩阵S:
计算偏最小二乘算法的增量形式的中间矩阵B、C、D:
C=uk-DUN-1(k-1) (13)
其他步骤及参数与具体实施方式一至六相同。
实施例
采用以下实施例验证本发明的有益效果:
本实施例为连续搅拌加热罐标准工业仿真过程,该被控对象是一个五输入三出输出系统。出于方便,本实施例中以对两个罐的热输入(单位为%)作为输入变量,以罐1的温度作为关键性能指标(输出变量)。被控对象的其他输入变量设置为标准工况的输入(分别为60%,55%和50%)
步骤一、初始化。将系统输入设为随即方波信号采集数据,进行控制算法初始化:
步骤A1、给定时间窗长度N=50,在系统约束范围内将系统的输入置为随机方波激励信号下的50组输入输出数据,构造输入矩阵U50(k-1)和输出矩阵Y50(k-1);
采用公式(1)计算输入矩阵的广义逆;
步骤A2、利用改进偏最小二乘算法计算回归模型初始值M(k-1);如公式(2)
步骤A3、利用被控对象的期望输出ysp(k)计算上一时刻跟踪误差E(k-1);
其中设定点设为:
利用公式(5)、(6)得到待求解二次型问题的相关参数Hk-1和fk-1;考虑到连续搅拌加热罐的热输入增量在一个采样间隔变动不超过0.2(Δuk≤0.2),热输入上限100下限为0(0<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。
本发明还可有其它多种实施例,在不背离本发明精神及其实质的情况下,本领域技术人员当可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (7)
1.一种数据驱动的系统关键性能指标自适应调节方法,其特征在于,所述方法具体包括以下步骤:
步骤一、将系统输入设为随机方波信号采集数据,进行所述方法的初始化:
步骤A1、设定时间窗长度,采集连续搅拌加热罐标准工业控制系统在随机方波激励信号下,过去时间窗长度的输入输出数据,构造截止上一时刻输入矩阵和截止上一时刻输出矩阵,并计算截止上一时刻输入矩阵的广义逆;
步骤A2、利用改进偏最小二乘法计算回归模型初始值进而得到预测模型;
步骤A3、利用连续搅拌加热罐标准工业控制系统的期望输出计算上一时刻跟踪误差,进而得到待求解二次型问题的相关参数;然后结合约束矩阵得到上一时刻控制序列的增量;
步骤A4、计算并储存当前时刻连续搅拌加热罐标准工业控制系统的输入;
步骤二、利用对连续搅拌加热罐标准工业控制系统的当前时刻测量数据进行所述方法的在线更新:
步骤B1、利用步骤A1中获得的截止上一时刻输入矩阵的广义逆和偏最小二乘算法计算当前时刻输入矩阵的广义逆;
步骤B2、采集连续搅拌加热罐标准工业控制系统的当前时刻输出的测量数据,并更新回归模型和预测模型;
步骤B3、利用当前时刻跟踪误差计算待求解二次型问题的相关参数;并结合约束矩阵得到当前时刻控制序列的增量;
步骤B4、计算并储存下一时刻连续搅拌加热罐标准工业控制系统的输入;进而由连续搅拌加热罐标准工业控制系统的输入确定系统关键性能指标;
步骤三、重复执行步骤二,直至系统运行结束。
6.根据权利要求5所述一种数据驱动的系统关键性能指标自适应调节方法,其特征在于,所述连续搅拌加热罐标准工业控制系统的输入的具体计算过程为:
uk+1=uk+Δuk
其中,uk+1为k+1时刻的连续搅拌加热罐标准工业控制系统的输入,uk为k时刻的连续搅拌加热罐标准工业控制系统的输入。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811519171.XA CN109358511B (zh) | 2018-12-12 | 2018-12-12 | 一种数据驱动的系统关键性能指标自适应调节方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811519171.XA CN109358511B (zh) | 2018-12-12 | 2018-12-12 | 一种数据驱动的系统关键性能指标自适应调节方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109358511A CN109358511A (zh) | 2019-02-19 |
CN109358511B true CN109358511B (zh) | 2022-04-19 |
Family
ID=65328641
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811519171.XA Active CN109358511B (zh) | 2018-12-12 | 2018-12-12 | 一种数据驱动的系统关键性能指标自适应调节方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109358511B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111158351B (zh) * | 2020-01-19 | 2021-06-01 | 哈尔滨工业大学 | 数据驱动的故障诊断与最优控制系统一体化设计方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100568144C (zh) * | 2008-09-04 | 2009-12-09 | 湖南大学 | 一种未知环境下移动机器人多行为融合自动导航方法 |
US8880321B2 (en) * | 2011-03-07 | 2014-11-04 | Toyota Motor Engineering & Manufacturing North America, Inc. | Adaptive air charge estimation based on support vector regression |
CN102298327B (zh) * | 2011-07-10 | 2013-10-30 | 东华大学 | 免疫数据驱动控制方法及解决涤纶工业丝卷绕系统控制问题的设备 |
CN104156560A (zh) * | 2014-07-12 | 2014-11-19 | 中国矿业大学 | 一种基于SaE-ELM的煤矿多等级突水预测方法 |
CN104991984A (zh) * | 2015-05-15 | 2015-10-21 | 广西大学 | 一种用于煮糖结晶的数据监测方法和系统 |
-
2018
- 2018-12-12 CN CN201811519171.XA patent/CN109358511B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109358511A (zh) | 2019-02-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103268082B (zh) | 一种基于灰色线性回归的热误差建模方法 | |
CN105589448A (zh) | 具有用以补偿模型失配的可调节积分分量的模型预测控制器 | |
CN107621269A (zh) | 光纤陀螺温度漂移误差补偿方法 | |
CN103309233A (zh) | 一种模糊pid控制器的设计方法 | |
CN108828954B (zh) | 气候风洞自适应预测控制系统及其控制方法 | |
CN110442911B (zh) | 一种基于统计机器学习的高维复杂系统不确定性分析方法 | |
CN112462599A (zh) | 一种高性能pid控制参数整定方法、装置及系统 | |
CN109358511B (zh) | 一种数据驱动的系统关键性能指标自适应调节方法 | |
CN106094526B (zh) | 一种广义预测工程化应用于脱硝控制系统的方法 | |
Stebel et al. | General tuning procedure for the nonlinear balance-based adaptive controller | |
CN104330972A (zh) | 一种基于模型自适应的综合预测迭代学习控制方法 | |
CN110687777A (zh) | 一种自耦pi协同控制理论新方法 | |
CN103279030B (zh) | 基于贝叶斯框架的动态软测量建模方法及装置 | |
CN109101759B (zh) | 一种基于正逆响应面法的参数识别方法 | |
CN110909492A (zh) | 一种基于极端梯度提升算法的污水处理过程软测量方法 | |
CN116382093A (zh) | 用于模型未知的非线性系统最优控制方法及设备 | |
CN111310110A (zh) | 一种高维耦合不确定系统混合状态估计方法 | |
Chu et al. | Integrated operation optimization strategy for batch process based on process transfer model under disturbance | |
CN114791334A (zh) | 一种压力传感器的标定简化方法 | |
CN112068423B (zh) | 有限值幂次吸引重复控制器设计方法 | |
CN104133371A (zh) | 一种室温控制对象的自抗扰比例积分控制算法 | |
CN106200619B (zh) | 附带控制器输出约束的pi控制回路性能评估方法 | |
Ruusu et al. | Sliding mode SISO control of model parameters for implicit dynamic feedback estimation of industrial tracking simulation systems | |
CN115455347A (zh) | 一种mems加速度计混合误差建模补偿方法 | |
JP5561519B2 (ja) | プラントシミュレータ |
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 |