CN113219823B - 一种基于约束跟踪的机械系统追踪-逃避-达点控制方法 - Google Patents

一种基于约束跟踪的机械系统追踪-逃避-达点控制方法 Download PDF

Info

Publication number
CN113219823B
CN113219823B CN202110090040.XA CN202110090040A CN113219823B CN 113219823 B CN113219823 B CN 113219823B CN 202110090040 A CN202110090040 A CN 202110090040A CN 113219823 B CN113219823 B CN 113219823B
Authority
CN
China
Prior art keywords
constraint
tracking
escape
arrival
servo
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
CN202110090040.XA
Other languages
English (en)
Other versions
CN113219823A (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.)
Nanjing University of Science and Technology
Original Assignee
Nanjing 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 Nanjing University of Science and Technology filed Critical Nanjing University of Science and Technology
Priority to CN202110090040.XA priority Critical patent/CN113219823B/zh
Publication of CN113219823A publication Critical patent/CN113219823A/zh
Application granted granted Critical
Publication of CN113219823B publication Critical patent/CN113219823B/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)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种基于约束跟踪的机械系统追踪‑逃避‑达点控制方法,基于拉格朗日建模方法,建立受控系统动力学模型;建立追逃约束和到达约束数学模型,分别将两种约束的控制目标转化成伺服约束,确定对应的约束矩阵、约束向量和约束跟踪误差;基于受控系统动力学模型,确定综合描述系统不确定性界值的函数,结合约束跟踪误差,构建自适应律;基于约束矩阵、约束向量和伺服约束跟踪误差,结合自适应律,构建自适应鲁棒控制器,对受控系统进行追踪‑逃避‑达点控制。本发明能够从约束跟踪的角度解决传统追逃问题,提高了控制精度和控制效率。

Description

