CN111025356A - 一种基于方差分量估计的gnss多系统自适应融合定位方法 - Google Patents

一种基于方差分量估计的gnss多系统自适应融合定位方法 Download PDF

Info

Publication number
CN111025356A
CN111025356A CN201911401417.8A CN201911401417A CN111025356A CN 111025356 A CN111025356 A CN 111025356A CN 201911401417 A CN201911401417 A CN 201911401417A CN 111025356 A CN111025356 A CN 111025356A
Authority
CN
China
Prior art keywords
gnss
variance
fusion positioning
component estimation
method based
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201911401417.8A
Other languages
English (en)
Inventor
许艳
许佳
蒋春华
江楠
刘广盛
王飞
孙张振
许国昌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Weihai Oreia Information Technology Co Ltd
Original Assignee
Weihai Oreia Information Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Weihai Oreia Information Technology Co Ltd filed Critical Weihai Oreia Information Technology Co Ltd
Priority to CN201911401417.8A priority Critical patent/CN111025356A/zh
Publication of CN111025356A publication Critical patent/CN111025356A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S19/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/42Determining position
    • G01S19/421Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system
    • G01S19/423Determining position by combining or switching between position solutions or signals derived from different satellite radio beacon positioning systems; by combining or switching between position solutions or signals derived from different modes of operation in a single system by combining or switching between position solutions derived from different satellite radio beacon positioning systems

Abstract

本发明提供一种基于方差分量估计的GNSS多系统自适应融合定位方法,其解决了现有GNSS多系统融合定位中不同系统观测值合理的权比问题,其包括以下步骤:对各个GNSS系统的观测值定初权;进行首次平差,得到各个GNSS系统观测值的残差信息;根据方差分量估计的方法计算各个GNSS系统方差因子并对其重新定权;确定多系统融合定位中各个GNSS系统观测值的权比,本发明可广泛应用于GNSS卫星精密定位领域。

Description

一种基于方差分量估计的GNSS多系统自适应融合定位方法
技术领域
本发明涉及GNSS卫星精密定位领域,尤其是涉及一种基于方差分量估计的GNSS多系统自适应融合定位方法。
背景技术
随着中国北斗二代卫星导航定位系统(BDS)的逐步实现,美国GPS系统的逐步现代化,俄罗斯GLONASS系统的逐步恢复完善,以及欧洲的Galileo系统逐步发展,越来越多的导航卫星能够为用户提供定位导航授时服务;单一导航系统服务存在可视卫星数目较少、星座几何分布较差以及卫星系统完好性薄弱等问题,因此更加稳健的多系统融合定位受到了广泛关注;多系统融合定位不仅比单系统具有更完善的星座分布,还具有更多的可视卫星数和更高的卫星可靠性等,然而由于不同系统之间存在星座差异,所以在考虑不同系统融合问题时,系统之间合理的权比分配就成为融合研究的核心问题。
在实际观测中,各类系统观测值的权会随时间、卫星升降变化,常数比并不能真实反映各类观测值实际精度,进而影响多系统融合定位的精度,因此,自适应地估计各个系统观测值的观测精度,合理调整各类观测值的权比,对GNSS多系统融合定位精度的提高有重要意义。
发明内容
本发明就是针对现有GNSS多系统融合定位中不同系统观测值合理的权比问题,提供一种基于方差分量估计的GNSS多系统自适应融合定位方法,可自适应调整不同GNSS系统观测值之间的融合比例因子,改善和提高融合定位精度。
为此,本发明中GNSS多系统自适应融合定位方法,包括以下步骤:
步骤1:对各个GNSS系统的观测值定初权;
步骤2:进行首次平差,得到各个GNSS系统观测值的残差信息;
步骤3:根据方差分量估计的方法计算各个GNSS系统方差因子并对其重新定权;
步骤4:确定多系统融合定位中各个GNSS系统观测值的权比。
优选的,对各个GNSS系统的观测值定初权,具体为:设GPS,GLONASS与BDS导航系统的观测值分别是L1,L2和L3,其权阵分别为P1,P2和P3,且观测值之间相互独立,相应的误差方程为:
Figure BDA0002347572150000021
Figure BDA0002347572150000022
Figure BDA0002347572150000023
则三类观测值的方差和单位权方差为:
Figure BDA0002347572150000024
Figure BDA0002347572150000025
可得法方程为:
Figure BDA0002347572150000026
根据二次型期望公式可得方差分量估计公式如下:
Figure BDA0002347572150000027
在GNSS融合定位中,将观测值分为GPS、GLONASS及BDS三类,并定各自的初始权分别为P1,P2和P3
优选的,进行首次平差,得到各个GNSS系统观测值的残差信息,具体为:
Vi TPiVi(i=1,2,3)
优选的,计算各个GNSS系统方差因子并对其重新定权,具体为:进行方差分量估计,求得各类观测值的单位权方差估值:
Figure BDA0002347572150000028
重新定权为:
Figure BDA0002347572150000031
优选的,确定多系统融合定位中各个GNSS系统观测值的权比,具体为:重复进行步骤2、3、4,直到
Figure BDA0002347572150000032
为止,满足条件后,计算结束输出参数估计结果。
本发明提出的基于方差分量估计的GNSS多系统自适应融合定位方法,将GNSS观测值根据各个系统进行分类,计算得到各个GNSS系统观测值的残差信息,根据方差分量估计的方法计算各个GNSS系统的方差因子并对其重新定权,从而自适应地合理确定多系统融合定位中各个GNSS系统观测值的权比,提高多系统融合定位的精度。
附图说明
图1为本发明实施例的计算流程示意图。
具体实施方式
下面参照附图详细说明本发明的具体实施方式。
如图1所示,本发明提供一种基于方差分量估计的GNSS多系统自适应融合定位方法,其包括以下步骤:
步骤1:对各个GNSS系统的观测值定初权,具体为:
设GPS,GLONASS与BDS导航系统的观测值分别是L1,L2和L3,其权阵分别为P1,P2和P3,且观测值之间相互独立。相应的误差方程为:
Figure BDA0002347572150000033
则三类观测值的方差和单位权方差为:
Figure BDA0002347572150000034
由式(1)可得法方程为:
Figure BDA0002347572150000035
式中,
Figure BDA0002347572150000041
Figure BDA0002347572150000042
根据二次型期望公式可得方差分量估计公式如下:
Figure BDA0002347572150000043
式中,
Figure BDA0002347572150000044
Figure BDA0002347572150000045
在GNSS融合定位中,将观测值分为GPS、GLONASS及BDS三类,并定各自的初始权分别为P1,P2和P3
步骤2:进行首次平差,得到各个GNSS系统观测值的残差信息,具体为:
Figure BDA0002347572150000046
步骤3:根据方差分量估计的方法计算各个GNSS系统方差因子并对其重新定权,具体为:根据式(4)进行方差分量估计,求得各类观测值的单位权方差估值:
Figure BDA0002347572150000047
再根据下式重新定权为:
Figure BDA0002347572150000048
其中c为任意常数,一般取
Figure BDA0002347572150000049
的某个值。
步骤4:确定多系统融合定位中各个GNSS系统观测值的权比,具体为:重复进行步骤2、3、4,即平差、方差分量估计、定权后再平差,直到
Figure BDA00023475721500000410
为止,满足各类观测值的方差因子相等的条件后,计算结束输出参数估计结果,比较计算结果得到,通过方差分量估计方法可以自适应地合理确定多系统融合定位中各个GNSS系统观测值的权比,提高多系统融合定位的精度。
惟以上所述者,仅为本发明的具体实施例而已,当不能以此限定本发明实施的范围,故其等同组件的置换,或依本发明专利保护范围所作的等同变化与修改,皆应仍属本发明权利要求书涵盖之范畴。

