一种用于北斗卫星导航系统共视数据的预处理方法
【专利摘要】本发明公开了一种用于北斗卫星导航系统共视数据的预处理方法，该方法步骤包括根据时频实验室的观测状态和卫星状态，选取合适的共视数据；对选取后的共视数据信号传播过程中的延迟进行修正；对修正后的共视数据进行粗大误差剔除处理；采用k阶Vondrak滤波器对共视数据进行平滑处理，通过最小二乘估计法迭代选取最佳平滑因子，消除随机误差。本发明所述技术方案计算量小，运算速度快，能够快速获得合适的平滑因子，能够有效的消除共视数据中的随机误差。
【专利说明】一种用于北斗卫星导航系统共视数据的预处理方法
【技术领域】
[0001]本发明涉及数据处理方法，特别是涉及一种用于北斗卫星导航系统共视时间频率传递的数据预处理方法。
【背景技术】
[0002]在北斗卫星导航系统共视时间频率传递过程中，参与共视的两个时频实验室在同一时刻观测同一颗北斗卫星，分别得到各自本地钟与该北斗卫星星载钟的时差数据，经通信链路将各自测得的数据相互交换，经数据处理去掉星载钟即可得到两时频实验室之间时差，从而实现北斗共视时间传递。
[0003]由于在接收北斗信号时存在干扰，为获得更精确的本地钟与星载钟时差，在进行数据交换之前，需要在本地进行数据预处理，包括:共视数据的选取，信号传播时延的修正，随机误差的削除。信号传播时延修正是指对时差数据中几何时延、空间传播介质时延、相对论时延、接收机硬件延迟等进行修正，可通过建立模型和校准的方式实现，而随机误差的削弱只能通过数据平滑的方式实现。
[0004]数据平滑是数据处理中一项很重要的工作，用户要从测量数据中得到有效的信息，必须对测量数据做预处理，尽可能减小随机误差的影响，同时还可以评估测量数据的质量。目前常用的平滑方法有曲线拟合，滑动平均方法和Vondrak平滑方法。曲线拟合方法的前提是给出函数形式，大多数为多项式，其阶数往往根据经验而定，人为因素影响较大；滑动平均方法要求数据必须是等间距而且采样点比较密集，具有很大的局限性。Vondrak平滑方法实质是通过选择平滑因子，在对测量序列的绝对拟合和绝对平滑之间寻求一个平衡，在未知拟合函数形式的情况下Vondrak方法也能对测量资料进行合理的平滑，其关键是选择合理的平滑因子。目前常用的选择平滑因子方法有:观测误差法、频率响应法、交叉认证法。观测误差法要求预先知道观测数据先验标准差，而通常我们并不能准确的知道标准差信息；频率响应法对于已经知道频率或周期的观测数据才能够得到较好的滤波效果；交叉认证方法是通过反复的试验选择平滑因子，运算量大，速度缓慢，选择的平滑因子随机性也比较大。
[0005]因此，需要提供一种共视数据的预处理方法，以快速选取平滑因子，降低运算量，提闻运行速度。
【发明内容】
[0006]本发明要解决的技术问题是提供一种用于北斗共视时间频率传递的数据预处理方法，解决目前共视数据处理中，常用平滑方法平滑效果不理想以及传统Vondrak平滑方法中平滑因子难以选择的问题。
[0007]为解决上述技术问题，本发明采用下述技术方案。
[0008]一种用于北斗卫星导航系统共视数据的预处理方法，该方法步骤包括
[0009]S1、根据时频实验室的观测状态和卫星状态，选取合适的共视数据；
[0010]S2、对选取后的共视数据信号传播过程中的延迟进行修正；
[0011]S3、对修正后的共视数据进行粗大误差剔除处理；
[0012]S4、采用k阶Vondrak滤波器对共视数据进行平滑处理，通过最小二乘估计法迭代选取最佳平滑因子，消除随机误差。
[0013]优选的，所述步骤SI中对共视数据的选取条件包括
[0014]保证两个时频实验室在同一时刻观测同一颗卫星；
[0015]被观测卫星的高度角不小于15° ;
[0016]根据北斗卫星导航系统特点延长全长跟踪时间。
[0017]优选的，所述步骤S2中需要修正的延迟项包括几何距离延迟、电离层延迟、对流层延迟，sagnac效应以及接收机延迟。
[0018]优选的，延迟修正方法包括
[0019]利用双频数据对电离层延迟进行修正；
[0020]根据北斗卫星导航系统接收机当行电文中提取的延迟参数和卫星星历建模，消除几何距离延迟、对流层延迟和sagnac效应；
[0021]通过校准接收机的方式克服接收机延迟。
[0022]优选的，所述步骤S3中采用3 σ准则剔除共视数据中的粗大误差。
[0023]优选的，所述步骤S4中采用最小二乘迭代法选取平滑因子，其步骤包括
[0024]S41、取平滑后的时差序列作为待估量，由Vondrak滤波基本原理，将时差序列和时差序列的k阶差分作为两类观测值，分别建立误差方程，从而将Vondrak滤波转化为最小二乘平差问题；
[0025]S42、给定初始平滑因子，利用最小二乘法求解观测值的误差方程，获得一次平滑结果;
[0026]S43、利用Helmert方差分量估计的方法，根据前次平滑得到的改正数向量V来估计方差分量，修正初始平滑因子；
[0027]S44、通过多次迭代，获得最佳平滑因子，同时迭代结束即完成消除随机误差。
[0028]本发明的有益效果如下:
[0029]本发明所述技术方案计算量小，运算速度快，能够快速获得合适的平滑因子，能够有效的消除共视数据中的随机误差。
【专利附图】
【附图说明】
[0030]下面结合附图对本发明的【具体实施方式】作进一步详细的说明；
[0031]图1示出本发明所述的一种用于北斗共视时间频率传递的数据预处理方法的示意图。
【具体实施方式】
[0032]下面结合一组实施例及附图对本发明做进一步描述。
[0033]本发明公开了一种用于北斗共视时间传递的数据预处理方法，该方法具体步骤包括:
[0034]第一步、共视数据的选取
[0035]在北斗共视时间传递技术中，选取共视数据时，除了接收的是同一颗卫星信号以外，还必须同时满足以下三个条件:
[0036]I)跟踪时刻相同，两个时频实验室须在同一时刻观测同一颗卫星，这样才可以完全消除卫星钟的误差；
[0037]2)高度角不小于15°，当卫星高度角过低时，电离层延迟误差难以改正，且多路径效应影响加剧；
[0038]3)跟踪长度，对于GPS来说，一次全长跟踪为13分钟，加上2分钟预置准备和I分钟数据处理，对于北斗共视来说，由于北斗卫星导航系统含GEO和IGSO卫星，可适当延长一次全长跟踪时间。
[0039]第二步、信号传播时延修正
[0040]选取共视数据后，应对卫星信号在传播过程中的各项延迟，包括几何距离延迟、电离层延迟、对流层延迟，sagnac效应以及接收机延迟等进行修正，并在共视数据中予以扣除。其中，电离层延迟可利用双频数据修正，几何距离延迟、对流层延迟、sagnac效应等可根据从北斗接收机导航电文中提取的延迟参数、卫星星历等建模消除，接收机延迟则通过校准的方式获得。
[0041]第三步粗大误差剔除
[0042]由于Vondrak平滑易受粗大误差的影响，在进行数据平滑之前首先采用3 σ准则剔除共视数据中的粗大误差。由于共视数据为每秒一次，前后两次测量相差不大，故可以直接利用粗差前后数据均值替代剔除掉的粗差。
[0043]第四步Vondrak平滑
[0044]采用Vondrak平滑的方法对经信号传播延迟修正及粗差剔除后的共视数据进行数据平滑，消除随机误差。
[0045]Vondrak方法平滑测量数据的基本原理是:
[0046]Q = F+λ 2S = min (I)
[0047]其中，F= Σρ^/ Ii)2,S = Σ (Akyi)^yiQ = I, 2，…，N)为时间引数为 Xi 的测量资料。Ii1为待求的平滑值，Pi为测量数据的权重。
[0048]本方案采用k阶Vondrak滤波器对共视平滑数据建模如下:
[0049]在和两点之间的间隔上用一个通过i，i+l，…，i+k共
k+1个点的k阶拉格朗日多项式来定义平滑结果y' =y' (t)，相应的，平滑值的k阶差分可表示为
【权利要求】
1.一种用于北斗卫星导航系统共视数据的预处理方法，其特征在于，该方法步骤包括 s1、根据时频实验室的观测状态和卫星状态，选取合适的共视数据； s2、对选取后的共视数据信号传播过程中的延迟进行修正； s3、对修正后的共视数据进行粗大误差剔除处理； s4、采用k阶Vondrak滤波器对共视数据进行平滑处理，通过最小二乘估计法迭代选取最佳平滑因子，消除随机误差。
2.根据权利要求1所述的数据预处理方法，其特征在于，所述步骤SI中对共视数据的选取条件包括 保证两个时频实验室在同一时刻观测同一颗卫星； 被观测卫星的高度角不小于15° ; 根据北斗卫星导航系统特点延长全长跟踪时间。
3.根据权利要求1所述的数据预处理方法，其特征在于，所述步骤S2中需要修正的延迟项包括几何距离延迟、电离层延迟、对流层延迟，sagnac效应以及接收机延迟。
4.根据权利要求3所述的数据预处理方法，其特征在于，延迟修正方法包括 利用双频数据对电离层延迟进行修正； 根据北斗卫星导航系统接收机当行电文中提取的延迟参数和卫星星历建模，消除几何距离延迟、对流层延迟和sagnac效应； 通过校准接收机的方式克服接收机延迟。
5.根据权利要求1所述的数据预处理方法，其特征在于，所述步骤S3中采用3σ准则剔除共视数据中的粗大误差。
6.根据权利要求1所述的数据预处理方法，其特征在于，所述步骤S4中采用最小二乘迭代法选取平滑因子，其步骤包括 s41、取平滑后的时差序列作为待估量，由Vondrak滤波基本原理,将时差序列和时差序列的k阶差分作为两类观测值，分别建立误差方程，从而将Vondrak滤波转化为最小二乘平差问题； s42、给定初始平滑因子，利用最小二乘法求解观测值的误差方程，获得一次平滑结果; s43、利用Helmert方差分量估计的方法，根据前次平滑得到的改正数向量V来估计方差分量，修正初始平滑因子； s44、通过多次迭代，获得最佳平滑因子，同时迭代结束即完成消除随机误差。
【文档编号】G01S19/23GK104199054SQ201410419574
【公开日】2014年12月10日 申请日期:2014年8月22日 优先权日:2014年8月22日
【发明者】易航, 王宏博, 张升康, 王海峰 申请人:北京无线电计量测试研究所