一种基于约束跟踪的机械系统追踪-逃避-达点控制方法
技术领域
本发明涉及机械动力学以及控制领域,具体涉及一种基于约束跟踪的机械系统追踪-逃避-达点控制方法。
背景技术
目前对于逃避问题,常常被表述成追逃问题,并且经常与避撞控制一起提及。而且过去的追逃问题解决方法比较单一,常常被当做一个单纯的追逐或者逃避问题而被解决,而没有让它们同时到达一个特定的目标地点。同时,目前对该类问题的控制方法大多只能选定单一的控制目标,一旦控制目标改变,就需要对整个控制架构进行重新架构,控制效率和控制灵活性不足。而且由于该系统也具有不确定干扰,因此有必要加入到达约束,开发一种针对机械系统追踪-逃避-到达问题的,灵活性更强、精准度更高、抗时变不确定性干扰的运动控制方法。
发明内容
本发明的目的在于提供一种基于约束跟踪的机械系统追踪-逃避-达点控制方法。
实现本发明目的的技术方案为:一种基于约束跟踪的机械系统追踪-逃避-达点控制方法,包括以下步骤:
步骤1,基于拉格朗日建模方法,建立受控系统动力学模型;
步骤2,建立追逃约束和到达约束数学模型,分别将两种约束的控制目标转化成伺服约束,确定对应的约束矩阵、约束向量和约束跟踪误差;
步骤3,基于步骤1的受控系统动力学模型,确定综合描述系统不确定性界值的函数,结合步骤2构建的约束跟踪误差,构建自适应律;
步骤4,基于步骤2构建的约束矩阵、约束向量和伺服约束跟踪误差,结合自适应律,构建自适应鲁棒控制器,对步骤1的受控系统进行追踪-逃避-达点控制。
进一步的,步骤1中,基于拉格朗日建模方法,建立受控系统动力学模型,具体方法为:
基于拉格朗日建模方法,建立受控系统动力学模型:
Figure BDA0002912107400000011
其中
Figure BDA0002912107400000015
是不确定性,t表示时间,q(t)∈Rn表示坐标,其为时间的函数,对应的
Figure BDA0002912107400000012
分别为速度与加速度,后面为了书写简便,统一省略t,M(q,σ,t)是惯性矩阵,
Figure BDA0002912107400000013
是科里奥利离心力,g(q,σ,t)是重力,
Figure BDA0002912107400000014
是摩擦力及其他外部的干扰,τ∈Rh是控制输入力矩;
针对系统的不确定性处理问题,对该模型进行分解,把其动力学模型分解成标称部分和不确定部分:
Figure BDA0002912107400000021
Figure BDA00029121074000000213
Figure BDA0002912107400000022
Figure BDA0002912107400000023
其中
Figure BDA0002912107400000024
是标称部分,ΔM(·),ΔC(·),Δg(·),ΔF(·)为不确定部分。
进一步的,步骤2中,建立追逃约束和到达约束数学模型,分别将两种约束的控制目标转化成伺服约束,确定对应的约束矩阵、约束向量和约束跟踪误差,具体方法为:
(一)对于追逃约束
定义函数:
Figure BDA0002912107400000025
其中:
Figure BDA0002912107400000026
e1(t)即为追逃约束的数学模型,
Figure BDA0002912107400000027
分别是追踪区域与逃避区域中的点,r为逃避区域的半径,R为追踪区域的半径;
求导,得到追逃一阶伺服约束的数学模型:
Figure BDA0002912107400000028
其中
Figure BDA0002912107400000029
c1(q,t)为关于q,t的函数,为简单表达,后面简写为c1
求二阶导数得到追逃二阶伺服约束的数学模型:
Figure BDA00029121074000000210
Figure BDA00029121074000000211
其中
Figure BDA00029121074000000212
从而得到追逃约束矩阵、约束向量;
追逃约束矩阵:
Figure BDA0002912107400000031
追逃约束向量:
Figure BDA0002912107400000032
Figure BDA0002912107400000033
Figure BDA0002912107400000034
从而得到追逃约束跟踪误差:
Figure BDA0002912107400000035
(二)对于到达约束
定义函数:
Figure BDA0002912107400000036
其中e2(t)为到达约束的数学模型,
Figure BDA0002912107400000037
为运动目标区域的参考点坐标,s为目标到达区域的半径;
求导,得到到达一阶伺服约束的数学模型:
Figure BDA0002912107400000038
其中c2(q,t)为q,t的函数,为方便书写,后面简写为c2
求二阶导数得到到达二阶伺服约束的数学模型:
Figure BDA0002912107400000039
从而得到到达约束矩阵、约束向量;
到达约束矩阵:
Figure BDA00029121074000000310
到达约束向量:
c2=lh(q)-ls2
Figure BDA0002912107400000041
从而得到到达约束跟踪误差:
Figure BDA0002912107400000042
其中l>0为一个常数,同时为了简化书写,后续书写中用β12代表
Figure BDA0002912107400000043
进一步的,步骤3中,基于步骤1的受控系统动力学模型,确定综合描述系统不确定性界值的函数,结合步骤2构建的约束跟踪误差,构建自适应律,具体方法为:
基于步骤1的受控系统动力学模型,分析其不确定性,确定不确定性参数σ的一般形式,并通过以下不等式进行放缩变换,得到综合描述系统不确定性界值的函数∏(·):
Figure BDA0002912107400000044
Figure BDA0002912107400000045
其中
Figure BDA0002912107400000046
Figure BDA0002912107400000047
其中α是函数∏(·)中的不确定性变量,κ是控制增益,P是单位矩阵,M(q,t)是惯性矩阵,
Figure BDA0002912107400000048
是科里奥利离心力,g(q,σ,t)是重力,
Figure BDA0002912107400000049
是摩擦力及其他外部的干扰,
Figure BDA00029121074000000410
对应标称部分,ΔM(·),ΔC(·),Δg(·),ΔF(·)对应不确定部分,ρE>-1为常数,A(q)是综合两种约束的约束矩阵,表示为[A1(q),A2(q)]T,c是综合两种约束的约束向量,表示为[c1,c2]T
Figure BDA00029121074000000411
(·)-1表示逆矩阵,I是单位矩阵;
基于步骤2构建的综合两种约束的伺服约束跟踪误差
Figure BDA00029121074000000412
即[β12]T,以及本步骤所构建的函数∏(·),构建能够自评估不确定性变量α的自适应律:
Figure BDA00029121074000000413
其中
Figure BDA00029121074000000414
是对α的估计值,
Figure BDA00029121074000000415
k1,k2∈R,k1,k2>0为设计参数。
进一步的,步骤4中,基于步骤2构建的约束矩阵、约束向量和伺服约束跟踪误差,结合自适应律,构建自适应鲁棒控制器,对步骤1的受控系统进行追踪-逃避-达点控制,具体方法为:
针对步骤1的机械系统,基于步骤3所构建的约束矩阵A(q)和约束向量c,以及构建的伺服约束跟踪误差β,结合自适应律
Figure BDA0002912107400000051
构建鲁棒控制器:
Figure BDA0002912107400000052
Figure BDA0002912107400000053
其中
Figure BDA0002912107400000054
Figure BDA0002912107400000055
Figure BDA0002912107400000056
Figure BDA0002912107400000057
Figure BDA0002912107400000058
其中
Figure BDA0002912107400000059
分别表示控制输入力矩τ的一部分,ε>0为常数,使其满足步骤2的伺服约束。
一种基于约束跟踪的机械系统追踪-逃避-达点控制系统,基于任一项所述的方法进行基于约束跟踪的机械系统追踪-逃避-达点控制。
一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现任一项所述的方法进行基于约束跟踪的机械系统追踪-逃避-达点控制。
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现任一项所述的方法进行基于约束跟踪的机械系统追踪-逃避-达点控制。
本发明与现有技术方案相比,其显著特点为:1)在传统追逃问题中加入到达约束,使得设计出的控制器更加符合实际要求;2)应用约束跟踪方法对机械系统追踪-逃避-到达问题进行运动控制,把预期运动特性转换成伺服约束,构建约束跟踪误差作为控制器设计的控制对象,并针对系统不确定性设计自适应律,由此构建自适应鲁棒控制策略,最终同时满足三个看似不同的约束,不仅可以使受控系统抵御时变不确定性干扰,同时当控制目标改变时只需重构伺服约束,无需重新设计控制器,具有较强的灵活性和精准性。
附图说明
图1为本发明基于约束跟踪的机械系统追踪-逃避-达点控制方法的控制跟踪误差的效果图。
图2为本发明基于约束跟踪的机械系统追踪-逃避-达点控制方法有无追逃约束的控制效果对比图。
图3为本发明基于约束跟踪的机械系统追踪-逃避-达点控制方法达点控制部分的控制效果图。
图4为本发明基于约束跟踪的机械系统追踪-逃避-达点控制方法的流程图。
具体实施方式
本发明提出一种基于约束跟踪的机械系统追踪-逃避-达点控制方法,包括以下步骤:
步骤1、基于拉格朗日建模方法,建立机械系统动力学模型:
Figure BDA0002912107400000061
其中
Figure BDA00029121074000000612
是不确定性,t表示时间,q(t)∈Rn表示坐标,其为时间的函数,对应的
Figure BDA0002912107400000062
分别为速度与加速度,后面为了书写简便,统一省略t,M(q,σ,t)是惯性矩阵,
Figure BDA0002912107400000063
是科里奥利离心力,g(q,σ,t)是重力,
Figure BDA0002912107400000064
是摩擦力及其他外部的干扰,τ∈Rh是控制输入力矩;
针对该系统的不确定性处理问题,对该模型进行分解,把其动力学模型分解成标称部分和不确定部分:
Figure BDA0002912107400000065
Figure BDA0002912107400000066
Figure BDA0002912107400000067
Figure BDA0002912107400000068
其中
Figure BDA0002912107400000069
是标称部分,ΔM(·),ΔC(·),Δg(·),ΔF(·)为不确定部分。
为了方便后续控制器的设计,做如下定义:
Figure BDA00029121074000000613
Figure BDA00029121074000000610
其中(·)-1表示逆矩阵,I是单位矩阵。
基于机械系统动力学模型,分析该系统的预期运动特性,如位置,速度,轨迹等,同时给之后的计算步骤提供标准形式,将该特性进行数学抽象,并构建可以描述该运动特性的一阶伺服约束,并写成矩阵形式:
Figure BDA00029121074000000611
其中A(q)=[Ali]m×n是约束矩阵,c=[c1c2…cm]T是约束向量,Ali(·)和cl(·)一阶可导,然后,对该伺服约束进行求导,得到二阶伺服约束为:
Figure BDA0002912107400000071
其中
Figure BDA0002912107400000072
是约束向量。
构建约束跟踪误差β,作为后续控制器设计的控制对象:
Figure BDA0002912107400000073
其中β=[β1β2…βm]T
步骤2、对于追逃约束,定义函数:
Figure BDA0002912107400000074
其中:
Figure BDA0002912107400000075
Figure BDA0002912107400000076
分别是追踪区域与逃避区域中的点,r为逃避区域的半径,R为追踪区域的半径。
对函数求导,得到追逃一阶伺服约束的数学模型:
Figure BDA0002912107400000077
其中
Figure BDA0002912107400000078
再对函数求二阶导数得到追逃二阶伺服约束的数学模型:
Figure BDA0002912107400000079
Figure BDA00029121074000000710
其中
Figure BDA00029121074000000711
从而可以得到
追逃约束矩阵
Figure BDA00029121074000000712
追逃约束向量c1=0
Figure BDA0002912107400000081
Figure BDA0002912107400000082
从而得到追逃约束跟踪误差:
Figure BDA0002912107400000083
对于到达约束,定义函数:
Figure BDA0002912107400000084
其中
Figure BDA0002912107400000085
为运动目标区域的参考点坐标,s为目标到达区域的半径。
对函数求导,得到到达一阶伺服约束的数学模型:
Figure BDA0002912107400000086
再对函数求二阶导数得到到达二阶伺服约束的数学模型:
Figure BDA0002912107400000087
从而可以得到:
到达约束矩阵
Figure BDA0002912107400000088
到达约束向量c2=lh(q)-ls2
Figure BDA0002912107400000089
从而可以得到,到达约束跟踪误差:
Figure BDA00029121074000000810
其中l>0为一个常数,
Figure BDA00029121074000000811
分别表示追逃约束与到达约束的伺服约束约束向量。
步骤3、基于步骤1的机械系统动力学模型,分析其不确定性(模型不确定性、外界干扰等),确定不确定性参数σ的一般形式,并通过以下不等式进行放缩变换,得到可以综合描述系统不确定性界值的函数∏(·):
Figure BDA0002912107400000091
Figure BDA0002912107400000092
其中
Figure BDA0002912107400000093
Figure BDA0002912107400000094
其中α是函数∏(·)中的不确定性变量,κ是控制增益,P是单位矩阵,ρE>-1为常数,A(q)是综合两种约束的约束矩阵,表示为[A1(q),A2(q)]T,c是综合两种约束的约束向量,表示为[c1,c2]T
基于步骤2所构建的综合两种约束的约束跟踪误差
Figure BDA0002912107400000095
即[β12]T,以及本步骤所构建的函数∏(·),构建可以自评估不确定性变量α的自适应律:
Figure BDA0002912107400000096
其中
Figure BDA0002912107400000097
是对α的估计值,
Figure BDA0002912107400000098
k1,k2∈R,k1,k2>0为设计参数。
针对步骤1的机械系统,基于步骤2所构建的约束矩阵A和约束向量c,包括以及所构建的伺服约束跟踪误差β,并结合自适应律
Figure BDA0002912107400000099
设计鲁棒控制器:
Figure BDA00029121074000000910
Figure BDA00029121074000000911
其中
Figure BDA00029121074000000912
Figure BDA00029121074000000913
Figure BDA00029121074000000914
Figure BDA00029121074000000915
Figure BDA00029121074000000916
在对以上公式计算时,式中的
Figure BDA00029121074000000917
应该将追逃约束与到达约束各个分量组合成步骤1当中的标准形式,之后分别带入计算。
其中ε>0为常数,使其满足步骤2的伺服约束。
本发明还提出一种基于约束跟踪的机械系统追踪-逃避-达点控制系统,基于所述的方法进行基于约束跟踪的机械系统追踪-逃避-达点控制。
一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现所述的方法进行基于约束跟踪的机械系统追踪-逃避-达点控制。
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现所述的方法进行基于约束跟踪的机械系统追踪-逃避-达点控制。
实施例
为验证本发明方案的有效性,对如下航空领域的追踪-逃避-到达问题进行自适应鲁棒控制。控制对象和目标分别为:
(1)控制对象为两架航空器,其运动方程为:
Figure BDA0002912107400000101
Figure BDA0002912107400000102
该系统写成步骤1中建立的拉格朗日动力学模型有:q=[x1,y1,x2,y2]T,τ=[τ1234]T
M=diag(m1,m1,m2,m2),C=0,g=0,F=[s1,s2,s3,s4]T
其中x1,y1,x2,y2分别为两架航天器起始的坐标,m1,2是航天器的质量,τ1,2,3,4为输入的控制力,s1,2,3,4是系统收到的外部干扰(包括风阻等等)。
(2)控制目标:a)使得两架航天器拥有最小安全距离r>0;
b)使得两架航天器拥有最大追踪距离R>r>0。
步骤1:受控的系统的动力学模型如上所述,考虑系统外部干扰F=s1,2,3,4的不确定性,即
Figure BDA0002912107400000103
步骤2:由控制目标可知,对于追逃约束,设想的目标区域为:
Figure BDA0002912107400000104
根据以上步骤2中的一阶伺服约束,对于追逃约束,得到:
Figure BDA0002912107400000105
Figure BDA0002912107400000106
再对一阶约束求导,得到:
Figure BDA0002912107400000111
Figure BDA0002912107400000112
Figure BDA0002912107400000113
对目标到达区域施加给追踪者的到达约束,认为:
追踪者初始的位置坐标
Figure BDA0002912107400000114
追踪者的目标到达区域的参考点坐标为:
Figure BDA0002912107400000115
从而获取追踪者的目标区域为:
Figure BDA0002912107400000116
其中r1为其目标区域的半径。
再结合步骤2中的一阶伺服约束,对于到达约束,得到:
Figure BDA0002912107400000117
Figure BDA0002912107400000118
对一阶约束求导,得到:
Figure BDA0002912107400000119
对目标到达区域施加给逃避者的到达约束,认为:
逃避者初始的位置坐标
Figure BDA00029121074000001110
逃避者的目标到达区域,参考点坐标为:
Figure BDA00029121074000001111
从而获取逃避者的目标区域为:
Figure BDA00029121074000001112
其中r2为其目标区域的半径。
再结合步骤2中的一阶伺服约束,对于到达约束,得到:
Figure BDA00029121074000001113
Figure BDA00029121074000001114
对一阶约束求导,得到:
Figure BDA00029121074000001115
因此,约束矩阵A(q)=[A1,A2,A3]T,约束向量c=[c1,c2,c3]T
Figure BDA00029121074000001116
相应的约束跟踪误差
Figure BDA00029121074000001117
步骤3:针对系统的不确定性,构建自适应律:
Figure BDA0002912107400000121
其中
Figure BDA0002912107400000122
于是有
Figure BDA0002912107400000123
将步骤1中的系统运动单元
Figure BDA0002912107400000124
步骤2中的伺服约束矩阵及向量
Figure BDA0002912107400000125
约束跟踪误差β,步骤3中的自适应律
Figure BDA0002912107400000126
参数γ,μ,∏等信息代入
Figure BDA0002912107400000127
Figure BDA0002912107400000128
Figure BDA0002912107400000129
得到自适应鲁棒控制策略:
Figure BDA00029121074000001210
步骤4、采用Matlab进行仿真,并与未加入追逃约束的控制效果对比,仿真结果如图1,2所示。
图1中a,b,c图分别表示应用该控制方法的控制跟踪误差β123,可以发现,本发明所提机械系统追踪-逃避-到达的自适应鲁棒控制方法可以使受控系统在很短的时间内呈现预期运动特性,且系统稳定后跟踪误差波动范围极小,具备较高的控制精度,图2中a,b图分别表示两架航天器的最小安全距离与最大追踪距离随时间的变化,可以看出,在加入追逃约束之后,系统的性能表现更加稳定,较没有追逃约束的系统性能更好。图3中的a,b图分别表示两架航天器与目标区域的距离,可以看出,随着时间的变化,航天器与目标区域越来越接近,直到距离为0,即到达目标区域。因此本发明能够迅速有效地解决复杂机械系统追踪-逃避-达点控制问题。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。

