一种刚性飞行器的非奇异固定时间神经网络控制方法与流程

本发明涉及一种刚性飞行器的非奇异固定时间神经网络控制方法，特别是存在外部干扰和转动惯量矩阵不确定的刚性飞行器姿态跟踪控制方法。

背景技术：

刚性飞行器姿态控制系统在刚性飞行器的健康，可靠的运动中扮演着重要的角色。在复杂的航天环境中，刚性飞行器姿态控制系统会受到各种外部干扰以及转动惯量矩阵不确定的影响。为了有效维持系统的性能，需要使其对外部干扰以及转动惯量矩阵不确定具有较强的鲁棒性。滑模变结构控制作为一种典型的非线性控制方法能够有效改善刚性飞行器的稳定性和操纵性，并且具有较强的鲁棒性，从而提高执行任务的能力。因此，研究刚性飞行器姿态控制系统的滑模变结构控制方法具有十分重要的意义。

滑模控制在解决系统不确定性和外部扰动方面被认为是一个有效的鲁棒控制方法。滑模控制方法具有算法简单、响应速度快、对外界噪声干扰和参数摄动鲁棒性强等优点。终端滑模控制是一种可以实现有限时间稳定性的传统滑模控制的改进方案。然而，现存的有限时间技术估计收敛时间需要知道系统的初始信息，这对于设计者是很难知道的。近年来，固定时间技术得到了广泛的应用，固定时间控制方法与现存的有限时间控制方法相比，具有无需知道系统的初始信息，也能保守估计系统的收敛时间的优越性。

神经网络是线性参数化近似方法的中一种，可以被任意的其他近似方法取代，比如rbf神经网络，模糊逻辑系统等等。利用神经网络逼近不确定的性质，有效的结合固定时间滑模控制技术，减少外部干扰及系统参数不确定性对系统控制性能的影响，实现刚性飞行器姿态的固定时间控制。

技术实现要素：

为了克服现有的刚性飞行器姿态控制系统存在的未知非线性问题，本发明提供一种刚性飞行器的非奇异固定时间神经网络控制方法，并且在系统存在外部干扰和转动惯量不确定的情况下，实现系统姿态跟踪误差和角速度误差的固定时间一致最终有界的控制方法。

为了解决上述技术问题提出的技术方案如下：

一种刚性飞行器的非奇异固定时间神经网络控制方法，包括以下步骤：

步骤1，建立刚性飞行器的运动学和动力学模型，初始化系统状态以及控制参数，过程如下：

1.1刚性飞行器系统的运动学方程为：

其中qv＝[q1,q2,q3]^t和q4分别为单位四元数的矢量部分和标量部分且满足分别为映射在空间直角坐标系x,y,z轴上的值；分别是qv和q4的导数；ω∈r³是刚性飞行器的角速度；i3是r^3×3单位矩阵；表示为：

1.2刚性飞行器系统的动力学方程为：

其中j∈r^3×3是飞行器的转动惯性矩阵；是飞行器的角加速度；u∈r³和d∈r³是控制力矩和外部扰动；ω^×表示为：

1.3刚性飞行器系统期望的运动学方程为：

其中qdv＝[qd1,qd2,qd3]^t和qd4分别为期望的单位四元数的矢量部分和标量部分且满足ωd∈r³为期望的角速度；分别为qdv,qd4的导数，为qdv的转置；表示为：

1.4由四元数描述的刚性飞行器相对姿态运动：

ωe＝ω-cωd(11)

其中ev＝[e1,e2,e3]^t和e4分别为姿态跟踪误差的矢量部分和标量部分；ωe＝[ωe1,ωe2,ωe3]^t∈r³为角速度误差；为相应的方向余弦矩阵并且满足||c||＝1和为c的导数；

根据式(1)-(11)，刚性飞行器姿态跟踪误差动力学和运动学方程为：

其中和分别为ev和e4的导数；为ev的转置；和分别为ωd和ωe的导数；(ωe+cωd)^×与ω^×等价；和分别表示为：

1.5转动惯性矩阵j满足j＝j0+δj，其中j0和δj分别表示j的标称部分和不确定部分，则式(14)重新写成：

进一步得到：

1.6对式(12)进行微分，得到：

其中为总不确定的集合；为ev的二阶导数；

步骤2，针对外部扰动和转动惯量不确定的刚性飞行器系统，设计所需的滑模面，过程如下：

选择非奇异固定时间滑模面为：

其中，λ1和λ2为正常数；m1,n1,p1,r1为正奇数，满足m1＞n1,p1＜r1＜2p1，和均为符号函数；i＝1,2,3；

