CN105843244A - 一种基于输出反馈的挠性航天器精细姿态控制方法 - Google Patents

Abstract

本发明提供了一种基于输出反馈的挠性航天器精细姿态控制方法，所述控制方法包括如下步骤：a)搭建挠性航天器动力学系统Σ₁，将挠性航天器动力学系统Σ₁变换为挠性航天器动力学系统Σ₂，所述挠性航天器动力学系统Σ₂引入航天器刚柔耦合动态干扰；b)构造外部系统Σ₃，所述刚柔耦合动态干扰通过外部系统Σ₃进行描述；c)设计干扰观测器，所述干扰观测器对所述刚柔耦合动态干扰进行估值；d)设计动态输出反馈H_∞控制器；e)将步骤c)中所诉的干扰观测器与步骤d)中所述的动态输出反馈H_∞控制器进行复合，得到挠性航天器精细姿态控制系统Σ₆；所述挠性航天器精细姿态控制系统Σ₆通过估值对刚柔耦合动态干扰进行补偿。

Description

一种基于输出反馈的挠性航天器精细姿态控制方法

技术领域

本发明涉及航空航天技术领域，特别涉及一种基于输出反馈的挠性航天器精细姿态控制方法。

背景技术

随着航天技术的发展，越来越多的航天器携带有大型的太阳帆板、天线等挠性机构，这些挠性机构使航天器的功能日益强大，但同时也给航天器的姿态控制问题带来了严峻的考验，尤其是当航天器的控制精度要求较高时，如对地观测卫星、太空望远镜等。因此，挠性航天器在天空运行时不仅会受到外部的环境干扰，还会受到刚柔耦合动态带来的严重影响。针对这类复杂问题，许多学者也提出了相应的控制方法，比较常见的有H_∞控制、滑模变结构控制等。但是H_∞控制为典型的干扰抑制方法，它把所有的干扰等价为单一的范数有界干扰，致使控制精度不高。而滑模变结构控制具有鲁棒性、快速性等优势，但是其控制率一般较为复杂，并且存在抖振现象，不利于实际的工程应用。另外，现阶段大多数姿态控制方法都通过假定系统的全部状态可测来设计控制器。但是，在实际的航天工程中难以满足上述假设，如在部分测量敏感器出现故障的情况下。而且由于角速度测量成本等原因，越来越多的学者在研究无角速度量测的情况下的姿态控制问题，实际上也是一种输出反馈的姿态控制问题。因此，基于输出反馈设计高精度的姿态控制器更具备理论及实用意义。早期的工作也有利用DOBC来估计刚柔耦合带来的干扰，但是前提都是假设系统的状态变量全部可测，并且假设干扰的变化率有界。

基于干扰观测器的控制(DOBC)充分利用了干扰的特性，实现了干扰的高精度估计与补偿，并且易于与其它控制方式相结合，因此，可以基于干扰观测器设计精细的复合姿态控制器来同时实现干扰的补偿与抑制，从而提升挠性航天器的姿态控制精度。

因此，需要一种能有效地补偿刚柔耦合干扰和抑制外部环境干扰的基于输出反馈的挠性航天器精细姿态控制方法。

发明内容

本发明的目的在于提供一种基于输出反馈的挠性航天器精细姿态控制方法，所述控制方法包括如下步骤：

a)搭建挠性航天器动力学系统Σ₁，将挠性航天器动力学系统Σ₁变换为挠性航天器动力学系统Σ₂，所述挠性航天器动力学系统Σ₂引入航天器刚柔耦合动态干扰；

b)构造外部系统Σ₃，所述刚柔耦合动态干扰通过外部系统Σ₃进行描述；

c)设计干扰观测器，所述干扰观测器对所述刚柔耦合动态干扰d₀进行估值

d)设计动态输出反馈H_∞控制器；

e)将步骤c)中所诉的干扰观测器与步骤d)中所述的动态输出反馈H_∞控制器进行复合，得到挠性航天器精细姿态控制系统Σ₆；

所述挠性航天器精细姿态控制系统Σ₆通过估值对刚柔耦合动态干扰d₀进行补偿。

优选地，所述控制方法，所述挠性航天器动力学系统Σ₁表述为：