Claims (4)

1.一种基于约束跟踪的机械系统追踪-逃避-达点控制方法,其特征在于,包括以下步骤:
步骤1,基于拉格朗日建模方法,建立受控系统动力学模型;
步骤2,建立追逃约束和到达约束数学模型,分别将两种约束的控制目标转化成伺服约束,确定对应的约束矩阵、约束向量和约束跟踪误差;
步骤3,基于步骤1的受控系统动力学模型,确定综合描述系统不确定性界值的函数,结合步骤2构建的约束跟踪误差,构建自适应律;
步骤4,基于步骤2构建的约束矩阵、约束向量和伺服约束跟踪误差,结合自适应律,构建自适应鲁棒控制器,对步骤1的受控系统进行追踪-逃避-达点控制;
步骤1中,基于拉格朗日建模方法,建立受控系统动力学模型,具体方法为:
基于拉格朗日建模方法,建立受控系统动力学模型:
Figure FDA0003843867060000011
其中
Figure FDA0003843867060000012
是不确定性,t表示时间,q(t)∈Rn表示坐标,其为时间的函数,对应的
Figure FDA0003843867060000013
分别为速度与加速度,后面为了书写简便,统一省略t,M(q,σ,t)是惯性矩阵,
Figure FDA0003843867060000014
是科里奥利离心力,g(q,σ,t)是重力,
Figure FDA0003843867060000015
是摩擦力及其他外部的干扰,τ∈Rh是控制输入力矩;
针对系统的不确定性处理问题,对该模型进行分解,把其动力学模型分解成标称部分和不确定部分:
Figure FDA0003843867060000016
Figure FDA0003843867060000017
Figure FDA0003843867060000018
Figure FDA0003843867060000019
其中
Figure FDA00038438670600000110
是标称部分,ΔM(·),ΔC(·),Δg(·),ΔF(·)为不确定部分;
步骤2中,建立追逃约束和到达约束数学模型,分别将两种约束的控制目标转化成伺服约束,确定对应的约束矩阵、约束向量和约束跟踪误差,具体方法为:
(一)对于追逃约束
定义函数:
Figure FDA00038438670600000111
其中:
Figure FDA0003843867060000021
e1(t)即为追逃约束的数学模型,
Figure FDA0003843867060000022
分别是追踪区域与逃避区域中的点,r为逃避区域的半径,R为追踪区域的半径;
求导,得到追逃一阶伺服约束的数学模型:
Figure FDA0003843867060000023
其中
Figure FDA0003843867060000024
c1(q,t)为关于q,t的函数,为简单表达,后面简写为c1
求二阶导数得到追逃二阶伺服约束的数学模型:
Figure FDA0003843867060000025
其中
Figure FDA0003843867060000026
从而得到追逃约束矩阵、约束向量;
追逃约束矩阵:
Figure FDA0003843867060000027
追逃约束向量:
c1=0
Figure FDA0003843867060000028
从而得到追逃约束跟踪误差:
Figure FDA0003843867060000029
(二)对于到达约束
定义函数:
Figure FDA0003843867060000031
其中e2(t)为到达约束的数学模型,
Figure FDA0003843867060000032
为运动目标区域的参考点坐标,s为目标到达区域的半径;
求导,得到到达一阶伺服约束的数学模型:
Figure FDA0003843867060000033
其中c2(q,t)为q,t的函数,为方便书写,后面简写为c2
求二阶导数得到到达二阶伺服约束的数学模型:
Figure FDA0003843867060000034
从而得到到达约束矩阵、约束向量;
到达约束矩阵:
Figure FDA0003843867060000035
到达约束向量:
c2=lh(q)-ls2
Figure FDA0003843867060000036
从而得到到达约束跟踪误差:
Figure FDA0003843867060000037
其中l>0为一个常数,同时为了简化书写,后续书写中用β12代表
Figure FDA0003843867060000038
步骤3中,基于步骤1的受控系统动力学模型,确定综合描述系统不确定性界值的函数,结合步骤2构建的约束跟踪误差,构建自适应律,具体方法为:
基于步骤1的受控系统动力学模型,分析其不确定性,确定不确定性参数σ的一般形式,并通过以下不等式进行放缩变换,得到综合描述系统不确定性界值的函数Π(·):
Figure FDA0003843867060000039
其中
Figure FDA00038438670600000310
Figure FDA0003843867060000041
其中α是函数Π(·)中的不确定性变量,κ是控制增益,P是单位矩阵,M(q,t)是惯性矩阵,
Figure FDA0003843867060000042
是科里奥利离心力,g(q,σ,t)是重力,
Figure FDA0003843867060000043
是摩擦力及其他外部的干扰,
Figure FDA0003843867060000044
对应标称部分,ΔM(·),ΔC(·),Δg(·),ΔF(·)对应不确定部分,ρE>-1为常数,A(q)是综合两种约束的约束矩阵,表示为[A1(q),A2(q)]T,c是综合两种约束的约束向量,表示为[c1,c2]T
Figure FDA0003843867060000045
(·)-1表示逆矩阵,I是单位矩阵;
基于步骤2构建的综合两种约束的伺服约束跟踪误差
Figure FDA0003843867060000046
即[β12]T,以及本步骤所构建的函数Π(·),构建能够自评估不确定性变量α的自适应律:
Figure FDA0003843867060000047
其中
Figure FDA0003843867060000048
是对α的估计值,
Figure FDA0003843867060000049
k1,k2∈R,k1,k2>0为设计参数;
步骤4中,基于步骤2构建的约束矩阵、约束向量和伺服约束跟踪误差,结合自适应律,构建自适应鲁棒控制器,对步骤1的受控系统进行追踪-逃避-达点控制,具体方法为:
针对步骤1的机械系统,基于步骤3所构建的约束矩阵A(q)和约束向量c,以及构建的伺服约束跟踪误差β,结合自适应律
Figure FDA00038438670600000410
构建鲁棒控制器:
Figure FDA00038438670600000411
Figure FDA00038438670600000412
其中
Figure FDA00038438670600000413
Figure FDA00038438670600000414
Figure FDA00038438670600000415
Figure FDA00038438670600000416
Figure FDA00038438670600000417
其中
Figure FDA00038438670600000418
分别表示控制输入力矩τ的一部分,ε>0为常数,使其满足步骤2的伺服约束。
2.一种基于约束跟踪的机械系统追踪-逃避-达点控制系统,其特征在于,基于权利要求1项所述的方法进行基于约束跟踪的机械系统追踪-逃避-达点控制。
3.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现权利要求1所述的方法进行基于约束跟踪的机械系统追踪-逃避-达点控制。
4.一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现权利要求1所述的方法进行基于约束跟踪的机械系统追踪-逃避-达点控制。
CN202110090040.XA 2021-01-22 2021-01-22 一种基于约束跟踪的机械系统追踪-逃避-达点控制方法 Active CN113219823B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110090040.XA CN113219823B (zh) 2021-01-22 2021-01-22 一种基于约束跟踪的机械系统追踪-逃避-达点控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110090040.XA CN113219823B (zh) 2021-01-22 2021-01-22 一种基于约束跟踪的机械系统追踪-逃避-达点控制方法

Publications (2)

Publication Number Publication Date
CN113219823A CN113219823A (zh) 2021-08-06
CN113219823B true CN113219823B (zh) 2022-11-04

Family

ID=77084475

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110090040.XA Active CN113219823B (zh) 2021-01-22 2021-01-22 一种基于约束跟踪的机械系统追踪-逃避-达点控制方法

Country Status (1)

Country Link
CN (1) CN113219823B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4860215A (en) * 1987-04-06 1989-08-22 California Institute Of Technology Method and apparatus for adaptive force and position control of manipulators
DE602007006748D1 (de) * 2007-01-23 2010-07-08 Gen Electric Multivariables Steuerungsdesignverfahren für mehrere Eingabe-/Ausgabe-Systeme mit mehreren Eingabe-/Ausgabe-Beschränkungen
CN108681257A (zh) * 2018-06-22 2018-10-19 合肥工业大学 一种主动防侧倾系统的控制器的设计方法
CN108803632A (zh) * 2018-09-19 2018-11-13 哈尔滨工程大学 基于饱和补偿技术的水面无人艇全状态约束轨迹跟踪控制方法
CN109100939A (zh) * 2018-09-19 2018-12-28 哈尔滨工程大学 考虑输入饱和的水面无人艇全状态约束轨迹跟踪控制方法
CN110046408A (zh) * 2019-03-29 2019-07-23 南京理工大学 基于bp神经网络微分与区间分析的区间不确定性优化方法
CN110376902A (zh) * 2019-08-19 2019-10-25 合肥工业大学 一种欠驱动机械系统伺服约束跟踪控制器的设计方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4860215A (en) * 1987-04-06 1989-08-22 California Institute Of Technology Method and apparatus for adaptive force and position control of manipulators
DE602007006748D1 (de) * 2007-01-23 2010-07-08 Gen Electric Multivariables Steuerungsdesignverfahren für mehrere Eingabe-/Ausgabe-Systeme mit mehreren Eingabe-/Ausgabe-Beschränkungen
CN108681257A (zh) * 2018-06-22 2018-10-19 合肥工业大学 一种主动防侧倾系统的控制器的设计方法
CN108803632A (zh) * 2018-09-19 2018-11-13 哈尔滨工程大学 基于饱和补偿技术的水面无人艇全状态约束轨迹跟踪控制方法
CN109100939A (zh) * 2018-09-19 2018-12-28 哈尔滨工程大学 考虑输入饱和的水面无人艇全状态约束轨迹跟踪控制方法
CN110046408A (zh) * 2019-03-29 2019-07-23 南京理工大学 基于bp神经网络微分与区间分析的区间不确定性优化方法
CN110376902A (zh) * 2019-08-19 2019-10-25 合肥工业大学 一种欠驱动机械系统伺服约束跟踪控制器的设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于干扰观测器的电液位置伺服系统跟踪控制;刘龙;《兵工学报》;20151130;第2053-2062页 *

