本发明涉及卫星定位导航中的时间传递技术，尤其涉及嵌入式远距离实时gps共视时间比对方面。
背景技术：
时间统一作为空间载体/分布式应用系统协同工作，高精度相对测量作为时间统一的基本条件，是现代国防高精度测量和全球一体化的重要研究内容，影响着航天、通信、交通、高速数字网同步等领域，如何实现大尺度范围内不同节点的实时时间同步，是现代时间基准统一和时间维持的重要保障。
gps共视时间比对作为时间基准统一的重要技术，建设的成本低，实现方法简单，同时随着gnss导航系统的完善，可用于共视观测卫星越来越多，同时gnss共视具备连续实时性特征，对于实现远距离实时时间同步具有无可比拟的优势。同时随着现在授时终端的模块化，利用多模多频板卡能够自主、方便地开发出各种满足特定需求的导航导航应用系统。将板卡与计算机、通讯技术相结合，开发能够实时gps共视时间比对算法。
技术实现要素：
本发明的目的是为了解决嵌入式远距离实时gps共视时间比对问题，目的在于解决相距几千公里的不同站间能够完成实时时间同步的问题。
本发明采用的技术方案为：
嵌入式远距离实时gps共视时间比对方法，是基于fpga原型的板卡，整个系统核心器件是dsp和fpga；gps共视时间比对部分是嵌入到dsp上，实现实时时间比对处理，主要完成任务包含gps共视与单机协同，在不增加运算复杂度和空间复杂度的前提下，保证钟差结果的稳定输出，包括以下步骤：
(1)实时数据解析：接收gps实时数据流，对不同测站的数据进行解析得到原始观测数据和电文数据，同时进行分类存储；
(2)数据预处理：对原始观测数据进行粗差剔除、周跳探测及修复和伪距平滑处理，得到高精度的平滑伪距值；
(3)单站钟差计算：根据高精度的平滑伪距值和电文数据，计算卫星位置和站星距离，同时对站星距进行误差项的改正，利用滤波的方式求解得到单站钟差；
(4)站间钟差求解：对不同站间的钟差数据进行交换，利用多星加权方法计算得到站间钟差结果。
其中，步骤(1)实时数据解析包括以下步骤：
(a)根据测站观测的gps实时数据流，首先解析测站编号标识；
(b)根据测站编号标识，分离和提取各个测站的原始观测数据和电文数据，并进行分类存储。
其中，其特征是步骤(2)数据预处理包括以下步骤：
(a)采用中值滤波和卡方检验的方法剔除原始观测数据中不好的数据；
(b)采用无几何距离和m-w组合来探测粗差剔除后的原始观测数据中包含的周跳，同时对周跳进行标识和修复；
(c)利用载波相位观测值对周跳探测及修复后的原始观测数据的伪距观测值进行平滑，得到高精度的平滑伪距值。
其中，步骤(3)单站钟差计算包括以下步骤：
(a)根据电文数据计算得到卫星的位置，根据卫星的位置、已知的测站坐标和平滑伪距值计算得到站星距；
(b)对站星距进行测站电离层、对流层、天线相位中心、地球自转以及dcb误差项的改正；
(c)利用经过误差改正后的平滑伪距值，通过滤波的方式计算单站钟差，并对钟差结果进行规范化处理。
本发明与现有技术相比优点为：
本发明提出的嵌入式远距离实时gps共视时间比对的方法，不仅解决了相距几千公里的不同站间时间同步问题，而且解决了远距离时间比对的模块化问题，具有很强的工程可实现性。
附图说明
图1是本发明的嵌入式远距离实时gps共视时间比对的方法中数据处理流程图。
具体实施方式
下面结合附图和实施例对本发明做进一步的描述。
如图1所示，嵌入式远距离实时gps共视时间比对方法，是基于fpga原型的板卡，整个系统核心器件是dsp和fpga；gps共视时间比对部分是嵌入到dsp上，实现实时时间比对处理，主要完成任务包含gps共视与单机协同，在不增加运算复杂度和空间复杂度的前提下，保证钟差结果的稳定输出，包括以下步骤：
(1)实时数据解析：接收gps实时数据流，对不同测站的数据进行解析得到原始观测数据和电文数据，同时进行分类存储；
(2)数据预处理：对原始观测数据进行粗差剔除、周跳探测及修复和伪距平滑处理，得到高精度的平滑伪距值；
(3)单站钟差计算：根据高精度的平滑伪距值和电文数据，计算卫星位置和站星距离，同时对站星距进行误差项的改正，利用滤波的方式求解得到单站钟差；
(4)站间钟差求解：对不同站间的钟差数据进行交换，利用多星加权方法计算得到站间钟差结果。
实时数据解析包括以下步骤：
(a)根据测站观测的gps实时数据流，首先解析测站编号标识；
(b)根据测站编号标识，分离和提取各个测站的原始观测数据和电文数据，并进行分类存储。
数据预处理包含粗差剔除、周跳探测及修复、相位平滑伪距。步骤如下：
(a)粗差剔除：主要是对实数采集的导航电文、原始观测数据和星历文件等进行有效性检验，剔除异常数据，检测当前共视时刻，整理数据格式为下一步的数据处理做准备。
(b)周跳探测及修复：周跳是发生在gps载波相位观测值中由于卫星信号被某障碍物阻挡、信号失锁等原因造成整周记数跳变的现象。本发明利用两种方式组合来探测周跳。一种是无几何距离组合观测值(geometry-free)参见公式(1),另一种是m-w组合观测值参见公式(2)。
式中，l3、p3表示观测值的无电离层组合的相位和伪距值，a、b表示组合系数，l1、l2表示gps双频载波相位观测值，p1、p2表示双频伪距观测值。f1、f2表示gpsl1、l2载波的频率。
mw的组合公式如下：
式中，l6表示观测值的melbourne-wubbena组合。
(c)相位平滑伪距：伪距的观测精度比较低，载波相位的观测精度比较高，利用载波相位的变化平滑伪距可以降低伪距中的观测噪声从而达到提高伪距的精度，本软件采用消除电离层影响的双频相位平滑伪距法，该方法首先要求出消除电离层以后的伪距和载波相位观测值，
单站钟差计算过程中关键的步骤是减弱或消除(电离层、对流层参见公式(3)、天线相位中心、地球自转、dcb)各项误差的影响。站间钟差计算通过交换共视站的数据，选取相同的观测卫星，选用多颗共视卫星采用加权平均的方法求取站间钟差。
单站钟差计算包括以下步骤：
(a)根据电文数据计算得到卫星的位置，根据卫星的位置、已知的测站坐标和平滑伪距值计算得到a站和b站的站星距；
(b)对a站和b站的站星距分别进行测站电离层、对流层、天线相位中心、地球自转以及dcb误差项的改正；
(c)利用经过误差改正后的平滑伪距值，通过滤波的方式分别计算a站和b站的钟差，并对钟差结果进行规范化处理。
对流层模型公式如下：
式中，z为卫星的天顶角，t为以开尔文为单位的地面站温度，p为大气压，单位为毫巴(mbar)，e为水气压，单位mbar，校正量b(h)和δr(h,z)可通过实际情况选择相应的数值。