CN104898418B - 一种挠性卫星自适应神经网络滑模姿态控制方法 - Google Patents

一种挠性卫星自适应神经网络滑模姿态控制方法 Download PDF

Info

Publication number
CN104898418B
CN104898418B CN201510172378.4A CN201510172378A CN104898418B CN 104898418 B CN104898418 B CN 104898418B CN 201510172378 A CN201510172378 A CN 201510172378A CN 104898418 B CN104898418 B CN 104898418B
Authority
CN
China
Prior art keywords
mrow
mover
msub
msubsup
msup
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
CN201510172378.4A
Other languages
English (en)
Other versions
CN104898418A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of 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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN201510172378.4A priority Critical patent/CN104898418B/zh
Publication of CN104898418A publication Critical patent/CN104898418A/zh
Application granted granted Critical
Publication of CN104898418B publication Critical patent/CN104898418B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)

Abstract

一种挠性卫星自适应神经网络滑模姿态控制方法，本发明涉及一种挠性卫星自适应神经网络滑模姿态控制方法。本发明是要解决挠性卫星由于帆板模态振动和天线转动造成姿态波动，降低***稳定性的问题而提出的一种挠性卫星自适应神经网络滑模姿态控制方法。该方法是通过步骤一、建立挠性卫星姿态动力学模型；步骤二、得到简化后的挠性卫星姿态动力学方程；步骤三、根据简化后的挠性卫星姿态动力学方程，利用RBF神经网络设计滑模姿态控制器；步骤四、进一步采用RBF神经网络逼近符号函数η′sgn(s)，削弱抖振对滑模姿态控制器的影响；得到削弱抖振后的滑模姿态控制器等步骤实现的。本发明应用于挠性卫星姿态控制领域。

Description