步骤3，设计固定时间神经网络控制器，过程如下：

3.1定义神经网络为：

其中为输入矢量，φi(xi)∈r⁴为神经网络基函数，为的转置，为理想的权值矢量，定义为：

其中wi∈r⁴为权值矢量，εi为近似误差，满足|εi|≤εn,i＝1,2,3，εn为很小的正常数；为取其最小值所有的集合；

3.2考虑固定时间控制器被设计为：

其中为3×3对称对角矩阵；为wi的估计值；φ(x)＝[φ(x1),φ(x2),φ(x3)]^t，γ＝diag(γ1,γ2,γ3)，k1＝diag(k11,k12,k13),k2＝diag(k21,k22,k23),k3＝diag(k31,k32,k33)，和均为3×3对称对角矩阵；为的导数，s＝[s1,s2,s3]^t；k11,k12,k13,k21,k22,k23,k31,k32,k33为正常数；和均为符号函数；sgn(s1),sgn(s2),sgn(s3)均为符号函数；

3.3设计更新律为：

其中γi＞0,τi＞0，为的估计，i＝1,2,3；φ(xi)选择为以下的sigmoid函数：

其中l1,l2,l3和l4为近似参数，φ(xi)满足0＜φ(xi)＜φ0，并且为两者中的最大值；

步骤4，固定时间稳定性证明，过程如下：

4.1证明刚性飞行器系统所有信号都是一致最终有界，设计李雅普诺夫函数为如下形式：

其中s^t是s的转置；为的转置；

对式(26)进行求导，得到：

其中为的二范数；为两者中的最小值；

则判定刚性飞行器系统所有信号都是一致最终有界的；

4.2证明固定时间收敛，设计李雅普诺夫函数为如下形式：

对式(28)进行求导，得到：

其中均为取其最小值；υ2为一个大于零的上界值；

基于以上分析，刚性飞行器系统的姿态跟踪误差和角速度误差在固定时间一致最终有界。

本发明在外界干扰和转动惯量不确定的因素下，运用一种刚性飞行器的非奇异固定时间神经网络控制方法，实现系统稳定跟踪，保证系统的姿态跟踪误差和角速度误差固定时间一致最终有界。本发明的技术构思为：针对含外部干扰和转动惯量不确定的刚性飞行器系统，利用滑模控制方法，再结合神经网络，设计了非奇异固定时间控制器。非奇异固定时间滑模面的设计不仅保证系统状态的固定时间收敛，而且解决了奇异值问题。本发明在系统存在外界干扰和转动惯量不确定的情况下，实现系统姿态跟踪误差和角速度误差的固定时间一致最终有界的控制方法。

本发明的有益效果为：所设计的固定时间滑模面有效的解决了奇异值问题；在系统存在外界干扰和转动惯量不确定的情况下，实现系统姿态跟踪误差和角速度误差的固定时间一致最终有界，并且收敛时间与系统的初始状态无关。

附图说明

图1为本发明的刚性飞行器姿态跟踪误差示意图；

图2为本发明的刚性飞行器角速度误差示意图；

图3为本发明的刚性飞行器滑模面示意图；

图4为本发明的刚性飞行器控制力矩示意图；

图5为本发明的刚性飞行器参数估计示意图；

图6为本发明的控制流程示意图。

具体实施方式

下面结合附图对本发明做进一步说明。

参照图1-图6，一种刚性飞行器的非奇异固定时间神经网络控制方法，所述控制方法包括以下步骤：

步骤1，建立刚性飞行器的运动学和动力学模型，初始化系统状态以及控制参数，过程如下：

1.5刚性飞行器系统的运动学方程为：

其中qv＝[q1,q2,q3]^t和q4分别为单位四元数的矢量部分和标量部分且满足分别为映射在空间直角坐标系x,y,z轴上的值；分别是qv和q4的导数；ω∈r³是刚性飞行器的角速度；i3是r^3×3单位矩阵；表示为：

1.6刚性飞行器系统的动力学方程为：

其中j∈r^3×3是飞行器的转动惯性矩阵；是飞行器的角加速度；u∈r³和d∈r³是控制力矩和外部扰动；ω^×表示为：

1.7刚性飞行器系统期望的运动学方程为：

其中qdv＝[qd1,qd2,qd3]^t和qd4分别为期望的单位四元数的矢量部分和标量部分且满足ωd∈r³为期望的角速度；分别为qdv,qd4的导数，为qdv的转置；表示为：

1.8由四元数描述的刚性飞行器相对姿态运动：

ωe＝ω-cωd(11)