Claims (5)

1.一种基于方差分量估计的GNSS多系统自适应融合定位方法,其特征在于,包括以下步骤:
步骤1:对各个GNSS系统的观测值定初权;
步骤2:进行首次平差,得到各个GNSS系统观测值的残差信息;
步骤3:根据方差分量估计的方法计算各个GNSS系统方差因子并对其重新定权;
步骤4:确定多系统融合定位中各个GNSS系统观测值的权比。
2.根据权利要求1所述的基于方差分量估计的GNSS多系统自适应融合定位方法,其特征在于所述步骤1通过以下方式实现,具体为:设GPS,GLONASS与BDS导航系统的观测值分别是L1,L2和L3,其权阵分别为P1,P2和P3,且观测值之间相互独立,相应的误差方程为:
Figure FDA0002347572140000011
Figure FDA0002347572140000012
Figure FDA0002347572140000013
则三类观测值的方差和单位权方差为:
Figure FDA0002347572140000014
Figure FDA0002347572140000015
可得法方程为:
Figure FDA0002347572140000016
根据二次型期望公式可得方差分量估计公式如下:
Figure FDA0002347572140000017
在GNSS融合定位中,将观测值分为GPS、GLONASS及BDS三类,并定各自的初始权分别为P1,P2和P3
3.根据权利要求2所述的基于方差分量估计的GNSS多系统自适应融合定位方法,其特征在于所述步骤2通过以下方式实现,具体为:
Vi TPiVi(i=1,2,3)。
4.根据权利要求3所述的基于方差分量估计的GNSS多系统自适应融合定位方法,其特征在于所述步骤3通过以下方式实现,具体为:进行方差分量估计,求得各类观测值的单位权方差估值:
Figure FDA0002347572140000021
重新定权为:
Figure FDA0002347572140000022
5.根据权利要求4所述的基于方差分量估计的GNSS多系统自适应融合定位方法,其特征在于所述步骤4通过以下方式实现,确定多系统融合定位中各个GNSS系统观测值的权比,具体为:重复进行步骤2、3、4,直到
Figure FDA0002347572140000023
为止,满足条件后,计算结束输出参数估计结果。
CN201911401417.8A 2019-12-30 2019-12-30 一种基于方差分量估计的gnss多系统自适应融合定位方法 Pending CN111025356A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911401417.8A CN111025356A (zh) 2019-12-30 2019-12-30 一种基于方差分量估计的gnss多系统自适应融合定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911401417.8A CN111025356A (zh) 2019-12-30 2019-12-30 一种基于方差分量估计的gnss多系统自适应融合定位方法

Publications (1)

Publication Number Publication Date
CN111025356A true CN111025356A (zh) 2020-04-17

Family

ID=70196512

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911401417.8A Pending CN111025356A (zh) 2019-12-30 2019-12-30 一种基于方差分量估计的gnss多系统自适应融合定位方法

Country Status (1)

Country Link
CN (1) CN111025356A (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1183478A (ja) * 1997-09-04 1999-03-26 Yasushi Hirozawa 地球情報供給システム
CN107229061A (zh) * 2017-07-18 2017-10-03 武汉大学 一种基于低轨卫星的星地差分实时精密定位方法
CN108415039A (zh) * 2018-01-25 2018-08-17 中国人民解放军63921部队 iGMAS多分析中心多卫星系统精密轨道产品综合方法
CN108828640A (zh) * 2018-07-02 2018-11-16 中国人民解放军战略支援部队信息工程大学 一种卫星导航定位观测值定权方法及装置
CN110018501A (zh) * 2019-05-09 2019-07-16 哈尔滨工程大学 一种基于系统间随机模型在线估计调整的多模精密单点定位方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH1183478A (ja) * 1997-09-04 1999-03-26 Yasushi Hirozawa 地球情報供給システム
CN107229061A (zh) * 2017-07-18 2017-10-03 武汉大学 一种基于低轨卫星的星地差分实时精密定位方法
CN108415039A (zh) * 2018-01-25 2018-08-17 中国人民解放军63921部队 iGMAS多分析中心多卫星系统精密轨道产品综合方法
CN108828640A (zh) * 2018-07-02 2018-11-16 中国人民解放军战略支援部队信息工程大学 一种卫星导航定位观测值定权方法及装置
CN110018501A (zh) * 2019-05-09 2019-07-16 哈尔滨工程大学 一种基于系统间随机模型在线估计调整的多模精密单点定位方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
何帆 等: "方差分量估计在多卫星导航系统联合定位中的应用" *

Similar Documents

Publication Publication Date Title
CN108196279B (zh) 基于实时数据流的卫星钟差解算及预报方法
CN107728180B (zh) 一种基于多维粒子滤波偏差估计的gnss精密定位方法
Zhang et al. An improved robust adaptive Kalman filter for GNSS precise point positioning
CN115144876B (zh) 低轨卫星导航增强终端测试评估方法及系统
CN108333604A (zh) 一种利用卫星定位的方法和装置、卫星授时方法和装置
CN107607969B (zh) 一种基于dcb改正的四系统伪距定位方法
CN108873029B (zh) 一种实现导航接收机钟差建模的方法
CN110346823B (zh) 可用于北斗精密单点定位的三频模糊度解算方法
CN110059361B (zh) 一种基于抗差卡尔曼滤波算法的实时区域对流层建模方法
KR102180302B1 (ko) 의사거리 측정치 및 이중차분 반송파 측정치를 이용한 전리층 보정정보 생성 시스템 및 방법
CN114966760B (zh) 一种电离层加权的非差非组合ppp-rtk技术实现方法
CN110208836B (zh) 基于卡尔曼滤波的gnss高适应性周跳探测与修复方法
CN110146908B (zh) 一种虚拟参考站观测数据的生成方法
CN114859389A (zh) 一种gnss多系统抗差自适应融合rtk解算方法
CN108196267A (zh) 一种基于gnss cp技术的不间断时间传递方法
CN107748377B (zh) 基于gnss的差分定位方法及其定位系统
CN110018501A (zh) 一种基于系统间随机模型在线估计调整的多模精密单点定位方法
CN114935770B (zh) 一种多历元加快精密单点定位收敛速度的方法及装置
CN111381253A (zh) 用于估计全球导航卫星系统接收器的位置的系统及其方法
CN113568020A (zh) 一种顾及硬件频间差的卫星导航定位误差修正方法和装置
CN111025356A (zh) 一种基于方差分量估计的gnss多系统自适应融合定位方法
CN111522032A (zh) 一种北斗三号系统用户完好性处理的优化方法及优化装置
CN110398762A (zh) 实时钟差估计中的模糊度固定方法、装置、设备及介质
CN115792979A (zh) 一种基于pdop贡献度的卫星分步选星方法
CN113703017B (zh) 一种卫星天线相位中心偏差计算方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
RJ01 Rejection of invention patent application after publication

Application publication date: 20200417