一种挠性卫星自适应神经网络滑模姿态控制方法
技术领域
本发明涉及姿态控制方法，特别涉及一种挠性卫星自适应神经网络滑模姿态控制方法。
背景技术
随着卫星功能和种类的增多，卫星的结构变得十分复杂、尺寸也变得更大。星体上通常需要安装用以实现各种功能的挠性附件，其中以太阳能帆板、运动天线的应用最为普遍，这类带有挠性附件的卫星统称为挠性卫星。这些挠性附件的存在，使得卫星的姿态控制***具有非线性、参数不确定性等特点，这就对卫星的姿态控制提出了更高的要求。如何能够在抑制模态振动的同时使卫星姿态控制***具有稳定性、快速性等优良性能，已经成为了航天工程的研究热点问题。
国内外许多学者对这挠性卫星的姿态控制进行了广泛的研究，目前的挠性卫星的姿态控制方法的自抗扰控制方法在非线性稳定性证明上缺少足够的理论依据，此外控制器参数多，实际应用时参数整定十分不便，而且参数的来源主要是凭借经验进行试凑，浪费时间。而PID/PD方法中***的响应速度和超调量相互制约，难以同时兼顾，***的鲁棒性较差，面对不确定性的干扰时，姿态控制效果不理想。
发明内容
本发明的目的是为了解决挠性卫星由于帆板模态振动和天线转动造成的姿态波动，降低***稳定性的问题而提出的一种挠性卫星自适应神经网络滑模姿态控制方法。
上述的发明目的是通过以下技术方案实现的：
步骤一、采用混合坐标法建立挠性卫星姿态动力学模型；
步骤二、根据挠性卫星姿态动力学模型，忽略挠性卫星姿态动力学方程中与模态相关的高阶耦合项，考虑卫星作惯性定向飞行，同时采用小角度假设，得到简化后的挠性卫星姿态动力学方程；
步骤三、根据简化后的挠性卫星姿态动力学方程，利用RBF神经网络设计滑模姿态控制器；
步骤四、进一步采用RBF神经网络逼近符号函数η′sgn(s)，削弱抖振对滑模姿态控制器的影响；得到削弱抖振后的滑模姿态控制器；其中，sgn(s)＝[sgn(s1),sgn(s2),sgn(s3)]T，η′为待设计参数且η′>0，s1、s2和s3分别为X、Y和Z三通道的滑模函数；s＝[s1,s2,s3]T；即完成了一种挠性卫星自适应神经网络滑模姿态控制方法。
发明效果
本发明针对挠性卫星由于帆板的模态振动和天线的转动扰动造成的姿态波动问题，采用一种改进的自适应神经网络滑模姿态控制方法，提出了一种姿态控制器的设计方案；有效抑制帆板振动和天线扰动，完成姿态稳定任务，同时削弱控制量抖振，提高姿态控制精度与稳定度。以滑模变结构控制为基础，运用了RBF神经网络的逼近原理，提出了基于RBF神经网络的滑模姿态控制器的设计方法。并且进一步采用RBF神经网络对符号函数进行逼近，其目的在于削弱抖振的影响并且使***具有良好的控制特性。如图2～5所示的仿真结果表明，该方法可以提高***的姿态控制精度、稳定度、鲁棒性等指标，对抖振的有明显的削弱作用，适于实际工程应用。
因此本发明在考虑了挠性影响的情况下进行设计，适合工程应用；采用滑模变结构控制设计姿态控制器，***抗干扰能力强，鲁棒性好；无需前馈补偿，避免了前馈补偿难以有效实现的问题；采用RBF神经网络对符号函数进行逼近，削弱了抖振的影响；***的姿态控制精度和稳定度较高。
附图说明
图1为具体实施方式一提出的一种挠性卫星自适应神经网络滑模姿态控制方法流程图；
图2为具体实施方式一提出的逼近sign函数前的卫星姿态角曲线；
图3为具体实施方式一提出的逼近sign函数前的卫星姿态角速度曲线；
图4为具体实施方式一提出的逼近sign函数后的卫星姿态角曲线；
图5为具体实施方式一提出的逼近sign函数后的卫星姿态角速度曲线；
图6为实施例提出的各坐标系位置关系示意图；
图7为实施例提出的天线在卫星本体安装示意图；
图8为实施例提出RBF神经网络结构图。
具体实施方式
具体实施方式一：本实施方式的一种挠性卫星自适应神经网络滑模姿态控制方法，具体是按照以下步骤制备的：
步骤一、建立带有运动天线的卫星姿态模型；采用混合坐标法建立挠性卫星姿态动力学模型；
步骤二、根据挠性卫星姿态动力学模型，忽略挠性卫星姿态动力学方程中与模态相关的高阶耦合项，考虑卫星作惯性定向飞行，同时采用小角度假设，得到简化后的挠性卫星姿态动力学方程；
步骤三、根据简化后的挠性卫星姿态动力学方程，利用RBF神经网络设计滑模姿态控制器；
步骤四、进一步采用RBF神经网络逼近符号函数η′sgn(s)，削弱抖振对滑模姿态控制器的影响；得到削弱抖振后的滑模姿态控制器；其中，sgn(s)＝[sgn(s1),sgn(s2),sgn(s3)]T，η′为待设计参数且η′>0，s1、s2和s3分别为X、Y和Z三通道的滑模函数；s＝[s1,s2,s3]T如图1即完成了一种挠性卫星自适应神经网络滑模姿态控制方法。
本实施方式效果：
本发明针对挠性卫星由于帆板的模态振动和天线的转动扰动造成的姿态波动问题，采用一种改进的自适应神经网络滑模姿态控制方法，提出了一种姿态控制器的设计方案；有效抑制帆板振动和天线扰动，完成姿态稳定任务，同时削弱控制量抖振，提高姿态控制精度与稳定度。以滑模变结构控制为基础，参考了RBF神经网络的逼近原理，提出了基于RBF神经网络的滑模姿态控制器的设计方法。并且采用RBF神经网络对符号函数进行逼近，其目的在于削弱抖振的影响并且使***具有良好的控制特性。如图2～5所示的仿真结果表明，该方法可以提高***的姿态控制精度、稳定度、鲁棒性等指标，对抖振的有明显的削弱作用，适于实际工程应用。
因此本实施方式在考虑了挠性影响的情况下进行设计，适合工程应用；采用滑模变结构控制设计姿态控制器，***抗干扰能力强，鲁棒性好；无需前馈补偿，避免了前馈补偿难以有效实现的问题；采用RBF神经网络对符号函数进行逼近，削弱了抖振的影响；***的姿态控制精度和稳定度较高。
具体实施方式二：本实施方式与具体实施方式一不同的是：步骤一中建立带有运动天线的卫星姿态模型；采用混合坐标法建立挠性卫星姿态动力学模型具体为：
(1)、含有两块帆板和一根运动天线的姿态动力学方程有以下形式(方程(2-1)中，第一个为卫星本体转动方程，第二个为天线转动方程)：
其中，Is∈R3×3为星体转动惯量阵；ωs＝[ω1,ω2,ω3]T∈R3为本体系相对于惯性系且投影分解在本体系中的姿态角速度矢量；ω1、ω2和ω3分别为本体系相对于惯性系且投影分解在本体系中X、Y和Z轴的姿态角速度；ηk∈Rn为挠性模态坐标，n为模态阶数，k为附件编号，k＝1、2时表示本体Y轴正和负方向帆板，k＝3表示天线；Fsk∈R3×n为附件振动与星体转动耦合系数；Rsa∈R3×3为天线与星体转动耦合系数；u∈R3是由执行机构(飞轮、动量轮、推力器等)提供的星体三个通道控制力矩矢量；d∈R3为卫星所受的干扰力矩，包括环境干扰力矩和部件安装误差所引起的干扰力矩等；Ia∈R3×3为天线转动惯量阵；ωa∈R3，ωa＝[ωa1,ωa2,ωa3]T为天线相对于天线支撑臂坐标系的角速度矢量；ωa1、ωa2和ωa3分别为天线相对于天线支撑臂坐标系X、Y和Z轴的角速度；Fa∈R3×n为天线转动与天线臂振动耦合系数；符号表示如下的反对称矩阵，类似的，表示如下的反对称矩阵Ta∈R3是天线转动驱动控制力矩；角标a表示天线；角标s表示星体；
(2)、建立附件模态方程(附件模态方程为天线和两块帆板的振动方程)为：
其中，ξk和Ωk为n维对角阵，ξk表示附件的阻尼比；Ωk表示附件的模态频率，k＝1或2时表示本体Y轴正和负方向帆板，k＝3时表示天线；
(3)、含有两块帆板和一根运动天线的姿态动力学方程与附件模态方程组成挠性卫星姿态动力学模型。其它步骤及参数与具体实施方式一相同。
具体实施方式三：本实施方式与具体实施方式一或二不同的是：步骤二中根据挠性卫星姿态动力学模型，忽略挠性卫星姿态动力学方程中与模态相关的高阶耦合项，考虑卫星作惯性定向飞行，同时采用小角度假设，得到简化后的挠性卫星姿态动力学方程具体过程为：
(1)、忽略挠性卫星姿态动力学方程中与模态相关的高阶耦合项，将挠性卫星姿态动力学模型简化为：
式中，
D为干扰及不确定性的总和，且D有界；F为未知的非线性项；
(2)、考虑卫星惯性定向飞行，且采用小角度假设，近似的认为卫星角速度ωs近似等于姿态角速度，则简化后的挠性卫星姿态动力学方程具体为：
其它步骤及参数与具体实施方式一或二相同。
具体实施方式四：本实施方式与具体实施方式一至三之一不同的是：步骤三中根据简化后的挠性卫星姿态动力学方程，利用RBF神经网络设计滑模姿态控制器具体过程为：
(1)、取滑模函数s为：
c为对角阵，是待设计参数；e是卫星姿态误差；卫星的期望姿态角Θd为0，误差e＝Θd-Θ＝-Θ，Θ为卫星的实际姿态角；
(2)、联立式(2-5)和式(2-6)得：
(3)、设计姿态控制器的基本形式为：
(4)、采用RBF神经网络对未知的非线性项F进行逼近：
取hj(x)为的高斯基函数形式，则：
F＝W*Th(x)+ε (2-9)
式中，x＝[x1,x2,…,xr]T为是神经网络的r维输入向量，cj为RBF神经网络中第j个节点的中心向量，cj＝[cj1,cj2,…,cjr]T；bj>0为节点j的基宽值，b＝[b1,b2,…,bm]T；h(x)＝[h1(x),h2(x),…,hj(x),…hm(x)]T为m×1维神经网络径向基向量，m为隐层数；为m×3理想网络权值阵；ε＝[ε1,ε2,ε3]T为神经网络逼近误差，且εi有界，i取1、2、3表示卫星的X、Y和Z三通道；将未知的非线性项F作为神经网络的输出；
(5)、将公式(2-9)的网络输入取为则网络输出的估计值为：
其中，是权值的估计值；或分别为卫星的X、Y和Z三通道权值的估计值；
(6)、将网络输出的估计值引入姿态控制器u即公式(2-8)中：
(7)、将控制器代入式(2-7)得
式中，
(8)、取第一个Lyapunov函数V1为：
式中，γ>0为常数；
将V1对时间求导得：
其中，si为s1、s2或s3；为或
(9)、取自适应律为：
则
如果令η′>|Di|+|εi|，则有***全局渐近稳定；其中，Di为D1、D2或D3；
于是，滑模姿态控制器表示为：
其它步骤及参数与具体实施方式一至三之一相同。
具体实施方式五：本实施方式与具体实施方式一至四之一不同的是：步骤四中进一步采用RBF神经网络逼近符号函数η′sgn(s)，削弱抖振对姿态控制器的影响具体过程为：
(1)、采用RBF神经网络逼近η′sgn(s)，令
H＝η′sgn(s) (2-18)
式中，H＝[H1,H2,H3]T，Hi＝η′sgn(si)；理想情况下有：
此时的为逼近符号函数时第i个通道的理想权值向量， 为或为逼近符号函数时的神经网络逼近误差，εf＝[εf1,εf2,εf3]T；εfi为εf1、εf2或εf3；xfi为逼近符号函数时的神经网络输入；xfi为xf1、xf2或xf3；i取1、2、3表示卫星的X、Y和Z三通道；Hi＝H1、H2或H3；
(2)、根据η′sgn(s)的形式，将网络的输入取为xfi＝si；H的估计值为：
式中，为理想权值的估计值；设计姿态控制器为：
式中，表示εf的估计值；
(3)、将上述控制器代入(2-7)得：
式中， 为或为或
(4)、取第二个Lyapunov函数V2为：
式中，γf和γc为常数且大于0；将V2对时间求导得
(5)、取自适应律为
则
令η′>|Di|+|εi|，有***全局渐近稳定；
得到削弱抖振后的滑模姿态控制器即最终完整的姿态控制器表示为：
式中，其它步骤及参数与具体实施方式一至四之一相同。
采用以下实施例验证本发明的有益效果：
实施例：
本实施例的一种挠性卫星自适应神经网络滑模姿态控制方法，具体是按照以下步骤制备的：
为了抑制挠性模态振动和天线转动扰动，维持卫星姿态***的长期稳定运行，采用变结构控制器设计的方法和RBF神经网络的基本原理，设计了RBF神经网络滑模姿态控制器。并在此基础上对符号函数进行逼近，用以削弱抖振的影响。
研究卫星姿态问题，采用如下坐标系：
(1)地心惯性坐标系OiXiYiZi
坐标系原点取在地心处，OiXi轴沿赤道平面与黄道平面的交线指向春分点，OiZi轴指向北极方向，与地球自转轴重合，OiYi轴与OiXi轴和OiZi轴构成的平面垂直并形成右手系。
(2)轨道坐标系OXoYoZo
原点取在卫星质心处，OXo轴沿轨道平面与水平面的交线指向卫星前进方向，OZo轴沿当地铅垂线指向地心，OYo轴与OZo轴和OXo轴组成的平面形成右手系。卫星对地定向时将轨道系作为参考系。
(3)卫星惯性参考坐标系OXrYrZr
坐标原点取在卫星质心处，OXr、OYr、OZr轴指向特定方向，形成惯性坐标系。卫星惯性定向时将此坐标系作为参考系。
(4)本体坐标系OXbYbZb
坐标原点同样取在卫星质心处，OXb轴指向飞行方向，OZb轴垂直指向星体对地安装面，OYb轴垂直于OZb轴与OXb轴构成的平面并形成右手系。将OXb、OYb、OZb分别称作卫星的滚动轴、俯仰轴、偏航轴。
根据上述定义给出如所示的坐标系位置关系图如图6所示。
采用混合坐标法建立挠性卫星姿态动力学模型，含有两块帆板和一根运动天线的姿态动力学方程有以下形式：
附件模态方程为：
其中，ωs＝[ω1,ω2,ω3]T∈R3为本体系相对于惯性系且投影分解在本体系中的姿态角速度矢量；Is∈R3×3为星体转动惯量阵；u∈R3是由执行机构提供的星体三个通道控制力矩矢量，可由飞轮、动量轮、推力器等提供；d∈R3为卫星所受的干扰力矩，包括环境干扰和部件安装误差所引起的力矩等；ωa∈R3为天线相对于其支撑臂坐标系的角速度即天线转动角速度；Ia∈R3×3为天线转动惯量阵；Ta∈R3是天线转动驱动控制力矩；ηk∈Rn为挠性模态坐标，n为模态阶数，k为附件编号，k＝1、2时表示两帆板，k＝3表示天线；ξk和Ωk为n维对角阵，分别表示附件的阻尼比和模态频率；Fsk∈R3×n为附件振动与星体转动耦合系数；Rsa∈R3×3为天线与星体转动耦合系数；Fa∈R3×n为天线转动与天线臂振动耦合系数；符号表示如下的反对称矩阵，有与类似的形式。
采用欧拉角描述卫星姿态，并考虑X-Y-Z转序，相应的转动姿态角分别为滚动角俯仰角θ和偏航角ψ。卫星作惯性定向飞行时，ωs可以表示为
得到卫星姿态运动学方程为
如图7所示为带有运动天线的卫星结构示意图。其中坐标系OXbYbZb为卫星本体坐标系，OXa1Ya1Za1为天线支撑臂坐标系，OXaYaZa为天线本体坐标系。假设天线安装在卫星本体偏航轴负方向，初始时刻天线面朝向本体俯仰方向，天线本体Xa轴与本体滚动轴方向相反，天线支撑臂坐标系与天线本体系初始时刻指向完全相同，中心在星体与支撑臂连接处。根据以上设定，Ia、ωa可以详细表述为：
Ia∈R3×3为天线相对于其自身本体坐标系的转动惯量阵；
ωa∈R3为天线相对于其支撑臂坐标系的角速度；
天线具有两个自由度分别为绕本体Xa轴的转动和绕Za轴的转动，转动的角度称为俯仰角α和方向角β。考虑天线运动方式为先作方位运动后作俯仰运动，则
于是，天线转动的方程为
(1)忽略卫星姿态动力学与模态附件方程中与模态相关的高阶耦合项并化简，可得
式中，
D为干扰及不确定性的总和，且D有界，F为未知的非线性项。考虑卫星惯性定向飞行，且采用小角度假设，则卫星角速度ωs近似等于姿态角速度有
期望姿态角Θd为0，误差e＝Θd-Θ＝-Θ，取滑模函数为：
式中，c>0，c为对角阵，是待设计参数；e是卫星姿态误差；卫星的期望姿态角Θd为0，误差e＝Θd-Θ＝-Θ，Θ为卫星的实际姿态角；则根据式(3-10)和(3-11)可得
设计姿态控制器为
式中，sgn(s)＝[sgn(s1),sgn(s2),sgn(s3)]T，η′为待设计参数且η′>0，F为未知的非线性项，需要采用RBF神经网络对其进行逼近。其中，径向基函数神经网络(Radial BasisFunction Neural Network)简称RBF神经网络，是一种具有单隐层的三层前馈网络，其结构如图8所示，
hj(x)取为的高斯基函数形式，则
F＝W*Th(x)+ε (3-14)
式中，x＝[x1,x2,…,xr]T为网络输入向量，cj为RBF神经网络中第j个节点的中心向量，cj＝[cj1,cj2,…,cjr]T；bj>0为节点j的基宽值，b＝[b1,b2,…,bm]T；h(x)＝[h1(x),h2(x),…,hj(x),…hm(x)]T为m×1维向量，m为隐层数，为m×3理想网络权值阵，ε＝[ε1,ε2,ε3]T为网络逼近误差，且εi有界，将未知的非线性项F作为神经网络的输出；网络输入取为则网络输出的估计值为
其中，是权值的估计值； 或分别为卫星的X、Y和Z三通道权值的估计值；
引入估计值的姿态控制器为
将控制器带入式(3-12)得
式中，
取如下的李雅普诺夫函数
式中，γ>0为常数。将上述李雅普诺夫函数对时间求导得
取自适应律为
则
令η′>|Di|+|εi|，其中i取1、2、3表示第i个通道，则有***全局渐近稳定。
完整姿态控制器表示为
(2)式(3-22)中的控制器有符号函数项η′sgn(s)，该符号函数将使***控制量产生抖振。在干扰及不确定性的上界|Di|，相应参数η′也较大，导致抖振更加严重。所以应该继续采用RBF神经网络逼近符号函数，使控制量连续化，从而削弱抖振的影响。
采用RBF神经网络逼近η′sgn(s)，令
H＝η′sgn(s) (3-23)
式中，H＝[H1,H2,H3]T，Hi＝η′sgn(si)。理想情况下有
此时的为逼近符号函数时第i个通道的理想权值向量， 为或为逼近符号函数时的神经网络逼近误差，εf＝[εf1,εf2,εf3]T；εfi为εf1、εf2或εf3；xfi为逼近符号函数时的神经网络输入；xfi为xf1、xf2或xf3；i取1、2、3表示卫星的X、Y和Z三通道；Hi＝H1、H2或H3；根据η′sgn(s)的形式，将网络的输入取为xfi＝si；H的估计值为：
式中，为理想权值的估计值；设计姿态控制器为
将上述控制器代入(3-12)得
式中， 为或为或
取第二个Lyapunov函数V2为：
式中，γs和γc为常数且大于0。将V2对时间求导得
取自适应律为
则
令η′>|Di|+|εi|，有***全局渐近稳定。
完整的姿态控制器表示为
式中，
为验证方法的有效性，进行仿真分析。仿真中考虑帆板和天线的前5阶模态，考虑环境干扰力矩，具体参数如下：
卫星主惯量：Ix＝15000,Iy＝6000,Iz＝13000(kg·m2)；
飞轮时间常数：τx＝0.1,τy＝0.1,τz＝0.1；
帆板模态频率：Ω＝diag(0.290；0.740；1.492；1.865；3.798)×2π(rad/s)；
阻尼比：ξ＝diag(0.0262 0.0267 0.0397 0.0259 0.0178)；
耦合系数：
初始姿态：
天线惯量：Ia＝diag[8；8；1](kg·m2)；
挠性耦合系数：
转动耦合系数：
天线转动耦合系数：
控制参数：逼近sign函数前c＝diag(0.05；0.08；0.06)，η＝0.00012；
逼近sign函数后c＝diag(0.05；0.08；0.06)。
RBF神经网络参数：两次估计分别采用6-5-1结构(如图8所示，其中r＝6、m＝5和1-5-1结构(如图8所示，其中r＝1、m＝5)，第一次估计(逼近sign函数前)高斯基函数参数取ci＝[-0.2,-0.02,0,0.02,0.2]T，bj＝0.5，γ＝30；第二次估计(逼近sign函数后)高斯基函数参数取ci＝[-0.05,-0.03,0,0.03,0.05]T，bj＝0.05，权初值为0，γ＝30,γs1＝γs2＝3,γs3＝4,γc＝30。
天线仅作俯仰运动，初始时刻俯仰角α为0°，并绕俯仰轴在±60°的范围内匀速扫描，速度为1°/s。
采用式(3-22)所示的逼近sign函数前的RBF神经网络滑模姿态控制器的仿真结果如图2和3所示。
采用式(3-32)所示的逼近sign函数后的RBF神经网络滑模姿态控制器的仿真结果如图4和5所示。
从仿真结果中可以看出，两种情况下姿态控制***最终均趋于稳定，并且有效抑制了天线转动所带来的扰动。采用RBF神经网络逼近sign后，有效削弱了抖振，且姿态精度与稳定度都达到了10-4量级，在50s之前***呈振荡趋势，是由于逼近sgn函数初始时段内存在一定的逼近误差产生的影响。
本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，本领域技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (2)

1.一种挠性卫星自适应神经网络滑模姿态控制方法，其特征在于一种挠性卫星自适应神经网络滑模姿态控制方法具体是按照以下步骤进行的：
步骤一、采用混合坐标法建立挠性卫星姿态动力学模型；
步骤二、根据挠性卫星姿态动力学模型，忽略挠性卫星姿态动力学方程中与模态相关的高阶耦合项，考虑卫星作惯性定向飞行，同时采用小角度假设，得到简化后的挠性卫星姿态动力学方程；具体过程为：
(1)、忽略挠性卫星姿态动力学方程中与模态相关的高阶耦合项，将挠性卫星姿态动力学模型简化为：
<mrow> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>s</mi> </msub> <mo>=</mo> <mi>F</mi> <mo>+</mo> <mi>G</mi> <mi>u</mi> <mo>+</mo> <mi>D</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式中，
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>G</mi> <mo>=</mo> <msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>I</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>F</mi> <mrow> <mi>s</mi> <mn>1</mn> </mrow> </msub> <msubsup> <mi>F</mi> <mrow> <mi>s</mi> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <mo>-</mo> <msub> <mi>F</mi> <mrow> <mi>s</mi> <mn>2</mn> </mrow> </msub> <msubsup> <mi>F</mi> <mrow> <mi>s</mi> <mn>2</mn> </mrow> <mi>T</mi> </msubsup> <mo>-</mo> <msub> <mi>F</mi> <mrow> <mi>s</mi> <mn>3</mn> </mrow> </msub> <msubsup> <mi>F</mi> <mrow> <mi>s</mi> <mn>3</mn> </mrow> <mi>T</mi> </msubsup> </mrow> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>F</mi> <mo>=</mo> <mi>G</mi> <mrow> <mo>(</mo> <mrow> <mo>-</mo> <msubsup> <mi>&amp;omega;</mi> <mi>s</mi> <mo>&amp;times;</mo> </msubsup> <msub> <mi>I</mi> <mi>s</mi> </msub> <msub> <mi>&amp;omega;</mi> <mi>s</mi> </msub> <mo>-</mo> <msubsup> <mi>&amp;omega;</mi> <mi>s</mi> <mo>&amp;times;</mo> </msubsup> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>a</mi> </mrow> </msub> <msub> <mi>&amp;omega;</mi> <mi>a</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>D</mi> <mo>=</mo> <mi>G</mi> <mo>&amp;lsqb;</mo> <mi>d</mi> <mo>-</mo> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>a</mi> </mrow> </msub> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>a</mi> </msub> <mo>+</mo> <msub> <mi>F</mi> <mrow> <mi>s</mi> <mn>3</mn> </mrow> </msub> <msubsup> <mi>F</mi> <mi>a</mi> <mi>T</mi> </msubsup> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>a</mi> </msub> <mo>+</mo> <msub> <mi>F</mi> <mrow> <mi>s</mi> <mn>1</mn> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mn>2</mn> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <msub> <mi>&amp;Omega;</mi> <mn>1</mn> </msub> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>+</mo> <msubsup> <mi>&amp;Omega;</mi> <mn>1</mn> <mn>2</mn> </msubsup> <msub> <mi>&amp;eta;</mi> <mn>1</mn> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <msub> <mi>F</mi> <mrow> <mi>s</mi> <mn>2</mn> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mn>2</mn> <msub> <mi>&amp;xi;</mi> <mn>2</mn> </msub> <msub> <mi>&amp;Omega;</mi> <mn>2</mn> </msub> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> <mo>+</mo> <msubsup> <mi>&amp;Omega;</mi> <mn>2</mn> <mn>2</mn> </msubsup> <msub> <mi>&amp;eta;</mi> <mn>2</mn> </msub> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>F</mi> <mrow> <mi>s</mi> <mn>3</mn> </mrow> </msub> <mrow> <mo>(</mo> <mrow> <mn>2</mn> <msub> <mi>&amp;xi;</mi> <mn>3</mn> </msub> <msub> <mi>&amp;Omega;</mi> <mn>3</mn> </msub> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>3</mn> </msub> <mo>+</mo> <msubsup> <mi>&amp;Omega;</mi> <mn>3</mn> <mn>2</mn> </msubsup> <msub> <mi>&amp;eta;</mi> <mn>3</mn> </msub> </mrow> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <mn>2</mn> <mo>-</mo> <mn>4</mn> </mrow> <mo>)</mo> </mrow> </mrow>
D为干扰及不确定性的总和，且D有界；F为未知的非线性项；u为姿态控制器；
(2)、考虑卫星惯性定向飞行，且采用小角度假设，近似的认为卫星角速度ωs近似等于姿态角速度则简化后的挠性卫星姿态动力学方程具体为：
<mrow> <mover> <mi>&amp;Theta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>=</mo> <mi>F</mi> <mo>+</mo> <mi>G</mi> <mi>u</mi> <mo>+</mo> <mi>D</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>5</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
步骤三、根据简化后的挠性卫星姿态动力学方程，利用RBF神经网络设计滑模姿态控制器；具体过程为：
(1)、取滑模函数s为：
<mrow> <mi>s</mi> <mo>=</mo> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>c</mi> <mi>e</mi> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>6</mn> <mo>)</mo> </mrow> </mrow>
卫星的期望姿态角Θd为0，误差e＝Θd-Θ＝-Θ，Θ为卫星的实际姿态角；
(2)、联立式(2-5)和式(2-6)得：
<mrow> <mover> <mi>s</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <mover> <mi>e</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>c</mi> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <mo>-</mo> <mover> <mi>&amp;Theta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>c</mi> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <mo>-</mo> <mi>F</mi> <mo>-</mo> <mi>G</mi> <mi>u</mi> <mo>-</mo> <mi>D</mi> <mo>+</mo> <mi>c</mi> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>7</mn> <mo>)</mo> </mrow> </mrow>
(3)、设计姿态控制器的基本形式为：
<mrow> <mi>u</mi> <mo>=</mo> <msup> <mi>G</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;lsqb;</mo> <mo>-</mo> <mi>F</mi> <mo>+</mo> <mi>c</mi> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>8</mn> <mo>)</mo> </mrow> </mrow>
(4)、采用RBF神经网络对未知的非线性项F进行逼近：
取hj(x)为的高斯基函数形式，则：
F＝W*Th(x)+ε (2-9)
式中，x＝[x1,x2,…,xr]T为是神经网络的r维输入向量，cj为RBF神经网络中第j个节点的中心向量，cj＝[cj1,cj2,…,cjr]T；bj＞0为节点j的基宽值，b＝[b1,b2,…,bm]T；h(x)＝[h1(x),h2(x),…,hj(x),…hm(x)]T为m×1维神经网络径向基向量，m为隐层数；为m×3理想网络权值阵；ε＝[ε1,ε2,ε3]T为神经网络逼近误差，且εi有界，i取1、2、3表示卫星的X、Y和Z三通道；将未知的非线性项F作为神经网络的输出；
(5)、将公式(2-9)的网络输入取为则网络输出的估计值为：
<mrow> <mover> <mi>F</mi> <mo>^</mo> </mover> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mover> <mi>W</mi> <mo>^</mo> </mover> <mi>T</mi> </msup> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>10</mn> <mo>)</mo> </mrow> </mrow>
其中，是权值的估计值； 或分别为卫星的X、Y和Z三通道权值的估计值；
(6)、将网络输出的估计值引入姿态控制器u即公式(2-8)中：
<mrow> <mi>u</mi> <mo>=</mo> <msup> <mi>G</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;lsqb;</mo> <mo>-</mo> <mover> <mi>F</mi> <mo>^</mo> </mover> <mo>+</mo> <mi>c</mi> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>11</mn> <mo>)</mo> </mrow> </mrow>
(7)、将控制器代入式(2-7)得
<mrow> <mtable> <mtr> <mtd> <mrow> <mover> <mi>s</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <mo>-</mo> <mi>F</mi> <mo>-</mo> <mi>G</mi> <mi>u</mi> <mo>-</mo> <mi>D</mi> <mo>+</mo> <mi>c</mi> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <mo>-</mo> <mi>F</mi> <mo>-</mo> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mo>-</mo> <mover> <mi>F</mi> <mo>^</mo> </mover> <mo>+</mo> <mi>c</mi> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>-</mo> <mi>D</mi> <mo>+</mo> <mi>c</mi> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <mi>F</mi> <mo>+</mo> <mover> <mi>F</mi> <mo>^</mo> </mover> <mo>-</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>D</mi> <mo>=</mo> <mo>-</mo> <mover> <mi>F</mi> <mo>~</mo> </mover> <mo>-</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>D</mi> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <mn>2</mn> <mo>-</mo> <mn>12</mn> </mrow> <mo>)</mo> </mrow> </mrow>
式中，
(8)、取第一个Lyapunov函数V1为：
<mrow> <msub> <mi>V</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>s</mi> <mi>T</mi> </msup> <mi>s</mi> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>&amp;gamma;</mi> <mi>t</mi> <mi>r</mi> <mrow> <mo>(</mo> <msup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>T</mi> </msup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>13</mn> <mo>)</mo> </mrow> </mrow>
将V1对时间求导得：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>=</mo> <msup> <mi>s</mi> <mi>T</mi> </msup> <mover> <mi>s</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>&amp;gamma;</mi> <mi>t</mi> <mi>r</mi> <mrow> <mo>(</mo> <mrow> <msup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>T</mi> </msup> <mover> <mover> <mi>W</mi> <mo>~</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> </mrow> <mo>)</mo> </mrow> <mo>=</mo> <msup> <mi>s</mi> <mi>T</mi> </msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mo>-</mo> <mover> <mi>F</mi> <mo>~</mo> </mover> <mo>-</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>D</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>-</mo> <mi>&amp;gamma;</mi> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>i</mi> <mi>T</mi> </msubsup> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <msup> <mi>s</mi> <mi>T</mi> </msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mo>-</mo> <msup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>T</mi> </msup> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>&amp;epsiv;</mi> <mo>-</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>D</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>-</mo> <mi>&amp;gamma;</mi> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>i</mi> <mi>T</mi> </msubsup> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msub> <mi>s</mi> <mi>i</mi> </msub> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>i</mi> <mi>T</mi> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>&amp;gamma;</mi> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>i</mi> <mi>T</mi> </msubsup> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <msup> <mi>s</mi> <mi>T</mi> </msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>&amp;epsiv;</mi> <mo>+</mo> <mi>D</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>i</mi> <mi>T</mi> </msubsup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msub> <mi>s</mi> <mi>i</mi> </msub> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>&amp;gamma;</mi> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>-</mo> <msup> <mi>s</mi> <mi>T</mi> </msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>&amp;epsiv;</mi> <mo>+</mo> <mi>D</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <mn>2</mn> <mo>-</mo> <mn>14</mn> </mrow> <mo>)</mo> </mrow> </mrow>
其中，si为s1、s2或s3；为或
(9)、取自适应律为：
<mrow> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mi>&amp;gamma;</mi> </mfrac> <msub> <mi>s</mi> <mi>i</mi> </msub> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>15</mn> <mo>)</mo> </mrow> </mrow>
则
<mrow> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>=</mo> <mo>-</mo> <msup> <mi>s</mi> <mi>T</mi> </msup> <mo>&amp;lsqb;</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>&amp;epsiv;</mi> <mo>+</mo> <mi>D</mi> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>16</mn> <mo>)</mo> </mrow> </mrow>
如果令η′＞|Di|+|εi|，则有***全局渐近稳定；其中，Di为D1、D2或D3；
于是，滑模姿态控制器表示为：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>u</mi> <mo>=</mo> <msup> <mi>G</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;lsqb;</mo> <mo>-</mo> <msup> <mover> <mi>W</mi> <mo>^</mo> </mover> <mi>T</mi> </msup> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>c</mi> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mi>&amp;gamma;</mi> </mfrac> <msub> <mi>s</mi> <mi>i</mi> </msub> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>17</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
步骤四、进一步采用RBF神经网络逼近符号函数η′sgn(s)，削弱抖振对滑模姿态控制器的影响，得到削弱抖振后的滑模姿态控制器；其中，sgn(s)＝[sgn(s1),sgn(s2),sgn(s3)]T，η′为待设计参数且η′＞0，s1、s2和s3分别为X、Y和Z三通道的滑模函数；s＝[s1,s2,s3]T；s为滑模函数；
步骤四中进一步采用RBF神经网络逼近符号函数η′sgn(s)，削弱抖振对姿态控制器的影响具体过程为：
(1)、采用RBF神经网络逼近η′sgn(s)，令
H＝η′sgn(s) (2-18)
式中，H＝[H1,H2,H3]T，Hi＝η′sgn(si)；理想情况下有：
<mrow> <msub> <mi>H</mi> <mi>i</mi> </msub> <mo>=</mo> <msubsup> <mi>W</mi> <mrow> <mi>f</mi> <mi>i</mi> </mrow> <mrow> <mo>*</mo> <mi>T</mi> </mrow> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&amp;epsiv;</mi> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>19</mn> <mo>)</mo> </mrow> </mrow>
此时的为逼近符号函数时第i个通道的理想权值向量， 为或εfi＞0为逼近符号函数时的神经网络逼近误差，εf＝[εf1,εf2,εf3]T；εfi为εf1、εf2或εf3；xfi为逼近符号函数时的神经网络输入；xfi为xf1、xf2或xf3；i取1、2、3表示卫星的X、Y和Z三通道；Hi＝H1、H2或H3；
(2)、根据η′sgn(s)的形式，将网络的输入取为xfi＝si；H的估计值为：
<mrow> <mover> <mi>H</mi> <mo>^</mo> </mover> <mo>=</mo> <msup> <mrow> <mo>&amp;lsqb;</mo> <msub> <mover> <mi>H</mi> <mo>^</mo> </mover> <mn>1</mn> </msub> <mo>,</mo> <msub> <mover> <mi>H</mi> <mo>^</mo> </mover> <mn>2</mn> </msub> <mo>,</mo> <msub> <mover> <mi>H</mi> <mo>^</mo> </mover> <mn>3</mn> </msub> <mo>&amp;rsqb;</mo> </mrow> <mi>T</mi> </msup> <mo>,</mo> <msub> <mover> <mi>H</mi> <mo>^</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <msubsup> <mover> <mi>W</mi> <mo>^</mo> </mover> <mrow> <mi>f</mi> <mi>i</mi> </mrow> <mi>T</mi> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>20</mn> <mo>)</mo> </mrow> </mrow>
式中，为理想权值的估计值；设计姿态控制器为：
<mrow> <mi>u</mi> <mo>=</mo> <msup> <mi>G</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;lsqb;</mo> <mo>-</mo> <msup> <mover> <mi>W</mi> <mo>^</mo> </mover> <mi>T</mi> </msup> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>c</mi> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <mover> <mi>H</mi> <mo>^</mo> </mover> <mo>+</mo> <msub> <mover> <mi>&amp;epsiv;</mi> <mo>^</mo> </mover> <mi>f</mi> </msub> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>21</mn> <mo>)</mo> </mrow> </mrow>
式中，表示εf的估计值；c为对角阵，是待设计参数；e是卫星姿态误差；
(3)、将上述控制器代入(2-7)得：
<mrow> <mtable> <mtr> <mtd> <mrow> <mover> <mi>s</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>=</mo> <mo>-</mo> <msup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>T</mi> </msup> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>&amp;epsiv;</mi> <mo>-</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mover> <mi>W</mi> <mo>^</mo> </mover> <mrow> <mi>f</mi> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mover> <mi>W</mi> <mo>^</mo> </mover> <mrow> <mi>f</mi> <mn>2</mn> </mrow> <mi>T</mi> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mover> <mi>W</mi> <mo>^</mo> </mover> <mrow> <mi>f</mi> <mn>3</mn> </mrow> <mi>T</mi> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mn>3</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <msub> <mover> <mi>&amp;epsiv;</mi> <mo>^</mo> </mover> <mi>f</mi> </msub> <mo>-</mo> <mi>D</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <msup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>T</mi> </msup> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>&amp;epsiv;</mi> <mo>+</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mrow> <mi>f</mi> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mrow> <mi>f</mi> <mn>2</mn> </mrow> <mi>T</mi> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mrow> <mi>f</mi> <mn>3</mn> </mrow> <mi>T</mi> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mn>3</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mover> <mi>&amp;epsiv;</mi> <mo>~</mo> </mover> <mi>f</mi> </msub> <mo>-</mo> <mi>D</mi> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <mn>2</mn> <mo>-</mo> <mn>22</mn> </mrow> <mo>)</mo> </mrow> </mrow>
式中， 为或 为或
(4)、取第二个Lyapunov函数V2为：
<mrow> <msub> <mi>V</mi> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msup> <mi>s</mi> <mi>T</mi> </msup> <mi>s</mi> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <mi>&amp;gamma;</mi> <mi>t</mi> <mi>r</mi> <mrow> <mo>(</mo> <msup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>T</mi> </msup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mo>)</mo> </mrow> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>&amp;gamma;</mi> <mi>f</mi> </msub> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mrow> <mi>f</mi> <mi>i</mi> </mrow> <mi>T</mi> </msubsup> <msub> <mover> <mi>W</mi> <mo>~</mo> </mover> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> <mo>+</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msub> <mi>&amp;gamma;</mi> <mi>c</mi> </msub> <msubsup> <mover> <mi>&amp;epsiv;</mi> <mo>~</mo> </mover> <mi>f</mi> <mi>T</mi> </msubsup> <msub> <mover> <mi>&amp;epsiv;</mi> <mo>~</mo> </mover> <mi>f</mi> </msub> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>23</mn> <mo>)</mo> </mrow> </mrow>
式中，γf和γc为常数且大于0；γ＞0为常数；将V2对时间求导得
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> <mo>=</mo> <msup> <mi>s</mi> <mi>T</mi> </msup> <mover> <mi>s</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <mi>&amp;gamma;</mi> <mi>t</mi> <mi>r</mi> <mrow> <mo>(</mo> <mrow> <msup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>T</mi> </msup> <mover> <mover> <mi>W</mi> <mo>~</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> </mrow> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>&amp;gamma;</mi> <mi>f</mi> </msub> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mrow> <mi>f</mi> <mi>i</mi> </mrow> <mi>T</mi> </msubsup> <msub> <mover> <mover> <mi>W</mi> <mo>~</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> <mo>+</mo> <msub> <mi>&amp;gamma;</mi> <mi>c</mi> </msub> <msubsup> <mover> <mi>&amp;epsiv;</mi> <mo>~</mo> </mover> <mi>f</mi> <mi>T</mi> </msubsup> <msub> <mover> <mover> <mi>&amp;epsiv;</mi> <mo>~</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>f</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <msup> <mi>s</mi> <mi>T</mi> </msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mo>-</mo> <msup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>T</mi> </msup> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>-</mo> <mi>&amp;epsiv;</mi> <mo>+</mo> <mfenced open = "[" close = "]"> <mtable> <mtr> <mtd> <mrow> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mrow> <mi>f</mi> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mn>1</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mrow> <mi>f</mi> <mn>2</mn> </mrow> <mi>T</mi> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mn>2</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mrow> <mi>f</mi> <mn>3</mn> </mrow> <mi>T</mi> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mn>3</mn> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mover> <mi>&amp;epsiv;</mi> <mo>~</mo> </mover> <mi>f</mi> </msub> <mo>-</mo> <mi>D</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <mi>&amp;gamma;</mi> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>i</mi> <mi>T</mi> </msubsup> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mi>f</mi> </msub> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mrow> <mi>f</mi> <mi>i</mi> </mrow> <mi>T</mi> </msubsup> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mi>c</mi> </msub> <msubsup> <mover> <mi>&amp;epsiv;</mi> <mo>~</mo> </mover> <mi>f</mi> <mi>T</mi> </msubsup> <msub> <mover> <mover> <mi>&amp;epsiv;</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>f</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>=</mo> <mo>-</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mi>i</mi> <mi>T</mi> </msubsup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msub> <mi>s</mi> <mi>i</mi> </msub> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>&amp;gamma;</mi> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>+</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>3</mn> </munderover> <msubsup> <mover> <mi>W</mi> <mo>~</mo> </mover> <mrow> <mi>f</mi> <mi>i</mi> </mrow> <mi>T</mi> </msubsup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msub> <mi>s</mi> <mi>i</mi> </msub> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mi>f</mi> </msub> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msup> <mi>s</mi> <mi>T</mi> </msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mi>&amp;epsiv;</mi> <mo>+</mo> <mi>D</mi> <mo>+</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mo>+</mo> <msubsup> <mover> <mi>&amp;epsiv;</mi> <mo>~</mo> </mover> <mi>f</mi> <mi>T</mi> </msubsup> <mrow> <mo>(</mo> <mrow> <mi>s</mi> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mi>c</mi> </msub> <msub> <mover> <mover> <mi>&amp;epsiv;</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>f</mi> </msub> </mrow> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <mn>2</mn> <mo>-</mo> <mn>24</mn> </mrow> <mo>)</mo> </mrow> </mrow>
(5)、取自适应律为
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mi>&amp;gamma;</mi> </mfrac> <msub> <mi>s</mi> <mi>i</mi> </msub> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>&amp;gamma;</mi> <mi>f</mi> </msub> </mfrac> <msub> <mi>s</mi> <mi>i</mi> </msub> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mover> <mi>&amp;epsiv;</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>f</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>&amp;gamma;</mi> <mi>c</mi> </msub> </mfrac> <mi>s</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <mn>2</mn> <mo>-</mo> <mn>25</mn> </mrow> <mo>)</mo> </mrow> </mrow>
则
<mrow> <msub> <mover> <mi>V</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> <mo>=</mo> <mo>-</mo> <msup> <mi>s</mi> <mi>T</mi> </msup> <mo>&amp;lsqb;</mo> <msup> <mi>&amp;eta;</mi> <mo>&amp;prime;</mo> </msup> <mi>sgn</mi> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>&amp;epsiv;</mi> <mo>+</mo> <mi>D</mi> <mo>&amp;rsqb;</mo> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>26</mn> <mo>)</mo> </mrow> </mrow>
令η′＞|Di|+|εi|，有***全局渐近稳定；
得到削弱抖振后的滑模姿态控制器即最终完整的姿态控制器表示为：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <mi>u</mi> <mo>=</mo> <msup> <mi>G</mi> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <mo>&amp;lsqb;</mo> <mo>-</mo> <msup> <mover> <mi>W</mi> <mo>^</mo> </mover> <mi>T</mi> </msup> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>+</mo> <mi>c</mi> <mover> <mi>e</mi> <mo>&amp;CenterDot;</mo> </mover> <mo>+</mo> <msubsup> <mover> <mi>W</mi> <mo>^</mo> </mover> <mi>f</mi> <mi>T</mi> </msubsup> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mover> <mi>&amp;epsiv;</mi> <mo>^</mo> </mover> <mi>f</mi> </msub> <mo>&amp;rsqb;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>i</mi> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mn>1</mn> <mi>&amp;gamma;</mi> </mfrac> <msub> <mi>s</mi> <mi>i</mi> </msub> <mi>h</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mover> <mi>W</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>&amp;gamma;</mi> <mi>f</mi> </msub> </mfrac> <msub> <mi>s</mi> <mi>i</mi> </msub> <mi>h</mi> <mrow> <mo>(</mo> <msub> <mi>x</mi> <mrow> <mi>f</mi> <mi>i</mi> </mrow> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mover> <mi>&amp;epsiv;</mi> <mo>^</mo> </mover> <mo>&amp;CenterDot;</mo> </mover> <mi>f</mi> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <msub> <mi>&amp;gamma;</mi> <mi>c</mi> </msub> </mfrac> <mi>s</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>27</mn> <mo>)</mo> </mrow> </mrow>
式中，即完成了一种挠性卫星自适应神经网络滑模姿态控制方法。
2.根据权利要求1所述一种挠性卫星自适应神经网络滑模姿态控制方法，其特征在于：步骤一中采用混合坐标法建立挠性卫星姿态动力学模型具体为：
(1)、含有两块帆板和一根运动天线的姿态动力学方程有以下形式：
<mrow> <mtable> <mtr> <mtd> <mrow> <msub> <mi>I</mi> <mi>s</mi> </msub> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>s</mi> </msub> <mo>+</mo> <msubsup> <mi>&amp;omega;</mi> <mi>s</mi> <mo>&amp;times;</mo> </msubsup> <msub> <mi>I</mi> <mi>s</mi> </msub> <msub> <mi>&amp;omega;</mi> <mi>s</mi> </msub> <mo>+</mo> <munder> <mi>&amp;Sigma;</mi> <mi>k</mi> </munder> <mrow> <mo>(</mo> <msub> <mi>F</mi> <mrow> <mi>s</mi> <mi>k</mi> </mrow> </msub> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>k</mi> </msub> <mo>+</mo> <msubsup> <mi>&amp;omega;</mi> <mi>s</mi> <mo>&amp;times;</mo> </msubsup> <msub> <mi>F</mi> <mrow> <mi>s</mi> <mi>k</mi> </mrow> </msub> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>k</mi> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>a</mi> </mrow> </msub> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>a</mi> </msub> <mo>+</mo> <msubsup> <mi>&amp;omega;</mi> <mi>s</mi> <mo>&amp;times;</mo> </msubsup> <msub> <mi>R</mi> <mrow> <mi>s</mi> <mi>a</mi> </mrow> </msub> <msub> <mi>&amp;omega;</mi> <mi>a</mi> </msub> <mo>=</mo> <mi>u</mi> <mo>+</mo> <mi>d</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>I</mi> <mi>a</mi> </msub> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>a</mi> </msub> <mo>+</mo> <msubsup> <mi>&amp;omega;</mi> <mi>a</mi> <mo>&amp;times;</mo> </msubsup> <msub> <mi>I</mi> <mi>a</mi> </msub> <msub> <mi>&amp;omega;</mi> <mi>a</mi> </msub> <mo>+</mo> <msub> <mi>F</mi> <mi>a</mi> </msub> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>3</mn> </msub> <mo>+</mo> <msubsup> <mi>&amp;omega;</mi> <mi>a</mi> <mo>&amp;times;</mo> </msubsup> <msub> <mi>F</mi> <mi>a</mi> </msub> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>3</mn> </msub> <mo>+</mo> <msubsup> <mi>R</mi> <mrow> <mi>s</mi> <mi>a</mi> </mrow> <mi>T</mi> </msubsup> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>s</mi> </msub> <mo>+</mo> <msubsup> <mi>&amp;omega;</mi> <mi>a</mi> <mo>&amp;times;</mo> </msubsup> <msubsup> <mi>R</mi> <mrow> <mi>s</mi> <mi>a</mi> </mrow> <mi>T</mi> </msubsup> <msub> <mi>&amp;omega;</mi> <mi>s</mi> </msub> <mo>=</mo> <msub> <mi>T</mi> <mi>a</mi> </msub> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mrow> <mn>2</mn> <mo>-</mo> <mn>1</mn> </mrow> <mo>)</mo> </mrow> </mrow>
其中，Is∈R3×3为星体转动惯量阵；ωs＝[ω1,ω2,ω3]T∈R3为本体系相对于惯性系且投影分解在本体系中的姿态角速度矢量；ω1、ω2和ω3分别为本体系相对于惯性系且投影分解在本体系中X、Y和Z轴的姿态角速度；ηk∈Rn为挠性模态坐标，n为模态阶数，k为附件编号，k＝1、2时表示本体Y轴正和负方向帆板，k＝3表示天线；Fsk∈R3×n为附件振动与星体转动耦合系数；Rsa∈R3×3为天线与星体转动耦合系数；u∈R3是由执行机构提供的星体三个通道控制力矩矢量；d∈R3为卫星所受的干扰力矩，包括环境干扰力矩和部件安装误差所引起的干扰力矩等；Ia∈R3×3为天线转动惯量阵；ωa∈R3，ωa＝[ωa1,ωa2,ωa3]T为天线相对于天线支撑臂坐标系的角速度矢量；ωa1、ωa2和ωa3分别为天线相对于天线支撑臂坐标系X、Y和Z轴的角速度；Fa∈R3×n为天线转动与天线臂振动耦合系数；符号表示如下的反对称矩阵，类似的，表示如下的反对称矩阵Ta∈R3是天线转动驱动控制力矩；角标a表示天线；角标s表示星体；
(2)、建立附件模态方程为：
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>+</mo> <mn>2</mn> <msub> <mi>&amp;xi;</mi> <mn>1</mn> </msub> <msub> <mi>&amp;Omega;</mi> <mn>1</mn> </msub> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>1</mn> </msub> <mo>+</mo> <msubsup> <mi>&amp;Omega;</mi> <mn>1</mn> <mn>2</mn> </msubsup> <msub> <mi>&amp;eta;</mi> <mn>1</mn> </msub> <mo>+</mo> <msubsup> <mi>F</mi> <mrow> <mi>s</mi> <mn>1</mn> </mrow> <mi>T</mi> </msubsup> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>s</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> <mo>+</mo> <mn>2</mn> <msub> <mi>&amp;xi;</mi> <mn>2</mn> </msub> <msub> <mi>&amp;Omega;</mi> <mn>2</mn> </msub> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>2</mn> </msub> <mo>+</mo> <msubsup> <mi>&amp;Omega;</mi> <mn>2</mn> <mn>2</mn> </msubsup> <msub> <mi>&amp;eta;</mi> <mn>2</mn> </msub> <mo>+</mo> <msubsup> <mi>F</mi> <mrow> <mi>s</mi> <mn>2</mn> </mrow> <mi>T</mi> </msubsup> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>s</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mn>3</mn> </msub> <mo>+</mo> <mn>2</mn> <msub> <mi>&amp;xi;</mi> <mn>3</mn> </msub> <msub> <mi>&amp;Omega;</mi> <mn>3</mn> </msub> <msub> <mover> <mi>&amp;eta;</mi> <mo>&amp;CenterDot;</mo> </mover> <mn>3</mn> </msub> <mo>+</mo> <msubsup> <mi>&amp;Omega;</mi> <mn>3</mn> <mn>2</mn> </msubsup> <msub> <mi>&amp;eta;</mi> <mn>3</mn> </msub> <mo>+</mo> <msubsup> <mi>F</mi> <mrow> <mi>s</mi> <mn>3</mn> </mrow> <mi>T</mi> </msubsup> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>s</mi> </msub> <mo>+</mo> <msubsup> <mi>F</mi> <mi>a</mi> <mi>T</mi> </msubsup> <msub> <mover> <mi>&amp;omega;</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>a</mi> </msub> <mo>=</mo> <mn>0</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>-</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
其中，ξk和Ωk为n维对角阵，ξk表示附件的阻尼比；Ωk表示附件的模态频率，k＝1或2时表示本体Y轴正和负方向帆板，k＝3时表示天线；
(3)、含有两块帆板和一根运动天线的姿态动力学方程与附件模态方程组成挠性卫星姿态动力学模型。
CN201510172378.4A 2015-04-13 2015-04-13 一种挠性卫星自适应神经网络滑模姿态控制方法 Active CN104898418B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510172378.4A CN104898418B (zh) 2015-04-13 2015-04-13 一种挠性卫星自适应神经网络滑模姿态控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510172378.4A CN104898418B (zh) 2015-04-13 2015-04-13 一种挠性卫星自适应神经网络滑模姿态控制方法