其中ev＝[e1,e2,e3]^t和e4分别为姿态跟踪误差的矢量部分和标量部分；ωe＝[ωe1,ωe2,ωe3]^t∈r³为角速度误差；为相应的方向余弦矩阵并且满足||c||＝1和为c的导数；

根据式(1)-(11)，刚性飞行器姿态跟踪误差动力学和运动学方程为：

其中和分别为ev和e4的导数；为ev的转置；和分别为ωd和ωe的导数；(ωe+cωd)^×与ω^×等价；和分别表示为：

1.5转动惯性矩阵j满足j＝j0+δj，其中j0和δj分别表示j的标称部分和不确定部分，则式(14)重新写成：

进一步得到：

1.6对式(12)进行微分，得到：

其中为总不确定的集合；为ev的二阶导数；

步骤2，针对外部扰动和转动惯量不确定的刚性飞行器系统，设计所需的滑模面，过程如下：

选择非奇异固定时间滑模面为：

其中，λ1和λ2为正常数；m1,n1,p1,r1为正奇数，满足m1＞n1,p1＜r1＜2p1，sgn(ei)和均为符号函数；i＝1,2,3；

步骤3，设计固定时间神经网络控制器，过程如下：

3.1定义神经网络为：

其中为输入矢量，φi(xi)∈r⁴为神经网络基函数，为的转置，为理想的权值矢量，定义为：

其中wi∈r⁴为权值矢量，εi为近似误差，满足|εi|≤εn,i＝1,2,3，εn为很小的正常数；为取其最小值所有的集合；

3.2考虑固定时间控制器被设计为：

其中为3×3对称对角矩阵；为wi的估计值；φ(x)＝[φ(x1),φ(x2),φ(x3)]^t，γ＝diag(γ1,γ2,γ3)，k1＝diag(k11,k12,k13),k2＝diag(k21,k22,k23),k3＝diag(k31,k32,k33)，和均为3×3对称对角矩阵；为的导数，s＝[s1,s2,s3]^t；k11,k12,k13,k21,k22,k23,k31,k32,k33为正常数；和均为符号函数；sgn(s1),sgn(s2),sgn(s3)均为符号函数；

3.3设计更新律为：

其中γi＞0,τi＞0，为的估计，i＝1,2,3；φ(xi)选择为以下的sigmoid函数：

其中l1,l2,l3和l4为近似参数，φ(xi)满足0＜φ(xi)＜φ0，并且为两者中的最大值；

步骤4，固定时间稳定性证明，过程如下：

4.1证明刚性飞行器系统所有信号都是一致最终有界，设计李雅普诺夫函数为如下形式：

其中s^t是s的转置；为的转置；

对式(26)进行求导，得到：

其中为的二范数；为两者中的最小值；

则判定刚性飞行器系统所有信号都是一致最终有界的；

4.2证明固定时间收敛，设计李雅普诺夫函数为如下形式：

对式(28)进行求导，得到：

其中均为取其最小值；υ2为一个大于零的上界值；

基于以上分析，刚性飞行器系统的姿态跟踪误差和角速度误差在固定时间一致最终有界。

为验证所提方法的有效性，本方法针对飞行器系统进行仿真验证。系统初始化参数设置如下：

系统的初始值：q(0)＝[0.3,-0.2,-0.3,0.8832]^t,ω(0)＝[1,0,-1]^t弧度/秒qd(0)＝[0,0,0,1]^t；期望角速度弧度/秒；转动惯性矩阵的标称部分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,m1＝9,n1＝5,p1＝3,r1＝5；控制器的参数如下：k1＝k2＝k3＝i3；更新律参数如下：sigmoid函数的参数选择如下：l1＝2,l2＝8,l3＝10,l4＝-0.5。

刚性飞行器的姿态跟踪误差和角速度误差的响应示意图分别如图1和图2所示，可以看出姿态跟踪误差和角速度误差都能在3.5秒左右收敛到平衡点的一个零域内；刚性飞行器的滑模面响应示意图如图3所示，可以看出滑模面能在2秒左右收敛到平衡点的一个零域内；刚性飞行器的控制力矩和参数估计响应示意图分别如图4和图5所示。

因此，本发明设计的固定时间滑模面有效的解决了奇异值问题；在系统存在外界干扰和转动惯量不确定的情况下，实现系统的姿态跟踪误差和角速度误差在固定时间一致最终有界，并且收敛时间与系统的初始状态无关。

以上阐述的是本发明给出的一个实施例表现出的优良优化效果，显然本发明不只是限于上述实施例，在不偏离本发明基本精神及不超出本发明实质内容所涉及范围的前提下对其可作种种变形加以实施。