1.一种北斗实时精密卫星钟差评估方法,其特征在于,包括以下步骤:
1)对实时精密卫星钟差和事后精密卫星钟差作差,得到产品间单差;对GEO、IGSO和MEO卫星分类定权,得到卫星分类权因子ls;
2)利用所述产品间单差和卫星分类权因子对时间尺度偏差进行迭代估计,得到时间尺度偏差的估值;
3)逐个历元计算扣除所述时间尺度偏差估值后的产品间单差,得到扣除时间尺度偏差后的产品间单差的时间序列,再计算该序列的标准差,将该标准差作为实时精密钟差的精度评定指标。
2.根据权利要求1所述的北斗实时精密卫星钟差评估方法,其特征在于,步骤1)中,实时精密卫星钟差和事后精密卫星钟差分别由以下关系式表示:其中,下标a和b代表提供实时或事后精密卫星钟差的分析中心;上标s代表卫星;t表示历元; 代表由分析中心a提供的实时精密卫星钟差, 代表由分析中心b提供的事后精密卫星钟差;Oa(t)、Ob(t)为时间尺度; 为初始钟误差;Cs(t)为相位观测值估计的钟差改正数;
则在各历元t中,对实时精密卫星钟差和事后精密卫星钟差作差,得到产品间单差,其表达式为:其中,Δ(·)表示差分算子; 为产品间单差;ΔOab(t)表示时间尺度偏差;
为初始钟误差偏差。
3.根据权利要求1所述的北斗实时精密卫星钟差评估方法,其特征在于,步骤1)中,卫星分类权因子的计算公式为:其中,ls为GEO、IGSO和MEO卫星的权因子;lGEO、lIGSO、lMEO均为常数。
4.根据权利要求3所述的北斗实时精密卫星钟差评估方法,其特征在于,lMEO=1、0<lGEO<1,0<lIGSO<1。
5.根据权利要求2所述的北斗实时精密卫星钟差评估方法,其特征在于,步骤2)中,对于历元t,设产品间单差 在第j轮迭代中的权为 则第j轮迭代中时间尺度偏差的估值 其中,[·]j为第j轮迭代标记符;m为该历元中卫星的总个数; 为第j轮迭代的等价权因子。
6.根据权利要求5所述的北斗实时精密卫星钟差评估方法,其特征在于,在第j轮迭代中,等价权因子 的计算公式为:其中, 为标准化的残差, 为残差;|·|为绝
对值符号; 为单位权中误差,m为历元t中卫星的总个数;
为 对应的协因数;k0、k1为经验值。
7.根据权利要求6所述的北斗实时精密卫星钟差评估方法,其特征在于,若 服从标准正态分布,k0取1~1.5,k1取2.5~6,k0=1.5,k1=3。
8.根据权利要求5所述的北斗实时精密卫星钟差评估方法,其特征在于,当时停止迭代;其中,|·|为绝对值符号;δ为阈值; 为第j+1轮迭代中时间尺度偏差的估值。
9.根据权利要求5~8之一所述的北斗实时精密卫星钟差评估方法,其特征在于,步骤
3) 中 ,精 度 评 定 指 标 其中 n 为 历 元 总 个 数 ;
为扣除时间尺度偏差后的产品间单差,其计算公式为