[发明专利]一种适用于室内定位导航的伪卫星导航电文模拟方法

权利要求书

1.一种适用于室内定位导航的伪卫星导航电文模拟方法，其特征在于，包括以下步骤：

(1)根据星历参考时间t_oe计算规化时间t_k,以及根据开普勒第三定律计算卫星的平均角速度n；

(2)根据规化时间t_k和平均角速度n，计算信号发射时刻的平近点角M_k；

(3)根据平近点角M_k和星历参数e_s，计算信号发射时刻的偏近点角E_k；

(4)根据偏近点角E_k与星历参数e_s，计算信号发射时刻的真近点角v_k；

(5)根据真近点角v_k，计算信号发射时刻的升交点角距Φ_k；

(6)根据升交点角距Φ_k，计算信号发射时刻的二次谐波摄动校正项δu_k,δr_k,δi_k；

(7)根据二次谐波摄动校正项δu_k,δr_k,δi_k，计算摄动校正后的升交点角距u_k、伪卫星矢径长度r_k和轨道倾角

(8)根据升交点角距u_k和伪卫星矢径长度r_k，计算信号发射时刻伪卫星在轨道平面的位置(x′_k,y′_k)；

(9)计算信号发射时刻的升交点赤经Ω_k，根据升交点赤经Ω_k和伪卫星在轨道平面的位置(x′_k,y′_k)，计算伪卫星在WGS-84地心地固坐标系中的坐标(x_k,y_k,z_k)，得到伪卫星坐标值，以模拟生成伪卫星导航电文。

2.根据权利要求1所述的适用于室内定位导航的伪卫星导航电文模拟方法，其特征在于，所述步骤(1)中的计算规化时间t_k，具体计算公式为，

t_k＝t-t_oe

其中，t_oe为星历参考时间，t_k为相对于t_oe的规化时间，t为计算伪卫星位置时刻，即当前时刻，是根据当前时刻t与星历参考时间t_oe来计算规化时间t_k。当计算得到的t_k大于302400s时，则t_k应该减去604800s，当t_k小于-302400s时，则t_k加上604800s。

所述步骤(1)中的计算卫星的平均角速度n，根据开普勒第三定律不同行星绕太阳运动的公转周期的平方分别与它们的轨道长半径的立方成正比，其公式为，

其中，T为公转周期，a_s为轨道长半径，M为太阳质量，G为引力常数；

令n₀代表卫星的平均角速度，则

由此可得

其中，u为常数，计算可知n₀为1.458555*10^(-4)

校正之后的卫星平均角速度n为

n＝n₀+Δn

将Δn置0处理，可得n＝n₀。

3.根据权利要求2所述的适用于室内定位导航的伪卫星导航电文模拟方法，其特征在于，所述步骤(2)中的根据规化时间t_k和平均角速度n，计算信号发射时刻的平近点角M_k，具体计算公式为，

M_k＝M₀+nt_k

令轨道长半径a_s等于ICD文件所要求的最大值的极限条件，则M_k＝M₀，

其中，M₀为星历给出的平近点角。