CN110488603B - 一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法 - Google Patents

一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法 Download PDF

Info

Publication number
CN110488603B
CN110488603B CN201910884253.2A CN201910884253A CN110488603B CN 110488603 B CN110488603 B CN 110488603B CN 201910884253 A CN201910884253 A CN 201910884253A CN 110488603 B CN110488603 B CN 110488603B
Authority
CN
China
Prior art keywords
rigid aircraft
fixed time
sat
neural network
diag
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
CN201910884253.2A
Other languages
English (en)
Other versions
CN110488603A (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.)
Zhejiang University of Technology ZJUT
Original Assignee
Zhejiang University of Technology ZJUT
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 Zhejiang University of Technology ZJUT filed Critical Zhejiang University of Technology ZJUT
Publication of CN110488603A publication Critical patent/CN110488603A/zh
Application granted granted Critical
Publication of CN110488603B publication Critical patent/CN110488603B/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/0265Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion
    • G05B13/027Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric the criterion being a learning criterion using neural networks only
    • 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
    • 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/047Adaptive 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 the criterion being a time optimal performance criterion

Abstract

一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法,针对具有集中不确定性的刚性飞行器姿态跟踪问题,设计了固定时间滑模面,保证了状态的固定时间收敛;引入神经网络逼近总不确定的函数,设计了神经网络固定时间控制器。本发明在外界干扰,转动惯量不确定,执行器饱和和故障的因素下,实现飞行器系统的姿态跟踪误差和角速度误差固定时间一致最终有界的控制。

Description

一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪 控制方法
技术领域
本发明涉及一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法,特别是存在外部干扰,转动惯量矩阵不确定,执行器饱和和故障的刚性飞行器姿态跟踪方法。
背景技术
刚性飞行器姿态控制系统在刚性飞行器的健康,可靠的运动中扮演着重要的角色。在复杂的航天环境中,刚性飞行器姿态控制系统会受到各种外部干扰以及刚性飞行器在长期不断任务时存在的老化和失效等故障等影响。为了有效维持系统的性能,需要使其对外部干扰以及执行器故障具有较强的鲁棒性;另外,刚性飞行器还存在转动惯量矩阵不确定,因此控制饱和也是飞行器经常出现的问题。综上所述,刚性飞行器在执行任务时,需要一种在短时间内使系统稳定收敛,高精度的容错控制方法。
滑模控制在解决系统不确定性和外部扰动方面被认为是一个有效的鲁棒控制方法。滑模控制方法具有算法简单、响应速度快、对外界噪声干扰和参数摄动鲁棒性强等优点。终端滑模控制是一种可以实现有限时间稳定性的传统滑模控制的改进方案。然而,现存的有限时间技术估计收敛时间需要知道系统的初始信息,这对于设计者是很难知道的。近年来,固定时间技术得到了广泛的应用,固定时间控制方法与现存的有限时间控制方法相比,具有无需知道系统的初始信息,也能保守估计系统的收敛时间的优越性。
神经网络是线性参数化近似方法的中一种,可以被任意的其他近似方法取代,比如RBF神经网络,模糊逻辑系统等等。利用神经网络逼近不确定的性质,有效的结合固定时间滑模控制技术,减少外部干扰及系统参数不确定性对系统控制性能的影响,实现刚性飞行器姿态的固定时间控制。
发明内容
为了克服现有的刚性飞行器姿态控制系统存在的未知非线性问题,本发明提供一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法,并且在系统存在外部干扰,转动惯量不确定,执行器饱和和故障的情况下,实现系统状态的固定时间一致最终有界的控制方法。
为了解决上述技术问题提出的技术方案如下:
一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法,包括以下步骤:
步骤1,建立刚性飞行器的运动学和动力学模型,初始化系统状态以及控制参数,过程如下:
1.1刚性飞行器系统的运动学方程为:
Figure BDA0002206812530000021
Figure BDA0002206812530000022
其中qv=[q1,q2,q3]T和q4分别为单位四元数的矢量部分和标量部分且满足
Figure BDA0002206812530000023
q1,q2,q3分别为映射在空间直角坐标系x,y,z轴上的值;
Figure BDA0002206812530000024
分别是qv和q4的导数;Ω∈R3是刚性飞行器的角速度;I3是R3×3单位矩阵;
Figure BDA0002206812530000025
表示为:
Figure BDA0002206812530000026
1.2刚性飞行器系统的动力学方程为:
Figure BDA0002206812530000027
其中J∈R3×3是刚性飞行器的转动惯性矩阵;
Figure BDA0002206812530000028
是刚性飞行器的角加速度;u∈R3和d∈R3是控制力矩和外部扰动;D=diag(D1,D2,D3)∈R3×3是3×3对称对角的执行器效率矩阵,满足0<Di(t)≤1,i=1,2,3;sat(u)=[sat(u1),sat(u2),sat(u3)]T为执行器产生的实际控制力矩,sat(ui)为带有饱和特性的执行器,表示为sat(ui)=sgn(ui)min{umi,|ui|},umi为最大提供的控制力矩,sgn(ui)为符号函数,min{umi,|ui|}为两者的最小值;为了表示控制约束,sat(u)表示为sat(u)=g(u)+ds(u),g(u)=[g1(u1),g2(u2),g3(u3)]T,gi(ui)为双曲正切函数
Figure BDA0002206812530000031
ds(u)=[ds1(u1),ds2(u2),ds3(u3)]T为近似误差矢量;根据中值定理,gi(ui)=miui,0<mi≤1;定义H=DM=diag(δ1m12m23m3)∈R3×3为3×3对称对角矩阵,M=diag(m1,m2,m3)∈R3×3为3×3对称对角矩阵;Dsat(u)重新表示为:Dsat(u)=Hu+Dds(u),满足0<h0≤Dimi≤1,i=1,2,3,h0为未知正常数;Ω×表示为:
Figure BDA0002206812530000032
1.3刚性飞行器系统期望的运动学方程为:
Figure BDA0002206812530000033
Figure BDA0002206812530000034
其中qdv=[qd1,qd2,qd3]T和qd4分别为期望的单位四元数的矢量部分和标量部分且满足
Figure BDA0002206812530000035
Ωd∈R3为期望的角速度;
Figure BDA0002206812530000036
分别为qdv,qd4的导数,
Figure BDA0002206812530000037
为qdv的转置;
Figure BDA0002206812530000038
表示为:
Figure BDA0002206812530000039
1.4由四元数描述的刚性飞行器相对姿态运动:
Figure BDA00022068125300000310
Figure BDA0002206812530000041
Ωe=Ω-CΩd (12)
其中ev=[e1,e2,e3]T和e4分别为姿态跟踪误差的矢量部分和标量部分;Ωe=[Ωe1e2e3]T∈R3为角速度误差;
Figure BDA0002206812530000042
为相应的方向余弦矩阵并且满足||C||=1和
Figure BDA0002206812530000043
Figure BDA0002206812530000044
为C的导数;
根据式(1)-(12),刚性飞行器姿态跟踪误差动力学和运动学方程为:
Figure BDA0002206812530000045
Figure BDA0002206812530000046
Figure BDA0002206812530000047
其中
Figure BDA0002206812530000048
Figure BDA0002206812530000049
分别为ev和e4的导数;
Figure BDA00022068125300000410
为ev的转置;
Figure BDA00022068125300000411
Figure BDA00022068125300000412
分别为Ωd和Ωe的导数;(Ωe+CΩd)×与Ω×等价;
Figure BDA00022068125300000413
Figure BDA00022068125300000414
分别表示为:
Figure BDA00022068125300000415
Figure BDA00022068125300000416
1.5转动惯性矩阵J满足J=J0+ΔJ,其中J0和ΔJ分别表示J的标称部分和不确定部分,则式(15)重新写成:
Figure BDA00022068125300000417
进一步得到:
Figure BDA0002206812530000051
1.6对式(13)进行微分,得到:
Figure BDA0002206812530000052
其中
Figure BDA0002206812530000053
为ev的二阶导数;
步骤2,针对外部扰动,转动惯量不确定,执行器饱和和故障的刚性飞行器系统,设计所需的滑模面,过程如下:
选择固定时间滑模面为:
Figure BDA0002206812530000054
其中
Figure BDA0002206812530000055
Figure BDA0002206812530000056
和sgn(ei)均为符号函数,λ1>0,λ2>0,a2>1,
Figure BDA0002206812530000057
Figure BDA0002206812530000058
为ei的导数,i=1,2,3;
定义S=[S1,S2,S3]T,对S求导,得到:
Figure BDA0002206812530000059
将式(20)代入(22),得到:
Figure BDA00022068125300000510
其中
Figure BDA00022068125300000511
Figure BDA00022068125300000512
为总确定;L=[L1,L2,L3]T,
Figure BDA00022068125300000513
Figure BDA0002206812530000061
Γ=diag(Γ123)∈R3×3为3×3对称的对角矩阵;i=1,2,3;
步骤3,设计神经网络固定时间控制器,其过程如下:
3.1定义神经网络为:
Gi(Xi)=Wi *TΦ(Xi)+εi (24)
其中
Figure BDA0002206812530000062
为输入矢量,Φi(Xi)∈R4为神经网络基函数,Wi *∈R4为理想的权值矢量,定义为:
Figure BDA0002206812530000063
其中Wi∈R4为权值矢量,εi为近似误差,满足|εi|≤εN,i=1,2,3,εN为很小的正常数;arg min{·}为Wi *取其最小值所有的集合;
3.2考虑固定时间控制器被设计为:
Figure BDA0002206812530000064
其中
Figure BDA0002206812530000065
为3×3对称的对角矩阵,
Figure BDA0002206812530000066
Figure BDA0002206812530000067
为Θi的估计值Φ(X)=[Φ(X1),Φ(X2),Φ(X3)]T;K1=diag(k11,k12,k13)∈R3×3为3×3对称的对角矩阵;K2=diag(k21,k22,k23)∈R3×3为3×3对称的对角矩阵;K3=diag(k31,k32,k33)∈R3×3为对称的对角矩阵;k11,k12,k13,k21,k22,k23,k31,k32,k33为正常数;0<r1<1,r2>1;||·||表示值的二范数;
3.2设计更新律为:
Figure BDA0002206812530000068
其中γi>0,pi>0,i=1,2,3,
Figure BDA0002206812530000069
Figure BDA00022068125300000610
的导数,Φ(Xi)选择为以下的sigmoid函数:
Figure BDA00022068125300000611
其中l1,l2,l3和l4为近似参数,Φ(Xi)满足0<Φ(Xi)<Φ0,并且
Figure BDA0002206812530000071
步骤4,固定时间稳定性证明,其过程如下:
4.1证明刚性飞行器系统所有信号都是一致最终有界,设计李雅普诺夫函数为如下形式:
Figure BDA0002206812530000072
其中
Figure BDA0002206812530000073
ST是S的转置;
Figure BDA0002206812530000074
Figure BDA0002206812530000075
的转置;
对式(29)进行微分,得到:
Figure BDA0002206812530000076
其中
Figure BDA0002206812530000077
min{·}表示最小值;
因此,刚性飞行器系统所有信号都是一致最终有界的;
4.2证明固定时间收敛,设计李雅普诺夫函数为如下形式:
Figure BDA0002206812530000078
对式(31)进行微分,得到:
Figure BDA0002206812530000079
其中
Figure BDA00022068125300000710
i=1,2,3;υ2为一个大于零的上界值;
基于以上分析,刚性飞行器系统的姿态跟踪误差和角速度误差在固定时间一致最终有界。
本发明在外界干扰,转动惯量不确定,执行器饱和和故障的因素下,运用自适应神经网络跟踪控制方法,实现系统稳定跟踪,保证系统状态实现固定时间一致最终有界。本发明的技术构思为:针对含外部干扰,转动惯量不确定,执行器饱和和故障的刚性飞行器系统,利用滑模控制方法,再结合神经网络,设计了神经网络固定时间控制器。固定时间滑模面的设计保证系统状态的固定时间收敛。本发明在系统存在外界干扰,转动惯量不确定,执行器饱和和故障的情况下,实现系统姿态跟踪误差和角速度误差的固定时间一致最终有界的控制方法。
本发明的有益效果为:在系统存在外界干扰,转动惯量不确定,执行器饱和和故障的情况下,实现系统姿态跟踪误差和角速度误差的固定时间一致最终有界,并且收敛时间与系统的初始状态无关。
附图说明
图1为本发明的刚性飞行器姿态跟踪误差示意图;
图2为本发明的刚性飞行器角速度误差示意图;
图3为本发明的刚性飞行器滑模面示意图;
图4为本发明的刚性飞行器控制力矩示意图;
图5为本发明的刚性飞行器参数估计示意图;
图6为本发明的控制流程示意图。
具体实施方式
下面结合附图对本发明做进一步说明。
参照图1-图6,一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法,所述控制方法包括以下步骤:
步骤1,建立刚性飞行器的运动学和动力学模型,初始化系统状态以及控制参数,过程如下:
1.1刚性飞行器系统的运动学方程为:
Figure BDA0002206812530000081
Figure BDA0002206812530000082
其中qv=[q1,q2,q3]T和q4分别为单位四元数的矢量部分和标量部分且满足
Figure BDA0002206812530000083
q1,q2,q3分别为映射在空间直角坐标系x,y,z轴上的值;
Figure BDA0002206812530000084
分别是qv和q4的导数;Ω∈R3是刚性飞行器的角速度;I3是R3×3单位矩阵;
Figure BDA0002206812530000085
表示为:
Figure BDA0002206812530000091
1.2刚性飞行器系统的动力学方程为:
Figure BDA0002206812530000092
其中J∈R3×3是刚性飞行器的转动惯性矩阵;
Figure BDA0002206812530000093
是刚性飞行器的角加速度;u∈R3和d∈R3是控制力矩和外部扰动;D=diag(D1,D2,D3)∈R3×3是3×3对称对角的执行器效率矩阵,满足0<Di(t)≤1,i=1,2,3;sat(u)=[sat(u1),sat(u2),sat(u3)]T为执行器产生的实际控制力矩,sat(ui)为带有饱和特性的执行器,表示为sat(ui)=sgn(ui)min{umi,|ui|},umi为最大提供的控制力矩,sgn(ui)为符号函数,min{umi,|ui|}为两者的最小值;为了更方便的表示控制约束,sat(u)表示为sat(u)=g(u)+ds(u),g(u)=[g1(u1),g2(u2),g3(u3)]T,gi(ui)为双曲正切函数
Figure BDA0002206812530000094
ds(u)=[ds1(u1),ds2(u2),ds3(u3)]T为近似误差矢量;根据中值定理,gi(ui)=miui,0<mi≤1;定义H=DM=diag(δ1m12m23m3)∈R3×3为3×3对称对角矩阵,M=diag(m1,m2,m3)∈R3×3为3×3对称对角矩阵;Dsat(u)重新表示为:Dsat(u)=Hu+Dds(u),满足0<h0≤Dimi≤1,i=1,2,3,h0为未知正常数;Ω×表示为:
Figure BDA0002206812530000095
1.3刚性飞行器系统期望的运动学方程为:
Figure BDA0002206812530000096
Figure BDA0002206812530000101
其中qdv=[qd1,qd2,qd3]T和qd4分别为期望的单位四元数的矢量部分和标量部分且满足
Figure BDA0002206812530000102
Ωd∈R3为期望的角速度;
Figure BDA0002206812530000103
分别为qdv,qd4的导数,
Figure BDA0002206812530000104
为qdv的转置;
Figure BDA0002206812530000105
表示为:
Figure BDA0002206812530000106
1.4由四元数描述的刚性飞行器相对姿态运动:
Figure BDA0002206812530000107
Figure BDA0002206812530000108
Ωe=Ω-CΩd (12)
其中ev=[e1,e2,e3]T和e4分别为姿态跟踪误差的矢量部分和标量部分;Ωe=[Ωe1e2e3]T∈R3为角速度误差;
Figure BDA0002206812530000109
为相应的方向余弦矩阵并且满足||C||=1和
Figure BDA00022068125300001010
Figure BDA00022068125300001011
为C的导数;
根据式(1)-(12),刚性飞行器姿态跟踪误差动力学和运动学方程为:
Figure BDA00022068125300001012
Figure BDA00022068125300001013
Figure BDA00022068125300001014
其中
Figure BDA00022068125300001015
Figure BDA00022068125300001016
分别为ev和e4的导数;
Figure BDA00022068125300001017
为ev的转置;
Figure BDA00022068125300001018
Figure BDA00022068125300001019
分别为Ωd和Ωe的导数;(Ωe+CΩd)×与Ω×等价;
Figure BDA00022068125300001020
Figure BDA00022068125300001021
分别表示为:
Figure BDA0002206812530000111
Figure BDA0002206812530000112
1.5转动惯性矩阵J满足J=J0+ΔJ,其中J0和ΔJ分别表示J的标称部分和不确定部分,则式(15)重新写成:
Figure BDA0002206812530000113
进一步得到:
Figure BDA0002206812530000114
1.6对式(13)进行微分,得到:
Figure BDA0002206812530000115
其中
Figure BDA0002206812530000116
为ev的二阶导数;
步骤2,针对外部扰动,转动惯量不确定,执行器饱和和故障的刚性飞行器系统,设计所需的滑模面,过程如下:
选择固定时间滑模面为:
Figure BDA0002206812530000117
其中
Figure BDA0002206812530000118
Figure BDA0002206812530000121
和sgn(ei)均为符号函数,λ1>0,λ2>0,a2>1,
Figure BDA0002206812530000122
Figure BDA0002206812530000123
为ei的导数,i=1,2,3;
定义S=[S1,S2,S3]T,对S求导,得到:
Figure BDA0002206812530000124
将式(20)代入(22),得到:
Figure BDA0002206812530000125
其中
Figure BDA0002206812530000126
Figure BDA0002206812530000127
为总确定;L=[L1,L2,L3]T,
Figure BDA0002206812530000128
Figure BDA0002206812530000129
Γ=diag(Γ123)∈R3×3为3×3对称的对角矩阵;i=1,2,3;
步骤3,设计神经网络固定时间控制器,其过程如下:
3.1定义神经网络为:
Gi(Xi)=Wi *TΦ(Xi)+εi (24)
其中
Figure BDA00022068125300001210
为输入矢量,Φi(Xi)∈R4为神经网络基函数,Wi *∈R4为理想的权值矢量,定义为:
Figure BDA00022068125300001211
其中Wi∈R4为权值矢量,εi为近似误差,满足|εi|≤εN,i=1,2,3,εN为很小的正常数;arg min{·}为Wi *取其最小值所有的集合;
3.2考虑固定时间控制器被设计为:
Figure BDA00022068125300001212
其中
Figure BDA00022068125300001213
为3×3对称的对角矩阵,
Figure BDA00022068125300001214
Figure BDA0002206812530000131
为Θi的估计值Φ(X)=[Φ(X1),Φ(X2),Φ(X3)]T;K1=diag(k11,k12,k13)∈R3×3为3×3对称的对角矩阵;K2=diag(k21,k22,k23)∈R3×3为3×3对称的对角矩阵;K3=diag(k31,k32,k33)∈R3×3为对称的对角矩阵;k11,k12,k13,k21,k22,k23,k31,k32,k33为正常数;0<r1<1,r2>1;||·||表示值的二范数;
3.2设计更新律为:
Figure BDA0002206812530000132
其中γi>0,pi>0,i=1,2,3,
Figure BDA0002206812530000133
Figure BDA0002206812530000134
的导数,Φ(Xi)选择为以下的sigmoid函数:
Figure BDA0002206812530000135
其中l1,l2,l3和l4为近似参数,Φ(Xi)满足0<Φ(Xi)<Φ0,并且
Figure BDA0002206812530000136
步骤4,固定时间稳定性证明,其过程如下:
4.1证明刚性飞行器系统所有信号都是一致最终有界,设计李雅普诺夫函数为如下形式:
Figure BDA0002206812530000137
其中
Figure BDA0002206812530000138
ST是S的转置;
Figure BDA0002206812530000139
Figure BDA00022068125300001310
的转置;
对式(29)进行微分,得到:
Figure BDA00022068125300001311
其中
Figure BDA00022068125300001312
min{·}表示最小值;
因此,刚性飞行器系统所有信号都是一致最终有界的;
4.2证明固定时间收敛,设计李雅普诺夫函数为如下形式:
Figure BDA00022068125300001313
对式(31)进行微分,得到:
Figure BDA0002206812530000141
其中
Figure BDA0002206812530000142
i=1,2,3;υ2为一个大于零的上界值;
基于以上分析,刚性飞行器系统的姿态跟踪误差和角速度误差在固定时间一致最终有界。
为验证所提方法的有效性,本方法针对飞行器系统进行仿真验证。系统初始化参数设置如下:
系统的初始值:q(0)=[0.3,-0.2,-0.3,0.8832]T,Ω(0)=[1,0,-1]T弧度/秒;qd(0)=[0,0,0,1]T;期望角速度
Figure BDA0002206812530000143
弧度/秒;转动惯性矩阵的标称部分J0=[40,1.2,0.9;1.2,17,1.4;0.9,1.4,15]千克*平方米,惯性矩阵的不确定部ΔJ=diag[sin(0.1t),2sin(0.2t),3sin(0.3t)];外部扰动d(t)=[0.2sin(0.1t),0.3sin(0.2t),0.5sin(0.2t)]T牛*米;滑模面的参数如下:λ1=1,λ2=1,a1=1.5,a2=1.5;控制器的参数如下:
Figure BDA0002206812530000144
K1=K2=K3=I3;更新律参数如下:ηi=2,εi=0.1,i=1,2,3,
Figure BDA0002206812530000145
sigmoid函数的参数选择如下:l1=2,l2=8,l3=10,l4=-0.5。最大的控制力矩umi=10牛*米,执行器效率值选择为:
Figure BDA0002206812530000146
刚性飞行器的姿态跟踪误差和角速度误差的响应示意图分别如图1和图2所示,可以看出跟踪姿态误差和角速度误差都能在5秒左右收敛到平衡点的一个零域内;刚性飞行器的滑模面响应示意图如图3所示,可以看出滑模面能在3秒左右收敛到平衡点的一个零域内;刚性飞行器的控制力矩如图4所示,可以看出控制力矩限幅在10牛*米内;参数估计响应示意图分别如图5所示。
因此,本发明在系统存在外界干扰,转动惯量不确定,执行器饱和和故障的情况下,实现系统姿态跟踪误差和角速度误差的固定时间一致最终有界,并且收敛时间与系统的初始状态无关。
以上阐述的是本发明给出的一个实施例表现出的优良优化效果,显然本发明不只是限于上述实施例,在不偏离本发明基本精神及不超出本发明实质内容所涉及范围的前提下对其可作种种变形加以实施。

