[发明专利]一种基于特征模型的挠性卫星控制方法

权利要求书

1.一种基于特征模型的挠性卫星控制方法，其特征在于通过以下步骤实现：

第一步，利用公式(2)确定挠性卫星动力学方程的时间尺度p，

p=min{1Mf·1,1Mf·,1Mf1,1Mf2,1Mu}---(2)]]>

其中挠性卫星动力学方程为公式组(1)，

φ·θ·ψ·T=C(x1)wxwywzT]]>

Isw·s+w~sIsws+Fslη··l+Fsrη··r=Tsη··l+2ξlwlη·l+wl2η·l+FslTw·s=0η··r+2ξrwrη·r+wr2η·r+FsrTw·s=0y=φθψ---(1)]]>

C(x1)=cosθcosψ0sinθcosψtanψcosθ1tanψsinθ-sinθ0cosθ,]]>w~s=0-wzwywz0-wx-wywx0]]>

φ，θ，ψ表示挠性卫星的俯仰、偏航和滚动姿态角，[w_x w_y w_z]^T表示卫星相对轨道坐标系的角速度在体坐标系中的坐标，w_s、分别表示挠性卫星中心体的角速度列阵和反对称阵，η_l、η_r分别为挠性卫星左、右太阳翼的模态坐标阵，ξ_l、ξ_r分别为挠性卫星左、右太阳翼的模态阻尼系数，F_sl、F_sr分别为挠性卫星左、右太阳翼与中心体的耦合系数，T_s表示作用在挠性卫星上的外力矩列阵，I_s表示挠性卫星惯量阵，x₁＝[φθψ]^T，y表示挠性卫星输出，w_l、w_r分别为挠性卫星左、右太阳翼的角速度，

f₁＝C(x₁)w_s，f2=-Is-1(w~sIsws+Fslη··l+Fsrη··r),]]>g=Is-1,]]>u＝T_s，f＝f₂+gu，

Mf·1=max|f·1|,]]>Mf·=max|f·|,]]>Mf1=max|f1|,]]>Mf2=max|f2|,]]>M_u＝max|gu|；

第二步，利用公式(3)和第一步得到的时间尺度p确定采样时间尺度h，

h=pd,d>5---(3);]]>

第三步，利用公式组(4)和第一步中确定的参数f₁、x₁、x₂、g，得到f_1i关于x_1j、x_2j的偏导数和g的上界M，

|∂f1i∂x1j|≤M]]>

|∂f1i∂x2j|≤M]]>

||g(k)||≤M                                (4)

其中x₂＝w_s，i，j＝1，2，3，k＝1，2，…，f_1i、x_1i、x_2i表示f₁、x₁、x₂的第i行；

第四步，利用公式(5)和第三步得到的M确定参数m，

m>lnd2ϵ6MlnNxCx2+Nx-1---(5)]]>

其中，N_x＞0表示挠性卫星输出|y_i|即y的上界，i＝1，2，3，ε为建模误差，C_x＞0；

第五步，根据第二步确定的采样时间尺度h、第三步确定的M和第四步确定的参数m，利用公式组(6)得到挠性卫星特征模型的系数范围，

|ai1(k)-2|≤Mh+6M(m+1)d2]]>

|ai2(k)+1|≤Mh+6M(m+1)d2---(6)]]>

|b_ij(k)|≤M²h²

其中a_i1、a_i2、b_ij表示特征模型的系数，b_ij∈R，j＝1，2，3，b_i＝[b_i1 b_i2 b_i3]；

第六步，利用梯度法对由第五步得到的特征模型的系数a_i1、a_i2、b_ij进行辨识得到辨识后的特征模型的系数

第七步，利用第六步得到的辨识后的特征模型的系数组成公式(7)的控制率，

u=b^1Tb^2Tb^3T-T×diag[u1,u2,u3]---(7)]]>

其中，u_i＝u_0i+u_Gi+u_Ii+u_Di，u_0i为维持/跟踪控制律，u_Gi为黄金分割控制率，u_Ii为逻辑积分控制率，u_Di为逻辑微分控制率；

第八步，将第七步确定的控制率代入公式组(1)的挠性卫星动力学方程中，控制挠性卫星的俯仰、偏航和滚动姿态角。

2.根据权利要求1所述的一种基于特征模型的挠性卫星控制方法，其特征在于：所述第四步中C_x的确定通过以下步骤完成，

A4.1、h²|h_i|为连续函数，且h²|h_i|≤ε，得到|x_1i|的范围，其中h为采样时间尺度，ε为建模误差；

A4.2、根据A4.1得到的|x_1i|的范围，根据公式|x_1i|≤C_x确定C_x的值。