[发明专利]一种基于区域多参考站联合解算的卫星钟差实时估计方法

权利要求书

1.一种基于区域多参考站联合解算的卫星钟差实时估计方法，其特征在于，包括以下步骤：

步骤(1)，对参考站获取的观测文件数据进行预处理；

步骤(2)，对实时获取的IGU预报轨道进行插值计算，以匹配所述观测文件的采样间隔；

步骤(3)，固定卫星轨道参数、参考站坐标和对流层参数，并对误差模型改正；

步骤(4)，采用测量的精度为0.2～0.3mm的载波相位观测值，通过非差无电离层组合，得到多个参考站的卫星无电离层组合相位观测方程，并分别得到所述多个参考站对应的卫星无电离层组合相位观测值L；所述无电离层组合相位观测方程为：

L＝ρ+c·δ_r-c·δ+N+Trop+d_M+ε_Φ   (1)

式中，ρ为卫星与接收机之间的几何距离；δ_r、δ分别为接收机和卫星钟差；c为光速；N为无电离层组合的模糊度参数；Trop为对流层延迟改正；d_M为相对论效应、潮汐改正和相位缠绕等误差项的模型改正；ε_Φ为相位的观测噪声；

步骤(5)，在同一参考站同一历元中，对所述卫星无电离层组合相位观测值进行星间差分，得到星间差分无电离层组合的相位观方程，并得到星间差分无电离层组合的相位观测值ΔL；所述星间差分无电离层组合的相位观方程为：

ΔL＝Δρ-c·Δδ+ΔN+ΔTrop+Δd_M+Δε_Φ   (2)

式中，Δρ为站、星几何距离星间差；Δδ为钟差星间差，即相对钟差；c为光速；ΔN为星间差分模糊度参数；ΔTrop为星间差分对流层延迟改正Δd_M为星间差分误差项模型改正；Δε_Φ为星间差分观测噪声；

遍历所述多个参考站，得到所述每个参考站对应的星间差分无电离层组合的相位观值；

步骤(6)，对所述多个参考站对应的星间差分无电离层组合的相位观测值进行历元间差分，得到每个参考站对应的星间差分观测方程，并得到每个参考站对应的星间差分观测值的历元间差▽ΔL；所述星间差分观测方程为：

▽ΔL＝▽Δρ-c·▽Δδ+▽ΔTrop+▽Δd_M+▽Δε_Φ   (3)

式中，▽Δρ为站、星几何距离星间差分的历元间差；c为光速；▽Δδ为相对钟差的历元间差；▽ΔTrop为对流层延迟改正星间差分的历元间差；▽Δd_M为误差项模型改正星间差分的历元间差；▽Δε_Φ为观测噪声星间差分的历元间差；

步骤(7)，所述多参考站对应的星间差分观测方程联立，采用最小二乘估计进行解算，得到相对于参考历元的相对卫星钟差：假设有m个参考站同步观测n颗GPS卫星，定义参考卫星号为r＝1，有如下观测方程：

AX＝B   (4)

其中，未知参数矩阵X：

X＝[Tw₁，...，Tw_m，▽Δδ^1，2，▽Δδ^1，3，...，▽Δδ^1，n]^T   (5)

未知参数中有m个观测站天顶对流层湿延迟Tw_m和n-1个相对卫星钟差历元间差▽Δδ^1，n；假定天顶对流层湿延迟不变；系数矩阵A为：

A=e1(n-1)×2mE·c(n-1)×(n-1)e2(n-1)×2mE·c(n-1)×(n-1)......em(n-1)×2mE·c(n-1)×(n-1)---(6)]]>

所述观测方程的观测向量B为：

B=u1(n-1)×1u2(n-1)×1...um(n-1)×1---(8)]]>

uk=▿Δρk1,2-▿ΔLk1,2+▿Δdryk1,2+▿ΔdMk1,2▿Δρk1,3-▿ΔLk1,3+▿Δdryk1,3+▿ΔdMk1,3...▿Δρk1,n-▿ΔLk1,n+▿Δdryk1,n+▿ΔdMk1,n---(9)]]>

其中，E为单位阵；c为光速；▽ΔM^r，j为投影函数星间差的历元间差；j为卫星号；k为观测站号；r为参考卫星号；为传播路径上的对流层干延迟星间差的历元间差；

步骤(8)，根据所述步骤(7)得到所述相对于参考历元的相对卫星钟差，利用IGU预报钟差及其内插值进行约束，得到卫星在任一历元的相对钟差。