Claims (1)

1.一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法,其特征在于:所述方法包括以下步骤:
步骤1,建立刚性飞行器的运动学和动力学模型,初始化系统状态以及控制参数,过程如下:
1.1刚性飞行器系统的运动学方程为:
Figure FDA0003645319090000011
Figure FDA0003645319090000012
其中qv=[q1,q2,q3]T和q4分别为单位四元数的矢量部分和标量部分且满足
Figure FDA0003645319090000013
q1,q2,q3分别为映射在空间直角坐标系x,y,z轴上的值;
Figure FDA0003645319090000014
分别是qv和q4的导数;Ω∈R3是刚性飞行器的角速度;I3是R3×3单位矩阵;
Figure FDA0003645319090000015
表示为:
Figure FDA0003645319090000016
1.2刚性飞行器系统的动力学方程为:
Figure FDA0003645319090000017
其中J∈R3×3是刚性飞行器的转动惯性矩阵;
Figure FDA0003645319090000018
是刚性飞行器的角加速度;u∈R3和d∈R3是控制力矩和外部扰动;D=diag(D1,D2,D3)∈R3×3是3×3对称对角的执行器效率矩阵,满足0<Di≤1,i=1,2,3;sat(u)=[sat(u1),sat(u2),sat(u3)]T为执行器产生的实际控制力矩,sat(ui)为带有饱和特性的执行器,表示为sat(ui)=sgn(ui)min{umi,|ui|},umi为最大提供的控制力矩,sgn(ui)为符号函数,min{umi,|ui|}为两者的最小值;为了表示控制约束,sat(u)表示为sat(u)=g(u)+ds(u),g(u)=[g1(u1),g2(u2),g3(u3)]T,gi(ui)为双曲正切函数
Figure FDA0003645319090000021
ds(u)=[ds1(u1),ds2(u2),ds3(u3)]T为近似误差矢量;根据中值定理,gi(ui)=miui,0<mi≤1;定义H=DM=diag(D1m1,D2m2,D3m3)∈R3×3为3×3对称对角矩阵,M=diag(m1,m2,m3)∈R3 ×3为3×3对称对角矩阵;Dsat(u)重新表示为:Dsat(u)=Hu+Dds(u),满足0<h0≤Dimi≤1,i=1,2,3,h0为未知正常数;Ω×表示为:
Figure FDA0003645319090000022
1.3刚性飞行器系统期望的运动学方程为:
Figure FDA0003645319090000023
Figure FDA0003645319090000024
其中qdv=[qd1,qd2,qd3]T和qd4分别为期望的单位四元数的矢量部分和标量部分且满足
Figure FDA0003645319090000025
Ωd∈R3为期望的角速度;
Figure FDA0003645319090000026
分别为qdv,qd4的导数,
Figure FDA0003645319090000027
为qdv的转置;
Figure FDA0003645319090000028
表示为:
Figure FDA0003645319090000029
1.4由四元数描述的刚性飞行器相对姿态运动:
Figure FDA00036453190900000210
Figure FDA00036453190900000211
Ωe=Ω-CΩd (12)
其中ev=[e1,e2,e3]T和e4分别为姿态跟踪误差的矢量部分和标量部分;Ωe=[Ωe1e2e3]T∈R3为角速度误差;
Figure FDA0003645319090000031
为相应的方向余弦矩阵并且满足||C||=1和
Figure FDA0003645319090000032
Figure FDA0003645319090000033
为C的导数;
根据式(1)-(12),刚性飞行器姿态跟踪误差动力学和运动学方程为:
Figure FDA0003645319090000034
Figure FDA0003645319090000035
Figure FDA0003645319090000036
其中
Figure FDA0003645319090000037
Figure FDA0003645319090000038
分别为ev和e4的导数;
Figure FDA0003645319090000039
为ev的转置;
Figure FDA00036453190900000310
Figure FDA00036453190900000311
分别为Ωd和Ωe的导数;(Ωe+CΩd)×与Ω×等价;
Figure FDA00036453190900000312
Figure FDA00036453190900000313
分别表示为:
Figure FDA00036453190900000314
Figure FDA00036453190900000315
1.5转动惯性矩阵J满足J=J0+ΔJ,其中J0和ΔJ分别表示J的标称部分和不确定部分,则式(15)重新写成:
Figure FDA00036453190900000316
进一步得到:
Figure FDA00036453190900000317
1.6对式(13)进行微分,得到:
Figure FDA0003645319090000041
其中
Figure FDA0003645319090000042
为ev的二阶导数;
步骤2,针对外部扰动,转动惯量不确定,执行器饱和故障的刚性飞行器系统,设计所需的滑模面,过程如下:
选择固定时间滑模面为:
Figure FDA0003645319090000043
其中
Figure FDA0003645319090000044
Figure FDA0003645319090000045
Figure FDA0003645319090000046
和sgn(ei)均为符号函数,λ1>0,λ2>0,a2>1,
Figure FDA0003645319090000047
Figure FDA0003645319090000048
为ei的导数,i=1,2,3;
定义S=[S1,S2,S3]T,对S求导,得到:
Figure FDA0003645319090000049
将式(20)代入(22),得到:
Figure FDA00036453190900000410
其中
Figure FDA00036453190900000411
Figure FDA00036453190900000412
为总确定;
Figure FDA00036453190900000413
Figure FDA00036453190900000414
Γ=diag(Γ123)∈R3×3为3×3对称的对角矩阵;i=1,2,3;
步骤3,设计神经网络固定时间控制器,其过程如下:
3.1定义神经网络为:
Gi(Xi)=Wi *TΦ(Xi)+εi (24)
其中
Figure FDA0003645319090000051
为输入矢量,Φ(Xi)∈R4为神经网络基函数,Wi *∈R4为理想的权值矢量,定义为:
Figure FDA0003645319090000052
其中Wi∈R4为权值矢量,εi为近似误差,满足|εi|≤εN,i=1,2,3,εN为很小的正常数;argmin{·}为Wi *取其最小值所有的集合;
3.2考虑固定时间控制器被设计为:
Figure FDA0003645319090000053
其中
Figure FDA0003645319090000054
为3×3对称的对角矩阵,
Figure FDA0003645319090000055
Figure FDA0003645319090000056
为Θi的估计值;Φ(X)=[Φ(X1),Φ(X2),Φ(X3)]T;K1=diag(k11,k12,k13)∈R3×3为3×3对称的对角矩阵;K2=diag(k21,k22,k23)∈R3×3为3×3对称的对角矩阵;K3=diag(k31,k32,k33)∈R3×3为对称的对角矩阵;k11,k12,k13,k21,k22,k23,k31,k32,k33为正常数;0<r1<1,r2>1;||·||表示值的二范数;
3.3设计更新律为:
Figure FDA0003645319090000057
其中γi>0,pi>0,i=1,2,3,
Figure FDA0003645319090000058
Figure FDA0003645319090000059
的导数,Φ(Xi)选择为以下的sigmoid函数:
Figure FDA00036453190900000510
其中l1,l2,l3和l4为近似参数,Φ(Xi)满足0<Φ(Xi)<Φ0,并且
Figure FDA00036453190900000511
步骤4,固定时间稳定性证明,其过程如下:
4.1证明刚性飞行器系统所有信号都是一致最终有界,设计李雅普诺夫函数为如下形式:
Figure FDA0003645319090000061
其中
Figure FDA0003645319090000062
ST是S的转置;
Figure FDA0003645319090000063
Figure FDA0003645319090000064
的转置;
对式(29)进行微分,得到:
Figure FDA0003645319090000065
其中
Figure FDA0003645319090000066
min{·}表示最小值;
因此,刚性飞行器系统所有信号都是一致最终有界的;
4.2证明固定时间收敛,设计李雅普诺夫函数为如下形式:
Figure FDA0003645319090000067
对式(31)进行微分,得到:
Figure FDA0003645319090000068
其中
Figure FDA0003645319090000069
Figure FDA00036453190900000610
υ2为一个大于零的上界值;
基于以上分析,刚性飞行器系统的姿态跟踪误差和角速度误差在固定时间一致最终有界。
CN201910884253.2A 2018-09-25 2019-09-19 一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法 Active CN110488603B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2018111145464 2018-09-25
CN201811114546.4A CN109143846A (zh) 2018-09-25 2018-09-25 一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法

