1.一种基于增强型幂次趋近律和快速终端滑模面的刚性航天飞行器有限时间自适应容错控制方法,其特征在于:所述控制方法包括以下步骤:步骤1,建立飞行器姿态容错控制系统的运动学和动力学模型,初始化系统状态以及控制参数,过程如下:
1.1飞行器姿态控制系统的动力学模型表达形式为:其中,ω, 分别是飞行器的角速度和角加速度;Ω∈Rn是反作用飞轮的角速度;×是运算符号,将运算符号×应用于a=[a1,a2,a3]T得a×=[0,-a3,a2;a3,0,-a1;-a2,a1,0];J∈R3×3是飞行器的转动惯性矩阵;Jω=diag([Jω1;Jω2,...,Jωn])∈Rn×n是反作用飞轮的转动惯性矩阵;D∈Rn是反作用飞轮控制力矩分配矩阵且行满秩;u∈R3和d(t)∈R3是控制输入和外部扰动;
1.2飞行器姿态控制系统的运动学模型表达形式为:其中,单位四元数 是飞行器姿态四元数且满足分别是q0和qv的导数;I3∈R3×3是3×3单位矩阵;
1.3假设转动惯性矩阵J=J0+ΔJ,其中J0和ΔJ分别表示J的标称部分和不确定部分;且控制输入表示为 其中E(t)=diag([e1(t),e2(t),...,en(t)])∈Rn×n是执行器控制效率矩阵;0≤ei(t)≤1是第i个反作用飞轮的效率因子;
是附加执行器故障矢量;uc=[uc1,uc2,...,ucn]T∈Rn是第n个执行器的控制力矩矢量;则式(1)重新写成:
1.4令 代入式(2),得到:
其中,
对式(5)进行微分,得到:
其中, 分别为P和qv的一阶导数和二阶导数;
将式(5)、式(6)代入式(4)后,在等式两边同时左乘PT得到:* T
其中,J =PJ0P且由于转动惯性矩阵J*是斜对称正定矩阵,则矩阵 满足以下斜对称关系:同时J*满足以下不等式:
其中,Jmin和Jmax是正常数,表示J*的下界和上界;
是干扰和不确定性的集合,满足||Td||≤υ0Φ,Φ=1+||ω||+||ω||2且υ0是正常数;
步骤2,在存在转动惯量不确定和外部扰动的情况下,基于飞行器的姿态控制系统,设计所需的滑模面,过程如下:
2.1选择快速终端滑模面s∈R3为:其中,λ1>0;λ2>0;a1>a2>1;函数sig(x)r=[|x1|rsign(x1),|x2|rsign(x2),|x3|rsign(x3)]T;
sign(·)为·符号函数;
对式(10)求导,得到:
其中,为s的一阶导数; |α(qv)|为α(qv)的绝对值;
如果α(qv)=0且β(qv)≠0,由于负分数幂 的存在会产生奇异性,为避免奇异性的产生,s的一阶导数改变为:其中,qve∈R3定义为:
其中,∈是很小的常数;|∈|是∈的绝对值;
然后,由式(7),式(10)和式(12)得到:其中,
步骤3,设计增强型幂次趋近律,过程如下:
3.1定义增强型幂次趋近律为:
其中,0<θ<1;K>0;0<p<1; θ>0;||s||为s的范数;
步骤4,设计有限时间自适应滑模控制器,过程如下:
4.1考虑有限时间自适应滑模控制器被设计为:其中,||P||为P的范数;||F||为F的范数;||Ps||为Ps的范数;||s||为s的范数;γ0=υ0; 正定矩阵DEDT满足:0<e0≤min{λmin(DEDT),1};λmin(·)表示矩阵的最小特征值;e0是一个正常数; 为γi的估计;i=0,1,2;
4.2设计自适应参数的更新律:
其中,ci和εi是正常数; 为 的一阶导数;i=0,1,2;
4.3设计李雅普诺夫函数:
其中, sT是s的转置;
对式(22)进行求导,如果将式(22)写成 的形式,则判定系统是有限时间一致最终有界;其中,
基于以上分析,滑模面s、飞行器姿态四元数qv和角速度ω是局部有限时间一致最终有界。