CN111585822A - 一种网络系统中数据丢包的提升模型预测补偿方法 - Google Patents

一种网络系统中数据丢包的提升模型预测补偿方法 Download PDF

Info

Publication number
CN111585822A
CN111585822A CN202010418993.XA CN202010418993A CN111585822A CN 111585822 A CN111585822 A CN 111585822A CN 202010418993 A CN202010418993 A CN 202010418993A CN 111585822 A CN111585822 A CN 111585822A
Authority
CN
China
Prior art keywords
data
output
model
linear
packet loss
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.)
Pending
Application number
CN202010418993.XA
Other languages
English (en)
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.)
Qingdao University of Science and Technology
Original Assignee
Qingdao University of Science and Technology
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 Qingdao University of Science and Technology filed Critical Qingdao University of Science and Technology
Priority to CN202010418993.XA priority Critical patent/CN111585822A/zh
Publication of CN111585822A publication Critical patent/CN111585822A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L41/00Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
    • H04L41/14Network analysis or design
    • H04L41/145Network analysis or design involving simulating, designing, planning or modelling of a network
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L41/00Arrangements for maintenance, administration or management of data switching networks, e.g. of packet switching networks
    • H04L41/14Network analysis or design
    • H04L41/147Network analysis or design for predicting network behaviour

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Data Exchanges In Wide-Area Networks (AREA)

Abstract

本发明公开一种网络系统中数据丢包的批量预测补偿方法,属于智能控制领域。本发明主要针对网络控制系统中输入侧和输出侧数据均有丢失的问题,所提出的控制方案如下:建立虚拟线性数据模型,在迭代域内描述多输入多输出线性重复系统的输入输出动态;利用线性数据模型建立一个基于提升技术的批量迭代预测模型,它实际上存在于计算机中,用于丢失数据的预测和补偿;提出一种网络系统中数据丢包的提升模型预测补偿方法;将所提方法推广到未知的非线性非仿射系统中。本发明公开的一种网络系统中数据丢包的提升模型预测补偿方法,主要针对网络控制系统中数据丢包问题,在数据驱动框架下,通过建立基于提升技术的批量迭代预测模型,对丢失数据进行批量预测并补偿,计算实时性高,能够有效避免网络系统由于数据丢失而导致控制性能的恶化,既可以应用于线性系统也可以应用于非线性系统。

Description

