[发明专利]一种基于增量式偏流角的动中成像卫星姿态控制方法

权利要求书

1.一种基于增量式偏流角的动中成像卫星姿态控制方法，其特征在于包括如下步骤：

（1）获取卫星当前时刻的目标滚动角φ_r、目标俯仰角θ_r和目标偏航角ψ_r，建立参考姿态坐标系；所述的参考姿态坐标系与轨道坐标系之间的转换矩阵用C_RO来表示，

CRO=cosψrcosθrcosψrsinθrsinφr+sinψrcosφr-cosψrsinθrcosφr+sinψrsinφr-sinψrcosθr-sinψrsinθrsinφr+cosψrcosφrsinψrsinθrcosφr+cosψrsinφrsinθr-cosθrsinφrcosθrcosφr]]>其中，轨道坐标系O_OX_OY_OZ_O的原点O_O在卫星质心，卫星的轨道平面是坐标平面，由卫星质心指向地心的坐标轴是Z_O轴，X_O轴在轨道平面上与Z_O轴垂直并指向卫星速度方向，Y_O轴与X_O轴、Z_O轴组成右手正交坐标系；

（2）根据卫星的轨道根数，获取轨道坐标系下卫星运行速度的前向分量径向分量由此得到参考姿态坐标系下卫星的运行速度vsR=CRO·vuo0vro;]]>其中μ为地球引力常数，a为卫星轨道半长轴，r为卫星的地心距，e为卫星轨道偏心率，f为真近点角；

（3）定义新的惯性坐标系I'，所述I'中的坐标(X',Y',Z')与J2000惯性坐标系I中的坐标(X,Y,Z)满足关系式

X′Y′Z′=CI′IXYZ]]>

其中，CI′I=10001000K,]]>K=1.0033633486；

（4）计算获得在新的惯性坐标系I'下，卫星上搭载的有效载荷光学系统的光轴矢量在J2000惯性坐标系I下的表达形式其中为J2000惯性坐标系到轨道坐标系的姿态转换矩阵的转置，s=CRO(3,1)CRO(3,2)CRO(3,3)]]>为参考姿态坐标系下的光轴指向单位矢量；R_e为地球半径，rsI=rcos(Ω)cos(u)-sin(Ω)cos(i)sin(u)sin(Ω)cos(u)+cos(Ω)cos(i)sin(u)sin(i)sin(u),]]>Ω为卫星轨道的升交点赤经，u为卫星轨道幅角，i为卫星轨道倾角，|·|表示取模，sin²γ'_m=1-cos²γ'_m；C_RO(3,1)、C_RO(3,2)、C_RO(3,3)分别表示矩阵C_RO的第3行第1个元素、第3行第2个元素、第3行第3个元素；

（5）根据步骤（4）的结果，计算得到由卫星惯性姿态参考角速度ω_RI引起的地面目标点相对于卫星的运行速度在参考姿态坐标系下的分量

vrsR=ωRI×00|rZbI|;]]>

（6）获取由于地球自转引起的卫星指向地面目标点的线速度在参考姿态坐标系下的分量verR=CROCOIverI,]]>其中verI=00ωe×ReI,]]>ReI=rsI+|rZbI′|·COIT·s|CI′I·COIT·s|,]]>ω_e表示地球自转角速度的大小；

（7）根据步骤（2）、步骤（5）、步骤（6）的结果，计算得到卫星光轴指向地面目标点相对于卫星的线速度在参考姿态坐标系下的表达式为

^Rv_es＝^Rv_er-^Rv_s-^Rv_re；

（8）根据步骤（7）的结果，计算得到所述有效载荷成像时的偏流角增量式中括号中的数字表示选取矢量^Rv_es中的第几个元素参与计算；

（9）取ψ_r=ψ_r+△ψ_p得到更新后的目标偏航角，并进一步得到更新后的参考姿态坐标系与轨道坐标系之间的转换矩阵C_RO新；

（10）利用C_RO新作为卫星姿态矩阵的目标值对卫星的三轴进行姿态控制。