[发明专利]一种高精度RDSS授时方法

权利要求书

1.一种RDSS授时方法，其特征在于，包括以下步骤：

步骤1、获得整分点时刻的中心站到卫星的上行时延值、整分点时刻卫星位置、卫星速度以及时延修正参数；

步骤2、利用步骤1获得的所述整分点时刻播报的中心站到卫星的上行时延值计算当前修正时刻对应的上行时延，具体步骤如下：

S201：累积步骤1所述的距离当前时刻最新的三个整分点时刻播发的中心站到卫星的上行时延值；假设当前授时卫星转发时刻为t时刻，距离当前时刻最新的三个整分点时刻t0、t1、t2，三个时刻离当前时刻由新及远，所播发的中心站到卫星时延值分别为CS0、CS1和CS2；

S202：根据S201所述的最新三个整分点时刻对应的中心站到卫星时延值推算上行时延的变化速度Vcs：

Vcs0＝(CS0-CS1)/(t0-t1)＝(CS0-CS1)/dT

Vcs1＝(CS1-CS2)/(t1-t2)＝(CS1-CS2)/dT

其中，dT是指RDSS授时播发时间间隔；

S203：根据S202所推算的中心站到卫星上行时延的变化速度Vcs，推算上行时延的加速度：

Acs＝(Vcs0-Vcs1)/dT

S204：根据S202所推算的中心站到卫星上行时延的变化速度Vcs，以及S203所推算的上行时延的加速度，计算当时授时卫星转发时刻的上行时延值CS：

CS＝CS0+(Vcs0+Acs×dt)×dt

CS＝CS0+((CS0-CS1)/dT+(Vcs0-Vcs1)×dt/dT)×dt

CS＝CS0+((CS0-CS1)/dT+((CS0-CS1)/dT-(CS1-CS2)/dT)×dt/dT)×dt

CS＝CS0+(CS0-CS1)×dt/dT+(CS0+CS2-2CS1)×dt²/dT²

其中，dt＝t-t0表示当前授时卫星转发时刻t时刻相对t0时刻的时间差；

步骤3、计算卫星到用户下行时延；

步骤4、利用用户机零值以及步骤1中的所述时延修正参数，计算包括电离层、对流层、Sagnac效应和用户机零值对时延的影响，即得到时延修正值；

步骤5、将RDSS授时信号发射时刻同时加上步骤2获得的上行时延、步骤3获得的下行时延以及步骤4获得的时延修正值，完成经修正的授时计算。

2.如权利要求1所述的一种RDSS授时方法，其特征在于，所述步骤3计算卫星到用户下行时延的步骤如下：

S301：累积步骤1所述的距离当前时刻最新的三个整分点时刻播发的卫星位置；假设当前授时卫星转发时刻为t时刻，距离当前时刻最新的三个整分点时刻t0、t1、t2，三个整分点时刻离当前时刻由新及远，所播发的整分点时刻卫星位置分别为S0(Xs0,Ys0,Zx0)、S1(Xs1,Ys1,Zx1)、S2(Xs2,Ys2,Zx2)；

S302：利用类似步骤2所述推算方法推算当前卫星转发时刻卫星位置：

Xs＝Xs0+(Xs0-Xs1)×dt/dT+(Xs0+Xs2-2 Xs1)×dt²/dT²

Ys＝Ys0+(Ys0-Ys1)×dt/dT+(Ys0+Ys2-2 Ys1)×dt²/dT²

Zs＝Zs0+(Zs0-Zs1)×dt/dT+(Zs0+Zs2-2 Zs1)×dt²/dT²

其中，dT是指RDSS授时播发时间间隔；dt为t时刻相对t0时刻的时间差dt＝t-t0；

S303：利用S302所推算的卫星转发时刻卫星位置的St(Xs,Ys,Zs)以及已知点用户位置U(Xu,Yu,Zu)计算卫星到用户时延SU：

SU＝[(Xs-Xu)²+(Ys-Yu)²+(Zs-Zu)²]^1/2。