一种网络系统中数据丢包的提升模型预测补偿方法
技术领域
本发明涉及网络控制技术领域,具体涉及一种网络系统中数据丢包的提升模型预测补偿方法。
背景技术
迭代学习控制(iterative learning control,ILC)最适合于具有可重复特性的过程,并且能够在有限的时间间隔内执行完美的跟踪。ILC使用以前操作的信息不断优化控制行为,随着迭代的增加更好地执行任务。由于其良好的控制性能和简单的控制结构,ILC得到了广泛的研究和应用。
近年来,网络控制系统以其成本低、安装简单、布线少、维护方便、可靠性高等优点,在各个领域得到了广泛的应用。然而,在执行重复任务的网络控制系统中ILC的实现仍然是一个开放的问题。此外,当网络控制系统受到时延、数据丢失、通信噪声等影响时,网络控制系统的性能会下降。因此,在网络控制系统中,如何提高数据噪声下所提控制方法的鲁棒性是迭代学习控制成功实现的关键。
在文献“数据丢失对迭代学习控制的影响分析”中,针对一类线性系统,提出一种存在数据丢失时的鲁棒迭代学习控制器设计方法,并将控制器设计问题转化为求取线性矩阵不等式的可行解;在文献“包含数据丢失的迭代学习控制系统收敛性分析”中,分析了数据丢失对测控信号无线传输迭代学习控制系统收敛性能的影响,并得出测量信号随机丢失会降低系统的收敛速度,而控制信号的随机丢失会显著影响系统的鲁棒收敛性;文献“网络控制系统中一种新形式的开闭环迭代学习控制”通过借鉴网络控制系统的分层结构形式,提出了一种远程开环本地闭环的迭代学习控制形式对控制系统进行前馈补偿同时降低了网络丢包对系统的实时影响。针对非线性仿射系统,文献“数据丢失下多智能体系统迭代学习跟踪控制”研究了据丢失下非线性多智能体系统的一致性跟踪问题;文献“测量数据丢失的一类非线性系统迭代学习控制”针对存在输出测量数据丢失的一类非线性系统,研究了P型迭代学习控制算法的收敛性问题。
综上所述,我们可以发现:大多数用于数据噪声的ILC方法仅限于线性系统和仿射非线性系统,换句话说,这些方法是基于模型的,需要被控对象的线性模型或仿射结构,由于未建模动态的存在,在实际系统中应用这些方法可能会遇到许多困难;现有的网络化ILC方案大多采用P型或D型学习律,只要选择得当,学习增益保持不变,难以应对外部干扰和测量噪声引起的被控对象的大变化;此外,当模型信息无法获得时,很难选择合适的学习增益;虽然已经提出了一些基于补偿的算法,但丢失的数据大多由以前的最新数据进行补偿;尽管也有一些算法是利用相应的预测模型对丢失的数据进行补偿的,但是预测数据的产生实时性差。
本发明考虑输入端和输出端的均有数据丢包的情况,提出一种基于提升模型的预测补偿算法,其中,对丢失数据的预测和补偿是批量进行的,可在控制器中直接调用使用,具有计算实时性高的特点,这是本发明的重要创新之处。
发明内容
本发明的目的是在数据驱动的框架下研究一种数据丢包的提升模型预测补偿方法,以解决网络控制系统中输入侧和输出侧数据均有丢失的问题。
为了实现上述目的,本发明采用如下技术方案,一种网络系统中数据丢包的提升模型预测补偿方法,包括以下步骤:
步骤1、建立一个虚拟线性数据模型,在迭代域内描述多输入多输出线性重复系统的输入输出动态;
步骤2、利用线性数据模型建立一个基于提升技术的批量迭代预测模型,它实际上存在于计算机中,用于预测和补偿输出数据丢失的;
步骤3、提出一种网络系统中数据丢包的提升模型预测补偿方法;
步骤4、将所提方法推广到未知的非线性非仿射系统中;
进一步地,步骤1中所述的建立一个虚拟线性数据模型,主要包括以下步骤:
步骤1.1、考虑一个重复地多输入多输出线性时变网络系统:
Figure BDA0002496201540000031
其中,yk(t)∈Rn表示系统的输出,uk(t)∈Rn表示系统的控制输入,xk(t)∈Rm是系统的状态,A(t)∈Rm×m,B(t)∈Rm×n和C(t)∈Rn×m表示未知的系统矩阵,t∈{0,1,…,N-1}代表时间,其中N表示终端时间,k∈{0,1,…}是迭代次数;
步骤1.2、由上述系统,我们可以得到下面的等式:
Figure BDA0002496201540000032
其中,
Figure BDA0002496201540000033
uk(t)=[u1,k(t),u2,k(t),…,un,k(t)]T∈Rn×1
Figure BDA0002496201540000034
i∈{0,…,t}。
该系统需要满足如下假设:
假设1:系统的初始状态xk(0)在每次迭代都是相同的,即xk(0)=x0,x0是一个常数向量。
步骤1.3、由上述等式利用微分中值定理,可建立如下的虚拟线性数据模型:
Figure BDA0002496201540000035
其中,
Figure BDA0002496201540000041
Figure BDA0002496201540000042
Figure BDA0002496201540000043
Δ是一个迭代差分算子,即ΔUk(N-1)=Uk(N-1)-Uk-1(N-1);
又进一步地,步骤2中所述的利用线性数据模型建立一个基于提升技术的批量迭代预测模型,主要包括以下步骤:
步骤2.1、考虑迭代估计算法:
Figure BDA0002496201540000044
其中,
Figure BDA0002496201540000045
表示
Figure BDA0002496201540000046
的估计,0<η<2和μ>0分别是步长因子和权重,αk是一个Bernoulli随机变量用来描述随机的数据丢失,其值随机取为0或者1,当αk=0时,表示数据丢失;当αk=1时,表示传输正常,数据未丢失;
为使算法具有更强的跟踪时变参数的能力,设计如下重置算法:
Figure BDA0002496201540000047
如果
Figure BDA0002496201540000048
Figure BDA0002496201540000049
其中,
Figure BDA00024962015400000410
表示第一次迭代的初始参数估计,ε是一个正的常数;
步骤2.2、建立迭代预测模型:
Figure BDA00024962015400000411
其中,
Figure BDA00024962015400000412
表示第k次迭代预测的输出;
又进一步地,步骤3中所述的提出一种网络系统中数据丢包的提升模型预测补偿方法,主要包括以下步骤:
步骤3.1、计算补偿输出:
Figure BDA0002496201540000051
其中,
Figure BDA0002496201540000052
表示第k次迭代补偿的输出;
步骤3.2、计算得到的控制输入:
Figure BDA0002496201540000053
其中,
Figure BDA0002496201540000054
代表计算得到的控制输入,
Figure BDA0002496201540000055
yd(t+1)∈Rn是期望的输出轨迹,
Figure BDA0002496201540000056
对应于
Figure BDA0002496201540000057
的每一个分量元素,ρ>0和λ>0分别是步长因子和权重;
步骤3.3、引入Bernoulli随机变量βk(t),提出控制律:
Figure BDA0002496201540000058
其中,βk(t)是一个Bernoulli随机变量用来描述随机的数据丢失,其值随机取为0或者1,当βk(t)=0时,表示数据丢失,当βk(t)=1时,表示数据未丢失;
又进一步地,步骤4中所述的将所提方法推广到未知的非线性非仿射系统中,主要包括以下步骤:
步骤4.1、考虑多输入多输出非线性非仿射离散时间系统:
yk(t+1)=f(yk(t),…,yk(t-ny),uk(t),…,uk(t-nu))
其中,yk(t)∈Rn和uk(t)∈Rn分别是系统的输出和输入,f(·)∈Rn代表未知的连续可微非线性函数,t∈{0,…,N-1},N是一个正整数,k是迭代次数,0<ny<N和0<nu<N代表系统的阶数;
上述多输入多输出非线性非仿射离散时间系统需要满足以下2个假设:
假设2:系统的初始输出保持不变,即满足yk(0)=y0,其中,y0是常数向量;
假设3:
Figure BDA0002496201540000059
是连续有界的,即满足
Figure BDA00024962015400000510
其中
Figure BDA00024962015400000511
是一个正的常数;
步骤4.2、上述多输入多输出非线性非仿射离散时间系统被等价地转化为如下的虚拟线性数据模型:
Figure BDA0002496201540000061
其中,
Figure BDA0002496201540000062
Figure BDA0002496201540000063
Figure BDA0002496201540000064
Δ是一个迭代差分算子,即ΔYk=Yk-Yk-1和ΔUk=Uk-Uk-1
步骤4.3、与步骤2–步骤3设计过程类似,针对非线性网络系统的数据丢包的提升模型预测补偿方法如下:
Figure BDA0002496201540000065
本发明在数据驱动框架下提出了一种网络系统中数据丢包的提升模型预测补偿方法,解决了输入和输出端数据均有丢失的问题,从而避免了由于数据丢失而导致控制性能的恶化。本发明的方法具有以下优点:
1、本发明的方法,可以批量地预测和补偿数据丢包情况,是一种前馈的方法,对丢失数据的预测和补偿可在控制器中直接调用使用,计算实时性高;
2、本发明的方法是基于数据驱动的,与精确的系统模型无关;
3、尽管本发明所提方法的设计和分析是从多输入多输出重复线性系统开始的,但是可以很容易地将结果推广到多输入多输出非线性非仿射系统;
结合附图阅读本发明的具体实施方式后,本发明的其他特点和优点将变得更加清楚。
附图说明
图1为本发明提出的网络系统中数据丢包的提升模型预测补偿方法在例1中不同情况下的跟踪误差比较图;
图2为本发明提出的网络系统中数据丢包的提升模型预测补偿方法在例1中不同情况下系统输出的比较图;
图3为本发明提出的网络系统中数据丢包的提升模型预测补偿方法在例1中不同丢包率下跟踪误差的收敛性比较图;
图4为本发明提出的网络系统中数据丢包的提升模型预测补偿方法在例2中不同情况下的跟踪误差比较图;
图5为本发明提出的网络系统中数据丢包的提升模型预测补偿方法在例2中不同情况下系统输出的比较图;
图6为本发明提出的网络系统中数据丢包的提升模型预测补偿方法在例2中不同丢包率下的跟踪误差收敛性比较图;
具体实施方式
为了更好地说明本发明的目的和优点,下面结合实施例和相应附图对发明内容做进一步详细说明。
考虑一个重复地多输入多输出线性时变网络系统:
Figure BDA0002496201540000071
其中,yk(t)∈Rn表示系统的输出,uk(t)∈Rn表示系统的控制输入,xk(t)∈Rm是系统的状态,A(t)∈Rm×m,B(t)∈Rm×n和C(t)∈Rn×m表示未知的系统矩阵,t∈{0,1,…,N}代表时间,其中N表示终端时间,k∈{0,1,…}是迭代次数;
由上述系统(a1),我们可以得到下面的等式:
Figure BDA0002496201540000081
其中,
Figure BDA0002496201540000082
uk(t)=[u1,k(t),u2,k(t),…,un,k(t)]T∈Rn×1
Figure BDA0002496201540000083
i∈{0,…,t};
该系统需要满足如下假设:
假设1:系统的初始状态xk(0)在每次迭代都是相同的,即xk(0)=x0,x0是一个常数向量;
由上述等式(a2)利用微分中值定理,可建立如下的虚拟线性数据模型:
Figure BDA0002496201540000084
其中,
Figure BDA0002496201540000085
Figure BDA0002496201540000086
Figure BDA0002496201540000087
Δ是一个迭代差分算子,即ΔUk=Uk-Uk-1
由于A(t)、B(t)和C(t)都是实矩阵,并且t∈{0,…,N}是有限的,因此,矩阵A(t)、B(t)和C(t)都是有界的;很明显,
Figure BDA0002496201540000088
也是有界的;又由于
Figure BDA0002496201540000089
的维数是有限的,所以
Figure BDA00024962015400000810
也是有界的;
虚拟线性数据模型是在迭代域中建立的,只利用输入输出数据,不依赖于被控系统的任何显式模型;因此,虚拟线性数据模型是有目的的,适用于非线性系统;而且,它没有特殊的物理意义,只描述了原被控对象在两次连续迭代之间的数据关系;此外,线性数据模型实际上存在于计算机中,并在随后的控制器设计中用作迭代预测模型;
由于(a1)式中所有的系数矩阵都是未知的,传统的基于模型的控制方案不能直接用于该系统;同时,从(a3)式可以看到,所有未知的系统矩阵都被压缩成一个参数矩阵
Figure BDA0002496201540000091
利用虚拟线性数据模型,
Figure BDA0002496201540000092
可由如下迭代估计算法进行估计:
Figure BDA0002496201540000093
其中,
Figure BDA0002496201540000094
表示第k次迭代
Figure BDA0002496201540000095
的估计,0<η<2和μ>0分别是步长因子和权重,αk是一个Bernoulli随机变量用来描述随机的数据丢失,其值随机取为0或者1,当αk=0时,表示数据丢失,当αk=1时,表示数据传输正常;
为使算法(a4)具有更强的跟踪时变参数的能力,设计如下重置算法:
Figure BDA0002496201540000096
如果
Figure BDA0002496201540000097
Figure BDA0002496201540000098
其中,
Figure BDA0002496201540000099
表示第一次迭代的初始参数估计,ε是一个正的常数;
当控制系统通过网络传输时,计算的输入和测量的输出通过网络传输,由于网络的拥塞和带宽的限制,I/O数据的传输常常出现延迟和丢失;因此,本发明是针对控制器端和执行器端都存在数据丢失的情况设计网络系统中数据丢包的提升模型预测补偿方法。
根据虚拟线性数据模型(a3)以及参数估计算法(a4)和(a5),建立如下的迭代预测模型:
Figure BDA00024962015400000910
其中,
Figure BDA00024962015400000911
表示第k次迭代预测的输出。
补偿的输出由下式计算:
Figure BDA0002496201540000101
其中,
Figure BDA0002496201540000102
表示第k次迭代补偿的输出;
计算得到的控制输入:
Figure BDA0002496201540000103
其中,
Figure BDA0002496201540000104
代表计算得到的控制输入,
Figure BDA0002496201540000105
yd(t+1)∈Rn是期望的输出轨迹,
Figure BDA0002496201540000106
Figure BDA0002496201540000107
的每一个分量元素,ρ>0和λ>0分别是步长因子和权重;
基于预测补偿的控制输入:
Figure BDA0002496201540000108
其中,βk(t)是一个Bernoulli随机变量用来描述随机的数据丢失,其值随机取为0或者1,当βk(t)=0时,表示数据丢失,当βk(t)=1时,表示数据未丢失;
从所提出的网络系统中数据丢包的预测补偿方法(a4)-(a9)可以看出,控制器(a8)和估计器(a4)位于同一网络侧,它们都使用从受控对象发送的I/O数据,这些数据通过使用Bernoulli随机变量αk反映出来;同时,在输出补偿器(a7)和输入补偿器(a9)中,两个Bernoulli随机变量αk和βk(t)用于处理输入和输出信号中的数据丢失;(a7)式是用于批量预测系统输出的迭代预测模型;
考虑在输入端和输出端都存在随机数据丢失的重复线性时变系统(a1)的ILC问题,如果线性时变系统(a1)满足假设1,且控制器参数选取正确,则所提出的网络系统中数据丢包的预测补偿方法(a4)-(a9)能够保证:
Figure BDA0002496201540000109
是有界的;当迭代次数k趋于无穷时,跟踪误差ek(t+1)收敛到一个很小的界。
定义参数估计误差
Figure BDA0002496201540000111
(a4)式等号两端同时减去
Figure BDA0002496201540000112
可得
Figure BDA0002496201540000113
其中,I∈RnN×nN是一个单位阵;
由(a10),很明显得到
Figure BDA0002496201540000114
考虑到下式是成立的
Figure BDA0002496201540000115
因为0<η<2,μ>0以及αk∈{0,1},所以必定存在一个常数0<d1<1,使得
Figure BDA0002496201540000116
因为
Figure BDA0002496201540000117
是给定有界的,因此由(a13)很容易得到
Figure BDA0002496201540000118
再由
Figure BDA0002496201540000119
的定义以及
Figure BDA00024962015400001110
的有界性可以得到
Figure BDA00024962015400001111
是有界的,还可得到
Figure BDA00024962015400001112
也是有界的;
定义输出估计误差
Figure BDA00024962015400001113
由(a3)和(a6),可得
Figure BDA00024962015400001114
通过解方程(a14),可得
Figure BDA00024962015400001115
其中,当k<0时,输出估计误差被设置为零;
根据(a13),已经得到了
Figure BDA00024962015400001116
因此,
Figure BDA00024962015400001117
再进一步,可得
Figure BDA0002496201540000121
是有界的;再由(a15)式,可得
Figure BDA0002496201540000122
是有界的,并且
Figure BDA0002496201540000123
的每一个元素σk(t+1)都是有界的;为了方便后续分析,记Lσ是σk(t+1)的上界,即||σk(t+1)||≤Lσ
因为
Figure BDA0002496201540000124
是有界的,并且βk(t)∈{0,1},存在一个合适的λ使得
Figure BDA0002496201540000125
Figure BDA0002496201540000126
Figure BDA0002496201540000127
根据(a8)和(a9),可得
Figure BDA0002496201540000128
Figure BDA0002496201540000129
C=[0,…,0,I2]T∈R(n*(t+1))×n,其中,0∈Rn×n是一个零矩阵,I2∈Rn×n是单位阵,则
Figure BDA00024962015400001210
其中,
Figure BDA00024962015400001211
如果t是小于零的,那么Δuk(t)=O,O=[0,…,0]T∈Rn
Ak(t)的特征方程为
Figure BDA0002496201540000131
根据(a18)式,通过选取合适的ρ,可得
Figure BDA0002496201540000132
注意如下引理:
引理1:记
Figure BDA0002496201540000133
如果
Figure BDA0002496201540000134
那么矩阵A的谱半径满足s(A)<1;
因此,可得
Figure BDA0002496201540000135
(a23)式意味着
Figure BDA0002496201540000136
进一步选择合适的ρ和λ使得
Figure BDA0002496201540000137
考虑到如下引理:
引理2:A∈Rt×t,对于
Figure BDA0002496201540000138
在赋范向量空间ν上存在一个合适的矩阵范数使得||Ak(t)||ν≤s(A)+δ。
可得
Figure BDA0002496201540000139
其中,||·||ν指的是合适的矩阵范数,0<d5,k(t)<0.5。
通过(a20)式,可得下式
Figure BDA00024962015400001310
然后,根据(a16),(a24)和(a25),可推得
Figure BDA0002496201540000141
Figure BDA0002496201540000142
的定义,则
Figure BDA0002496201540000143
其中,
Figure BDA0002496201540000144
由于
Figure BDA0002496201540000145
是有界的,且t是有界的,因此
Figure BDA0002496201540000146
是有界的,为方便后续讨论,记Lu
Figure BDA0002496201540000147
的上界,即
Figure BDA0002496201540000148
根据(a20),重写(a27),可得
Figure BDA0002496201540000149
根据方程(a28),可得
Figure BDA00024962015400001410
考虑到下式是成立的
Figure BDA00024962015400001411
根据(a17)和(a30),可得
Figure BDA00024962015400001412
通过选取合适的ρ和λ使得
Figure BDA0002496201540000151
其中,0<d6,k(t)=1-ρd3,k(t)<1是一个常数。
对于任意的t和k,||σk(t+1)||≤Lσ是有界的。又αk∈{0,1},所以||αk-1σk-1(t+1)||+||αkσk(t+1)||≤2Lσ是有界的。则由(a26)、(a29)以及(a32)式可得
Figure BDA0002496201540000152
定义||emax,k||=max{||ek(t+1)||}。假设||emax,k||=||ek(τ)||,其中τ是{1,…,N}上的时刻,根据(a33),则
Figure BDA0002496201540000153
定义d7,k(t)=ρd3,k(t),重写(a34),则
Figure BDA0002496201540000154
由(a32)式,可知0<d7,k(t)=1-d6,k(t)<1,则
Figure BDA0002496201540000155
根据(a24)式,很明显0<d5,k(t)<0.5,并且d7,k(t)是在(0,1)域内的,则
Figure BDA0002496201540000156
因此,(a35)式可重写为
Figure BDA0002496201540000157
因此,
Figure BDA0002496201540000158
Figure BDA0002496201540000159
都是有界的,很明显
Figure BDA00024962015400001510
是收敛于一个小的界的。而且由
Figure BDA00024962015400001511
的定义,则
Figure BDA0002496201540000161
(a39)式可被进一写为
Figure BDA0002496201540000162
在(a40)式等号两端取范数,可得
Figure BDA0002496201540000163
因为
Figure BDA0002496201540000164
Figure BDA0002496201540000165
都是有界的,所以当k→∞时,||ek(t+1)||可以收敛于一个有界的范围,即
Figure BDA0002496201540000166
将所提方法推广到未知的非线性非仿射系统中,考虑多输入多输出非线性非仿射离散时间系统
yk(t+1)=f(yk(t),…,yk(t-ny),uk(t),…,uk(t-nu)) (a42)
其中,yk(t)∈Rn和uk(t)∈Rn分别是系统的输出和输入;f(·)∈Rn代表未知的连续可微非线性函数;t∈{0,…,N},N是一个正整数;k是迭代次数;0<ny<N和0<nu<N代表系统的阶数。
上述多输入多输出非线性非仿射离散时间系统需要满足以下2个假设:
假设2:系统的初始输出保持不变,即满足yk(0)=y0,其中,y0是常数向量。
假设3:
Figure BDA0002496201540000167
是连续有界的,即满足
Figure BDA0002496201540000168
其中
Figure BDA0002496201540000169
是一个正的常数。
上述多输入多输出非线性非仿射离散时间系统(a42)可被表述成
yk(t+1)=gt(yk(0),uk(0),uk(1),…,uk(t)) (a43)
其中,gi(·),i={0,…,t}是f(·)的复合函数。
定义
Figure BDA0002496201540000171
uk(i)=[u1,k(i),u2,k(i),…,un,k(i)]T∈Rn,i∈{0,…,t}。则(a43)可重写为
进一步,利用微分中值定理以及假设2,上述多输入多输出非线性非仿射离散时间系统(a42)被等价地转化为如下的虚拟线性数据模型:
Figure BDA0002496201540000173
其中,
Figure BDA0002496201540000174
Figure BDA0002496201540000175
Figure BDA0002496201540000176
Figure BDA0002496201540000177
Δ是一个迭代差分算子,即ΔYk=Yk-Yk-1和ΔUk=Uk-Uk-1
针对非线性网络系统的数据丢包的提升模型预测补偿方法如下:
Figure BDA0002496201540000178
考虑如下假设:
假设4:
Figure BDA0002496201540000179
的元素是非零的,并且元素的符号都是不变的,即
Figure BDA00024962015400001710
Figure BDA0002496201540000181
i={0,…,t},l={1,…,n},其中
Figure BDA0002496201540000182
是一个正的常数。在本文中,不失一般性的假设
Figure BDA0002496201540000183
针对多输入多输出非线性系统(a42),在假设2-4满足的条件下,其中随机数据丢包发生在输入和输出端;通过选择合适的控制器参数,所提出的非线性网络系统的数据丢包的提升模型预测补偿方法(a46)能够保证
Figure BDA0002496201540000184
在所有迭代和时间上都是有界的,并且当k趋于无穷大时,跟踪误差ek(t+1)收敛到有界的范围。
实际上,
Figure BDA0002496201540000185
表示控制系统在迭代域中的微分信号,它在迭代方向上变化缓慢,特别是当Δuk(t)不太大时,换言之,
Figure BDA0002496201540000186
行为的数值变化对迭代变化因子不敏感;因此,在收敛性分析中可以把
Figure BDA0002496201540000187
看作是迭代不变量,也就是说,
Figure BDA0002496201540000188
可以看作是迭代不变量。后续分析与线性系统的情况类似,不再赘述。
为了验证本发明方法的正确性,对本发明的方法进行了以下仿真:
实施例子1:考虑一个离散时间线性网络控制系统
Figure BDA0002496201540000189
其中,
Figure BDA00024962015400001810
C=[1 0.5],t∈{0,…,20}。期望轨迹为yd(t)=10-6(t-1)3(4-0.03(t-1))
在本例中,我们分三种情况对本发明进行验证。
情况一:测量输出和计算输入的数据损失率均为10%。
初始值设置为yk(0)=0,
Figure BDA00024962015400001811
u0(t)=0,
Figure BDA00024962015400001812
控制器参数设置为λ=0.01,ρ=0.003,η=0.05,μ=0.2,ε=0.00001。应用所提的网络系统的数据丢包的提升模型预测补偿方法(a4)–(a9),跟踪误差的收敛性如图1所示,其中,纵轴代表最大的绝对跟踪误差,其定义为
Figure BDA00024962015400001813
第500次迭代时系统输出的跟踪性能如图2所示。从图1和图2可以看出,在数据丢失率为10%的影响下,所提出的数据丢包的提升模型预测补偿方法能够实现迭代方向上的渐近收敛。
情况二:测量输出和计算输入的数据损失率为0%,即在输入和输出侧均没有数据丢失。
仿真结果如图1和图2中的黑色实线所示。从图中可以发现,emax(k)随着迭代的增加逐渐收敛,系统输出可以完美地跟踪期望的轨迹,在数据传输正常的情况下,具有良好的收敛速度。
通过比较,可以得出结论,所提的网络系统的数据丢包的提升模型预测补偿方法可以补偿丢包带来的不利影响,从而在略微降低收敛速度的情况下实现良好的跟踪性能。
情况三:对于测量输出和计算输入,考虑不同的数据损失率,10%、15%和20%。
这种情况的目的是验证所提出的数据丢包的预测补偿方法在数据丢失率不同的情况下的有效性。
为了比较的公平,设置相同的控制器参数以及初始值,然后应用所提出的提升模型预测补偿方法(a4)–(a9),仿真结果如图3所示。从图3可以看出,所提出的预测补偿方法在处理不同的数据丢失率时是有效的,但是当数据丢失率增加时,控制性能会变差。
实施例子2:考虑一个具有如下形式的非线性网络系统
Figure BDA0002496201540000191
该系统的期望轨迹为yd(t+1)=5sin(2tπ/50)+cos(2tπ/10),t∈{0,1,…,99};
在本例中分三种情况对本发明方法进行验证。
情况一:测量输出和计算输入的数据损失率都是10%。
初始值被设置为yk(0)=0,
Figure BDA0002496201540000201
u0(t)=1,
Figure BDA0002496201540000202
控制器参数被选为λ=4,ρ=0.02,η=1.5,μ=0.1,ε=0.00001。最大跟踪误差曲线和第400次迭代时的系统输出分别如图4和图5所示。可以看出,所提出的预测补偿方法实现了跟踪误差沿迭代方向的渐近收敛,系统输出对期望轨迹的跟踪效果令人满意。
情况二:测量输出和计算输入的数据丢失率均为0%,即输入侧和输出侧均不发生数据丢失。
在相同的仿真条件下,应用所提的预测补偿方法,跟踪误差为图4中的黑色实线所示,系统输出如图5中的虚线所示。可以看出,该系统可以达到完美收敛,并具有良好的速度。
相比之下,说明所提的预测补偿方法可以在较小的数据丢失率下补偿数据丢失对收敛速度的不利影响。
情况三:对于测量输出和计算输入,考虑不同的数据包丢失率,分别为10%、15%和20%;
本案例情况的目的是测试设计的针对非线性系统的提升模型预测补偿方法在不同数据丢失率下的广泛有效性。仿真结果如图6所示。
从图6可以得出结论,所提出的网络系统中数据丢包的提升模型预测补偿方法在处理不同的数据丢失率时是有效的,但是当数据丢失率增加时,收敛性能变差。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行了进一步的详细说明,所应理解的是,实施例仅用于解释本发明,并不用于限定本发明的保护范围,凡与本发明在相同原理和构思条件下修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (5)

1.一种网络系统中数据丢包的提升模型预测补偿方法,其特征在于:所述方法包括以下步骤:
步骤1、建立一个虚拟线性数据模型,在迭代域内描述多输入多输出线性重复系统的输入输出动态;
步骤2、利用线性数据模型建立一个基于提升技术的批量迭代预测模型,用于预测和补偿输出数据的丢失;
步骤3、建立网络系统中数据丢包的提升模型预测补偿方法;
步骤4、将所提方法推广到未知的非线性非仿射系统中。
2.根据权利要求1所述方法,其特征在于:步骤1中所述的建立一个虚拟线性数据模型主要包括以下步骤:
步骤1.1、考虑一个重复的多输入多输出线性时变网络系统:
Figure FDA0002496201530000011
其中,yk(t)∈Rn表示系统的输出,uk(t)∈Rn表示系统的控制输入,xk(t)∈Rm是系统的状态,A(t)∈Rm×m,B(t)∈Rm×n和C(t)∈Rn×m表示未知的系统矩阵;t∈{0,1,…,N}代表时间,其中N表示终端时间;k∈{0,1,…}是迭代次数;
步骤1.2、由上述系统(1),我们可以得到下面的等式:
Figure FDA0002496201530000012
其中,
Figure FDA0002496201530000013
uk(t)=[u1,k(t),u2,k(t),…,un,k(t)]T∈Rn×1
Figure FDA0002496201530000014
该系统需要满足如下假设:
假设1:系统的初始状态xk(0)在每次迭代都是相同的,即xk(0)=x0,x0是一个常数向量;
步骤1.3、由上述等式(2)利用微分中值定理,可建立如下的虚拟线性数据模型:
Figure FDA0002496201530000021
其中,
Figure FDA0002496201530000022
Figure FDA0002496201530000023
Figure FDA0002496201530000024
Δ是一个迭代差分算子,即ΔUk=Uk-Uk-1
3.根据权利要求1所述方法,其特征在于:步骤2中所述的利用线性数据模型建立一个基于提升技术的批量迭代预测模型,主要包括以下步骤:
步骤2.1、提出迭代估计算法:
Figure FDA0002496201530000025
其中,
Figure FDA0002496201530000026
表示
Figure FDA0002496201530000027
的估计,0<η<2和μ>0分别是步长因子和权重,αk是一个Bernoulli随机变量用来描述随机的数据丢失,其值随机取为0或者1,当αk=0时,表示数据丢失;当αk=1时,表示传输正常,数据未丢失;
为使算法(4)具有更强的跟踪时变参数的能力,设计如下重置算法:
Figure FDA0002496201530000028
如果
Figure FDA0002496201530000029
Figure FDA00024962015300000210
其中,
Figure FDA00024962015300000211
表示第一次迭代的初始参数估计,ε是一个正的常数;
步骤2.2、建立迭代预测模型:
Figure FDA00024962015300000212
其中,
Figure FDA00024962015300000213
表示第k次迭代预测的输出。
4.根据权利要求1所述方法,其特征在于:步骤3中所述的提出一种网络系统中数据丢包的提升模型预测补偿方法,主要包括以下步骤:
步骤3.1、计算补偿输出:
Figure FDA0002496201530000031
其中,
Figure FDA0002496201530000032
表示第k次迭代补偿的输出;
步骤3.2、计算得到的控制输入:
Figure FDA0002496201530000033
其中,
Figure FDA0002496201530000034
代表计算得到的控制输入;
Figure FDA0002496201530000035
yd(t+1)∈Rn是期望的输出轨迹,
Figure FDA0002496201530000036
对应于
Figure FDA0002496201530000037
的每一个分量元素,ρ>0和λ>0分别是步长因子和权重;
步骤3.3、引入Bernoulli随机变量βk(t),提出控制律:
Figure FDA0002496201530000038
其中,βk(t)是一个Bernoulli随机变量用来描述随机的数据丢失,其值随机取为0或者1,当βk(t)=0时,表示数据丢失,当βk(t)=1时,表示数据未丢失。
5.根据权利要求1所述方法,其特征在于:步骤4中所述的将所提方法推广到未知的非线性非仿射系统中,主要包括以下步骤:
步骤4.1、考虑多输入多输出非线性非仿射离散时间系统:
yk(t+1)=f(yk(t),…,yk(t-ny),uk(t),…,uk(t-nu)) (10)
其中,yk(t)∈Rn和uk(t)∈Rn分别是系统的输出和输入;f(·)∈Rn代表未知的连续可微非线性函数;t∈{0,…,N},N是一个正整数;k是迭代次数;0<ny<N和0<nu<N代表系统的阶数;
上述多输入多输出非线性非仿射离散时间系统(10)需要满足以下2个假设:
假设2:系统的初始输出保持不变,即满足yk(0)=y0,其中,y0是常数向量;
假设3:
Figure FDA0002496201530000041
是连续有界的,即满足
Figure FDA0002496201530000042
其中
Figure FDA0002496201530000043
是一个正的常数;
步骤4.2、上述多输入多输出非线性非仿射离散时间系统(10)被等价地转化为如下的虚拟线性数据模型:
Figure FDA0002496201530000044
其中,
Figure FDA0002496201530000045
Figure FDA0002496201530000046
Figure FDA0002496201530000047
Δ是一个迭代差分算子,即ΔYk=Yk-Yk-1和ΔUk=Uk-Uk-1
步骤4.3、与步骤2–步骤3设计过程类似,提出的针对非线性网络系统的数据丢包的提升模型预测补偿方法如下:
Figure FDA0002496201530000048
CN202010418993.XA 2020-05-18 2020-05-18 一种网络系统中数据丢包的提升模型预测补偿方法 Pending CN111585822A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010418993.XA CN111585822A (zh) 2020-05-18 2020-05-18 一种网络系统中数据丢包的提升模型预测补偿方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010418993.XA CN111585822A (zh) 2020-05-18 2020-05-18 一种网络系统中数据丢包的提升模型预测补偿方法

Publications (1)

Publication Number Publication Date
CN111585822A true CN111585822A (zh) 2020-08-25

Family

ID=72115540

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010418993.XA Pending CN111585822A (zh) 2020-05-18 2020-05-18 一种网络系统中数据丢包的提升模型预测补偿方法

Country Status (1)

Country Link
CN (1) CN111585822A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113985737A (zh) * 2021-10-27 2022-01-28 湘潭大学 具有时间延迟和数据包丢失的网络化控制系统的研究
CN115576206A (zh) * 2022-11-07 2023-01-06 中国矿业大学 一种多速率重介分选信息物理系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104639293A (zh) * 2015-01-30 2015-05-20 北方工业大学 网络化控制系统丢包补偿装置及方法
CN110083064A (zh) * 2019-04-29 2019-08-02 辽宁石油化工大学 一种基于非策略q-学习的网络最优跟踪控制方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104639293A (zh) * 2015-01-30 2015-05-20 北方工业大学 网络化控制系统丢包补偿装置及方法
CN110083064A (zh) * 2019-04-29 2019-08-02 辽宁石油化工大学 一种基于非策略q-学习的网络最优跟踪控制方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
吕云凯: "数据驱动的鲁棒学习控制系统综合与分析", 《中国优秀博硕士学位论文全文数据库(硕士)基础科学辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113985737A (zh) * 2021-10-27 2022-01-28 湘潭大学 具有时间延迟和数据包丢失的网络化控制系统的研究
CN115576206A (zh) * 2022-11-07 2023-01-06 中国矿业大学 一种多速率重介分选信息物理系统

Similar Documents

Publication Publication Date Title
CN110221542B (zh) 一种二阶非线性多智能体系统固定时间协同跟踪控制方法
Wang et al. Consensus tracking for nonlinear multi-agent systems with unknown disturbance by using model free adaptive iterative learning control
Chu et al. H∞ observer-based event-triggered sliding mode control for a class of discrete-time nonlinear networked systems with quantizations
Wei et al. Leader‐following consensus of second‐order multi‐agent systems with arbitrarily appointed‐time prescribed performance
CN110209148B (zh) 一种基于描述系统观测器的网络化系统的故障估计方法
Yoo et al. Guaranteed performance design for distributed bounded containment control of networked uncertain underactuated surface vessels
Liu et al. Adaptive fixed-time hierarchical sliding mode control for switched under-actuated systems with dead-zone constraints via event-triggered strategy
CN110119087B (zh) 一种有向通信下二阶多智能体系统固定时一致性跟踪方法
Sathishkumar et al. Mixed H∞ and passivity-based resilient controller for nonhomogeneous Markov jump systems
Chu et al. H∞ non-fragile observer-based dynamic event-triggered sliding mode control for nonlinear networked systems with sensor saturation and dead-zone input
CN111585822A (zh) 一种网络系统中数据丢包的提升模型预测补偿方法
Wang et al. Finite-time adaptive neural control for nonstrict-feedback stochastic nonlinear systems with input delay and output constraints
Wang et al. Fixed-time containment control for nonlinear multi-agent systems with external disturbances
CN113406886B (zh) 单连杆机械臂的模糊自适应控制方法、系统及存储介质
Wang et al. Event-triggered decentralized output-feedback control for interconnected nonlinear systems with input quantization
CN111382499A (zh) 一种化学循环反应器系统故障与扰动的联合估计方法
Chi et al. Compensation‐based data‐driven ILC with input and output package dropouts
Ren et al. Event-triggered non-fragile H∞ fault detection for discrete time-delayed nonlinear systems with channel fadings
CN113589689A (zh) 一种基于多参数自适应神经网络的滑模控制器设计方法
Wang et al. Command filter based globally stable adaptive neural control for cooperative path following of multiple underactuated autonomous underwater vehicles with partial knowledge of the reference speed
Yang et al. Adaptive predictive control of periodic non‐linear auto‐regressive moving average systems using nearest‐neighbour compensation
Si et al. Finite-time adaptive neural control for uncertain nonlinear time-delay systems with actuator delay and full-state constraints
Ao et al. Finite-time leaderless consensus control of a group of Euler-Lagrangian systems with backlash nonlinearities
Esmaeili et al. Data-driven observer-based model-free adaptive discrete-time terminal sliding mode control of rigid robot manipulators
CN113433825B (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
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20200825

WD01 Invention patent application deemed withdrawn after publication