一种bds三频伪距相位组合的周跳探测修复方法及装置的制造方法
【技术领域】
[0001]本发明涉及卫星定位技术领域，尤其涉及一种BDS三频伪距相位组合的周跳探测 修复方法及装置。
【背景技术】
[0002] 接收机在接收卫星信号的过程中，由于障碍物遮挡、接收机故障和低信噪比等因 素造成卫星信号暂时中断，致使载波相位观测值整周计数发生突变，而不足一周的部分依 然正确的现象称为周跳，周跳的存在严重影响了载波观测值，因此，在高精度导航和定位中 必须予以消除。周跳的确定是载波相位测量中特有的问题，要想获得高精度的载波观测值， 就需要首先确定周跳，可知BDS((BeiDou Navigation Satellite System，中国北斗卫星导 航系统)播发三个频率的载波信号，如何利用三频组合的优势探测周跳成为研究的热点。
[0003] 现有技术中的BDS三频伪距相位组合的周跳探测修复方法虽然探测精度高，通过 选取合适的组合系数可以探测大部分周跳，但是在低采样率的情况下，由于采样间隔大，历 元间站星距离变化大，噪声残差和电离层残差影响变大，使得周跳检验量变大，导致会有周 跳解算错误的情况出现，降低了周跳探测和修复的精度。
【发明内容】
[0004] 有鉴于此，本发明提供一种BDS三频伪距相位组合的周跳探测修复方法及装置，用 于避免周跳解算错误的情况出现，提高周跳探测和修复的精度。
[0005] 本发明通过以下技术手段解决上述问题：
[0006] 本发明提供一种BDS三频伪距相位组合的周跳探测修复方法，包括:获取多项式拟 合方程;根据所述多项式拟合方程和所述多项式拟合方程的残差确定所述多项式拟合方程 的中误差，并根据所述多项式拟合方程确定当前历元的组合载波观测值的预测值;确定所 述当前历元的组合载波观测值的预测值与所述当前历元的组合载波观测值的实测值的差 值的绝对值是否小于k倍的所述多项式拟合方程的中误差;其中，所述k为正整数;在确定所 述当前历元的载波预测值与所述当前历元的组合载波实测值的差值的绝对值大于k倍的所 述多项式拟合方程的中误差时，确定当前历元发生周跳;获取当前历元各频点的周跳值，并 修复当前历元各频点的周跳值。
[0007] 进一步，所述获取当前历元各频点的周跳值包括:根据方程
解算当前历元各频点的周跳值，得到ΛΝιΑ^^Λη;其中，A为根据三频伪距相位组合系数 xi，yi，zi得到的探测矩阵，i = 1，2，3，ΛΝ为各频点的周跳值集合，Λη为当前历元的三个三 频伪距相位组合的周跳值的集合。
[0008] 进一步，在所述修复当前历元各频点的周跳值后，还包括:根据所述探测矩阵Α和 修复后的当前历元各频点的周跳值确定修复后的组合周跳值ΛΝβ;确定是否大 于零;其中，为所述组合观测值的探测阈值;在确定I△乂大于零时，将所述当 前历元的组合载波观测值的预测值与所述当前历元的组合载波观测值的实测值的差值取 整，并根据所述探测矩阵A和取整后得到的周跳值再次解算各频点的周跳值，修复再次解算 后得到的各频点的周跳值，并根据所述探测矩阵A和修复后的所述再次解算后得到的各频 点的周跳值确定再次修复后的组合周跳值，直至确定所述再次修复后的组合周跳值的绝对 值与组合观测值的探测阈值的差值小于等于零。
[0009] 进一步，所述获取多项式拟合方程包括:获取m个无整周跳变的载波观测值;其中， 所述m为正整数;采用最小二乘法，并根据所述m个无整周跳变的载波观测值得到多项式拟 合方S
其中，所述m个无整周跳变的载波观 测值为所述当前历元之前获取的值，i = l，2,......m，m>k+l，k为正整数，ao，ai,…，ak为多 项式拟合方程的系数;根据所述多项式拟合方程和所述多项式拟合方程的残差确定所述多 项式拟合方程的中误差包括:根据公式
确定所述多项式拟合方程的中误 差，其中，σ为多项式拟合方程的中误差，Vi为多项式拟合方程的残差，m为无整周跳变的载 波观测值的个数，η为拟合多项式的阶数。
[0010] 进一步，所述k等于3。
[0011]本发明提供一种BDS三频伪距相位组合的周跳探测修复装置，包括:获取单元，用 于获取多项式拟合方程;确定单元，用于根据所述多项式拟合方程和所述多项式拟合方程 的残差确定所述多项式拟合方程的中误差，并根据所述多项式拟合方程确定当前历元的组 合载波观测值的预测值;所述确定单元，还用于确定所述当前历元的组合载波观测值的预 测值与所述当前历元的组合载波观测值的实测值的差值的绝对值是否小于k倍的所述多项 式拟合方程的中误差;其中，所述k为正整数;所述确定单元，还用于在确定所述当前历元的 载波预测值与所述当前历元的组合载波实测值的差值的绝对值大于k倍的所述多项式拟合 方程的中误差时，确定当前历元发生周跳;所述获取单元，还用于获取当前历元各频点的周 跳值;处理单元，用于修复当前历元各频点的周跳值。
[0012] 进一步，所述获取单元，具体用于根据方茜
挥算当前历 元各频点的周跳值，得到ΛΝ=Α，Λη;其中，A为根据三频伪距相位组合系数^，71，21得到 的探测矩阵，i = 1，2，3，ΛΝ为各频点的周跳值集合，Λη为当前历元的三个三频伪距相位组 合的周跳值的集合。
[0013] 进一步，所述确定单元，还用于根据所述探测矩阵Α和修复后的当前历元各频点的 周跳值确定修复后的组合周跳值ΛΝβ;所述确定单元，还用于确定|Δ# 5| - n&Vs是否大于零； 其中，为所述组合观测值的探测阈值；所述处理单元，还用于在所述确定单元确定 ΙΔΛ^Ι-πΑ%大于零时，将所述当前历元的组合载波观测值的预测值与所述当前历元的组 合载波观测值的实测值的差值取整，并根据所述探测矩阵Α和取整后得到的周跳值再次解 算各频点的周跳值，修复再次解算后得到的各频点的周跳值，并根据所述探测矩阵A和修复 后的所述再次解算后得到的各频点的周跳值确定再次修复后的组合周跳值，直至确定所述 再次修复后的组合周跳值的绝对值与组合观测值的探测阈值的差值小于等于零。
[0014] 进一步，所述获取单元，具体用于获取m个无整周跳变的载波观测值;其中，所述m 为正整数;所述获取单元，具体用于采用最小二乘法，并根据所述m个无整周跳变的载波观 测值得到多项式拟合方S
;其中，所述m个无 整周跳变的载波观测值为所述当前历元之前获取的值，i = l，2,......m，m>k+l，k为正整 数，，…，ak为多项式拟合方程的系数；所述确定单元，具体用于根据公式
商定所述多项式拟合方程的中误差，其中，σ为多项式拟合方程的中误差， I为多项式拟合方程的残差，m为无整周跳变的载波观测值的个数，η为拟合多项式的阶数。
[0015] 进一步，所述k等于3。
[0016]本发明提供的一种BDS三频伪距相位组合的周跳探测修复方法及装置，根据多项 式拟合方程获取到当前历元的组合载波观测值的预测值，在确定当前历元的载波预测值与 所述当前历元的组合载波实测值的差值的绝对值大于k倍的多项式拟合方程的中误差时， 确定当前历元发生周跳，此时，获取当前历元各频点的周跳值，并修复当前历元各频点的周 跳值，这样，在低采样率的情况下，能够避免周跳解算错误的情况出现，提高周跳探测和修 复的精度。
【附图说明】
[0017]图1是本发明提供的一种BDS三频伪距相位组合的周跳探须彳修复方法的流程示意 图；
[0018] 图2是本发明提供的另一种BDS三频伪距相位组合的周跳探测修复方法的流程示 意图；
[0019] 图3是本发明提供的一种BDS三频伪距相位组合的周跳探测修复装置的功能示意 图。
【具体实施方式】
[0020] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完 整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于 本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他 实施例，都属于本发明保护的范围。
[0021] 参见图1，是本发明提供的一种BDS三频伪距相位组合的周跳探测修复方法的流程 示意