Also Published As

Publication number Publication date
CN113219823A (zh) 2021-08-06

Similar Documents

Publication Publication Date Title
Chai et al. Six-DOF spacecraft optimal trajectory planning and real-time attitude control: a deep neural network-based approach
Shin Adaptive dynamic surface control for a hypersonic aircraft using neural networks
Hofer et al. Application of an approximate model predictive control scheme on an unmanned aerial vehicle
Grande et al. Experimental validation of Bayesian nonparametric adaptive control using Gaussian processes
Gandolfo et al. Trajectory tracking control of a pvtol aircraft based on linear algebra theory
Trapiello et al. Position‐heading quadrotor control using LPV techniques
Dai et al. Path planning of spatial rigid motion with constrained attitude
Wang et al. A new redundancy resolution for underwater vehicle–manipulator system considering payload
Zong et al. Nonlinear adaptive filter backstepping flight control for reentry vehicle with input constraint and external disturbances
Kim et al. Numerical time-scale separation for rotorcraft nonlinear optimal control analyses
Nivison et al. Development of a robust deep recurrent neural network controller for flight applications
Venkatraman et al. Online instrumental variable regression with applications to online linear system identification
Shahrooei et al. Falsification of learning-based controllers through multi-fidelity Bayesian optimization
Glida et al. Trajectory tracking control of a coaxial rotor drone: Time-delay estimation-based optimal model-free fuzzy logic approach
Hu et al. A robust adaptive nonlinear fault-tolerant controller via norm estimation for reusable launch vehicles
CN113219823B (zh) 一种基于约束跟踪的机械系统追踪-逃避-达点控制方法
Heidrich et al. Optimal information filtering for robust aerocapture trajectory generation and guidance
CN115047769A (zh) 一种基于约束跟随的无人作战平台避障-到达控制方法
Sharma et al. Near optimal finite-time terminal controllers for space trajectories via SDRE-based approach using dynamic programming
Shamsudin et al. Recursive Gauss-Newton based training algorithm for neural network modelling of an unmanned rotorcraft dynamics
Kang et al. Two-stage iterative finite-memory neural network identification for unmanned aerial vehicles
Ma et al. An improved central difference Kalman filter for satellite attitude estimation with state mutation
Yao et al. Koopman-operator-based safe learning control for spacecraft attitude reorientation with angular velocity constraints
Lotufo et al. Identification and control of a quadrotor from experimental data
Lal et al. Formally Verified Switching Logic for Recoverability of Aircraft Controller

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