CN108845493B - 带有输出约束的机械臂系统的固定时间跟踪控制方法 - Google Patents

带有输出约束的机械臂系统的固定时间跟踪控制方法 Download PDF

Info

Publication number
CN108845493B
CN108845493B CN201810955717.XA CN201810955717A CN108845493B CN 108845493 B CN108845493 B CN 108845493B CN 201810955717 A CN201810955717 A CN 201810955717A CN 108845493 B CN108845493 B CN 108845493B
Authority
CN
China
Prior art keywords
fixed time
equal
tracking
mechanical arm
trajectory
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
Application number
CN201810955717.XA
Other languages
English (en)
Other versions
CN108845493A (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.)
Qufu Normal University
Original Assignee
Qufu Normal University
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 Qufu Normal University filed Critical Qufu Normal University
Priority to CN201810955717.XA priority Critical patent/CN108845493B/zh
Publication of CN108845493A publication Critical patent/CN108845493A/zh
Application granted granted Critical
Publication of CN108845493B publication Critical patent/CN108845493B/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/0205Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system
    • G05B13/024Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric not using a model or a simulator of the controlled system 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

带有输出约束的机械臂系统的固定时间跟踪控制方法
技术领域
本发明涉及一种带有输出约束的机械臂系统的固定时间跟踪控制方法。
背景技术
机械臂系统虽然目前还不如人手那样灵活,但它具有能不断重复工作和劳动、不知疲劳、不怕危险、抓举重物的力量比人手大等特点,因此,机械臂已受到许多部门的重视,并越来越广泛地得到了应用。例如,机床加工工件的装卸,特别是在自动化车床、组合机床上的使用较为普遍;在装配作业中应用广泛,在电子行业中它可以用来装配印制电路板,在机械行业中它可以用来组装零部件;它可以在劳动条件差,单调重复易疲劳的工作环境工作,以代替人的劳动;它可以在危险的场合下工作,如军用品的装卸、危险品及有害物质的搬运等;还可用于宇宙及海洋的开发及军事工程和生物医学方面的研究和实验等。
由于机械臂系统可化为标准的二阶链式系统,其一直被科学界作为基准系统来研究。最近几十年,机械臂系统的轨迹跟踪控制得到了很大关注。现有很多技术对机械臂系统的轨迹跟踪进行了研究,滑模控制技术是主要的研究手段。针对机械臂系统,现有技术提出了全局非奇异终端滑模控制器以及有限时间终端滑模控制器。同时研究了多网络机械臂系统的分布式协同跟踪控制问题。值得注意的是现有技术只能得到有限时间稳定的结果,这意味着在有限时间控制中,停息时间的上界依赖于系统初始值。实际应用中,最好能够预先知道停息时间的上限,而且这个上限应该是一个正常数,与系统初始条件无关,即为固定时间稳定问题。相比于有限时间稳定,固定时间要求控制系统是全局有限时间稳定的并且停息时间的上界是一个与系统初始条件无关的常数。在实际应用中,固定时间控制比有限时间控制更可取,因为固定时间方法可以产生一个控制律,规定一个与操作域无关的过渡时间。为此,需要解决机械臂系统的固定时间轨迹跟踪控制问题。
针对输入通道带有扰动和不确定的二阶链式非线性系统,现有技术提出了固定时间终端滑模控制器,并将该思想成功应用于二阶多智能网络系统的一致跟踪控制问题中。针对带有不匹配不确定性的高阶严格反馈非线性系统,现有技术研究了固定时间状态反馈控制器设计问题。注意到实际中的很多系统可以转化为多输入系统,现有技术考虑了干扰不确定非完整系统的固定时间镇定控制器设计问题。另一方面,为满足实际的需求,约束非线性系统控制问题显得格外重要,例如桥式起重机在作业过程中,台车必须在一定的运行范围内运行,一旦台车超出运行界,便会与边界发生激烈的碰撞,造成安全事故。最近几年,带有状态约束或输出约束的非线性系统的控制设计成为了新的研究热点,其中障碍李雅普诺夫函数法主要用来处理链式非线性系统的约束控制问题。现有技术研究了带有输出约束的不确定欠驱动水面船舶的自适应神经网络轨迹跟踪控制问题。同时现有技术研究了带有输出约束和输入死区的机械臂系统的自适应神经网络跟踪控制。然而现有技术还不能解决带有输出约束的机械臂系统的固定时间轨迹跟踪控制问题,新的技术还尚待提出。
发明内容
本发明为了解决上述问题,提出了一种带有输出约束的机械臂系统的固定时间跟踪控制方法,本发明利用障碍李雅普诺夫函数法和增加幂积分器技术,有效的解决了带有输出约束的机械臂系统的固定时间轨迹跟踪控制问题。
为了实现上述目的,本发明采用如下技术方案:
一种带有输出约束的机械臂系统的固定时间跟踪控制方法,包括以下步骤:
描述刚体结构机械臂系统的动态方程;
不考虑外部干扰时,利用障碍李雅普诺夫函数法和有限时间控制算法,设立连续的固定时间轨迹跟踪控制律;
根据得到的控制律设计固定时间轨迹跟踪控制器,使得机械臂系统的动态方程的状态轨迹可以在固定时间内跟踪到预设的期望轨迹,同时在控制过程中跟踪误差不会超出预先设定的界;
依靠得到的固定时间轨迹跟踪控制器进行跟踪控制。
进一步的,刚体结构的机械臂系统,其动态方程描述为:
Figure BDA0001772629530000031
其中,系统广义坐标
Figure BDA0001772629530000032
是广义控制输入力,
Figure BDA0001772629530000033
为对称正定矩阵,
Figure BDA0001772629530000034
为中心力矩和哥氏力矩向量,
Figure BDA0001772629530000035
是引力矩向量,
Figure BDA0001772629530000036
表示外部干扰,
Figure BDA0001772629530000037
为n维欧氏空间。
进一步的,预设的期望轨迹qd
Figure BDA0001772629530000038
qdi,1≤i≤n为qd的第i个分量。
进一步的,预先设定的界,即|qi-qdi|<bi,1≤i≤n其中bi>0,qi,1≤i≤n为广义坐标
Figure BDA0001772629530000039
的第i个分量。
进一步的,当考虑外部干扰时,跟踪误差收敛到原点附近的有界域内。
进一步的,引入跟踪误差变量e1=q-qd
Figure BDA00017726295300000310
得到轨迹跟踪误差系统为:
Figure BDA0001772629530000041
Figure BDA0001772629530000042
其中e1=(e11,e12,...,e1n)T,e2=(e21,e22,...,e2n)T,设计控制器使得e1,e2的固定时间内收敛到零。
更进一步的,不考虑外部干扰时的固定时间轨迹跟踪控制器的设计过程包括:将e2看作虚拟控制输入,设计虚拟控制律
Figure BDA0001772629530000043
使得e1固定时间内收敛到零;设计实际的控制输入τ使得e2能在固定时间内跟踪到其虚拟控制律
Figure BDA0001772629530000044
更进一步的,虚拟控制输入的确定过程包括:选取李雅普诺夫函数,其沿轨迹跟踪误差系统关于时间的导数满足设定条件,选取虚拟控制输入,并代入设定条件,求解参数。
更进一步的,选取实际控制输入τ为:
Figure BDA0001772629530000045
其中x1=q为广义坐标,
Figure BDA0001772629530000047
为广义坐标导数,M(x1)为对称正定矩阵,C(x1,x2)为中心力矩和哥氏力矩阵,G(x1)是引力矩向量,
Figure BDA0001772629530000048
为期望轨迹的二阶导数,中间变量
Figure BDA0001772629530000049
虚拟控制输入
Figure BDA00017726295300000410
跟踪误差变量e1=[e11,e12,…,e1n]T=q-qd
Figure BDA00017726295300000411
函数
Figure BDA0001772629530000051
函数
Figure BDA0001772629530000052
常数bi>0,1≤i≤n,函数Πi和Ξi分别由(19)和(21)定义,
Figure BDA0001772629530000053
是定义在区间(-bi,bi)的一阶连续可微的函数,常数k2>0,γ>0,ρ=p1/p2>1,d=1+1/p,,1<p=p3/p4<2,pi,1≤i≤4是奇整数。
一种带有输出约束的机械臂系统的固定时间跟踪控制器,根据描述刚体结构机械臂系统的动态方程,在不考虑外部干扰时,利用障碍李雅普诺夫函数法和有限时间控制算法,设立连续的固定时间轨迹跟踪控制律;根据得到的控制律确定虚拟控制输入和实际输入,使得机械臂系统的动态方程的状态轨迹可以在固定时间内跟踪到预设的期望轨迹,同时在控制过程中跟踪误差不会超出预先设定的界。
与现有技术相比,本发明的有益效果为:
1、针对带有输出约束的机械臂系统,研究了固定时间轨迹跟踪控制器设计问题。当不考虑外部干扰时,基于障碍李雅普诺夫函数技术,设计了固定时间跟踪控制器,使得轨迹跟踪误差固定时间内收敛到零,同时在控制过程中不会超出预先设定的界,具有较好的收敛性;
2、当考虑外部干扰时,本发明的控制器可使轨迹跟踪误差有限时间内到达一个原点附近的有界域内,具有较好的抗干扰性和有效性。
附图说明
构成本申请的一部分的说明书附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。
图1为两关节刚体机械臂系统示意图;
图2为无干扰作用下的闭环系统响应效果图;
图3为干扰作用下的闭环系统响应效果图;
图4为本实施例的流程图;
具体实施方式:
下面结合附图与实施例对本发明作进一步说明。
应该指出,以下详细说明都是例示性的,旨在对本申请提供进一步的说明。除非另有指明,本实施例使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
在本发明中,术语如“上”、“下”、“左”、“右”、“前”、“后”、“竖直”、“水平”、“侧”、“底”等指示的方位或位置关系为基于附图所示的方位或位置关系,只是为了便于叙述本发明各部件或元件结构关系而确定的关系词,并非特指本发明中任一部件或元件,不能理解为对本发明的限制。
本发明中,术语如“固接”、“相连”、“连接”等应做广义理解,表示可以是固定连接,也可以是一体地连接或可拆卸连接;可以是直接相连,也可以通过中间媒介间接相连。对于本领域的相关科研或技术人员,可以根据具体情况确定上述术语在本发明中的具体含义,不能理解为对本发明的限制。
本实施例利用障碍李雅普诺夫函数法和增加幂积分器技术,研究了带有输出约束的机械臂系统的固定时间轨迹跟踪控制问题。本实施例首次考虑了机械臂系统的固定时间轨迹跟踪控制问题。首先,当不考虑外部干扰时,利用障碍李雅普诺夫函数法和有限时间控制方法,设立了连续的固定时间轨迹跟踪控制律。通过选取合适的李雅普诺夫函数证明了跟踪误差系统的固定时间收敛性。同时所设计的控制律可使得误差信号在控制过程中不会超出预先设定的界。其次,当考虑外部干扰时,本实施例给出了上述控制律的详细的抗干扰性能理论分析。建立了跟踪误差信号与设计参数之间的密切关系。最后,仿真结果说明了控制算法的有效性。
为方便理解,先介绍一些算法预备知识:
给定非线性系统
Figure BDA0001772629530000071
其中状态
Figure BDA0001772629530000072
函数
Figure BDA0001772629530000078
在包含原点x=0的开领域D内是连续的并且满足f(0)=0,为简化表示,系统的初始值表示为x(0)=x0
定义1:如果存在包含原点的开领域U,函数T(x0):U\{0}→(0,∞),使得系统(1)的解x(t,x0),x0∈U对任意的t∈[0,T(x0)]是有意义的并且
Figure BDA0001772629530000073
其中
Figure BDA0001772629530000074
则称系统(1)是有限时间收敛的。
如果系统(1)是李雅普诺夫稳定的并且是有限时间收敛的,则称其是有限时间稳定的。进一步地,如果
Figure BDA0001772629530000075
则称系统(1)是全局有限时间稳定的。
引理1:针对系统(1),如果存在连续正定的函数
Figure BDA0001772629530000076
使得
Figure BDA0001772629530000077
其中ε>0,λ∈(0,1),D0∈D是一个包含原点的开领域,则称系统(1)是有限时间稳定的。如果
Figure BDA0001772629530000081
则称系统(1)是全局有限时间稳定的。另外,停息时间T(x0)满足
Figure BDA0001772629530000082
很显然,引理1中的停息时间T(x0)的上界依赖于系统初始条件x0。但是一旦系统初始条件不可知,就不能预知停息时间T(x0)的上界。且通过式(3)可知,如果初始值x0充分大,停息时间T(x0)也可能变的非常大。这种现象限制了很多实际应用,比如在切换系统中,人们更期望事先知道下一个准确的切换时间。为此,提出了固定时间稳定的概念。
定义2:如果系统(1)是全局有限时间稳定的并且停息时间T(x0)是有界的,即存在一个正常数使得
Figure BDA0001772629530000083
则称系统(1)的原点是固定时间稳定的。
引理2:考虑系统(1),如果存在正定函数
Figure BDA0001772629530000084
使得
Figure BDA0001772629530000085
并且
Figure BDA0001772629530000086
其中α>0,β>0,p>0,q>0,pk<1,qk>1,则称系统(1)的原点是固定时间稳定的,并且
Figure BDA0001772629530000087
引理3:令x1,...,xn,x,y是实数,且0<b≤1,则(|x1|+…+|xn|)b≤|x1|b+…+|xn|b≤n1-p(|x1|+…+|xn|)b。当b=p/q≤1时,其中p>0,q>0是奇数,则|xb-yb|≤21-b|x-y|b
引理4:令c,d是正实数,则对任意的实值函数
Figure BDA0001772629530000091
成立不等式
Figure BDA0001772629530000092
引理5:对任意的正常数向量
Figure BDA0001772629530000093
则对满足|x|<|b|向量的
Figure BDA0001772629530000094
成立不等式
Figure BDA0001772629530000095
引理6:令x1,x2,...,xn≥0并且p>1。那么
Figure BDA0001772629530000096
本实施例研究刚体结构的机械臂系统,其动态方程描述为:
Figure BDA0001772629530000097
其中系统广义坐标
Figure BDA0001772629530000098
是广义控制输入力,
Figure BDA0001772629530000099
为对称正定矩阵,
Figure BDA00017726295300000910
为中心力矩和哥氏力矩向量,
Figure BDA00017726295300000911
是引力矩向量,
Figure BDA00017726295300000912
表示外部干扰。
给定期望跟踪轨迹
Figure BDA00017726295300000913
本实施例的控制目标是:设计固定时间轨迹跟踪控制器τ使得系统(6)的状态轨迹q可以在固定时间内跟踪到期望轨迹qd同时在控制过程中跟踪误差不会超出预先设定的界,即|qi-qdi|<bi,1≤i≤n其中bi>0。当考虑外部干扰时,跟踪误差收敛到原点附近的有界域内。为实现上述目标,本实施例需要的假设条件如下:
假设1:矩阵M-1(q)满足||M-1(q)||≤β1,其中常数β1>0。
假设2:参考轨迹qd及其一阶和二阶导数是有界的。
假设3:外部干扰满足|di|≤l<+∞,1≤i≤n,其中l是已知常数。
通过定义
Figure BDA00017726295300000914
系统(6)可重新改写为
Figure BDA0001772629530000101
Figure BDA0001772629530000102
下面给出可实现上述控制目标的控制器设计过程。
引入跟踪误差变量e1=q-qd
Figure BDA0001772629530000103
可得轨迹跟踪误差系统为
Figure BDA0001772629530000104
Figure BDA0001772629530000105
其中e1=(e11,e12,...,e1n)T,e2=(e21,e22,...,e2n)T。接下来,只需设计控制器τ使得e1,e2固定时间内收敛到零。
情况一:不考虑外部干扰时的固定时间轨迹跟踪控制器设计
设计过程可以分成两步。首先,将e2看作虚拟控制输入,设计虚拟控制律
Figure BDA0001772629530000106
使得e1固定时间内收敛到零。然后设计实际的控制输入τ使得e2能在固定时间内跟踪到其虚拟控制律
Figure BDA0001772629530000107
第一步:虚拟控制输入e2设计。
选取李雅普诺夫函数
Figure BDA0001772629530000108
其沿系统(8)关于时间的导数满足
Figure BDA0001772629530000109
选取虚拟控制输入
Figure BDA00017726295300001010
从而将(11)代入到(10),可知
Figure BDA0001772629530000111
其中
Figure BDA0001772629530000112
1<p=p3/p4<2,pi,1≤i≤4是奇整数。
第二步:控制输入τ设计。
引入
Figure BDA0001772629530000113
并选取李雅普诺夫函数
Figure BDA0001772629530000114
其中函数Vi的定义如下
Figure BDA0001772629530000115
现有已知
Figure BDA0001772629530000116
是正定,径向无界并且可微。对(12)中的最后一项,由引理3,可推得
Figure BDA0001772629530000117
进而,由引理4可知
Figure BDA0001772629530000118
其中
Figure BDA0001772629530000119
是定义在区间(-bi,bi)的一阶连续可微的函数,即βi∈C1。将(16)代入到(12)中,可得
Figure BDA00017726295300001110
另外,函数Vi沿系统(8)关于时间的导数可计算为
Figure BDA0001772629530000121
其中
Figure BDA0001772629530000122
为简化表达,定义
Figure BDA0001772629530000123
进而,由(11)可知
Figure BDA0001772629530000124
至此,
Figure BDA0001772629530000125
关于e1i的导数可表示为
Figure BDA0001772629530000126
从而可得
Figure BDA0001772629530000127
其中
Figure BDA0001772629530000128
的表达式为
Figure BDA0001772629530000129
基于式(22)中的结果,由引理4,可知(18)式右端的第一项满足
Figure BDA0001772629530000131
其中
Figure BDA0001772629530000132
的表达式为
Figure BDA0001772629530000133
将(24)代入到(18),可知
Figure BDA0001772629530000134
联立(17)和(25),可得
Figure BDA0001772629530000135
参照(26)的结构,选取实际控制输入τ为
Figure BDA0001772629530000136
因此,可得
Figure BDA0001772629530000141
进一步可得
Figure BDA0001772629530000142
结合(13),(29),引理5,进一步可推知
Figure BDA0001772629530000143
其中c=max{1/2,21-1/p/(2-1/p)}。再考虑到引理3,可知
Figure BDA0001772629530000144
紧接着,由引理6,可算得
Figure BDA0001772629530000145
因此,由(31)和(32)可得
Figure BDA0001772629530000146
其中c1=γ61-(1+ρ)/2c-(1+ρ)/2,0<k1<k2-λcd/2
该部分的主要结果归纳在下面的定理1中。
定理1:考虑满足假设条件1-2的系统(6)。如果控制输入τ由(27)定义,则系统轨迹q固定时间内跟踪到期望轨迹qd
证明:由于d/2<1,(1+ρ)/2>1,所以利用引理2,可知(6)与(27)构成的闭环系统是固定时间稳定的,并且停息时间的上界t1可根据引理2计算其详细值。
情况二:考虑外部干扰时的抗干扰性能分析
在这部分,我们考虑外部干扰的影响,主要结果归纳在下面的定理2中。
定理2:考虑满足假设条件1-3的系统(6)。如果控制输入选取为(27),则轨迹跟踪误差有限时间内收敛到有界域Δ内,其中
Figure BDA0001772629530000151
其中正常数χ,χ12i的具体定义在下面给出。
证明:尽管本部分考虑了外部干扰,证明过程与定理1的思路基本一致。将(27)代入到(26),并利用假设条件3可知
Figure BDA0001772629530000152
其中μ=lβ/(2-1/ρ)。下面的证明分成三步。首先,定义两个领域
Figure BDA0001772629530000161
Figure BDA0001772629530000162
其中ξ=[ξ12,...,ξn]。整体的分析过程可描述为:首先证明一旦
Figure BDA0001772629530000163
则有
Figure BDA0001772629530000164
但是这并不能保证Δ1就是吸引域,因为在领域Δ1内部,不能保证
Figure BDA0001772629530000165
所以状态进入领域Δ1后,又可能逃离该领域。然而,我们可以寻找一个更大的包含Δ1的吸引域Δ2。接着我们说明对任意的初始状态
Figure BDA0001772629530000166
Δ2是一个吸引域并且系统状态可在有限时间内到达该领域。最后,给出跟踪误差的收敛估计。
第一步:假定
Figure BDA0001772629530000167
Figure BDA0001772629530000168
接下来分两种情况讨论。
情况1:
Figure BDA0001772629530000169
由于1<p<2,d=1+1/p,所以(2-1/p)/d<1。利用引理3,可推得
Figure BDA00017726295300001610
进而,可知
Figure BDA00017726295300001611
因此,
Figure BDA0001772629530000171
另外,注意到
Figure BDA0001772629530000172
由引理3,推得
Figure BDA0001772629530000173
从而,
Figure BDA0001772629530000174
基于以上结果,可知
Figure BDA0001772629530000175
因此由(35)可知
Figure BDA0001772629530000176
情况2:
Figure BDA0001772629530000177
如果
Figure BDA0001772629530000178
Figure BDA0001772629530000179
意味着
Figure BDA00017726295300001710
另一方面,基于(38)中的结果,可得到
Figure BDA00017726295300001711
至此,根据以上结果,可推得
Figure BDA00017726295300001712
因此,综合情况1和2中的结果可知,一旦
Figure BDA00017726295300001713
则有
Figure BDA00017726295300001714
第二步:下面将要说明一旦
Figure BDA00017726295300001715
则有
Figure BDA00017726295300001716
事实上,由第一步中的分析可知,只需说明
Figure BDA00017726295300001717
接下来,给出分析。对任意的
Figure BDA00017726295300001718
基于(31),可知
Figure BDA0001772629530000181
因此,(44)意味着
Figure BDA0001772629530000182
Figure BDA0001772629530000183
从而,只有
Figure BDA0001772629530000184
就有
Figure BDA0001772629530000185
由Δ2的定义可知,系统状态
Figure BDA0001772629530000186
能在有限时刻t1内到达领域Δ2内,即对任意的t>t1,有V(t)<χ。接下来,给出跟踪误差的收敛估计。
第三步:由函数V的定义可知,当t>t1时,有
Figure BDA0001772629530000187
其中
Figure BDA0001772629530000188
由(45),进一步可推知
Figure BDA0001772629530000189
进而,可得
Figure BDA00017726295300001810
现在给出e2i的收敛估计。由引理3,可知
Figure BDA00017726295300001811
进而成立
Figure BDA00017726295300001812
如果
Figure BDA00017726295300001813
由上面的不等式以及Vi的定义,可推得
Figure BDA00017726295300001814
如果
Figure BDA00017726295300001815
也可得到同样的不等式。由于当t≥t1时,Vi(t)≤V(t)≤χ,所以当t≥t1时有
Figure BDA00017726295300001816
定义
Figure BDA00017726295300001817
基于上述不等式以及e1i(t)的收敛估计,可知当t>t1时有
Figure BDA0001772629530000191
为验证上述系统有效,进行验证仿真实验。
在本部分,为验证控制算法的有效性,我们考虑下面的两关节刚体机械臂系统。
Figure BDA0001772629530000192
其中:
Figure BDA0001772629530000193
Figure BDA0001772629530000194
m22=m2r2 2+J2,C12(q2)=m2r1r2sin(q2),
G1=[(m1+m2)r1cos(q2)+m2r2cos(q1+q2)],
G2=m2r2cos(q1+q2),g=9.8m/s。
在仿真中,选取物理参数r1=1m,r2=0.8m,J1=5kg·m,J2=5kg·m,m1=0.5kg,m2=1.5kg。控制参数选取为p=5/3,k2=1,γ=2,ρ=3,b1=1,b2=1。系统初始值设置为q1(0)=0.1deg,q2(0)=-0.1deg,
Figure BDA0001772629530000195
期望位置轨迹选取为q1d=0.6deg,q2d=-0.5deg。仿真结果如图2所示。不难发现机械臂的角度位置q很快地跟踪到期望轨迹qd=(0.6,-0.5)T,同时在控制过程中不会超出预先设定的界|q1-0.6|<1,|q2+0.5|<1。与此同时,由图2可知,当角度位置q到达期望位置qd后,角速度
Figure BDA0001772629530000201
很快地收敛到零。仿真结果说明了控制策略的有效性。
最后,验证本实施例提出的固定时间控制器对外部干扰的鲁棒性。为此,假定输入通道的外部干扰为d1(t)=10/(1+t),d2(t)=20/(1+t)。仿真结果如图3所示。从图3可知,即使有外部干扰的存在,本实施例设计的固定时间跟踪控制算法仍然能够得到满意的控制效果。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。

