一种用于卫星双向时间传递的卫星轨道摄动补偿方法

摘要

本发明公开了一种用于卫星双向时间传递的卫星轨道摄动补偿方法，通过搭建轨道摄动补偿卫星双向时间传递系统、距离计算模块（7）确定地球站到卫星的距离、控制量计算模块（8）确定1PPS相位控制字、脉冲调整模块（9）调节秒脉冲和数据补偿模块（10）补偿1PPS调节操作产生的测量偏差，完成对卫星双向时间传递系统中卫星轨道摄动的补偿。本发明抵消了卫星轨道摄动在正反两个方向上引起的时间信号传输时延的变化，极大地降低卫星轨道摄动的影响，提高了卫星双向实时时间传递的精度。

主权项

一种用于卫星双向时间传递的卫星轨道摄动补偿方法，其特征在于该方法的具体步骤为：第一步  搭建轨道摄动补偿卫星双向时间传递系统轨道摄动补偿卫星双向时间传递系统，包括：卫星双向时间传递系统(1)和轨道摄动补偿装置(6)；其中，卫星双向时间传递系统(1)，包括：主原子钟、主地球站(2)、从原子钟和从地球站(3)；轨道摄动补偿装置(6)，包括：距离计算模块(7)、控制量计算模块(8)、脉冲调整模块(9)和数据补偿模块(10)；轨道摄动补偿卫星双向时间传递系统是在卫星双向时间传递系统(1)从原子钟和从地球站(3)间串联轨道摄动补偿装置(6)来实现的；轨道摄动补偿卫星双向时间传递系统中，主原子钟输出的10MHz和1PPS信号直接和主地球站(2)的参考输入端相连；从原子钟输出的10MHz和1PPS信号则与轨道摄动补偿装置(6)参考输入端相连，轨道摄动补偿装置(6)的1PPS和10MHz输出端口再与从地球站(3)的参考输入端相连，轨道摄动补偿装置(6)数据采集端口与从地球站(3)数据输出端口相连；距离计算模块(7)的功能为：根据本地地球站地理坐标和通信卫星的卫星定点坐标，分别计算主地球站(2)、从地球站(3)到卫星的距离，并将计算结果输入至控制量计算模块(8)；控制量计算模块(8)的功能为：将主地球站(2)、从地球站(3)到卫星的距离求差并除以光速，得到两站触发调制脉冲到达卫星转发器的时间差，减去两站的粗略钟差，得到相位超前滞后控制量，上报给脉冲调整模块(9)；脉冲调整模块(9)的功能为：根据输入的控制量，对本地秒信号进行超前滞后调节，控制量为正则表示滞后调节，控制量为负表示超前调节；数据补偿模块(10)的功能为：根据相位控制量对最终的双向数据进行补偿，得到最终的两站秒脉冲时差数据；第二步  距离计算模块(7)确定地球站到卫星的距离根据实际时间传递需要，将卫星双向时间传递系统主地球站(2)和从地球站(3)布置在适当的位置；利用GPS接收机测得的或测绘部门给出的地球站点大地坐标，通过转换公式(1)得到地球站在WGS84坐标系下的直角坐标；$\left\{\begin{array}{c}x=\left(N+h\right)\cos \phi \cos \lambda \\ y=\left(N+h\right)\cos \phi \sin \lambda \\ z=\left[N\left(1-e^2\right)+h\right]\sin \phi \end{array}\right.---\left(1\right)$公式(1)中，(x,y,z)为站点的直角坐标；e为椭球偏心率；N为卯酉圆曲率半径；φ为纬度；λ为经度；h为海拔；椭球偏心率和卯酉圆曲率半径由椭球体长半轴半径a和短半轴半径b计算得到，公式如下：$\left\{\begin{array}{c}{e}^2=\frac{a^2-b^2}{a^2}\\ N=\frac{a}{\sqrt{1-e^2{\sin }^2\phi }}\end{array}\right.---\left(2\right)$为了表示方便，卫星双向时间传递主、从地球站(3)分别用A站和B站来表示，地球站坐标分别用(x_A,y_A,z_A)，(x_B,y_B,z_B)表示，根据卫星的经度、纬度、高度，可得卫星的直角坐标为(x_S,y_S,z_S)，则A站、B站两地球站到卫星的距离分别为：$\left\{\begin{array}{c}{r}_A=\sqrt{{(x_A-x_S)}^2+{(y_A-y_S)}^2+{(z_A-z_S)}^2}\\ r_B=\sqrt{{(x_B-x_S)}^2+{(y_B-y_S)}^2+{(z_B-z_S)}^2}\end{array}\right.---\left(3\right)$公式(3)中，r_A为A站到卫星的距离；r_B为B站到卫星的距离；第三步  控制量计算模块(8)确定1PPS相位控制字根据两站到卫星的距离差，再确定相位调节控制量；首先测量A站、B站的粗略钟差，采用外接GPS定时接收机与时间间隔计数器进行测量，或者直接建立卫星双向链路进行测量，测得的粗略钟差为Δτ_BA，B站超前时Δτ_BA为正值；考虑脉冲调整模块(9)安装在B站的情况，相位调节量为：${\tau }_B=\frac{r_B-r_A}{c}+{\mathrm{\Delta \tau }}_{BA}---\left(4\right)$公式(4)中，c为光速；当τ_B为正数时，则将1PPS滞后调节τ_B s；当τ_B为负数时，则将1PPS滞后调节1‑τ_B s；脉冲调整模块(9)是以100ns为最小步进进行相位调节，所以实际输送给脉冲调整模块(9)的相位控制字为：公式(5)中，Phase_adjust为相位控制字；[·]表示取整数；第四步  脉冲调整模块(9)调节秒脉冲脉冲调整模块(9)接收从站原子钟(5)输入的10MHz和1PPS信号，根据相位控制字对输入的1PPS进行相位调节，并将调节后的1PPS信号和对应的相参10MHz信号输送给从地球站(3)中；第五步  数据补偿模块(10)补偿1PPS调节操作产生的测量偏差完成相位调整后，将调整后秒脉冲信号输送至卫星双向时间传递系统从地球站(3)，从地球站(3)将测得本地调整脉冲与远程站秒脉冲之时间差，差值为正，则表示从站秒脉冲超前于主站秒脉冲；差值为负值，则表示从站秒脉冲滞后于主站秒脉冲；此时，从地球站(3)给出的测量值为主站原子钟(4)1PPS与经过调节后轨道摄动补偿装置(6)输出1PPS之间的时差，而非主站原子钟(4)和从站原子钟(5)的秒脉冲之间的时差，为了得到主站原子钟(4)、从站原子钟(5)原始秒脉冲之差，需要对测量数据进行补偿，具体补偿方式是在从站测量数据上加上τ_B分量；${\mathrm{\Delta T}}_{BA}=\Delta {\hat{T}}_{BA}+{\tau }_B---\left(6\right)$公式(6)中，ΔT_BA为最终的主从站钟差测量数据，为从站卫星双向时间传递系统测量数据；至此，完成对卫星双向时间传递系统中卫星轨道摄动的补偿。