CN113341707A - 一种多时间尺度多电机系统的事件触发协同控制方法 - Google Patents

一种多时间尺度多电机系统的事件触发协同控制方法 Download PDF

Info

Publication number
CN113341707A
CN113341707A CN202110521127.8A CN202110521127A CN113341707A CN 113341707 A CN113341707 A CN 113341707A CN 202110521127 A CN202110521127 A CN 202110521127A CN 113341707 A CN113341707 A CN 113341707A
Authority
CN
China
Prior art keywords
motor
event
triggered
time
sampling
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.)
Granted
Application number
CN202110521127.8A
Other languages
English (en)
Other versions
CN113341707B (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.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
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 China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN202110521127.8A priority Critical patent/CN113341707B/zh
Publication of CN113341707A publication Critical patent/CN113341707A/zh
Application granted granted Critical
Publication of CN113341707B publication Critical patent/CN113341707B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive 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/042Adaptive 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)
  • Control Of Multiple Motors (AREA)

Abstract

本发明涉及一种多时间尺度多电机系统的事件触发协同控制方法,属于多电机系统控制技术领域。其步骤包括:步骤一:给出关于多电机奇异摄动系统网络结构的假设,描述每个电机系统的具有多时间尺度的奇异摄动系统模型;步骤二:设计一个基于时间采样的分布式事件触发控制策略,提出控制输入的分布式协议和事件触发函数以及最终多电机趋同达到一致性的数学定义;步骤三:根据以上假设和问题陈述,详细推导为实现控制目标需要满足的不等式条件;步骤四:计算控制器增益。本发明平衡了系统的控制性能与通信资源尽可能节约之间的矛盾关系。

Description