Claims (7)

1.一种带有输出约束的机械臂系统的固定时间跟踪控制方法,其特征是:包括以下步骤:
描述刚体结构机械臂系统的动态方程;
不考虑外部干扰时,利用障碍李雅普诺夫函数法和有限时间控制算法,设立连续的固定时间轨迹跟踪控制律;e1=q-qd
Figure FDA0002973601150000011
将e2看作虚拟控制输入,设计虚拟控制律
Figure FDA0002973601150000012
使得e1固定时间内收敛到零;设计实际的控制输入τ使得e2能在固定时间内跟踪到其虚拟控制律
Figure FDA0002973601150000013
虚拟控制输入的确定过程包括:选取障碍李雅普诺夫函数,
Figure FDA0002973601150000014
其沿轨迹跟踪误差系统关于时间的导数满足设定条件,
Figure FDA0002973601150000015
选取虚拟控制输入,
Figure FDA0002973601150000016
并代入设定条件,求解参数,
Figure FDA0002973601150000017
其中k2>0,γ>0,ρ=p1/p2>1,d=1+1/p,
Figure FDA0002973601150000018
pi是奇整数,1≤i≤4;选取实际控制输入τ为:
Figure FDA0002973601150000019
其中x1=q为广义坐标,
Figure FDA00029736011500000110
为广义坐标导数,M(x1)为对称正定矩阵,C(x1,x2)为中心力矩和哥氏力矩阵,G(x1)是引力矩向量,
Figure FDA00029736011500000111
为期望轨迹的二阶导数,中间变量
Figure FDA0002973601150000021
其中ρ=p1/p2>1,1<p=p3/p4<2,e2i=qi-qdi,虚拟控制输入
Figure FDA0002973601150000022
跟踪误差变量e1=[e11,e12,…,e1n]T=q-qd
Figure FDA0002973601150000023
其中qd为预设的期望轨迹,
Figure FDA0002973601150000024
函数
Figure FDA0002973601150000025
函数
Figure FDA0002973601150000026
常数bi>0,1≤i≤n,函数Πi和Ξi分别由
Figure FDA0002973601150000027
Figure FDA0002973601150000028
定义,
Figure FDA0002973601150000029
是定义在区间(-bi,bi)的一阶连续可微的函数,常数k2>0,γ>0,ρ=p1/p2>1,d=1+1/p,1<p=p3/p4<2,pi是奇整数,1≤i≤4;
根据得到的固定时间轨迹跟踪控制律设计固定时间轨迹跟踪控制器,使得机械臂系统的动态方程的状态轨迹可以在固定时间内跟踪到预设的期望轨迹,同时在控制过程中跟踪误差不会超出预先设定的界;
依靠得到的固定时间轨迹跟踪控制器进行跟踪控制。
2.如权利要求1所述的一种带有输出约束的机械臂系统的固定时间跟踪控制方法,其特征是:刚体结构的机械臂系统,其动态方程描述为:
Figure FDA0002973601150000031
其中,系统广义坐标
Figure FDA0002973601150000032
是广义控制输入力,
Figure FDA0002973601150000033
为对称正定矩阵,
Figure FDA0002973601150000034
为中心力矩和哥氏力矩向量,
Figure FDA0002973601150000035
是引力矩向量,
Figure FDA0002973601150000036
表示外部干扰,
Figure FDA0002973601150000037
为n维欧氏空间。
3.如权利要求2所述的一种带有输出约束的机械臂系统的固定时间跟踪控制方法,其特征是:预设的期望轨迹qd
Figure FDA0002973601150000038
qdi,1≤i≤n为qd的第i个分量。
4.如权利要求3所述的一种带有输出约束的机械臂系统的固定时间跟踪控制方法,其特征是:预先设定的界,即|qi-qdi|<bi,1≤i≤n其中bi>0,qi,1≤i≤n为广义坐标
Figure FDA0002973601150000039
的第i个分量。
5.如权利要求1所述的一种带有输出约束的机械臂系统的固定时间跟踪控制方法,其特征是:当考虑外部干扰时,跟踪误差收敛到原点附近的有界域内。
6.如权利要求3所述的一种带有输出约束的机械臂系统的固定时间跟踪控制方法,其特征是:引入跟踪误差变量e1=q-qd
Figure FDA00029736011500000310
得到轨迹跟踪误差系统为:
Figure FDA00029736011500000312
Figure FDA00029736011500000311
其中e1=(e11,e12,...,e1n)T,e2=(e21,e22,...,e2n)T,设计控制器使得e1,e2的固定时间内收敛到零。
7.一种带有输出约束的机械臂系统的固定时间跟踪控制器,其特征是:根据描述刚体结构机械臂系统的动态方程,在不考虑外部干扰时,利用障碍李雅普诺夫函数法和有限时间控制算法,设立连续的固定时间轨迹跟踪控制律;
e1=q-qd
Figure FDA0002973601150000041
将e2看作虚拟控制输入,设计虚拟控制律
Figure FDA0002973601150000042
使得e1固定时间内收敛到零;设计实际的控制输入τ使得e2能在固定时间内跟踪到其虚拟控制律
Figure FDA0002973601150000043
虚拟控制输入的确定过程包括:选取障碍李雅普诺夫函数,
Figure FDA0002973601150000044
其沿轨迹跟踪误差系统关于时间的导数满足设定条件,
Figure FDA0002973601150000045
选取虚拟控制输入,
Figure FDA0002973601150000046
并代入设定条件,求解参数,
Figure FDA0002973601150000047
其中k2>0,γ>0,ρ=p1/p2>1,d=1+1/p,
Figure FDA0002973601150000048
pi是奇整数,1≤i≤4;选取实际控制输入τ为:
Figure FDA0002973601150000049
其中x1=q为广义坐标,
Figure FDA00029736011500000410
为广义坐标导数,M(x1)为对称正定矩阵,C(x1,x2)为中心力矩和哥氏力矩阵,G(x1)是引力矩向量,
Figure FDA00029736011500000411
为期望轨迹的二阶导数,中间变量
Figure FDA0002973601150000051
其中ρ=p1/p2>1,1<p=p3/p4<2,e2i=qi-qdi,虚拟控制输入
Figure FDA0002973601150000052
跟踪误差变量e1=[e11,e12,…,e1n]T=q-qd
Figure FDA0002973601150000053
其中qd为预设的期望轨迹,
Figure FDA0002973601150000054
函数
Figure FDA0002973601150000055
函数
Figure FDA0002973601150000056
常数bi>0,1≤i≤n,函数Πi和Ξi分别由
Figure FDA0002973601150000057
Figure FDA0002973601150000058
定义,
Figure FDA0002973601150000059
是定义在区间(-bi,bi)的一阶连续可微的函数,常数k2>0,γ>0,ρ=p1/p2>1,d=1+1/p,1<p=p3/p4<2,pi是奇整数,1≤i≤4;
根据得到的控制律确定虚拟控制输入和实际输入,使得机械臂系统的动态方程的状态轨迹可以在固定时间内跟踪到预设的期望轨迹,同时在控制过程中跟踪误差不会超出预先设定的界。
CN201810955717.XA 2018-08-21 2018-08-21 带有输出约束的机械臂系统的固定时间跟踪控制方法 Active CN108845493B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810955717.XA CN108845493B (zh) 2018-08-21 2018-08-21 带有输出约束的机械臂系统的固定时间跟踪控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810955717.XA CN108845493B (zh) 2018-08-21 2018-08-21 带有输出约束的机械臂系统的固定时间跟踪控制方法