${\mathrm{\Σ}}_1:\left\{\begin{array}{c}J\stackrel{\·\·}{\mathrm{\θ}}+F\stackrel{\·\·}{\mathrm{\η}}=u+d_1\\ \stackrel{\·\·}{\mathrm{\η}}+C_d\stackrel{\·}{\mathrm{\η}}+\mathrm{\Λ}\mathrm{\η}+F^T\stackrel{\·\·}{\mathrm{\θ}}=0\end{array}\right.$

其中，θ为航天器姿态角，为航天器姿态角θ的二阶导数，J为航天器转动惯量，F为挠性附件与本体的刚柔耦合矩阵，F^T为所述刚柔耦合矩阵的转置矩阵，u为控制输入，d₁为环境干扰力矩，η为挠性附件的模态，为挠性附件模态η的导数，为挠性附件模态η的二阶导数，C_d为模态阻尼矩阵，Λ为刚度矩阵。

优选地，所述控制方法，所述模态阻尼矩阵C_d表示为C_d＝diag{2ζ_iω_i}(i＝1,2,...N)，其中N为模态阶数，ζ_i为模态阻尼，ω_i为模态频率，所述刚度矩阵Λ表示为

优选地，所述控制方法，所述的系统Σ₂表述为：

${\mathrm{\Σ}}_2:(J-{\mathrm{FF}}^T)\stackrel{\·\·}{\mathrm{\θ}}=F(C_d\stackrel{\·}{\mathrm{\η}}+\mathrm{\Λ}\mathrm{\η})+u+d_1$

其中，为航天器姿态角θ的二阶导数，J为航天器转动惯量，F为挠性附件与本体的刚柔耦合矩阵，F^T为所述刚柔耦合矩阵的装置矩阵，u为控制输入，d₁为环境干扰力矩，η为挠性附件的模态，为挠性附件模态η的导数，C_d为模态阻尼矩阵，Λ为刚度矩阵。

优选地，所述控制方法，步骤b)中所述刚柔耦合动态干扰d₀表示为其中，F为挠性附件与本体的刚柔耦合矩阵，C_d为模态阻尼矩阵，Λ为刚度矩阵，η为挠性附件的模态，为挠性附件模态η的导数；

所述外部系统Σ₃对刚柔耦合动态干扰d₀描述为：

${\mathrm{\Σ}}_3:\left\{\begin{array}{c}\stackrel{\·}{w}=Ww+H(u+d_1)\\ d_0=Vw\end{array}\right.$

其中，V＝[FΛ FC_d]w为外部系统Σ₃的干扰状态变量，为w的导数，W、H和V为定义的系数矩阵，I为单位矩阵。

优选地，所述控制方法，所述系数矩阵W中，矩阵M表示为M＝I-F^TJ^-1F，其中，I为单位矩阵。

优选地，所述控制方法，步骤c)中所述干扰观测器的设计步骤如下：

(1)构造航天器姿态角输入矩阵

(2)将所述系统Σ₂转化为状态空间系统Σ₄，所述状态空间系统Σ₄表述为：

${\mathrm{\Σ}}_4:\left\{\begin{array}{c}\stackrel{\·}{x}=Ax+B(u+d_0+d_1)\\ y=Cx\end{array}\right.$

其中，A和B为系数矩阵，y为量测输出，C为量测矩阵；

(3)借助量测输出y设计干扰观测器，所述干扰观测器表述为：

$\left\{\begin{array}{c}\widehat{d_0}=V\hat{w}\\ \hat{w}=v-Ly\\ \stackrel{\·}{v}=(W+LVBV)\hat{w}+(H+LCB)u\end{array}\right.$

其中，是刚柔耦合动态干扰d₀的估计值，v为辅助变量，为辅助变量v的导数，y为量测输出，L为干扰观测器增益矩阵；

所述干扰控制器的观测误差动态Σ₅表述为：

${\mathrm{\Σ}}_5:\stackrel{\·}{e_w}=(W+LCBV)e_w+LVAx+(LVB+H)d_1$

其中， 为干扰状态变量w的估计值，为e_w的导数。

优选地，所述控制方法，步骤d)中所述动态输出反馈H_∞控制器表述为：

$\left\{\begin{array}{c}{\stackrel{\·}{x}}_k=A_k{x}_k+B_{k}y\\ u_1=C_k{x}_k+D_{k}y\end{array}\right.$

其中，u₁为动态输出反馈H_∞控制器输入，x_k为控制器状态，A_x、B_x、C_x、D_x为待定的控制器参数矩阵，所述动态输出反馈H_∞控制器对环境干扰进行抑制。

优选地，所述控制方法，步骤e)中所述的挠性航天器精细姿态控制系统Σ₆表述为：

${\mathrm{\Σ}}_6:\left\{\begin{array}{c}{\stackrel{\·}{x}}_k=A_k{x}_k+B_{k}y\\ u=C_k{x}_k+D_{k}y-\widehat{d_0}\end{array}\right.$

其中，u为控制输入，为所述干扰观测器对所述刚柔耦合动态干扰d₀的估计值；

所述精细姿态控制系统Σ₆通过估计值对刚柔耦合动态干扰d₀进行补偿。

优选地，所述控制方法，所述干扰控制器增益矩阵L和所述待定控制器参数矩阵A_x、B_x、C_x、D_x通过凸优化算法求解，所述求解过程如下：

联立系统Σ₄、Σ₅和Σ₆，得到：

$\left[\begin{array}{c}\stackrel{\·}{x}\\ \stackrel{\·}{x_k}\\ \stackrel{\·}{e_w}\end{array}\right]=\left[\begin{array}{ccc}A+{\mathrm{BD}}_{k}C& {\mathrm{BC}}_k& BV\\ B_k& A_k& 0\\ LCA& 0& W+LCBV\end{array}\right]\left[\begin{array}{c}x\\ x_k\\ e_w\end{array}\right]+\left[\begin{array}{c}B\\ 0\\ LCB+H\end{array}\right]d_1$

求解如下凸优化问题：

$\begin{array}{cc}\left[\begin{array}{cccccc}{\mathrm{\&Phi;}}_{11}& {\mathrm{\&Phi;}}_{12}& {\mathrm{\&Phi;}}_{13}& {\mathrm{\&Phi;}}_{14}& {\mathrm{\&Phi;}}_{15}& {\mathrm{\&Phi;}}_{16}\\ *& {\mathrm{\&Phi;}}_{22}& {\mathrm{\&Phi;}}_{23}& 0& 0& 0\\ *& *& -{\mathrm{\&gamma;}}^{2}I& 0& 0& 0\\ *& *& *& -I& 0& 0\\ *& *& *& *& -I& 0\\ *& *& *& *& *& -I\end{array}\right]<0& \left[\begin{array}{cc}{P}_2& S\\ S& S\end{array}\right]>0\end{array}$

其中，

${\mathrm{\&Phi;}}_{11}=\left[\begin{array}{cc}{P}_{2}A+B_{c}C& A_c\\ SA+B_{c}C& A_c\end{array}\right]+{\left[\begin{array}{cc}{P}_{2}A+B_{c}C& A_c\\ SA+B_{c}C& A_c\end{array}\right]}^T,{\mathrm{\&Phi;}}_{12}=\left[\begin{array}{c}{P}_{2}BV+{\left(YCA\right)}^T\\ SBV\end{array}\right],$

${\mathrm{\&Phi;}}_{13}=\left[\begin{array}{c}{P}_{2}B\\ SB\end{array}\right],{\mathrm{\&Phi;}}_{14}=\left[\begin{array}{c}{C}_1^T\\ 0\end{array}\right],{\mathrm{\&Phi;}}_{15}=\left[\begin{array}{cc}{P}_2& 0\\ 0& S\end{array}\right],{\mathrm{\&Phi;}}_{16}=\left[\begin{array}{cc}{\mathrm{BD}}_{c}C& {\mathrm{BC}}_c\\ {\mathrm{BD}}_{c}C& {\mathrm{BC}}_c\end{array}\right],$

Φ₂₂＝Q₁W+YCBV+(Q₁W+YCBV)^T，Φ₂₃＝(YCB+Q₁H)，

得到未知参数：

$L=Q_1^{-1}Y,\left[\begin{array}{cc}{A}_k& B_k\\ C_k& D_k\end{array}\right]=\left[\begin{array}{cc}{S}^{-1}& 0\\ 0& I\end{array}\right]\left[\begin{array}{cc}{A}_c& B_c\\ C_c& D_c\end{array}\right].$

本发明提供的一种基于输出反馈的挠性航天器精细姿态控制方法，针对刚柔耦合动态建立了干扰的模型，充分利用了干扰的信息，采用输出反馈来进行控制器及干扰观测器的设计，相比于状态反馈的控制方法，更具备工程价值。

应当理解，前述大体的描述和后续详尽的描述均为示例性说明和解释，并不应当用作对本发明所要求保护内容的限制。

附图说明

参考随附的附图，本发明更多的目的、功能和优点将通过本发明实施方式的如下描述得以阐明，其中：

图1示意性示出本发明一种基于输出反馈的挠性航天器精细姿态控制方法的设计流程图；

图2示出了本发明实施例中挠性航天器精细姿态控制的模块框图。

具体实施方式

通过参考示范性实施例，本发明的目的和功能以及用于实现这些目的和功能的方法将得以阐明。然而，本发明并不受限于以下所公开的示范性实施例；可以通过不同形式来对其加以实现。说明书的实质仅仅是帮助相关领域技术人员综合理解本发明的具体细节。

在下文中，将参考附图描述本发明的实施例。在附图中，相同的附图标记代表相同或类似的部件，或者相同或类似的步骤。

本发明提供了一种基于输出反馈的挠性航天器精细姿态控制方法，如图1所示本发明基于输出反馈的挠性航天器精细姿态控制方法的设计流程图，本实施例中航天器精细姿态控制方法100，通过将航天器刚柔耦合动态干扰引入到挠性航天器动力学系统，借组外部系统对航天器刚柔耦合动态干扰进行描述。设计干扰观测器对刚柔耦合动态干扰进行估值，设计动态输出反馈H_∞控制器对环境干扰进行抑制，并将干扰观测器和动态输出反馈H_∞控制器进行复合通过刚柔耦合动态干扰的估计值对刚柔耦合动态干扰进行补偿。

出于说明性的目的，本发明所提供的基于输出反馈的挠性航天器精细姿态控制方法通过不同模块实现，如图2所示本发明实施例中挠性航天器精细姿态控制的模块框图200，具体地，所述的模块包括但不限于：那外部系统模块201、干扰观测模块202、动态输出反馈模203块和精细姿态控制模块204。

所述外部系统模块201中，外部系统对航天器刚柔耦合动态干扰进行描述。

所述干扰观测模块202中，干扰观测器对航天器刚柔耦合动态干扰进行估值。

所述动态输出反馈模块203中，动态输出反馈H_∞控制器对环境干扰进行抑制。

所述精细姿态控制模块204中，精细姿态控系统通过航天器刚柔耦合动态干扰的估计值完成对航天器刚柔耦合动态干扰进行补偿。

下面结合图1具体描述本实施例中航天器姿态控制方法：

步骤S101：将航天器刚柔耦合干扰引入到挠性航天器动力学系统：

搭建挠性航天器动力学系统Σ₁，所述系统Σ₁表述为：

其中，θ为航天器姿态角，为航天器姿态角θ的二阶导数，J为航天器转动惯量，F为挠性附件与本体的刚柔耦合矩阵，F^T为所述刚柔耦合矩阵的转置矩阵，u为控制输入，d₁为环境干扰力矩，η为挠性附件的模态，为挠性附件模态η的导数，为挠性附件模态η的二阶导数，C_d为模态阻尼矩阵，Λ为刚度矩阵。模态阻尼矩阵C_d表示为C_d＝diag{2ζ_iω_i}(i＝1,2,...N)，其中N为模态阶数，ζ_i为模态阻尼，ω_i为模态频率；刚度矩阵Λ表示为

对挠性航天器动力学系统Σ₁通过数学变换转变为系统Σ₂，系统Σ₂表述为：系统Σ₂中为挠性航天器刚柔耦合动态干扰。

步骤S102：将刚柔耦合动态干扰通过外部系统Σ ₃ 进行描述

令用d₀表示刚柔耦合动态干扰，利用外部系统Σ₃对刚柔耦合动态干扰d₀进行描述：

${\mathrm{\&Sigma;}}_3:\left\{\begin{array}{c}\stackrel{\&CenterDot;}{w}=Ww+H(u+d_1)\\ d_0=Vw\end{array}\right.$

其中，V＝[FΛ FC_d]w为外部系统Σ₃的干扰状态变量，为w的导数，W、H和V为定义的系数矩阵，数矩阵W中，矩阵M表示为M＝I-F^TJ^-1F，I为单位矩阵。

步骤S103：设计干扰观测器

本实施例步骤S101中在航天器动力学系统中引入了刚柔耦合动态干扰，需要对引入的刚柔耦合动态干扰进行估值。本发明采用干扰观测器对刚柔耦合动态干扰进行估值，具体的干扰参测器的设计在下面具体说明。

干扰观测器的设计步骤具体如下：

(1)构造航天器姿态角输入矩阵

(2)将系统Σ₂转化为状态空间系统Σ₄，所述状态空间系统Σ₄表述为：

${\mathrm{\&Sigma;}}_4:\left\{\begin{array}{c}\stackrel{\&CenterDot;}{x}=Ax+B(u+d_0+d_1)\\ y=Cx\end{array}\right.$

其中，A和B为系数矩阵，y为量测输出，C为量测矩阵；

(3)借助量测输出y设计干扰观测器，干扰观测器具体通过下式表述：

$\left\{\begin{array}{c}\widehat{d_0}=V\hat{w}\\ \hat{w}=v-Ly\\ \stackrel{\&CenterDot;}{v}=(W+LVBV)\hat{w}+(H+LCB)u\end{array}\right.---\left(1\right)$

其中，是刚柔耦合动态干扰d₀的估计值，v为辅助变量，为辅助变量v的导数，y为量测输出，L为干扰观测器增益矩阵；干扰控制器的观测误差动态Σ₅表述为：

${\mathrm{\&Sigma;}}_5:\stackrel{\&CenterDot;}{e_w}=(W+LCBV)e_w+LVAx+(LVB+H)d_1$

其中， 为干扰状态变量w的估计值，为e_w的导数。通过式(1)设计的干扰观测器对刚柔耦合动态干扰d₀进行估计值

步骤S104：设计动态输出反馈H _∞ 控制器：

动态输出反馈H_∞控制器具体通过下式表述：

$\left\{\begin{array}{c}{\stackrel{\&CenterDot;}{x}}_k=A_k{x}_k+B_{k}y\\ u_1=C_k{x}_k+D_{k}y\end{array}\right.---\left(2\right)$

其中，u₁为动态输出反馈H_∞控制器输入，x_k为控制器状态，A_x、B_x、C_x、D_x为待定的控制器参数矩阵，所述动态输出反馈H_∞控制器对环境干扰进行抑制。

步骤S105：干扰观测器与动态输出反馈H _∞ 控制器进行复合

为了清楚的说明航天器的精细姿态控制，在航天器动力学系统中引入的刚柔耦合动态干扰需要通过估计值进行补偿。本发明中采用精细姿态控系统将干扰观测器与动态输出反馈H_∞控制器进行复合对刚柔耦合动态干扰的估计值进行补偿。

本实施例中，具体地，将干扰观测器的表达式(1)与动态输出反馈H_∞控制器的表达式(2)进行复合，得到挠性航天器精细姿态控制系统Σ₆，挠性航天器精细姿态控制系统Σ₆具体表述为：

${\&Sigma;}_6:\left\{\begin{array}{c}{\stackrel{\&CenterDot;}{x}}_k=A_k{x}_k+B_{k}y\\ u=C_k{x}_k+D_{k}y-\widehat{d_0}\end{array}\right.$

其中，u为控制输入，为干扰观测器对刚柔耦合动态干扰d₀的估计值。

上述挠性航天器精细姿态控制系统Σ₆中，控制输入u在动态输出反馈H_∞控制器输入u₁的基础上减掉干扰观测器对刚柔耦合动态干扰d₀的估计值干扰观测器对刚柔耦合动态干扰d₀的估计值近似于刚柔耦合动态干扰d₀时，即实现挠性航天器精细姿态控制系统Σ₆对刚柔耦合动态干扰通过刚柔耦合动态干扰d₀的估计值进行补偿。

干扰控制器增益矩阵L和所述待定控制器参数矩阵A_x、B_x、C_x、D_x通过凸优化算法求解，求解过程如下：

联立系统Σ₄、Σ₅和Σ₆，得到：

$\left[\begin{array}{c}\stackrel{\&CenterDot;}{x}\\ \stackrel{\&CenterDot;}{x_k}\\ \stackrel{\&CenterDot;}{e_w}\end{array}\right]=\left[\begin{array}{ccc}A+{\mathrm{BD}}_{k}C& {\mathrm{BC}}_k& BV\\ B_k& A_k& 0\\ LCA& 0& W+LCBV\end{array}\right]\left[\begin{array}{c}x\\ x_k\\ e_w\end{array}\right]+\left[\begin{array}{c}B\\ 0\\ LCB+H\end{array}\right]d_1$

求解如下凸优化问题：

$\begin{array}{cc}\left[\begin{array}{cccccc}{\mathrm{\&Phi;}}_{11}& {\mathrm{\&Phi;}}_{12}& {\mathrm{\&Phi;}}_{13}& {\mathrm{\&Phi;}}_{14}& {\mathrm{\&Phi;}}_{15}& {\mathrm{\&Phi;}}_{16}\\ *& {\mathrm{\&Phi;}}_{22}& {\mathrm{\&Phi;}}_{23}& 0& 0& 0\\ *& *& -{\mathrm{\&gamma;}}^{2}I& 0& 0& 0\\ *& *& *& -I& 0& 0\\ *& *& *& *& -I& 0\\ *& *& *& *& *& -I\end{array}\right]<0& \left[\begin{array}{cc}{P}_2& S\\ S& S\end{array}\right]>0\end{array}$

其中，

${\mathrm{\&Phi;}}_{11}=\left[\begin{array}{cc}{P}_{2}A+B_{c}C& A_c\\ SA+B_{c}C& A_c\end{array}\right]+{\left[\begin{array}{cc}{P}_{2}A+B_{c}C& A_c\\ SA+B_{c}C& A_c\end{array}\right]}^T,{\mathrm{\&Phi;}}_{12}=\left[\begin{array}{c}{P}_{2}BV+{\left(YCA\right)}^T\\ SBV\end{array}\right],$

${\mathrm{\&Phi;}}_{13}=\left[\begin{array}{c}{P}_{2}B\\ SB\end{array}\right],{\mathrm{\&Phi;}}_{14}=\left[\begin{array}{c}{C}_1^T\\ 0\end{array}\right],{\mathrm{\&Phi;}}_{15}=\left[\begin{array}{cc}{P}_2& 0\\ 0& S\end{array}\right],{\mathrm{\&Phi;}}_{16}=\left[\begin{array}{cc}{\mathrm{BD}}_{c}C& {\mathrm{BC}}_c\\ {\mathrm{BD}}_{c}C& {\mathrm{BC}}_c\end{array}\right],$

Φ₂₂＝Q₁W+YCBV+(Q₁W+YCBV)^T，Φ₂₃＝(YCB+Q₁H)，

得到未知参数：

$L=Q_1^{-1}Y,\left[\begin{array}{cc}{A}_k& B_k\\ C_k& D_k\end{array}\right]=\left[\begin{array}{cc}{S}^{-1}& 0\\ 0& I\end{array}\right]\left[\begin{array}{cc}{A}_c& B_c\\ C_c& D_c\end{array}\right].$

结合这里披露的本发明的说明和实践，本发明的其他实施例对于本领域技术人员都是易于想到和理解的。说明和实施例仅被认为是示例性的，本发明的真正范围和主旨均由权利要求所限定。

Claims (10)

1.一种基于输出反馈的挠性航天器精细姿态控制方法，其特征在于，所述控制方法包括如下步骤：

a)搭建挠性航天器动力学系统Σ₁，将挠性航天器动力学系统Σ₁变换为挠性航天器动力学系统Σ₂，所述挠性航天器动力学系统Σ₂引入航天器刚柔耦合动态干扰；

b)构造外部系统Σ₃，所述刚柔耦合动态干扰通过外部系统Σ₃进行描述；

c)设计干扰观测器，所述干扰观测器对所述刚柔耦合动态干扰进行估值；

d)设计动态输出反馈H_∞控制器；

e)将步骤c)中所诉的干扰观测器与步骤d)中所述的动态输出反馈H_∞控制器进行复合，得到挠性航天器精细姿态控制系统Σ₆；

所述挠性航天器精细姿态控制系统Σ₆通过估值对刚柔耦合动态干扰进行补偿。

2.根据权利要求1所述的控制方法，其特征在于，所述挠性航天器动力学系统Σ₁表述为：

${\mathrm{\&Sigma;}}_1:\left\{\begin{array}{c}J\stackrel{\&CenterDot;\&CenterDot;}{\mathrm{\&theta;}}+F\stackrel{\&CenterDot;\&CenterDot;}{\mathrm{\&eta;}}=u+d_1\\ \stackrel{\&CenterDot;\&CenterDot;}{\mathrm{\&eta;}}+C_d\stackrel{\&CenterDot;}{\mathrm{\&eta;}}+\mathrm{\&Lambda;}\mathrm{\&eta;}+F^T\stackrel{\&CenterDot;\&CenterDot;}{\mathrm{\&theta;}}=0\end{array}\right.$

其中，θ为航天器姿态角，为航天器姿态角θ的二阶导数，J为航天器转动惯量，F为挠性附件与本体的刚柔耦合矩阵，F^T为所述刚柔耦合矩阵的转置矩阵，u为控制输入，d₁为环境干扰力矩，η为挠性附件的模态，为挠性附件模态η的导数，为挠性附件模态η的二阶导数，C_d为模态阻尼矩阵，Λ为刚度矩阵。

3.根据权利要求2所述的控制方法，其特征在于，所述模态阻尼矩阵C_d表示为C_d＝diag{2ζ_iω_i}(i＝1,2,...N)，其中N为模态阶数，ζ_i为模态阻尼，ω_i为模态频率，所述刚度矩阵Λ表示为

4.根据权利要求1所述的控制方法，其特征在于，所述的系统Σ₂表述为：

${\mathrm{\&Sigma;}}_2:(J-{\mathrm{FF}}^T)\stackrel{\&CenterDot;\&CenterDot;}{\mathrm{\&theta;}}=F(C_d\stackrel{\&CenterDot;}{\mathrm{\&eta;}}+\mathrm{\&Lambda;}\mathrm{\&eta;})+u+d_1$

其中，为航天器姿态角θ的二阶导数，J为航天器转动惯量，F为挠性附件与本体的刚柔耦合矩阵，F^T为所述刚柔耦合矩阵的装置矩阵，u为控制输入，d₁为环境干扰力矩，η为挠性附件的模态，为挠性附件模态η的导数，C_d为模态阻尼矩阵，Λ为刚度矩阵。

5.根据权利要求1所述的控制方法，其特征在于，步骤b)中所述刚柔耦合动态干扰d₀表示为其中，F为挠性附件与本体的刚柔耦合矩阵，C_d为模态阻尼矩阵，Λ为刚度矩阵，η为挠性附件的模态，为挠性附件模态η的导数；