一种多时间尺度多电机系统的事件触发协同控制方法
技术领域
本发明涉及一种多时间尺度多电机系统的事件触发协同控制方法,属于多电机系统控制技术领域。
背景技术
在现代大型生产生活中,关键输送设备直接影响了企业安全生产。在煤炭生产线中,带式输送机是关键的运输设备,其具有成本低、运输量大、结构简单等优点。带式输送机作业效率高被广泛应用在矿山开采、金属冶炼、发电厂等各种行业。带式输送机的稳定运行是工业安全生产的重要前提,其中带式输送机运行多依赖于多电机驱动。在现代工业的生产环节中,常常会利用多台电机的协同控制以实现产品的传输和卷曲。例如在冶金、建材、纺织、造纸等行业中,为了避免货物发生堆积、防止物件被拉断,多个电机之间必须保持协调运行,进而才能提高产品的生产质量和运输效率。因此研究多电机转速、转矩等输出参数的一致性控制具有重要的价值,对生产效率的提高具有深远的意义。
值得注意的是,在实际应用中很多多电机系统的个体具有明显的多时间尺度特性,例如无人机系统、电力系统等。随着现代工业水平的不断提高和自然科学关系网的错综繁杂,系统的结构与规模也日趋精细、庞大,从而对系统的控制质量也提出了更高的要求,因此研究多时间尺度系统的控制算法具有重大意义。传统多电机系统的设计方法不能直接应用,会带来病态数值问题,使得系统的分析与设计变得更为复杂。另外,由于多电机之间的通信一般都是在网络环境下进行的,而实际的网络通信资源非常有限,限制了网络通信协议的达成。
发明内容
本发明的目的在于克服现有技术的不足,提供了一种平衡系统的控制性能尽可能良好与通信资源尽可能节约之间的矛盾关系的多时间尺度多电机系统的事件触发协同控制方法,以解决病态数值问题。
本发明的技术方案如下:
本发明针对带式输送机工程应用中经常遇到的具有多时间尺度的奇异摄动系统模型,设计了在有向生成树拓扑结构下的多电机系统的一致性控制算法。首先设计了基于采样和事件触发策略的多电机系统的控制器模型,然后利用矩阵变换、Lyapunov函数、LMI等方法推导了使系统误差达到稳定的充分条件。这样可以通过求解LMIs就可以得到使得系统稳定的控制器增益,从而完成控制器的设计。其次,为了平衡系统的控制性能尽可能良好与通信资源尽可能节约之间的矛盾关系,本发明提出了以预期传输速率为指标的控制器增益选择算法,使得所设计的控制器既能实现多电机系统的一致性控制目标,同时多电机之间的平均信息传输速率为预期值。
具体地,本发明一种多时间尺度多电机系统的事件触发协同控制方法,其步骤包括:
步骤一:给出关于多电机奇异摄动系统网络结构的假设,描述每个电机系统的具有多时间尺度的奇异摄动系统模型;
步骤二:设计基于时间采样的分布式事件触发控制策略,提出控制输入的事件触发控制协议和事件触发函数以及最终多电机趋同达到一致性的数学定义;
步骤三:根据以上假设和问题陈述,详细推导为实现控制目标需要满足的不等式条件;
步骤四:计算控制器增益。
进一步地,步骤一中,假设存在具有N个相同电机的奇异摄动系统网络。第i个电机的动力学由以下微分方程描述:
Figure BDA0003064013560000021
其中,i=1,2,...,N,
Figure BDA0003064013560000022
Figure BDA0003064013560000023
分别是第i个节点的状态向量和输入控制,
Figure BDA0003064013560000024
Figure BDA0003064013560000025
都是常数矩阵,矩阵Eμ具有以下结构:
Figure BDA0003064013560000026
其中0<μ<<1是奇异摄动参数,n1+n2=n。式(1)的初始状态给定为
Figure BDA0003064013560000027
为了减少电机之间的通信损耗,步骤二中,本发明设计了一个基于时间采样的分布式事件触发控制策略。每个电机i通过采样器在固定时刻(t=kh)对自身状态信息进行采样,(
Figure BDA0003064013560000028
h>0为采样周期),并通过零阶保持器保证在下一个采样时刻(k+1)h到来之前采样数据不变。如果电机i在t=kh的采样数据满足事件触发条件,那么电机i将向其邻居广播自身当前的状态信息,同时其邻居的存储器也立刻得到更新。
事件触发控制协议设计表述为:
Figure BDA0003064013560000031
其中K是后面需要设计的反馈增益矩阵,邻接矩阵
Figure BDA0003064013560000032
表示通信网络的结构。
Figure BDA0003064013560000033
代表的是电机i的事件触发时刻,
Figure BDA0003064013560000034
表示事件触发时刻电机i的状态信息。考虑到电机i在
Figure BDA0003064013560000035
内,邻居电机也可能会触发事件,所以用
Figure BDA0003064013560000036
表示邻居电机j最新广播的状态信息。
每个电机i的事件触发函数由下式给出:
Figure BDA0003064013560000037
其中Φ>0是加权矩阵,σi>0是阈值参数,且:
Figure BDA0003064013560000038
Figure BDA0003064013560000039
所以,电机i下一次触发时刻
Figure BDA00030640135600000310
为:
Figure BDA00030640135600000311
可以直观地看出,电机i的下一个广播时刻
Figure BDA00030640135600000312
不仅与其自身的采样数据有关,而且还与其邻居的采样数据有关。当每次采样时刻到来时,电机i首先根据采样数据计算出自身的状态误差,再判断当前的状态误差是否超过了给定的阀值,即判断事件触发函数fi(kh)是否大于0。如果发现fi(kh)>0,那么就记电机i触发了一次事件。当一次事件触发后,电机i就立刻将自身当前的采样数据存储下来,并迅速更新自身控制器的输入ui和阈值参数,同时把自身当前的状态信息广播给它的出邻居电机。当电机j接收到来自其入邻居电机i的最新触发时刻的状态信息时,该电机j也立即更新它的控制器输入和阈值参数。
如果fi(kh)<0时,则意味着采样数据的误差仍然在允许的范围内,电机i就不需要向邻居发送自身当前的状态信息,也不需要改变自身控制器的输入信号,直到下一次采样时刻到来后再进行判断。从这种意义上来讲,事件触发传输策略可以显著降低采样数据传输和控制信号更新的频率,这意味着可以节省大量的通信资源和计算资源。
此外,事件驱动控制的一个重要问题是排除Zeno现象的发生,由于本发明采用了固定周期的采样策略,所以两次事件触发的间隔至少大于采样周期h,而h>0,因此Zeno现象是显然不存在的。
将第k次采样时的测量误差定义为:
Figure BDA0003064013560000041
由于事件触发间隔
Figure BDA0003064013560000042
包含多个采样区间,所以将
Figure BDA0003064013560000043
分成
Figure BDA0003064013560000044
个采样区间[kh,(k+1)h),得到在t∈[kh,(k+1)h)上的闭环系统:
Figure BDA0003064013560000045
Figure BDA0003064013560000046
可得到如下多电机闭环系统:
Figure BDA0003064013560000047
其中L是有向图的拉普拉斯矩阵。
zi(t)=x1(t)-xi+1(t),εi(kh)=e1(kh)-ei+1(kh),
Figure BDA0003064013560000048
Figure BDA0003064013560000049
不难发现
Figure BDA00030640135600000410
Figure BDA00030640135600000411
其中T1=[1,-IN-1],T2=[0,-IN-1]T。代入多电机闭环系统方程中可得:
Figure BDA00030640135600000412
其中
Figure BDA0003064013560000051
再构造时滞函数τ(t)=t-kh,kh≤t<(k+1)h,显然τ(t)是分段线性的,当t≠kh时其导数
Figure BDA0003064013560000052
当t=kh时
Figure BDA0003064013560000053
是不连续的。由于kh≤t<(k+1)h,所以0≤τ(t)<h。因此可以得到:
Figure BDA0003064013560000054
其中初始状态给定为z(θ)=φ(θ),-h≤θ≤0,
Figure BDA0003064013560000055
很明显,在事件触发控制策略与事件触发控制协议下,当且仅当系统渐近稳定时,多电机系统可以渐近达成一致性。因此,多电机系统的一致性可以转化为系统的稳定性问题。
如果存在标量μ0和适当的事件触发控制协议,使得
Figure BDA0003064013560000056
Figure BDA0003064013560000057
对任意μ∈(0,μ0]都成立,则称电机系统达到了一致性。
步骤三中,用于实现控制目标需要满足的不等式条件如下:
在事件触发策略下,对任意
Figure BDA0003064013560000058
有下面不等式成立:
Figure BDA0003064013560000059
其中
Figure BDA00030640135600000510
Figure BDA00030640135600000511
的定义域为
Figure BDA00030640135600000512
给定采样周期h>0,事件触发参数Λ=diag{σ12,…,σN},以及奇异摄动参数上界μ0,如果有向图
Figure BDA00030640135600000515
具有一棵生成树,若存在矩阵Pi(i=1,2,3,4,5),Xi(i=1,2,3,4,5),且Pi=Pi T(i=1,2,3,4),
Figure BDA00030640135600000513
存在对称正定矩阵Q>0,R>0,Y>0,使得下列不等式成立:
X1>0
Figure BDA00030640135600000514
Figure BDA0003064013560000061
Figure BDA0003064013560000062
Figure BDA0003064013560000063
其中sym(·)为对称函数,如sym(A)=A+AT,*代表大矩阵中相应的对称矩阵块,
Figure BDA0003064013560000064
Figure BDA0003064013560000065
λ=λmax(T1 TT1)
Figure BDA0003064013560000066
则对任意μ∈(0,μ0],在事件触发策略下闭环系统是全局渐近稳定的。
分布式控制器存在的充分条件以及控制器增益K的显示表达式为:
给定h>0,ρ>0,δ>0,Λ=diag{σ12,…,σN},以及奇异摄动参数上界μ0,若存在矩阵Pi(i=1,2,3,4,5)满足Pi=Pi T(i=1,2,3,4),存在对称正定矩阵
Figure BDA0003064013560000071
使得下列线性矩阵不等式LMIs成立:
Figure BDA0003064013560000072
Figure BDA0003064013560000073
其中
Figure BDA0003064013560000074
Figure BDA0003064013560000075
λ=λmax(T1 TT1),
Figure BDA0003064013560000076
则对任意μ∈(0,μ0],多电机系统在事件触发策略下的渐近一致性可以实现,并且控制器增益
Figure BDA0003064013560000081
然而在实际应用中,由于电机i的事件触发传输策略高度依赖于其邻居的采样数据和自身的采样数据,因此参数σi的选择不仅会影响电机i的传输次数,而且还会影响其邻居的传输次数。此外,分布式多电机系统的一致性控制性能与各电机的采样数据传输次数密切相关,事件触发的频率越高,那么系统的控制性能就能得到更好的保证,但相反通信资源的损耗和浪费也会大大增加。因此,就需要在系统性能和节约资源两方面作以权衡,但同时获得事件触发传输策略的一致控制器K和所有参数σi和Φ是一项困难的任务。
为此本发明设计了以下算法来获得较好的控制器增益,见步骤四。首先引入一个评价系统通信传输快慢的指标,在时间段[0,Th)上系统网络的平均传输速率定义为:
Figure BDA0003064013560000082
其中
Figure BDA0003064013560000083
代表某一采样时刻所有电机的采样数据量,
Figure BDA0003064013560000084
是数据传输信号函数,
Figure BDA0003064013560000085
表示在时间段[0,Th)内整个系统中所有电机触发事件向外传输的数据量。所以平均传输速率S的物理意义为:有限时间内系统总传输数据量与总采样数据量之比。
由定义可知:0<S<1,并且S越小表明系统越节能。极端地,当S=1时系统从事件触发控制转变成常规采样控制系统。那么对于一个事件触发控制系统,S取值为多少时,系统才叫好呢,或是系统的控制性能与低耗性能才会达到平衡呢。为此本发明引入预期传输速率S*,S*是人为预先规定的,并认为S与S*越接近则系统的性能越良好。定义实际速率与预期速率的误差为ε=|S-S*|,从而评价系统好坏的标准变成了系统传输率是否和预期情况相一致。故多电机系统一致性问题转变成:以满足控制器存在的充分条件为约束,以使ε达到最小值为目标的参数优化问题。为了解决这一问题,本发明提出了以下算法。
在确保采用事件触发策略和分布式控制协议的多电机系统能达到一致性的前提下,控制器增益的优化选择可以通过以下步骤来实现。
算法:协同设计算法
第一步:选择选择参数σi的增长步长κi,采样周期h,采样次数T,预期传输速率S*,奇异摄动参数上界μ0,以及较大的初始误差ε0和事件触发参数Λ=diag{σ12,…,σN}的初始值Λ0=0。其中令σi的初值为0是为了让事件更容易触发,以及方便下一步σi从小到大以步长κi不断增加。初始ε0选择较大是由于算法的目标是寻找尽量小的ε,这样方便ε下一步的更新;
第二步:将以上参数带入控制器存在的充分条件的LMIs,求解能使系统稳定的σi的最大值σim
第三步:
Figure BDA0003064013560000092
第四步:输出最优的参数ε000,K0
有益效果
发明研究了具有多时间尺度的多电机系统基于事件触发策略的分布式控制协议,并通过Lyapunov稳定性理论严格推导了使多电机系统达到一致性的充分条件。其次,为了平衡系统的控制性能和系统的节能性,提出了以为达到预定传输速率为目标的控制器增益的选择算法;
本发明针对带式输送机工程应用中经常遇到的具有多时间尺度的奇异摄动系统模型,设计了在有向生成树拓扑结构下的多电机系统的一致性控制算法。首先设计了基于采样和事件触发策略的多电机系统的控制器模型,然后利用矩阵变换、Lyapunov函数、LMI等方法推导了使系统误差达到稳定的充分条件。这样可以通过求解LMIs就可以得到使得系统稳定的控制器增益,从而完成控制器的设计;
为了平衡系统的控制性能尽可能良好与通信资源尽可能节约之间的矛盾关系,本发明提出了以预期传输速率为指标的控制器增益选择算法,使得所设计的控制器既能实现多电机系统的一致性控制目标,同时多电机之间的平均信息传输速率为预期值。
附图说明
图1是本发明一种多时间尺度多电机系统的事件触发协同控制方法的流程图;
图2是多电机网络拓扑结构图;
图3是各电机的角速度(转速)误差轨迹;
图4是各电机的角加速度误差轨迹;
图5是各电机的转矩误差轨迹;
图6是各电机的事件触发时刻。
具体实施方式
下面将结合附图对本发明技术方案的实施例进行详细地描述。
本发明一种多时间尺度多电机系统的事件触发协同控制方法,如图1所示,其步骤包括:
步骤一:给出关于多电机奇异摄动系统网络结构,描述每个电机系统的具有多时间尺度的奇异摄动系统模型;
步骤二:提出控制输入的事件触发控制协议和事件触发函数以及最终多电机趋同达到一致性的数学定义;
步骤三:根据以上假设和问题陈述,详细推导为实现控制目标需要满足的不等式条件;
步骤四:计算控制器增益。
实施例:
1、给出关于多电机奇异摄动系统网络结构,描述每个电机系统的具有多时间尺度的奇异摄动系统模型:
假设存在具有N个相同电机的奇异摄动系统网络。第i个电机的动力学由以下微分方程描述:
Figure BDA0003064013560000111
其中,i=1,2,...,N,
Figure BDA0003064013560000112
Figure BDA0003064013560000113
分别是第i个节点的状态向量和输入控制,
Figure BDA0003064013560000114
Figure BDA0003064013560000115
都是常数矩阵,矩阵Eμ具有以下结构:
Figure BDA0003064013560000116
其中0<μ<<1是奇异摄动参数,n1+n2=n。式(1)的初始状态给定为
Figure BDA0003064013560000117
2、提出控制输入的事件触发控制协议和事件触发函数以及最终多电机趋同达到一致性的数学定义:
为了减少电机之间的通信损耗,本发明设计了一个基于时间采样的分布式事件触发控制策略。每个电机i通过采样器在固定时刻(t=kh)对自身状态信息进行采样,(
Figure BDA0003064013560000118
h>0为采样周期),并通过零阶保持器保证在下一个采样时刻(k+1)h到来之前采样数据不变。如果电机i在t=kh的采样数据满足事件触发条件,那么电机i将向其邻居广播自身当前的状态信息,同时其邻居的存储器也立刻得到更新。事件触发控制协议设计表述为:
Figure BDA0003064013560000119
其中K是后面需要设计的反馈增益矩阵,邻接矩阵
Figure BDA00030640135600001110
表示通信网络的结构。
Figure BDA00030640135600001111
代表的是电机i的事件触发时刻,
Figure BDA00030640135600001112
表示事件触发时刻电机i的状态信息。考虑到电机i在
Figure BDA0003064013560000121
内,邻居电机也可能会触发事件,所以用
Figure BDA0003064013560000122
表示邻居电机j最新广播的状态信息。
每个电机i的事件触发函数由下式给出:
Figure BDA0003064013560000123
其中Φ>0是加权矩阵,σi>0是阈值参数,且:
Figure BDA0003064013560000124
Figure BDA0003064013560000125
所以,电机i下一次触发时刻
Figure BDA0003064013560000126
为:
Figure BDA0003064013560000127
可以直观地看出,电机i的下一个广播时刻
Figure BDA0003064013560000128
不仅与其自身的采样数据有关,而且还与其邻居的采样数据有关。当每次采样时刻到来时,电机i首先根据采样数据计算出自身的状态误差,再判断当前的状态误差是否超过了给定的阀值,即判断事件触发函数fi(kh)是否大于0。如果发现fi(kh)>0,那么就记电机i触发了一次事件。当一次事件触发后,电机i就立刻将自身当前的采样数据存储下来,并迅速更新自身控制器的输入ui和阈值参数,同时把自身当前的状态信息广播给它的出邻居电机。当电机j接收到来自其入邻居电机i的最新触发时刻的状态信息时,该电机j也立即更新它的控制器输入和阈值参数。
如果fi(kh)<0时,则意味着采样数据的误差仍然在允许的范围内,电机i就不需要向邻居发送自身当前的状态信息,也不需要改变自身控制器的输入信号,直到下一次采样时刻到来后再进行判断。从这种意义上来讲,事件触发传输策略可以显著降低采样数据传输和控制信号更新的频率,这意味着可以节省大量的通信资源和计算资源。
此外,事件驱动控制的一个重要问题是排除Zeno现象的发生,由于本发明采用了固定周期的采样策略,所以两次事件触发的间隔至少大于采样周期h,而h>0,因此Zeno现象是显然不存在的。
将第k次采样时的测量误差定义为:
Figure BDA0003064013560000131
其中,kh代表的是电机i的采样时刻,xi(kh)代表的是采样时刻电机i的状态信息,
Figure BDA0003064013560000132
代表的是电机i的事件触发时刻,
Figure BDA0003064013560000133
表示事件触发时刻电机i的状态信息。
由于事件触发间隔
Figure BDA0003064013560000134
包含多个采样区间,所以将
Figure BDA0003064013560000135
分成
Figure BDA0003064013560000136
个采样区间[kh,(k+1)h),得到在t∈[kh,(k+1)h)上的闭环系统:
Figure BDA0003064013560000137
Figure BDA0003064013560000138
可得到如下多电机闭环系统:
Figure BDA0003064013560000139
其中L是有向图的拉普拉斯矩阵。
令zi(t)=x1(t)-xi+1(t),εi(kh)=e1(kh)-ei+1(kh),
Figure BDA00030640135600001310
Figure BDA00030640135600001311
不难发现
Figure BDA00030640135600001312
Figure BDA00030640135600001313
其中T1=[1,-IN-1],T2=[0,-IN-1]T。代入多电机闭环系统方程中可得:
Figure BDA00030640135600001314
其中
Figure BDA00030640135600001315
再构造时滞函数τ(t)=t-kh,kh≤t<(k+1)h,显然τ(t)是分段线性的,当t≠kh时其导数
Figure BDA00030640135600001316
当t=kh时
Figure BDA00030640135600001317
是不连续的。由于kh≤t<(k+1)h,所以0≤τ(t)<h。因此可以得到:
Figure BDA00030640135600001318
Figure BDA0003064013560000141
其中初始状态给定为z(θ)=φ(θ),-h≤θ≤0,
Figure BDA0003064013560000142
很明显,在事件触发控制策略与事件触发控制协议下,当且仅当系统渐近稳定时,多电机系统可以渐近达成一致性。因此,多电机系统的一致性可以转化为系统的稳定性问题。
如果存在标量μ0和适当的事件触发控制协议,使得
Figure BDA0003064013560000143
对任意μ∈(0,μ0]都成立,则称电机系统达到了一致性。
3、根据以上假设和问题陈述,推导为实现控制目标需要满足的不等式条件:
用于实现控制目标需要满足的不等式条件如下:
在事件触发策略下,对任意
Figure BDA0003064013560000144
有下面不等式成立:
Figure BDA0003064013560000145
其中
Figure BDA0003064013560000146
Figure BDA0003064013560000147
的定义域为
Figure BDA0003064013560000148
给定采样周期h>0,Λ=diag{σ12,…,σN},以及奇异摄动参数上界μ0,如果有向图
Figure BDA00030640135600001412
具有一棵生成树,若存在矩阵Pi(i=1,2,3,4,5),Xi(i=1,2,3,4,5),且Pi=Pi T(i=1,2,3,4),
Figure BDA0003064013560000149
存在对称正定矩阵Q>0,R>0,Y>0,使得下列不等式成立:
X1>0 (12)
Figure BDA00030640135600001410
Figure BDA00030640135600001411
Figure BDA0003064013560000151
Figure BDA0003064013560000152
其中sym(·)为对称函数,如sym(A)=A+AT,*代表大矩阵中相应的对称矩阵块,
Figure BDA0003064013560000153
Figure BDA0003064013560000154
λ=λmax(T1 TT1)
Figure BDA0003064013560000155
则对任意μ∈(0,μ0],在事件触发策略下闭环系统是全局渐近稳定的。
4、计算控制器增益:
分布式控制器存在的充分条件以及控制器增益K的显示表达式为:
给定采样周期h>0,ρ>0,δ>0,Λ=diag{σ12,…,σN},以及奇异摄动参数上界μ0,若存在矩阵Pi(i=1,2,3,4,5)满足Pi=Pi T(i=1,2,3,4),存在对称正定矩阵
Figure BDA0003064013560000161
使得下列不等式成立:
Figure BDA0003064013560000162
Figure BDA0003064013560000163
其中
Figure BDA0003064013560000164
Figure BDA0003064013560000165
λ=λmax(T1 TT1),
Figure BDA0003064013560000166
则对任意μ∈(0,μ0],多电机系统在事件触发策略下的渐近一致性可以实现,并且控制器增益的表达式为:
Figure BDA0003064013560000167
然而在实际应用中,由于电机i的事件触发传输策略高度依赖于其邻居的采样数据和自身的采样数据,因此参数σi的选择不仅会影响电机i的传输次数,而且还会影响其邻居的传输次数。此外,分布式多电机系统的一致性控制性能与各电机的采样数据传输次数密切相关,事件触发的频率越高,那么系统的控制性能就能得到更好的保证,但相反通信资源的损耗和浪费也会大大增加。因此,就需要在系统性能和节约资源两方面作以权衡,但同时获得事件触发传输策略的一致控制器K和所有参数σi和Φ是一项困难的任务。
为此本发明设计了以下算法来获得较好的控制器增益。首先引入一个评价系统通信传输快慢的指标,在时间段[0,Th)上系统网络的平均传输速率定义为:
Figure BDA0003064013560000171
其中
Figure BDA0003064013560000172
代表某一采样时刻所有电机的采样数据量,
Figure BDA0003064013560000173
是数据传输信号函数,
Figure BDA0003064013560000174
表示在时间段[0,Th)内整个系统中所有电机触发事件向外传输的数据量。所以平均传输速率S的物理意义为:有限时间内系统总传输数据量与总采样数据量之比。
由定义可知:0<S<1,并且S越小表明系统越节能。极端地,当S=1时系统从事件触发控制转变成常规采样控制系统。那么对于一个事件触发控制系统,S取值为多少时,系统才叫好呢,或是系统的控制性能与低耗性能才会达到平衡呢。为此本发明引入预期传输速率S*,S*是人为预先规定的,并认为S与S*越接近则系统的性能越良好。定义实际速率与预期速率的误差为ε=|S-S*|,从而评价系统好坏的标准变成了系统传输率是否和预期情况相一致。故多电机系统一致性问题转变成:以满足控制器存在的充分条件为约束,以使ε达到最小值为目标的参数优化问题。为了解决这一问题,本发明提出了以下算法。
在确保采用事件触发策略和分布式控制协议的多电机系统能达到一致性的前提下,控制器增益的优化选择可以通过以下步骤来实现。
算法:协同设计算法
第一步:选择选择参数σi的增长步长κi,采样周期h,采样次数T,预期传输速率S*,奇异摄动参数上界μ0,以及较大的初始误差ε0和事件触发参数Λ=diag{σ12,…,σN}的初始值Λ0=0。其中令σi的初值为0是为了让事件更容易触发,以及方便下一步σi从小到大以步长κi不断增加。初始ε0选择较大是由于算法的目标是寻找尽量小的ε,这样方便ε下一步的更新;
第二步:将以上参数带入控制器存在的充分条件的LMIs,求解能使系统稳定的σi的最大值σim
第三步:
Figure BDA0003064013560000182
第四步:输出最优的参数ε000,K0
5、仿真算例分析
利用4台电机组成的多电机系统作为仿真对象,各智能个体的电机电阻R、电感L、等效粘滞摩擦系数beq、等效电机转动惯量Jeq、电机转矩系数km、反电动势常数ke、齿轮箱传动比kt等参数设置情况如表1所示。
表1各电机参数设置表
Figure BDA0003064013560000191
各电机的转速、角加速度、转矩初始值为x1(0)=[27,12,30,15]T,x2(0)=[2.1,1,0.3,3.8]T,x3(0)=[2.3,4.6,3,5.2]T。4台电机的网络拓扑结构如图2所示,可以看出它是具有生成树的有向图。不失一般性,假设所有的权重均为1。
则以上拓扑相对应的邻接矩阵和Laplacian矩阵为:
Figure BDA0003064013560000192
接下来,选择采样周期h=0.03s,事件触发函数的阈值参数σ1=0.021,σ2=0.073,σ3=0.064,σ4=0.039,通过求解LMIS,可以得到以下可行解:
触发函数的常量矩阵φ=[0.0015 0.0017 0.0018;0.0017 0.0067 0.002;0.0021 0.002 0.035],控制器增益K=[0.842 2.135 0.786]。从而保证了在给定初始条件下的4台电机经过有限的时间后最终都能实现转速、转矩等参数的一致性输出。通过Matlab仿真可知4台电机的转速,角加速度,以及转矩误差轨迹分别如图3至图5所示,同时每台电机的事件触发时间序列如图6所示。
从图3至图5可以看出4台电机在不同初始条件下,经过控制器的作用最终能实现转速、角加速度等状态的一致性同步,这说明了所提出的设计方法是有效的,能够实现预期的控制目标。这在实际生产应用中具有重大的意义,比如在生产实际中要考虑电气设备的安全运行以及保证产品的生产质量,这时就需要让多台电机共同协作,统一输出来完成某项作业,因此多电机的一致性控制就起到了关键作用。另外从图6可知,在采用事件触发传输策略后,网络通信的次数明显减少。多台电机不用在每个采样时刻都进行互相通信,而只在事件触发时刻间歇地进行数据传输。而且通过程序计算可得出每个电机的事件触发频率如表2所示。
表2多电机事件触发次数统计表
Figure BDA0003064013560000201
根据上表可知,每个电机在[0,30)s内的触发次数分别为158,44,65,74次,也就是说在总体4000次的采样次数上只需要在其中的341次进行电机之间的信息传输,就能达到预期的控制目标。这无疑是很大程度上减少了网络带宽负担和通信资源的浪费,从而说明了事件触发策略的经济性和实用性。
以上所述的具体实施方式,对本发明的目的、技术方案和有益效果进行了进一步地详细说明,所应理解的是,以上所述仅为本发明的具体实施方式而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种多时间尺度多电机系统的事件触发协同控制方法,其特征在于,其步骤包括:
步骤一:给出关于多电机奇异摄动系统网络结构的假设,描述每个电机系统的具有多时间尺度的奇异摄动系统模型;
步骤二:设计基于时间采样的分布式事件触发控制策略,提出控制输入的事件触发控制协议和事件触发函数以及最终多电机趋同达到一致性的数学定义;
步骤三:根据以上假设和问题陈述,详细推导为实现控制目标需要满足的不等式条件;
步骤四:计算控制器增益。
2.根据权利要求1所述的事件触发协同控制方法,其特征在于,步骤一中,假设存在具有N个相同电机的奇异摄动系统网络,第i个电机的动力学由以下微分方程描述:
Figure FDA0003064013550000011
其中,i=1,2,...,N,
Figure FDA0003064013550000012
Figure FDA0003064013550000013
分别是第i个节点的状态向量和输入控制,
Figure FDA0003064013550000014
Figure FDA0003064013550000015
都是常数矩阵,矩阵Eμ具有以下结构:
Figure FDA0003064013550000016
其中0<μ<<1是奇异摄动参数,n1+n2=n,式中的初始状态给定为
Figure FDA0003064013550000017
3.根据权利要求1所述的事件触发协同控制方法,其特征在于,步骤二中,所述基于时间采样的分布式事件触发控制策略指:每个电机i通过采样器在固定时刻t=kh对自身状态信息进行采样,
Figure FDA0003064013550000018
h>0为采样周期,并通过零阶保持器保证在下一个采样时刻(k+1)h到来之前采样数据不变;如果电机i在t=kh的采样数据满足事件触发条件,那么电机i将向其邻居广播自身当前的状态信息,同时其邻居的存储器也立刻得到更新。
4.根据权利要求3所述的事件触发协同控制方法,其特征在于,所述事件触发条件的事件触发控制协议设计表述为:
Figure FDA0003064013550000021
其中K是后面需要设计的反馈增益矩阵,邻接矩阵
Figure FDA0003064013550000022
表示通信网络的结构;
Figure FDA0003064013550000023
代表的是电机i的事件触发时刻,
Figure FDA0003064013550000024
表示事件触发时刻电机i的状态信息;考虑到电机i在
Figure FDA0003064013550000025
内,邻居电机也可能会触发事件,所以用
Figure FDA0003064013550000026
表示邻居电机j最新广播的状态信息;
每个电机i的事件触发函数由下式给出:
Figure FDA0003064013550000027
其中Φ>0是加权矩阵,σi>0是阈值参数,且:
Figure FDA0003064013550000028
Figure FDA0003064013550000029
所以,电机i下一次触发时刻
Figure FDA00030640135500000210
为:
Figure FDA00030640135500000211
其中,电机i的下一个广播时刻
Figure FDA00030640135500000212
不仅与其自身的采样数据有关,而且还与其邻居的采样数据有关;当每次采样时刻到来时,电机i首先根据采样数据计算出自身的状态误差,再判断当前的状态误差是否超过了给定的阀值,即判断事件触发函数fi(kh)是否大于0;如果发现fi(kh)>0,那么就记电机i触发了一次事件;当一次事件触发后,电机i就立刻将自身当前的采样数据存储下来,并迅速更新自身控制器的输入ui和阈值参数,同时把自身当前的状态信息广播给它的出邻居电机;当电机j接收到来自其入邻居电机i的最新触发时刻的状态信息时,该电机j也立即更新它的控制器输入和阈值参数;
如果fi(kh)<0时,则意味着采样数据的误差仍然在允许的范围内,电机i就不需要向邻居发送自身当前的状态信息,也不需要改变自身控制器的输入信号,直到下一次采样时刻到来后再进行判断;
将第k次采样时的测量误差定义为:
Figure FDA0003064013550000031
其中,kh代表的是电机i的采样时刻,xi(kh)代表的是采样时刻电机i的状态信息,
Figure FDA0003064013550000032
代表的是电机i的事件触发时刻,
Figure FDA0003064013550000033
表示事件触发时刻电机i的状态信息。
5.根据权利要求4所述的事件触发协同控制方法,其特征在于,由于事件触发间隔
Figure FDA0003064013550000034
包含多个采样区间,所以将
Figure FDA0003064013550000035
分成
Figure FDA0003064013550000036
个采样区间[kh,(k+1)h),得到在t∈[kh,(k+1)h)上的闭环系统:
Figure FDA0003064013550000037
Figure FDA0003064013550000038
可得到如下多电机闭环系统:
Figure FDA0003064013550000039
其中L是有向图的拉普拉斯矩阵;
令zi(t)=x1(t)-xi+1(t),εi(kh)=e1(kh)-ei+1(kh),
Figure FDA00030640135500000310
Figure FDA00030640135500000311
不难发现
Figure FDA00030640135500000312
Figure FDA00030640135500000313
其中T1=[1,-IN-1],T2=[0,-IN-1]T,代入多电机闭环系统方程中可得:
Figure FDA00030640135500000314
其中
Figure FDA00030640135500000315
再构造时滞函数τ(t)=t-kh,kh≤t<(k+1)h,显然τ(t)是分段线性的,当t≠kh时其导数
Figure FDA00030640135500000316
当t=kh时
Figure FDA00030640135500000317
是不连续的;由于kh≤t<(k+1)h,所以0≤τ(t)<h,因此可以得到:
Figure FDA00030640135500000318
其中初始状态给定为z(θ)=φ(θ),-h≤θ≤0,
Figure FDA0003064013550000041
如果存在标量μ0和适当的事件触发控制协议,使得
Figure FDA0003064013550000042
Figure FDA0003064013550000043
对任意μ∈(0,μ0]都成立,则称电机系统达到了一致性。
6.根据权利要求1所述的事件触发协同控制方法,其特征在于,步骤三中,用于实现控制目标需要满足的不等式条件如下:
在事件触发策略下,对任意
Figure FDA0003064013550000044
有下面不等式成立:
Figure FDA0003064013550000045
其中
Figure FDA0003064013550000046
Figure FDA0003064013550000047
的定义域为
Figure FDA0003064013550000048
给定采样周期h>0,事件触发参数Λ=diag{σ12,…,σN},以及奇异摄动参数上界μ0,如果有向图
Figure FDA0003064013550000049
具有一棵生成树,若存在矩阵Pi(i=1,2,3,4,5),Xi(i=1,2,3,4,5),且Pi=Pi T(i=1,2,3,4),
Figure FDA00030640135500000410
存在对称正定矩阵Q>0,R>0,Y>0,使得下列不等式成立:
X1>0
Figure FDA00030640135500000411
Figure FDA00030640135500000412
Figure FDA0003064013550000051
Figure FDA0003064013550000052
Figure FDA0003064013550000053
Figure FDA0003064013550000054
其中sym(·)为对称函数,如sym(A)=A+AT,*代表大矩阵中相应的对称矩阵块,
Figure FDA0003064013550000055
Figure FDA0003064013550000056
Figure FDA0003064013550000057
则对任意μ∈(0,μ0],在事件触发策略下闭环系统是全局渐近稳定的。
7.根据权利要求1所述的事件触发协同控制方法,其特征在于,分布式控制器存在的充分条件以及控制器增益K的显示表达式为:
给定采样周期h>0,ρ>0,δ>0,Λ=diag{σ1,σ2,…,σN},以及奇异摄动参数上界μ0,若存在矩阵Pi(i=1,2,3,4,5)满足Pi=Pi T(i=1,2,3,4),存在对称正定矩阵
Figure FDA0003064013550000061
使得下列不等式成立:
Figure FDA0003064013550000062
Figure FDA0003064013550000063
Figure FDA0003064013550000064
Figure FDA0003064013550000065
其中
Figure FDA0003064013550000066
Figure FDA0003064013550000067
λ=λmax(T1 TT1),
Figure FDA0003064013550000068
则对任意μ∈(0,μ0],多电机系统在事件触发策略下的渐近一致性可以实现,并且控制器增益的表达式为:
Figure FDA0003064013550000071
8.根据权利要求1所述的事件触发协同控制方法,其特征在于,步骤四中,所述控制器增益的算法如下:首先引入一个评价系统通信传输快慢的指标,在时间段[0,Th)上系统网络的平均传输速率定义为:
Figure FDA0003064013550000072
其中
Figure FDA0003064013550000073
代表某一采样时刻所有电机的采样数据量,
Figure FDA0003064013550000074
是数据传输信号函数,
Figure FDA0003064013550000075
表示在时间段[0,Th)内整个系统中所有电机触发事件向外传输的数据量,由定义可知:0<S<1,并且S越小表明系统越节能。
9.根据权利要求8所述的事件触发协同控制方法,其特征在于,当S=1时系统从事件触发控制转变成常规采样控制系统。
10.根据权利要求8所述的事件触发协同控制方法,其特征在于,所述控制器增益的优化选择通过协同设计算法来实现,其步骤如下:
第一步:选择选择参数σi的增长步长κi,采样周期h,采样次数T,预期传输速率S*,奇异摄动参数上界μ0,以及较大的初始误差ε0和事件触发参数Λ=diag{σ1,σ2,…,σN}的初始值Λ0=0,其中令σi的初值为0是为了让事件更容易触发,以及方便下一步σi从小到大以步长κi不断增加,初始ε0选择较大是由于算法的目标是寻找尽量小的ε,这样方便ε下一步的更新;
第二步:将以上参数带入控制器存在的充分条件的LMIs的公式(17)和公式(18),求解能使系统稳定的σi的最大值σim
第三步:
Figure FDA0003064013550000077
Figure FDA0003064013550000082
第四步:输出最优的参数ε0,Λ0,Φ0,K0
CN202110521127.8A 2021-05-13 2021-05-13 一种多时间尺度多电机系统的事件触发协同控制方法 Active CN113341707B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110521127.8A CN113341707B (zh) 2021-05-13 2021-05-13 一种多时间尺度多电机系统的事件触发协同控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110521127.8A CN113341707B (zh) 2021-05-13 2021-05-13 一种多时间尺度多电机系统的事件触发协同控制方法

