利用抗差多项式拟合法探测修复GNSS周跳

《利用抗差多项式拟合法探测修复GNSS周跳》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、第３３卷第３期大地测量与地球动力学Ｖｏｌ．３３Ｎｏ．３２０１３年６月ＪＯＵＲＮＡＬＯＦＧＥＯＤＥＳＹＡＮＤＧＥＯＤＹＮＡＭＩＣＳＪｕｎ．，２０１３文章编号：１６７１５９４２（２０１３）０３０１２９０５?利用抗差多项式拟合法探测修复ＧＮＳＳ周跳１，２）２）１）２）王福丽成英燕韦铖王晓明(１）山东科技大学测绘学院，青岛２６６５９０)２）中国测绘科学研究院，北京１００８３０摘要多项式拟合法探测周跳的能力受载波相位观测值的观测误差和采样率的较大限制。将观测数据分成若干段处理，虽然对结果有所改善，但

2、降低了工作效率。采用历元间求差与抗差估计相结合的方法对多项式进行改进。并通过ＧＰＳ／ＧＬＯＮＡＳＳ实测数据进行实验，结果表明，利用抗差多项式拟合探测修复周跳的结果更加稳定可靠。关键词多项式拟合；历元间求差；抗差估计；ＧＰＳ／ＧＬＯＮＡＳＳ；探测修复周跳中图分类号：Ｐ２０７文献标识码：ＡＧＮＳＳＣＹＣＬＥＳＬＩＰＤＥＴＥＣＴＩＯＮＡＮＤＣＯＲＲＥＣＴＩＯＮＵＳＩＮＧＲＯＢＵＳＴＰＯＬＹＮＯＭＩＡＬＦＩＴＴＩＮＧ１，２）２）１）２）ＷａｎｇＦｕｌｉ，ＣｈｅｎｇＹｉｎｇｙａｎ，ＷｅｉＣｈｅｎｇａｎｄＷａ

3、ｎｇＸｉａｏｍｉｎｇ(１）ＣｏｌｌｅｇｅｏｆＧｅｏｍａｔｉｃｓ，ＳｈａｎｄｏｎｇＵｎｉｖｅｒｓｉｔｙｏｆＳｃｉｅｎｃｅａｎｄＴｅｃｈｎｏｌｏｇｙ，Ｑｉｎｇｄａｏ２６６５９０)２）ＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｕｒｖｅｙｉｎｇａｎｄＭａｐｐｉｎｇ，Ｂｅｉｊｉｎｇ１００８３０ＡｂｓｔｒａｃｔＴｈｅａｂｉｌｉｔｙｏｆｐｏｌｙｎｏｍｉａｌｆｉｔｔｉｎｇｆｏｒｄｅｔｅｃｔｉｎｇａｎｄｃｏｒｒｅｃｔｉｎｇｃｙｃｌｅｓｌｉｐｓｉｓｓｔｒｏｎｇｌｙａｆｆｅｃｔｅｄｂｙｔｈｅｅｒｒｏｒｓｏｆｔｈｅｃａｒｒｉｅ

4、ｒｐｈａｓｅｏｂｓｅｒｖａｔｉｏｎｓａｎｄｔｈｅｉｎｔｅｒｖａｌｓｏｆｔｈｅｏｂｓｅｒｖａｔｉｏｎｄａｔａ．Ｔｈｅｏｂｓｅｒｖａｔｉｏｎｄａｔａｂｅｉｎｇｄｉｖｉｄｅｄｉｎｔｏｓｅｖｅｒａｌｐｉｅｃｅｓｗｉｌｌｉｍｐｒｏｖｅｔｈｅｒｅｓｕｌｔｓ，ｈｏｗｅｖｅｒ，ｔｈｅｅｆｆｉｃｉｅｎｃｙｗｉｌｌｂｅｄｅｃｒｅａｓｅｄ．Ａｍｅｔｈｏｄｆｏｒｃｏｍｂｉｎｉｎｇｅｐｏｃｈｄｉｆｆｅｒｅｎｃｅｄｗｉｔｈｒｏｂｕｓｔｅｓｔｉｍａｔｉｏｎｔｏｉｍｐｒｏｖｅｔｈｅｍｅｔｈｏｄｏｆｐｏｌｙｎｏｍｉａｌｆｉｔｔ