Publications (2)

Publication Number Publication Date
CN104898418A CN104898418A (zh) 2015-09-09
CN104898418B true CN104898418B (zh) 2017-08-25

Family

ID=54031137

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510172378.4A Active CN104898418B (zh) 2015-04-13 2015-04-13 一种挠性卫星自适应神经网络滑模姿态控制方法

Country Status (1)

Country Link
CN (1) CN104898418B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105468007B (zh) * 2015-12-02 2018-07-06 哈尔滨工业大学 一种基于干扰观测器的挠性卫星轨迹线性化姿态控制方法
CN105607485B (zh) * 2016-02-04 2018-08-14 河北科技师范学院 基于故障特征模型的挠性充液卫星姿态自适应容错控制方法
CN106777469B (zh) * 2016-11-15 2020-07-03 上海航天控制技术研究所 用于挠性卫星挠性振动评估的地面物理仿真验证***和验证方法
CN107831775B (zh) * 2017-11-14 2021-06-08 哈尔滨工业大学深圳研究生院 基于挠性航天器无角速度测量的姿态控制方法
CN108319135B (zh) * 2018-01-18 2019-12-27 西北工业大学 基于神经网络鲁棒自适应的空间柔性***展开控制方法
CN108398885B (zh) * 2018-03-29 2019-03-26 湖南大学 旋翼飞行机械臂自适应RBFNNs测噪自抗扰控制方法
CN118068694B (zh) * 2024-04-19 2024-07-05 南京信息工程大学 一种自适应神经网络模糊积分滑模控制方法及装置

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102073276B (zh) * 2011-02-21 2012-06-27 哈尔滨工业大学 采用rbf神经网络的桡性结构自适应变结构控制方法
CN103116275B (zh) * 2013-03-01 2016-04-06 河海大学常州校区 基于滑模补偿的微陀螺仪鲁棒神经网络控制***及方法
CN103324087B (zh) * 2013-06-19 2015-10-07 河海大学常州校区 基于神经网络的微陀螺仪的自适应反演控制***及方法

Also Published As

Publication number Publication date
CN104898418A (zh) 2015-09-09

Similar Documents

Publication Publication Date Title
CN104898418B (zh) 一种挠性卫星自适应神经网络滑模姿态控制方法
CN104527994B (zh) 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法
Jia et al. Integral backstepping sliding mode control for quadrotor helicopter under external uncertain disturbances
Sun et al. Adaptive backstepping control of spacecraft rendezvous and proximity operations with input saturation and full-state constraint
CN106292681B (zh) 一种基于观测器和在线控制分配的卫星主动容错控制方法
CN104898683B (zh) 一种挠性卫星神经网络反步滑模姿态控制方法
Xin et al. Indirect robust control of spacecraft via optimal control solution
CN105785762B (zh) 一种基于自适应反步滑模的两轴惯性稳定平台高精度控制方法
CN105353763B (zh) 一种非合作目标航天器相对轨道姿态有限时间控制方法
CN106933241B (zh) 基于故障解耦的单框架控制力矩陀螺航天器容错控制方法
CN104570742B (zh) 基于前馈pid控制的异面交叉快变轨道快速高精度相对指向控制方法
CN104267732A (zh) 基于频域分析的挠性卫星高稳定度姿态控制方法
Wu et al. Sliding-mode control for staring-mode spacecraft using a disturbance observer
CN112572835B (zh) 一种具有姿态切换的卫星在轨角动量管理及控制方法
CN104656447A (zh) 一种航天器抗干扰姿态跟踪的微分几何非线性控制方法
Shi et al. A robust attitude controller for a spacecraft equipped with a robotic manipulator
CN103863578B (zh) 火星着陆器喷气推力器和控制力矩陀螺复合控制***
Sun Saturated adaptive output-constrained control of cooperative spacecraft rendezvous and docking
Wu et al. Robust sliding mode control for spacecraft global fast-tracking manoeuvre
CN110502038B (zh) 一种机动过程中天线预置的高稳定度控制方法
Chen et al. Decoupling attitude control of a hypersonic glide vehicle based on a nonlinear extended state observer
Dell’Elce et al. Comparison between analytical and optimal control techniques in the differential drag based rendez-vous
CN103869823B (zh) 火星着陆器喷气推力器和质量矩复合控制***
CN106005483A (zh) 一种模块化手机星的主动姿态控制方法
Khelfi et al. Robust control with sliding mode for a quadrotor unmanned aerial vehicle

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Sun Yanchao

Inventor after: Ma Guangfu

Inventor after: Li Chuanjiang

Inventor after: Zhang Chao

Inventor after: Wang Pengyu

Inventor after: Jiang Lisong

Inventor before: Zhang Chao

Inventor before: Sun Yanchao

Inventor before: Wang Pengyu

Inventor before: Jiang Lisong

Inventor before: Ma Guangfu

Inventor before: Li Chuanjiang

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant