一种高精度RDSS授时方法与流程
本发明属于rdss授时技术领域，具体涉及一种高精度rdss授时方法。
背景技术：
rdss(radiodeterminationsatelliteservice)提供定位授时服务10余年，广泛应用于通信、电力、商务及国防建设等领域，其授时服务精度为50ns。
rdss授时工作原理如下：不同于gps等卫星导航系统基于星载原子钟时间基准的设计，rdss时间基准采用基于地面的高精度原子钟，这种设计具有更高的时间基准精度。此外，rdss星座采用geo卫星，服务于中国及周边区域。rdss授时信号的传输路径是从地面中心站经卫星转发至地面用户
目前rdss授时接收机利用当前时刻整分点以及前一整分点时刻播报的中心站到卫星时延推算当前时刻的上行时延，同时利用整分点播报的卫星位置和速度推算转发时刻卫星位置，进而计算下行时延；在此基础之上进行大气时延、硬件零值、地球自转修正等误差修正，授时精度约为50ns，随着用户对精度要求的提高，现有的授时精度的劣势也越发明显。
技术实现要素：
有鉴于此，本发明的目的是提供一种高精度rdss授时方法，能够使rdss授时精度得到大幅度提升。
一种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)×dt2/dt2
其中，dt＝t-t0表示当前授时卫星转发时刻t时刻相对t0时刻的时间差；
步骤3、计算卫星到用户下行时延；
步骤4、利用用户机零值以及步骤1中的所述时延修正参数，计算包括电离层、对流层、sagnac效应和用户机零值对时延的影响，即得到时延修正值；
步骤5、将rdss授时信号发射时刻同时加上步骤2获得的上行时延、步骤3获得的下行时延以及步骤4获得的时延修正值，完成经修正的授时计算。
进一步的，所述步骤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-2xs1)×dt2/dt2
ys＝ys0+(ys0-ys1)×dt/dt+(ys0+ys2-2ys1)×dt2/dt2
zs＝zs0+(zs0-zs1)×dt/dt+(zs0+zs2-2zs1)×dt2/dt2
其中，dt是指rdss授时播发时间间隔；dt为t时刻相对t0时刻的时间差dt＝t-t0；
s303：利用s302所推算的卫星转发时刻卫星位置的st(xs,ys,zs)以及已知点用户位置u(xu,yu,zu)计算卫星到用户时延su：
su＝[(xs-xu)2+(ys-yu)2+(zs-zu)2]1/2。
进一步的，所述步骤3计算卫星到用户下行时延的步骤如下：
s311：累积步骤1所述的距离当前时刻最新的三个整分点时刻播发的卫星位置；假设当前授时卫星转发时刻为t时刻，距离当前时刻最新的三个整分点时刻t0、t1、t2，三个整分点时刻离当前时刻由新及远，所播发的整分点时刻卫星位置分别为s0(xs0,ys0,zx0)、s1(xs1,ys1,zx1)、s2(xs2,ys2,zx2)；
s312：利用s311所述的整分点所播发的卫星位置以及已知点用户位置u(xu,yu,zu)计算t0、t1、t2整分点时刻对应的卫星到用户时延：
su0＝[(xs0-xu)2+(ys0-yu)2+(zs0-zu)2]1/2
su1＝[(xs1-xu)2+(ys1-yu)2+(zs1-zu)2]1/2
su2＝[(xs2-xu)2+(ys2-yu)2+(zs2-zu)2]1/2
s313：利用类似步骤2所述推算方法推算当前卫星转发时刻对应的卫星到用户的下行时延su：
su＝su0+(su0-su1)×dt/dt+(su0+su2-2su1)×dt2/dt2
其中，dt是指rdss授时播发时间间隔；dt为t时刻相对t0时刻的时间差dt＝t-t0。
较佳的，当针对北斗二号的授时时，dt＝1min。
较佳的，当针对北斗三号的授时时，dt＝6s。
本发明具有如下有益效果：
本发明的一种高精度rdss授时方法，改进了rdss授时上行时延的推算方法，简单有效，从而使授时服务精度有大幅度的提升；同样方法也可以用于推算卫星转发时刻对应的卫星位置，或者说用于下行时延的推算。通过我们的试验分析及验证，采用本发明的授时方法，可以从当前传统rdss授时接收机的50ns提高到约10ns的授时精度(见附图6)，攻克了rdss授时服务精度的瓶颈，取得了突破性的进展，具有很大的应用潜力和经济效益。
附图说明
图1为本发明的rdss授时信号传输路径图；
图2为本发明的高精度rdss授时方法流程图；
图3为本发明的高精度rdss授时方法中上行时延求解流程图；
图4为本发明的高精度rdss授时方法中下行时延求解方法一流程图；
图5为本发明的高精度rdss授时方法中下行时延求解方法二流程图；
图6为采用本发明方法后的授时误差与当前方法误差结果比较图。
具体实施方式
下面结合附图并举实施例，对本发明进行详细描述。
10年来，传统rdss授时用户机虽然满足了系统指标要求，但是存在很明显的周期性误差，峰峰值大约在正负50ns以内，一直没有突破性发展。本实验室多年来致力于rdss定位授时技术的研究发展，经过数次分析、研究、评估、论证，最终取得rdss授时技术10几年来的突破性进展，在已授权的申请号为cn201510537516.4《一种改进的高精度rdss授时方法》中提出两点改进意见：一是研究发现geo卫星运动产生的多普勒频移对rdss授时具有很大的影响；二是利用卫星速度推算卫星位置也会带来很大的误差，需要提出新的改进策略。本发明在上述基础之上，通过系统与接收机端的联调测试分析，进一步取得突破进展：一是改进了传统授时接收机计算上行时延的算法，并在接收机硬件算法可实现考虑的基础之上提出一种简单的上行时延推算方法；二是将类似的推导算法用于卫星位置的推导，或者说是卫星到用户下行时延的推导。
综上，国内近十余年对rdss授时技术的研究一直未有中断，特别时对rdss授时周期性缓变误差的分析，主要关注点在电离层时延误差、对流层时延误差、卫星星历误差、地球自转引起的sagnac效应、以及接收机零值几个方面，相关研究报道不少，但是一直没有找到rdss授时周期性缓变误差的根本原因。国外相关研究无，因为rdss是我国卫星导航系统的一个专有特色。本发明在研究rdss定位授时技术及其应用10多年的基础上，通过系统端与接收机端的联试联调，针对rdss授时改进在分别北京、三亚、喀什测试终端的测试分析研究，通过长期的数据采集、试验分析、质量监测、改进研究及评估验证，最终定位授时缓变误差原因，并研究发现通过改进接收机端上行时延的计算算法，可以使rdss授时服务精度大幅度提升。
当前传统rdss接收机在计算中心站到卫星上行时，利用距离当前卫星转发时刻t最近的两个整分点时刻t0，t1所播发的上行时延cs0,cs1计算当前卫星转发时刻的上行时延：假设授时卫星转发时刻t相对最近整分点t0时刻的时间差dt＝t-t0，中心站到卫星上行时延cs＝cs0+(vcs0+acs×dt)×dt。随着卫星运动的周期性变化，一直以来授时服务存在周期性缓变误差；本发明在此基础之上经过多次研究论证，提出一种新的改进算法，利用距离授时卫星转发时刻t最近的三个整分点时刻播报的上行时延来推算t时刻的上行时延：cs＝cs0+((cs0-cs1)/dt+(vcs0-vcs1)×dt/dt)×dt，经过三亚、喀什、北京三站授时服务试验验证，采用该方法可以使rdss授时服务精度从当前传统接收机的50ns提高到约10ns的授时服务精度(见附图6)。
有鉴于此，本发明提供了一种改进的高精度rdss授时方法，能够使rdss授时精度得到大幅度提升。
一种高精度rdss授时方法，如图2所示，包括以下步骤：
步骤1、获得整分点时刻的中心站到卫星的上行时延值、整分点时刻卫星位置、卫星速度以及时延修正参数；
步骤2、利用步骤1获得的所述整分点时刻播报的中心站到卫星的上行时延值计算当前修正时刻对应的上行时延，如图3所示，具体步骤如下：
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授时播发时间间隔，为1分钟。
s203：根据s202所推算的两分钟的中心站到卫星上行时延的变化速度vcs，推算上行时延的加速度：
acs＝(vcs0-vcs1)/dt。
s204：根据s202所推算的两分钟的中心站到卫星上行时延的变化速度vcs，以及s203所推算的上行时延的加速度，计算当时授时卫星转发时刻的上行时延值cs：
根据步骤1假设当前授时卫星转发时刻为t时刻，则t时刻相对t0时刻的时间差dt＝t-t0
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)×dt2/dt2
即，最终上行时延的计算仅用到距离当前时刻最新的三个整分点时刻播发的中心站到卫星的上行时延值。
步骤3、求解当前rdss授时信号卫星转发时刻对应的卫星位置，从而得到卫星到用户下行时延：
传统方法利用卫星位置和速度推算当前时刻卫星位置，本方法仅利用整分点时刻播报的卫星位置计算当前rdss授时信号卫星转发时刻对应的卫星位置：
如图4所示，方法一的步骤如下：
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-2xs1)×dt2/dt2
ys＝ys0+(ys0-ys1)×dt/dt+(ys0+ys2-2ys1)×dt2/dt2
zs＝zs0+(zs0-zs1)×dt/dt+(zs0+zs2-2zs1)×dt2/dt2
其中，dt是指rdss授时播发时间间隔，为1分钟；dt为t时刻相对t0时刻的时间差dt＝t-t0。
s303：利用s302所推算的卫星转发时刻卫星位置的st(xs,ys,zs)以及已知点用户位置u(xu,yu,zu)计算卫星到用户时延su：
su＝[(xs-xu)2+(ys-yu)2+(zs-zu)2]1/2。
如图5，方法二具体步骤如下：
s311：累积步骤1所述的距离当前时刻最新的三个整分点时刻播发的卫星位置；假设当前授时卫星转发时刻为t时刻，距离当前时刻最新的三个整分点时刻t0、t1、t2，三个整分点时刻离当前时刻由新及远，所播发的整分点时刻卫星位置分别为s0(xs0,ys0,zx0)、s1(xs1,ys1,zx1)、s2(xs2,ys2,zx2)。
s312：利用s311所述的整分点所播发的卫星位置以及已知点用户位置u(xu,yu,zu)计算t0、t1、t2整分点时刻对应的卫星到用户下行时延：
su0＝[(xs0-xu)2+(ys0-yu)2+(zs0-zu)2]1/2
su1＝[(xs1-xu)2+(ys1-yu)2+(zs1-zu)2]1/2
su2＝[(xs2-xu)2+(ys2-yu)2+(zs2-zu)2]1/2
s313：利用类似步骤2所述推算方法推算当前卫星转发时刻对应的卫星到用户的下行时延su：
su＝su0+(su0-su1)×dt/dt+(su0+su2-2su1)×dt2/dt2
其中，dt是指rdss授时播发时间间隔，为1分钟；dt为t时刻相对t0时刻的时间差dt＝t-t0。
步骤4、利用用户机零值以及步骤1中的所述时延修正参数，计算包括电离层、对流层、sagnac效应和用户机零值对时延的影响，即得到时延修正值；
步骤5、将rdss授时信号发射时刻同时加上步骤2获得的上行时延、步骤3获得的下行时延以及步骤4获得的时延修正值，完成经修正的授时计算。
此外，也可以利用超过3个点以上的短时间内整分点时刻中心站到卫星的上行数据，通过拟合、差分或拟合差分的方式计算当前时刻的上行时延，通过这种方法服务效果会更好，主要考虑接收端的可实现性。
其次，本发明所述方法同样可以用到北斗三号rdss授时，只是与北斗二号rdss授时不同的是，北斗三号rdss授时的播发频率为6s，即dt＝6s，本发明所用数据及时间对应的改为最近的三个6s周期播发数据，计算相应的当前时刻上行时延和下行时延。
综上所述，以上仅为本发明的较佳实施例而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
- 还没有人留言评论。精彩留言会获得点赞！