Publications (2)

Publication Number Publication Date
CN108845493A CN108845493A (zh) 2018-11-20
CN108845493B true CN108845493B (zh) 2021-04-27

Family

ID=64188206

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810955717.XA Active CN108845493B (zh) 2018-08-21 2018-08-21 带有输出约束的机械臂系统的固定时间跟踪控制方法

Country Status (1)

Country Link
CN (1) CN108845493B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109581873A (zh) * 2018-12-26 2019-04-05 河海大学 未知执行器死区切换系统的有限时间指定性能控制算法
CN109782601B (zh) * 2019-01-31 2021-09-14 浙江大学 一种协调机械臂自适应神经网络同步鲁棒控制器设计方法
CN110095989B (zh) * 2019-05-29 2023-02-14 南京邮电大学 一种基于反步法的分布式多拉格朗日系统跟踪控制策略
CN110456641B (zh) * 2019-07-23 2022-03-11 西北工业大学 一种固定时预定性能循环神经网络机械臂控制方法
CN112148036B (zh) * 2020-09-11 2021-08-03 中国地质大学(武汉) 网络化机器人系统的固定时间估计器的双边跟踪控制方法
CN112987575B (zh) * 2021-03-05 2023-09-15 中国矿业大学 一种电液伺服系统位置闭环跟踪误差限定控制方法
CN113534666B (zh) * 2021-07-29 2023-03-03 河南科技大学 多目标约束下单关节机械臂系统的轨迹跟踪控制方法
CN116141339B (zh) * 2023-04-19 2023-07-21 珞石(北京)科技有限公司 一种七自由度机械臂预定时间轨迹跟踪控制方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5959861A (en) * 1996-10-31 1999-09-28 Fujitsu Limited Adaptive robust control device
CN104216284A (zh) * 2014-08-14 2014-12-17 浙江工业大学 机械臂伺服系统的有限时间协同控制方法
CN104570740A (zh) * 2015-01-21 2015-04-29 江南大学 一种输入饱和机械臂系统的周期自适应学习控制方法
CN105093934A (zh) * 2015-08-17 2015-11-25 哈尔滨工业大学 考虑干扰与模型不确定性的多机器人系统分布式有限时间跟踪控制方法
CN105319972A (zh) * 2015-11-27 2016-02-10 燕山大学 基于快速终端滑模的遥操作机器人固定时间控制方法
CN106249602A (zh) * 2016-09-30 2016-12-21 山东大学 桥式吊车有限时间轨迹跟踪控制器及其设计方法
CN107662208A (zh) * 2017-08-24 2018-02-06 浙江工业大学 一种基于神经网络的柔性关节机械臂有限时间自适应反步控制方法
CN107817818A (zh) * 2017-12-07 2018-03-20 中国人民解放军国防科技大学 一种模型不确定飞艇航迹跟踪有限时间控制方法
CN107942672A (zh) * 2017-12-06 2018-04-20 浙江工业大学 一种基于对称时不变障碍李雅普诺夫函数的四旋翼飞行器输出受限反步控制方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5959861A (en) * 1996-10-31 1999-09-28 Fujitsu Limited Adaptive robust control device
CN104216284A (zh) * 2014-08-14 2014-12-17 浙江工业大学 机械臂伺服系统的有限时间协同控制方法
CN104570740A (zh) * 2015-01-21 2015-04-29 江南大学 一种输入饱和机械臂系统的周期自适应学习控制方法
CN105093934A (zh) * 2015-08-17 2015-11-25 哈尔滨工业大学 考虑干扰与模型不确定性的多机器人系统分布式有限时间跟踪控制方法
CN105319972A (zh) * 2015-11-27 2016-02-10 燕山大学 基于快速终端滑模的遥操作机器人固定时间控制方法
CN106249602A (zh) * 2016-09-30 2016-12-21 山东大学 桥式吊车有限时间轨迹跟踪控制器及其设计方法
CN107662208A (zh) * 2017-08-24 2018-02-06 浙江工业大学 一种基于神经网络的柔性关节机械臂有限时间自适应反步控制方法
CN107942672A (zh) * 2017-12-06 2018-04-20 浙江工业大学 一种基于对称时不变障碍李雅普诺夫函数的四旋翼飞行器输出受限反步控制方法
CN107817818A (zh) * 2017-12-07 2018-03-20 中国人民解放军国防科技大学 一种模型不确定飞艇航迹跟踪有限时间控制方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Barrier Lyapunov functions for the control of output-constrained nonlinear systems;Keng Peng Tee,等;《Automatica》;20091231;第45卷;第918-927页 *
Global Adaptive Tracking control of robot manipulators Using Neural Networks with Finite-time Learning Convergence;Chenguang Yang,等;《International Journal of Control,Automation and Systems》;20171231;第15卷(第4期);第1916-1924页 *
刚性机械手有限时间自适应跟踪控制;陈刚,等;《Proceedings of the 32nd Chinese Control Conference》;20130728;第3003-3008页 *
控制方向未知的全状态约束非线性系统的鲁棒自适应跟踪控制;王春晓,等;《控制理论与应用》;20180228;第35卷(第2期);第153-161页 *

Also Published As

Publication number Publication date
CN108845493A (zh) 2018-11-20

Similar Documents

Publication Publication Date Title
CN108845493B (zh) 带有输出约束的机械臂系统的固定时间跟踪控制方法
CN108828955B (zh) 基于有限时间扩张状态观测器的精准航迹跟踪控制方法
CN108983606B (zh) 一种机械臂系统的鲁棒滑模自适应控制方法
CN109343350A (zh) 一种基于模型预测控制的水下机器人路径跟踪控制方法
CN109514564B (zh) 一种复合二次型多关节机械臂最优控制方法
Chwa Robust nonlinear disturbance observer based adaptive guidance law against uncertainties in missile dynamics and target maneuver
Zhao et al. Motion trajectory tracking and sway reduction for double-pendulum overhead cranes using improved adaptive control without velocity feedback
Lin et al. Guidance law design by adaptive fuzzy sliding-mode control
Xia et al. Modeling and compound control for unmanned turret system with coupling
Liu et al. Finite-time attitude control: a finite-time passivity approach
Lamraoui et al. Path following control of fully actuated autonomous underwater vehicle based on LADRC
Truong et al. An adaptive terminal sliding mode control scheme via neural network approach for path-following control of uncertain nonlinear systems
CN110908389A (zh) 一种针对不确定水下机器人的自适应快速速度跟踪控制方法
Wang et al. Adaptive Sliding Mode Fixed‐Time Tracking Control Based on Fixed‐Time Sliding Mode Disturbance Observer with Dead‐Zone Input
Adelhedi et al. Adaptive sliding mode control design of a SCARA robot manipulator system under parametric variations
Chaudhary et al. Fractional order fast terminal sliding mode control scheme for tracking control of robot manipulators
Cengiz et al. Optimal controller design for autonomous quadrotor landing on moving platform
Tianyi et al. A fractional order sliding mode controller design for spacecraft attitude control system
Fang et al. Trajectory tracking control for manned submersible system with disturbances via disturbance characterization index approach
Vizer et al. Gray-box LPV model identification of a 2-DoF surgical robotic manipulator by using an H∞-norm-based local approach
Wang et al. A novel prescribed performance controller with unknown dead-zone and impactive disturbance
Zhang et al. High-gain observer-based model predictive control for cross tracking of underactuated autonomous underwater vehicles
Yazdanzad et al. Combined neural network feedforward and RISE feedback control structure for a 5 DOF Upper-limb exoskeleton robot with asymptotic tracking
Gao et al. Tracking controllers of nonlinear output‐constrained surface ships subjected to external disturbances
Loueipour Design of dynamic positioning control system for an ROV with unknown dynamics using modified time delay estimation

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