Publications (2)

Publication Number Publication Date
CN110488603A CN110488603A (zh) 2019-11-22
CN110488603B true CN110488603B (zh) 2022-10-28

Family

ID=64823432

Family Applications (2)

Application Number Title Priority Date Filing Date
CN201811114546.4A Withdrawn CN109143846A (zh) 2018-09-25 2018-09-25 一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法
CN201910884253.2A Active CN110488603B (zh) 2018-09-25 2019-09-19 一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法

Family Applications Before (1)

Application Number Title Priority Date Filing Date
CN201811114546.4A Withdrawn CN109143846A (zh) 2018-09-25 2018-09-25 一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法

Country Status (1)

Country Link
CN (2) CN109143846A (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110989337B (zh) * 2019-12-10 2020-12-01 北京理工大学 考虑舵机失效及气动扰动的飞行器旋转控制系统及方法
CN110989338B (zh) * 2019-12-10 2020-12-01 北京理工大学 考虑气动非线性的飞行器旋转稳定控制系统及方法
CN111498148A (zh) * 2020-04-23 2020-08-07 北京航空航天大学 基于fdnn的智能航天器控制器和控制方法
CN111781942B (zh) * 2020-06-23 2021-10-12 南京航空航天大学 一种基于自构造模糊神经网络的容错飞行控制方法
CN113985905A (zh) * 2021-09-29 2022-01-28 西安工业大学 一种考虑执行器饱和及故障的无退绕姿态控制方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103760906B (zh) * 2014-01-29 2016-06-01 天津大学 神经网络与非线性连续无人直升机姿态控制方法
CN104267596B (zh) * 2014-08-15 2017-06-13 浙江工业大学 一种小车倒立摆系统的有限时间解耦控制方法
CN104698847B (zh) * 2015-02-10 2017-04-12 浙江工业大学 一种转台伺服系统的非奇异终端滑模指定性能控制方法
CN106843254B (zh) * 2017-03-08 2019-08-09 北京航天自动控制研究所 一种实时主动重构容错控制方法
CN107662208B (zh) * 2017-08-24 2020-07-31 浙江工业大学 一种基于神经网络的柔性关节机械臂有限时间自适应反步控制方法

Also Published As

Publication number Publication date
CN110488603A (zh) 2019-11-22
CN109143846A (zh) 2019-01-04

Similar Documents

Publication Publication Date Title
CN110488603B (zh) 一种考虑执行器受限问题的刚性飞行器自适应神经网络跟踪控制方法
CN110543183B (zh) 一种考虑执行器受限问题的刚体飞行器固定时间姿态跟踪控制方法
CN110543184B (zh) 一种刚性飞行器的固定时间神经网络控制方法
CN109188910B (zh) 一种刚性飞行器的自适应神经网络容错跟踪控制方法
CN107703952B (zh) 一种刚性飞行器的非奇异固定时间自适应姿态控制方法
CN109062240B (zh) 一种基于神经网络估计的刚性飞行器固定时间自适应姿态跟踪控制方法
Zou et al. Fixed-time attitude tracking control for rigid spacecraft without angular velocity measurements
CN108490783B (zh) 基于增强型双幂次趋近律和快速终端滑模面的刚性航天飞行器有限时间自适应容错控制方法
CN107577145B (zh) 编队飞行航天器反步滑模控制方法
CN110471438B (zh) 一种刚性飞行器的固定时间自适应姿态跟踪控制方法
CN107422741B (zh) 基于学习的保预设性能集群飞行分布式姿态追踪控制方法
CN110488854B (zh) 一种基于神经网络估计的刚性飞行器固定时间姿态跟踪控制方法
Yang et al. Formation control of multiple underwater vehicles subject to communication faults and uncertainties
CN107703753B (zh) 一种空间机械臂的容错控制方法
Zhang et al. Decentralized coordinated control for multiple spacecraft formation maneuvers
CN113589689B (zh) 一种基于多参数自适应神经网络的滑模控制器设计方法
CN108958043B (zh) 一种刚性飞行器的固定时间自适应姿态容错控制方法
CN110501911A (zh) 一种考虑执行器受限问题的刚性飞行器自适应固定时间姿态跟踪控制方法
CN110488855B (zh) 一种基于神经网络估计的刚性飞行器自适应固定时间姿态容错控制方法
CN113220007B (zh) 执行机构故障的挠性航天器有限时间姿态协同控制方法
CN110515389B (zh) 一种考虑执行器受限问题的刚性飞行器自适应固定时间姿态镇定方法
CN113568423A (zh) 一种考虑电机故障的四旋翼无人机智能容错控制方法
CN111427264B (zh) 一种复杂遥操作技术的神经自适应固定时间控制方法
CN116923730A (zh) 一种具有自调节预设性能约束的航天器姿态主动容错控制方法
CN108762065B (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