Publications (2)

Publication Number Publication Date
CN113341707A true CN113341707A (zh) 2021-09-03
CN113341707B CN113341707B (zh) 2021-12-17

Family

ID=77469526

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110521127.8A Active CN113341707B (zh) 2021-05-13 2021-05-13 一种多时间尺度多电机系统的事件触发协同控制方法

Country Status (1)

Country Link
CN (1) CN113341707B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113885313A (zh) * 2021-09-29 2022-01-04 杭州电子科技大学 一种多智能车系统周期采样事件触发协同控制方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120062239A1 (en) * 2010-09-10 2012-03-15 Rostron Joseph R Directional fault sectionalizing system
CN102749845A (zh) * 2012-06-15 2012-10-24 华中科技大学 基于事件触发机制的电力系统状态反馈控制器构建方法
CN105334888A (zh) * 2015-11-16 2016-02-17 江苏科技大学 一种基于触发机制的网络串级温度控制系统及其方法
CN108448576A (zh) * 2018-04-09 2018-08-24 南京邮电大学 弹性事件触发机制下的多机电力系统广域阻尼控制方法
CN110011583A (zh) * 2019-04-12 2019-07-12 中国矿业大学 基于奇异摄动理论的永磁同步电机滑模控制系统及建模方法
CN110109351A (zh) * 2019-04-08 2019-08-09 广东工业大学 一种基于指定性能的多智能体一致性控制方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120062239A1 (en) * 2010-09-10 2012-03-15 Rostron Joseph R Directional fault sectionalizing system
CN102749845A (zh) * 2012-06-15 2012-10-24 华中科技大学 基于事件触发机制的电力系统状态反馈控制器构建方法
CN105334888A (zh) * 2015-11-16 2016-02-17 江苏科技大学 一种基于触发机制的网络串级温度控制系统及其方法
CN108448576A (zh) * 2018-04-09 2018-08-24 南京邮电大学 弹性事件触发机制下的多机电力系统广域阻尼控制方法
CN110109351A (zh) * 2019-04-08 2019-08-09 广东工业大学 一种基于指定性能的多智能体一致性控制方法
CN110011583A (zh) * 2019-04-12 2019-07-12 中国矿业大学 基于奇异摄动理论的永磁同步电机滑模控制系统及建模方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
MANISHA BHANDARI等: "Event-Triggered Composite Control of a Two Time Scale System", 《IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS—II: EXPRESS BRIEFS》 *
XIANG REN,等: "Event-triggered stability of singularly perturbed system with time delay", 《2020 39TH CHINESE CONTROL CONFERENCE (CCC)》 *
王燕舞等: "混杂奇异摄动系统的研究综述", 《控制与决策》 *
苏枝轩,等: "多时间尺度多智能体系统的事件触发协同控制", 《第31届中国过程控制会议(CPCC 2020)摘要集》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113885313A (zh) * 2021-09-29 2022-01-04 杭州电子科技大学 一种多智能车系统周期采样事件触发协同控制方法

