1.一种适用于伪卫星位置的计算方法,其特征在于,该方法包括以下步骤:判断是否为伪星系统;若是,
读取所述伪星系统的坐标常数;
计算归一化时间;
计算卫星的平均角速度;
分别计算信号发射时刻的平近点角、偏近点角、真近点角、升交点角距、摄动校正项;
计算摄动校正后的升交点角距、卫星矢量长度和轨道倾角;
计算伪卫星在轨道平面的位置;
计算信号发射时刻的升交点赤经;
计算伪卫星在地固坐标系下的位置。
2.根据权利要求1所述的一种适用于伪卫星位置的计算方法,其特征在于,归一化时间tk通过以下公式计算得到:tk=t-toe,toe为参数时间。
3.根据权利要求2所述的一种适用于伪卫星位置的计算方法,其特征在于,计算卫星的平均角速度n通过以下公式计算得到:n=n0+Δn,其中, A=(as)2,A表示物体与参考地固坐标系原点之间的距离,μ表示伪卫星系统引力常数,Δn表示平均运动角速度校正值。
4.根据权利要求3所述的一种适用于伪卫星位置的计算方法,其特征在于,信号发射时刻的平近点角Mk通过以下方法计算:Mk=M0+n·tk,M0为星历参考时间的平近点角。
5.根据权利要求4所述的一种适用于伪卫星位置的计算方法,其特征在于,信号发射时刻的偏近点角Ek通过以下方法计算得到:Mk=Ek-es·sin Ek。
6.根据权利要求5所述的一种适用于伪卫星位置的计算方法,其特征在于,信号发射时刻的真近点角vk通过以下方法计算得到:
7.根据权利要求6所述的一种适用于伪卫星位置的计算方法,其特征在于,信号发射时刻的升交点角距φk通过以下方法计算得到:φk=vk+ω,ω为轨道近地角距。
8.根据权利要求7所述的一种适用于伪卫星位置的计算方法,其特征在于,信号发射时刻的摄动校正项δuk,δrk,δik通过以下方法计算得到:其中,Cus,Cuc,Crs,Crc,Cis,Cic表示星历参数。
9.根据权利要求8所述的一种适用于伪卫星位置的计算方法,其特征在于,所述摄动校正后的升交点角距uk通过以下方法计算:uk=φk+δuk
所述卫星矢量长度rk通过以下方法计算:rk=A(1-es·cosEk)+δrk所述轨道倾角ik通过以下方法计算:
表示轨道倾角对时间的变化率。
10.根据权利要求9所述的一种适用于伪卫星位置的计算方法,其特征在于,信号发射时刻伪卫星在轨道平面的位置(xk,yk)通过以下方法计算:信号发射时刻的升交点赤经Ωk通过以下方法计算:Ω0表示表示周内秒等于0时的轨道升交点赤经;
伪卫星在地心地固坐标直角坐标系中的坐标通过以下方法计算: