[发明专利]基于GNSS绝对定位数据的卫星姿态导引角计算方法及系统

权利要求书

1.一种基于GNSS绝对定位数据的卫星姿态导引角计算方法，，其特征在于，包括：

步骤S1：根据GNSS实时绝对定位数据中的时间、位置、速度转换到惯性坐标系下，得到惯性坐标系下的位置、速度；

步骤S2：将GNSS实时绝对定位的时间下位置速度数值积分递推至当前星时，得到当前星时时刻的位置和速度；

步骤S3：将当前时刻的位置速度信息转换得到对应时刻的轨道参数；

步骤S4：根据当前星时时刻轨道参数，计算当前时刻偏航导引角和俯仰导引角。

2.根据权利要求1所述的基于GNSS绝对定位数据的卫星姿态导引角计算方法，其特征在于，在所述步骤S1中：

步骤S1.1：

根据时间T₀计算J2000地心惯性系到地固坐标系的坐标变换矩阵：

M_ECI2ECF＝EP·ER·NR·PR

其中，EP为极移矩阵，简化为3×3单位阵；ER为地球自转矩阵，NR为章动矩阵，PR为岁差矩阵；

步骤S1.2：

将GNSS绝对定位数据位置速度信息由WGS84坐标系转换至J2000.0惯性坐标系，得到位置r_ECI0、速度v_ECI0为：

r_ECI0＝(M_ECI2ECF)^-1·r_ECF0

v_ECI0＝(M_ECI2ECF)^-1v_ECF0+ω_e^×r_ECF0

其中，r_ECFO为T₀时刻WGS84坐标系下的位置，v_ECFO为T₀时刻WGS84坐标系下的速度；上表T表示转置，上表·表示对时间求导。

3.根据权利要求1所述的基于GNSS绝对定位数据的卫星姿态导引角计算方法，其特征在于，在所述步骤S2中：

根据T₀时刻位置速度，采用四阶龙格库塔公式，数值积分递推至当前星时T，得到T时刻的位置r_ECI、速度v_ECI。