5、ｉｎｇｉｓｐｒｏｐｏｓｅｄ．ＴｈｅｎｕｍｅｒｉｃａｌｅｘｐｅｒｉｍｅｎｔｉｓｐｅｒｆｏｒｍｅｄｂｙｕｓｉｎｇｔｈｅｒｅａｌｄａｔａｏｆＧＰＳ／ＧＬＯＮＡＳＳ．Ｔｈｅｒｅｓｕｌｔｓｉｎｄｉｃａｔｅｔｈａｔｔｈｅｍｅｔｈｏｄｏｆｒｏｂｕｓｔｐｏｌｙｎｏｍｉａｌｆｉｔｔｉｎｇｉｓｍｏｒｅｒｉｇｏｒｏｕｓｉｎｔｈｅｏｒｙａｎｄｍｏｒｅｐｏｗｅｒｆｕｌｉｎｐｒａｃｔｉｃｅｆｏｒｄｅｔｅｃｔｉｎｇａｎｄｃｏｒｒｅｃｔｉｎｇｃｙｃｌｅｓｌｉｐｓｃｏｍｐａｒｉｎｇｗｉｔｈｔｈｅｎｏｒｍａｌｐｏｌｙｎｏｍｉａ

6、ｌｆｉｔｔｉｎｇｍｅｔｈｏｄ．Ｋｅｙｗｏｒｄｓ：ｐｏｌｙｎｏｍｉａｌｆｉｔｔｉｎｇ；ｅｐｏｃｈｄｉｆｆｅｒｅｎｃｅｄ；ｒｏｂｕｓｔｅｓｔｉｍａｔｉｏｎ；ＧＰＳ／ＧＬＯＮＡＳＳ；ｄｅｔｅｃｔｉｏｎａｎｄｃｏｒｒｅｃｔｉｏｎｆｏｒｃｙｃｌｅｓｌｉｐｓ．［１－５］法、电离层残差法、相位减伪距法等。这些方法１引言对载波相位观测值周跳探测和修复都有其局限性，在多卫星系统组合定位中，ＧＮＳＳ载波相位整如高次差法受数据采样间隔影响显著，且只能探测周模糊度和周跳探测与修复是必须要解决的问题。大周跳，对小周跳无能

7、为力，当连续发生周跳时该方目前，探测周跳的方法主要有多项式拟合法、高次差法失效；电离层残差法受电离层延迟变化的影收稿日期：２０１２１１１９基金项目：地心坐标系维护与推广应用（Ａ１２０１）；卫星连续运行站数据处理（Ａ１９２８）作者简介：王福丽，女，１９８６年生，硕士，主要从事多卫星系统数据质量控制方面研究．Ｅ－ｍａｉｌ：ｗａｎｇｆｕｌｉ１９８８＿＠１２６．ｃｏｍ１３０大地测量与地球动力学３３卷［６－８］响，并存在多值性；相位减伪距法只适合双频接机钟偏差。收机，且可能出现将粗差误作周跳的情形。鉴

8、于此，利用式（６）对式（４）进行时间单位的标准化处本文将以多项式拟合法为基础，采用历元间求差与理，抗差估计相结合的方法对传统算法中存在问题进行ｔｉ－ｔ１Ｔｉ＝２×－１（６）改进，并利用ＧＰＳ／ＧＬＯＮＡＳＳ实测数据验证改进后ｔｉ－ｔ１的算法的有效性。这样就把ｔｉ∈［ｔ１，ｔｎ］一一映射为Ｔｉ∈［－１，１］，则式（４）变为等价的时间多项式：２传统的多项式拟合法２２Δφｉ＝ａ１（Ｔｉ－Ｔｉ－１）＋ａ２（Ｔｉ－Ｔｉ－１）＋…＋ｎｎ传统的多项式拟合法的数学模型为［９］：ａｎ（Ｔｉ－Ｔｉ－１