所述外部系统Σ₃对刚柔耦合动态干扰d₀描述为：

${\mathrm{\&Sigma;}}_3:\left\{\begin{array}{c}\stackrel{\&CenterDot;}{w}=Ww+H(u+d_1)\\ d_0=Vw\end{array}\right.$

其中，V＝[FΛ FC_d]w为外部系统Σ₃的干扰状态变量，为w的导数，W、H和V为定义的系数矩阵，I为单位矩阵。

6.根据权利要求5所述的控制方法，其特征在于，所述系数矩阵W中，矩阵M表示为M＝I-F^TJ^-1F，其中，I为单位矩阵。

7.根据权利要求1所述的控制方法，其特征在于，步骤c)中所述干扰观测器的设计步骤如下：

(1)构造航天器姿态角输入矩阵

(2)将所述系统Σ₂转化为状态空间系统Σ₄，所述状态空间系统Σ₄表述为：

${\mathrm{\&Sigma;}}_4:\left\{\begin{array}{c}\stackrel{\&CenterDot;}{x}=Ax+B(u+d_0+d_1)\\ y=Cx\end{array}\right.$

其中，A和B为系数矩阵，y为量测输出，C为量测矩阵；

(3)借助量测输出y设计干扰观测器，所述干扰观测器表述为：

$\left\{\begin{array}{c}\widehat{d_0}=V\hat{w}\\ \hat{w}=v-Ly\\ \stackrel{\&CenterDot;}{v}=(W+LVBV)\hat{w}+(H+LCB)u\end{array}\right.$

其中，是刚柔耦合动态干扰d₀的估计值，v为辅助变量，为辅助变量v的导数，y为量测输出，L为干扰观测器增益矩阵；

所述干扰控制器的观测误差动态Σ₅表述为：

${\mathrm{\&Sigma;}}_5:\stackrel{\&CenterDot;}{e_w}=(W+LCBV)e_w+LVAx+(LVB+H)d_1$

其中， 为干扰状态变量w的估计值，为e_w的导数。

8.根据权利要求1所述的控制方法，其特征在于，步骤d)中所述动态输出反馈H_∞控制器表述为：

$\left\{\begin{array}{c}{\stackrel{\&CenterDot;}{x}}_k=A_k{x}_k+B_{k}y\\ u_1=C_k{x}_k+D_{k}y\end{array}\right.$

其中，u₁为动态输出反馈H_∞控制器输入，x_k为控制器状态，A_x、B_x、C_x、D_x为待定的控制器参数矩阵，所述动态输出反馈H_∞控制器对环境干扰进行抑制。

9.根据权利要求1所述的控制方法，其特征在于，步骤e)中所述的挠性航天器精细姿态控制系统Σ₆表述为：

${\mathrm{\&Sigma;}}_6:\left\{\begin{array}{c}{\stackrel{\&CenterDot;}{x}}_k=A_k{x}_k+B_{k}y\\ u=C_k{x}_k+D_{k}y-\widehat{d_0}\end{array}\right.$

其中，u为控制输入，为所述干扰观测器对所述刚柔耦合动态干扰d₀的估计值。

10.根据权利要求1、5、7、8或9所述的控制方法，其特征在于，所述干扰控制器增益矩阵L和所述待定控制器参数矩阵A_x、B_x、C_x、D_x通过凸优化算法求解，所述求解过程如下：

联立系统Σ₄、Σ₅和Σ₆，得到：

$\left[\begin{array}{c}\stackrel{\&CenterDot;}{x}\\ \stackrel{\&CenterDot;}{x_k}\\ \stackrel{\&CenterDot;}{e_w}\end{array}\right]=\left[\begin{array}{ccc}A+{\mathrm{BD}}_{k}C& {\mathrm{BC}}_k& BV\\ B_{k}C& A_k& 0\\ LCA& 0& W+LCBV\end{array}\right]\left[\begin{array}{c}x\\ x_k\\ e_w\end{array}\right]+\left[\begin{array}{c}B\\ 0\\ LCB+H\end{array}\right]d_1$

求解如下凸优化问题：

$\begin{array}{cc}\left[\begin{array}{cccccc}{\mathrm{\&Phi;}}_{11}& {\mathrm{\&Phi;}}_{12}& {\mathrm{\&Phi;}}_{13}& {\mathrm{\&Phi;}}_{14}& {\mathrm{\&Phi;}}_{15}& {\mathrm{\&Phi;}}_{16}\\ *& {\mathrm{\&Phi;}}_{22}& {\mathrm{\&Phi;}}_{23}& 0& 0& 0\\ *& *& -{\mathrm{\&gamma;}}^{2}I& 0& 0& 0\\ *& *& *& -I& 0& 0\\ *& *& *& *& -I& 0\\ *& *& *& *& *& -I\end{array}\right]<0& \left[\begin{array}{cc}{P}_2& S\\ S& S\end{array}\right]>0\end{array}$

其中，

${\mathrm{\&Phi;}}_{11}=\left[\begin{array}{cc}{P}_{2}A+B_{c}C& A_c\\ SA+B_{c}C& A_c\end{array}\right]+{\left[\begin{array}{cc}{P}_{2}A+B_{c}C& A_c\\ SA+B_{c}C& A_c\end{array}\right]}^T,{\mathrm{\&Phi;}}_{12}=\left[\begin{array}{c}{P}_{2}BV+{\left(YCA\right)}^T\\ SBV\end{array}\right],$

${\mathrm{\&Phi;}}_{13}=\left[\begin{array}{c}{P}_{2}B\\ SB\end{array}\right],{\mathrm{\&Phi;}}_{14}=\left[\begin{array}{c}{C}_1^T\\ 0\end{array}\right],{\mathrm{\&Phi;}}_{15}=\left[\begin{array}{cc}{P}_2& 0\\ 0& S\end{array}\right],{\mathrm{\&Phi;}}_{16}=\left[\begin{array}{cc}{\mathrm{BD}}_{c}C& {\mathrm{BC}}_c\\ {\mathrm{BD}}_{c}C& {\mathrm{BC}}_c\end{array}\right],$

Φ₂₂＝Q₁W+YCBV+(Q₁W+YCBV)^T，Φ₂₃＝(YCB+Q₁H)，

得到未知参数：

$L=Q_1^{-1}Y,\left[\begin{array}{cc}{A}_k& B_k\\ C_k& D_k\end{array}\right]=\left[\begin{array}{cc}{S}^{-1}& 0\\ 0& I\end{array}\right]\left[\begin{array}{cc}{A}_c& B_c\\ C_c& D_c\end{array}\right].$