Also Published As

Publication number Publication date
CN113341707B (zh) 2021-12-17

Similar Documents

Publication Publication Date Title
Cao et al. Active disturbance rejection‐based event‐triggered bipartite consensus control for nonaffine nonlinear multiagent systems
CN110705034B (zh) 一种基于事件触发的永磁同步电机位置跟踪控制方法
CN112327633A (zh) 具有时滞和干扰的领导跟随多智能体系统一致性的方法
CN113110059B (zh) 基于事件触发的单连杆机械臂系统实际跟踪的控制方法
CN107463095B (zh) 一种具有时变采样周期的输出反馈控制器设计方法
CN113341707B (zh) 一种多时间尺度多电机系统的事件触发协同控制方法
CN112379626A (zh) 具有外部干扰的多智能体系统分组领导跟随一致的方法
Liu Tracking control of multi-agent systems using a networked predictive PID tracking scheme
CN109818792B (zh) 一种基于二阶线性系统时变耦合复杂动态网络模型的控制器
CN114114904B (zh) 一种有限时间和固定时间分布式事件触发一致性方法
CN112936266A (zh) 网络化机械臂系统自适应一致性切换控制方法及系统
Liu et al. Event-based consensus control of multi-agent systems by L∞ theory
CN117762014A (zh) 一种采用分布式事件触发控制的有限时间协同控制方法
CN108614432A (zh) 一种基于粒子群算法的网络环境电机控制器设计算法
Ren et al. Event-Based Predefined-Time Fuzzy Formation Control for Nonlinear Multi-Agent Systems with Unknown Disturbances
CN114879492A (zh) 一种完全数据驱动的多智能体系统事件触发控制方法
CN109600077B (zh) 一种多永磁同步电机的分布式自适应协调控制方法
CN108958022B (zh) 非线性主-从时延系统的时滞相关有限时间同步控制方法
CN109039166B (zh) 一种永磁同步直线伺服系统速度环pi-ip控制参数自校正方法
CN106873374A (zh) 基于调度与控制器参数动态重构的大型电子束焊接装置网络化控制方法
CN113031644A (zh) 一种面向通信时滞的飞行器编队控制系统事件触发方法、装置及介质
CN113703315B (zh) 一种非线性混合阶多智能体机电系统事件触发协同控制方法
Beschi et al. A new two degree-of-freedom event-based PI control strategy
Hanchevici et al. A Hybrid PID‐Fuzzy Control for Linear SISO Systems with Variant Communication Delays
CN115616985B (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
GR01 Patent grant
GR01 Patent grant