[发明专利]一种导航卫星姿态测量方法

权利要求书

1、一种导航卫星姿态测量方法，其特征在于，包括以下步骤：

步骤(1)，接收机观测卫星信号的载波相位、观测点的大地纬度、经度和导航电文；

步骤(2)，解析出载波相位φ，并根据观测点的大地纬度、经度和导航电文解析卫星的高度角α、方位角Ω；

步骤(3)，采用蚁群算法搜索卫星的基线姿态角和β，具体为根据公式：

cos[2π(Nij+δij)]]]>

计算基线姿态角和β；其中，b为基线长度，λ为射频载波的波长，αⁱ和α^j分别是卫星i和j对基线的载波平面的高度角，Ωⁱ和Ω^j分别是卫星i和卫星j对基线的载波平面的方位角，和β则为基线矢量的航向角和俯仰角；N^ij为相位双差整周模糊度，δ^ij为观测噪声，是均值为零的高斯白噪声；

步骤(4)，计算基线矢量b；

步骤(5)，反算双差整周模糊度矩阵N。

2、根据权利要求1所述的一种导航卫星姿态测量方法，其特征在于，步骤(4)中，基线矢量b的计算方法如下：

根据极坐标和直角坐标的关系得出基线在当地水平直角坐标系下的三维分量[X_zs Y_zs Z_zs]：

计算基线在地心坐标系下的三维分量[ΔX ΔY ΔZ]，即基线矢量b：

ΔXΔYΔZ=T×XzsYzsZzs;]]>

T=-sinBpcosLp-sinLpcosBpcosLp-sinBpsinLpcosLpcosBpsinLpcosBp0sinBp;]]>

其中，B_p，L_p为观测点的大地纬度和经度，由接收机直接得到。

3、根据权利要求3所述的一种导航卫星姿态测量方法，其特征在于，步骤(5)中，根据公式N=N^-Qb^N^-1·QN^·(b^-b)]]>计算双差整周模糊度矩阵N；其中b为基线矢量，为双差整周模糊度的实数估计值，和是协方差矩阵中的元素，是基线矢量的实数估计值。