CN113970772B - 一种面向城市环境的多频bds-2/bds-3/ins车载组合定位方法 - Google Patents

一种面向城市环境的多频bds-2/bds-3/ins车载组合定位方法 Download PDF

Info

Publication number
CN113970772B
CN113970772B CN202111257863.3A CN202111257863A CN113970772B CN 113970772 B CN113970772 B CN 113970772B CN 202111257863 A CN202111257863 A CN 202111257863A CN 113970772 B CN113970772 B CN 113970772B
Authority
CN
China
Prior art keywords
bds
ambiguity
frequency
observation
model
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.)
Active
Application number
CN202111257863.3A
Other languages
English (en)
Other versions
CN113970772A (zh
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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN202111257863.3A priority Critical patent/CN113970772B/zh
Publication of CN113970772A publication Critical patent/CN113970772A/zh
Application granted granted Critical
Publication of CN113970772B publication Critical patent/CN113970772B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining 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/396—Determining accuracy or reliability of position or pseudorange measurements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining 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/42—Determining position
    • G01S19/421—Determining 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/426—Determining 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 or signals derived from different modes of operation in a single system
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining 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/42—Determining position
    • G01S19/43—Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44—Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO 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/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
    • G01S19/39—Determining 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/42—Determining position
    • G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
    • G01S19/47—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement the supplementary measurement being an inertial measurement, e.g. tightly coupled inertial
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00—Road transport of goods or passengers
    • Y02T10/10—Internal combustion engine [ICE] based vehicles
    • Y02T10/40—Engine management systems

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种面向城市环境的多频BDS‑2/BDS‑3/INS车载组合定位方法，该发明利用BDS多频观测值信息（BDS‑2三频和BDS‑3四频）进行车载组合定位。面向BDS多频模糊度固定，根据城市观测环境的复杂程度，使用不同的模糊度固定策略：首先基于GF模型单历元取整固定到多频宽巷模糊度，如果观测环境较为复杂，则直接使用已经固定的宽巷模糊度参与定位解算，如果观测环境良好，允许对信号连续观测，则在已经固定的宽巷模糊度基础上再结合GB模型继续固定到窄巷模糊度后参与定位解算。最后根据不同的模糊度固定策略结果，构建相应的组合定位量测模型，既可以快速固定模糊度，又能根据观测条件达到最佳解算精度，进一步提高组合定位的稳定性和实时性。

Description

一种面向城市环境的多频BDS-2/BDS-3/INS车载组合定位 方法
技术领域
本发明属于北斗(BDS)导航定位技术领域，具体涉及一种面向城市车载动态环境的多频BDS/INS组合定位技术。
背景技术
随着自动驾驶技术的发展，城市复杂环境中的车载定位成为近年来研究的热点。由于存在复杂的信号遮挡，仅依靠全球导航卫星系统(GNSS)很难在诸如城市峡谷等复杂环境中获得稳定的高精度定位结果。惯性导航(INS)作为无源自主的导航定位系统，能够有效弥补卫星导航信号易受遮挡或干扰、导航数据更新率低以及缺乏姿态信息等方面的不足。因此，GNSS/INS组合定位可以克服各自的缺点，有效地提高组合系统的稳定性和连续性。
传统的GNSS/INS紧组合定位使用双频载波观测值和伪距观测值构建量测模型，在面对城市环境中卫星信号的复杂遮挡时，仍存在适应性不足，无法兼得定位精度和连续性的问题，从而影响组合定位的稳定性。此外，采用LAMBDA算法固定模糊度，必须基于浮点模糊度解及其协方差矩阵进行搜索，以获得正确的整数解，从而带来计算复杂性。在城市复杂环境中，GNSS信号易受干扰，一旦信号失锁或发生周跳，模糊度需要重新固定，导致定位可靠性不足，实时性难以保证。
随着北斗三号系统(BDS-3)的正式开通，可以向外公开播发四频信号的数据，多频卫星观测信息为定位性能的提升提供了新的机遇及挑战。通过引入多频GNSS观测信息，获取连续单历元高精度GNSS观测值，有望提升紧组合定位系统整体输出稳定性和可靠性。目前，国内外大部分关于车载组合定位的研究主要集中在GNSS双频观测值与INS的组合，对于多频GNSS与INS的研究甚少，尤其是城市动态环境下。因此，有必要研究面向城市环境下多频BDS-2/BDS-3与INS的车载组合定位，来进一步提升组合定位的可靠性和实时性。
发明内容
为了克服传统紧组合定位在城市环境下易受环境制约和使用LAMBDA方法实时性不高等缺陷，提高组合定位的可靠性和实时性，本发明针对北斗导航系统的特点，提出一种面向城市环境的多频BDS-2/BDS-3/INS车载组合定位方法。
为达到上述目的，本发明的技术方案如下：
一种面向城市环境的多频BDS-2/BDS-3/INS车载组合定位方法，包括以下步骤：
步骤(1)构建BDS多频线性组合载波相位和伪距双差模型，获取多频线性组合载波相位和伪距观测值；
步骤(2)基于BDS-2三频信号及BDS-3四频信号的特点和差异，根据波长、电离层和噪声因素的综合影响选择合适的观测值线性组合，解算BDS多频载波相位单历元模糊度；
步骤(3)根据城市观测环境的复杂程度，使用不同的模糊度固定策略：首先基于无几何(GF)模型单历元取整固定到多频宽巷模糊度，如果观测环境较为复杂，则直接使用已经固定的宽巷模糊度参与定位解算，如果观测环境良好，允许对信号连续观测，则在已经固定的宽巷模糊度基础上再结合几何(GB)模型继续固定到窄巷模糊度后参与定位解算；
步骤(4)针对不同的模糊度固定策略的结果建立相应的观测向量，构建组合定位系统量测模型和状态模型；
步骤(5)在步骤(4)的基础上，使用抗差滤波方法剔除城市环境下观测值出现的异常值，以消除和削弱粗差对组合定位结果的影响。
所述步骤(1)中，构建BDS多频线性组合载波相位和伪距双差模型，获取多频线性组合载波相位和伪距观测值，具体包括：
BDS多频观测值的线性组合频率、波长和双差整周模糊度分别表示为：
f(k)＝i1·f1+i2·f2+…+ik·fk (1)
其中，c表示光速；f表示频率；k表示组合系数的个数，例如三频观测值，k＝3。因此，相应的多频双差载波相位和伪距线性组合观测值分别表示为：
其中，和/>表示频率fi对应的双差载波相位和伪距观测值。在短基线下，通过双差可以消除大气层延迟误差、卫星钟差、接收机钟差，上述多频双差伪距和相位线性组合观测值相对应的观测方程，可表示为：
其中，表示双差卫地距；/>和/>分别表示/>和/>的观测值噪声。
不失一般性，设各频率的相位观测值噪声相等且独立，相应的多频双差载波相位和伪距线性组合观测值噪声精度满足：
其中，和/>分别表示等价双差载波相位和伪距观测值噪声精度，/>表示噪声放大因子，表示为：
所述步骤(2)中，基于BDS-2三频信号及BDS-3四频信号的特点和差异，根据波长、电离层和噪声因素的综合影响选择合适的观测值线性组合，解算BDS多频载波相位单历元模糊度，具体包括：
BDS多频观测值线性组合包括BDS-2三频和BDS-3四频的线性组合，其中，BDS-2三频信号分别为B1、B2和B3，BDS-3四频信号分别为B1C、B1I、B3I和B2a。通过不同形式的组合，可以有无穷多个选择。因此，需要根据一定的准则选取最优组合作为多频相位单历元模糊度固定模型，通常需要满足波长较长，电离层延迟尺度因子较小以及观测值噪声放大因子较小等条件。
通过上述的准则进行组合系数计算后，确定超宽巷、宽巷和窄巷的线性组合，采用分步无几何(GF)模型按照模糊度固定的难度依次固定超宽巷(EWL)、宽巷(WL)和窄巷(NL)模糊度。因此，基于GF模型的超宽巷、宽巷和窄巷组合模糊度的解算如下式所述：
其中，表示浮点模糊度，下标ek、wk和nk分别表示多频超宽巷、宽巷和窄巷组合；/>表示模糊度固定的相位组合观测值。需要注意的是，BDS-3四频相对于BDS-2三频，线性组合更多，在计算超宽巷/宽巷模糊度时，可以考虑计算两个超宽巷来固定宽巷模糊度。
所述步骤(3)中，根据城市观测环境的复杂程度，使用不同的模糊度固定策略：首先基于GF模型单历元取整固定到多频宽巷模糊度，如果观测环境较为复杂，则直接使用已经固定的宽巷模糊度参与定位解算，如果观测环境良好，允许对信号连续观测，则在已经固定的宽巷模糊度基础上再结合GB模型继续固定到窄巷模糊度后参与定位解算，具体包括：
在步骤(2)中，利用GF-MCAR方法依次固定超宽巷、宽巷和窄巷模糊度，但在城市环境下，受到卫星信号的频繁遮挡和多径影响，导致观测噪声很大且窄巷模糊度固定困难，可靠性难以保证。因此，要根据城市环境的复杂情况，确定不同的模糊度固定策略。当城市观测环境复杂时，仅采用GF-MCAR方法固定到宽巷模糊度，对步骤(2)的浮点模糊度舍入取整即可固定。当城市观测良好时，依然使用GF-MCAR方法固定到宽巷模糊度，使用取整的方法固定浮点模糊度，再结合GB模型继续固定到窄巷模糊度后参与定位解算，利用LAMBDA方法固定窄巷模糊度，模型表示如下：
其中，A表示系数矩阵，X表示待估位置参数。利用上述模型，最终求得原始频点的模糊度。
上述针对城市不同的观测环境，使用不同的模糊度固定策略，在观测环境复杂时，窄巷模糊度的难以固定，仅固定到宽巷，依然可以获得较高精度的观测值；在观测环境良好时，为了避免窄巷模糊度的固定不可靠，使用GF和GB组合模型进行解算，利用取整的方法固定超宽巷/宽巷，再利用LAMBDA方法固定原始频点观测值。通过上述策略，既可以获取到精度相对较高的观测值，保证定位结果的可靠性，又可以避免每一步都使用LAMBDA方法带来的计算复杂性，从而提高组合定位的实时性。
所述步骤(4)中，针对不同的模糊度固定策略的结果建立相应的观测向量，构建组合定位系统量测模型和状态模型，具体包括：
若将模糊度固定到NL，则NL模糊度即为各频点的模糊度整数解，将其代入Kalman滤波的观测方程，观测方程表示为：
Z＝HX+V (15)
其中，E是接收机到卫星的方向余弦构成的矢量矩阵；ε表示各观测量的噪声。
若将模糊度固定到WL，则观测方程的观测量为WL固定后的双差载波组合观测量与惯导预测的双差几何距离之差，因此观测方程为：
系统模型由INS的动力学模型决定，其系统连续状态方程表示如下：
其中，F表示系统的状态转移矩阵；W是系统的噪声向量；G是系统的动态噪声矩阵；X表示状态参数向量，其可以表示为：
X＝[δr,δv,φ,bg,ba] (21)
其中，δr，δv和φ分别表示位置误差、速度误差和姿态误差；bg和ba分别表示陀螺仪和加速度计的零偏。
所述步骤(5)中，在步骤(4)的基础上，使用抗差滤波方法剔除城市环境下观测值出现的异常值，以消除和削弱粗差对组合定位结果的影响，具体包括：
在城市环境下，信号受到遮挡或干扰时，容易导致信号中断，卫星观测值易出现异常值。因此，在组合定位中，将抗差模型引入到EKF中，进一步提高组合定位的定位结果和可靠性。
本发明的有益效果是：
本发明方法利用BDS多频观测值信息进行车载组合定位，有利于提高组合定位的精度。针对BDS-2三频和BDS-3四频的特点和差异，根据一定的准则选取合适的线性组合观测值。面向BDS多频模糊度固定，根据城市观测环境的复杂程度，使用不同的模糊度固定策略：首先基于GF模型单历元取整固定到多频宽巷模糊度，如果观测环境较为复杂，则直接使用已经固定的宽巷模糊度参与定位解算，如果观测环境良好，允许对信号连续观测，则在已经固定的宽巷模糊度基础上再结合GB模型继续固定到窄巷模糊度后参与定位解算，这样既可以快速固定模糊度，又能根据观测条件达到最佳解算精度。根据不同的模糊度固定策略结果，构建相应的组合定位量测模型，结合了单历元固定的实时性和多历元固定的高精度，既克服了传统基于LAMBDA方法易受环境制约导致频繁重新固定的问题，又相较全程宽巷模糊度固定提高了定位精度，进一步提高组合定位的稳定性和实时性。
附图说明
图1是本发明涉及的多频BDS-2/BDS-3/INS车载组合定位流程图。
具体实施方式
下面结合附图和具体实施方式，进一步阐明本发明，应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。
图1为本发明的流程图，本发明实施例公开的是一种面向城市环境的多频BDS-2/BDS-3/INS车载组合定位方法。
本发明的具体实施步骤如下：
步骤(1)中，构建BDS多频线性组合载波相位和伪距双差模型，获取多频线性组合载波相位和伪距观测值，具体包括：
BDS多频双差载波相位和伪距线性组合观测值分别表示为：
其中，和/>表示频率fi对应的双差载波相位和伪距观测值。在短基线下，通过双差可以消除大气层延迟误差、卫星钟差、接收机钟差，上述多频双差伪距和相位线性组合观测值相对应的观测方程，可表示为：
其中，表示双差卫地距；/>和/>分别表示/>和/>的观测值噪声。
步骤(2)中，基于BDS-2三频信号及BDS-3四频信号的特点和差异，根据波长、电离层和噪声因素的综合影响选择合适的观测值线性组合，解算BDS多频载波相位单历元模糊度，具体包括：
BDS多频观测值线性组合包括BDS-2三频和BDS-3四频的线性组合，其中，BDS-2三频信号分别为B1、B2和B3，BDS-3四频信号分别为B1C、B1I、B3I和B2a。根据一定的准则选取最优组合作为多频相位单历元模糊度固定模型，通常需要满足波长较长，电离层延迟尺度因子较小以及观测值噪声放大因子较小等条件。因此，对于BDS-2三频，可以选择和/>依次进行固定，其中伪距线性组合可以选择/>对于BDS-3四频，可以选择/>和/>依次进行固定，其中伪距线性组合可以选择/>
确定超宽巷、宽巷和窄巷的线性组合，采用分步无几何(GF)模型按照模糊度固定的难度依次固定超宽巷(EWL)、宽巷(WL)和窄巷(NL)模糊度。因此，基于GF模型的超宽巷、宽巷和窄巷组合模糊度的解算如下式所述：
其中，表示浮点模糊度，下标ek、wk和nk分别表示多频超宽巷、宽巷和窄巷组合；/>表示模糊度固定的相位组合观测值。需要注意的是，BDS-3四频相对于BDS-2三频，线性组合更多，在计算超宽巷/宽巷模糊度时，可以考虑计算两个超宽巷来固定宽巷模糊度。
步骤(3)中，根据城市环境的复杂程度，确定不同的模糊度固定策略：首先基于GF模型单历元取整固定到多频宽巷模糊度，如果观测环境较为复杂，则直接使用已经固定的宽巷模糊度参与定位解算，如果观测环境良好，允许对信号连续观测，则在已经固定的宽巷模糊度基础上再结合GB模型继续固定到窄巷模糊度后参与定位解算；具体包括：
在上述步骤(2)中，利用GF-MCAR方法依次固定超宽巷、宽巷和窄巷模糊度，但在城市环境下，受到卫星信号的频繁遮挡和多径影响，导致观测噪声很大且窄巷模糊度固定困难，可靠性难以保证。因此，要根据城市环境的复杂情况，确定不同的模糊度固定策略。当城市观测环境复杂时，仅采用GF-MCAR方法固定到宽巷模糊度，对步骤(2)的浮点模糊度舍入取整即可固定。当城市观测良好时，依然使用GF-MCAR方法固定到宽巷模糊度，使用取整的方法固定浮点模糊度，再结合GB模型继续固定到窄巷模糊度后参与定位解算，利用LAMBDA方法固定模糊度，模型表示如下：
其中，A表示系数矩阵，X表示待估位置参数。利用上述模型，最终求得原始频点的模糊度。
步骤(4)中，针对不同的模糊度固定策略的结果建立相应的观测向量，构建组合定位系统量测模型和状态模型，具体包括：
若将模糊度固定到NL，则NL模糊度即为各频点的模糊度整数解，将其代入Kalman滤波的观测方程，观测方程表示为：
Z＝HX+V (9)
其中，E是接收机到卫星的方向余弦构成的矢量矩阵；ε表示各观测量的噪声。
若将模糊度固定到WL，则观测方程的观测量为WL固定后的双差载波组合观测量与惯导预测的双差几何距离之差，因此观测方程为：
系统模型由INS的动力学模型决定，其系统连续状态方程表示如下：
其中，F表示系统的状态转移矩阵；W是系统的噪声向量；G是系统的动态噪声矩阵；X表示状态参数向量，其可以表示为：
X＝[δr,δv,φ,bg,ba] (15)
其中，δr，δv和φ分别表示位置误差、速度误差和姿态误差；bg和ba分别表示陀螺仪和加速度计的零偏。
步骤(5)中，在上述步骤(4)的基础上，使用抗差滤波方法剔除城市环境下观测值出现的异常值，以消除和削弱粗差对组合定位结果的影响，具体包括：
在城市环境下，信号受到遮挡或干扰时，容易导致信号中断，卫星观测值易出现异常值。因此，在组合定位中，将抗差模型引入到EKF中，进一步提高组合定位的定位结果和可靠性。
需要说明的是，以上内容仅仅说明了本发明的技术思想，不能以此限定本发明的保护范围，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰均落入本发明权利要求书的保护范围之内。

Claims (7)

1.一种面向城市环境的多频BDS-2/BDS-3/INS车载组合定位方法，其特征在于，包括如下步骤：
步骤(1)构建BDS多频线性组合载波相位和伪距双差模型，获取多频线性组合载波相位和伪距观测值；
步骤(2)基于BDS-2三频信号及BDS-3四频信号的特点和差异，根据波长、电离层和噪声因素的综合影响选择合适的观测值线性组合，解算BDS多频载波相位单历元模糊度；
步骤(3)根据城市观测环境的复杂程度，使用不同的模糊度固定策略：首先基于GF模型单历元取整固定到多频宽巷模糊度，如果观测环境较为复杂，则直接使用已经固定的宽巷模糊度参与定位解算，如果观测环境良好，允许对信号连续观测，则在已经固定的宽巷模糊度基础上再结合GB模型继续固定到窄巷模糊度后参与定位解算；
当城市观测环境较为复杂时，采用GF-MCAR方法固定到宽巷模糊度，对步骤(2)的浮点模糊度舍入取整即可固定；当城市观测良好时，依然使用GF-MCAR方法固定到宽巷模糊度，使用取整的方法固定浮点模糊度；为了提高定位的可靠性，对于窄巷模糊度使用GB模型进行解算，利用LAMBDA方法固定模糊度，模型表示如下：
其中，A表示系数矩阵，X表示待估位置参数；利用上述模型，最终求得原始频点的模糊度；
步骤(4)在步骤(3)的基础上针对不同的模糊度固定策略的结果建立相应的观测向量，构建BDS-2/BDS-3/INS组合定位系统量测模型和状态模型；
步骤(5)在步骤(4)的基础上，使用抗差滤波方法剔除城市环境下观测值出现的异常值，以消除和削弱粗差对组合定位结果的影响。
2.根据权利要求1所述的一种面向城市环境的多频BDS-2/BDS-3/INS车载组合定位方法，其特征在于：步骤(1)BDS多频双差载波相位和伪距线性组合观测值分别表示为：
其中，和/>表示频率fi对应的双差载波相位和伪距观测值；在短基线下，通过双差消除大气层延迟误差、卫星钟差、接收机钟差，上述多频双差载波相位和伪距线性组合观测值相对应的观测方程，表示为：
其中，表示双差卫地距；/>和/>分别表示/>和/>的观测值噪声。
3.根据权利要求1所述的一种面向城市环境的多频BDS-2/BDS-3/INS车载组合定位方法，其特征在于：步骤(2)基于BDS-2三频信号及BDS-3四频信号的特点和差异，根据波长、电离层和噪声因素的综合影响选择合适的观测值线性组合，按照一定的准则选取最优组合，通常需要满足波长较长，电离层延迟尺度因子较小以及观测值噪声放大因子较小条件。
4.根据权利要求2所述的一种面向城市环境的多频BDS-2/BDS-3/INS车载组合定位方法，其特征在于：步骤(2)解算BDS多频载波相位单历元模糊度，具体包括：
采用分步无几何(GF)模型按照模糊度固定的难度依次固定超宽巷(EWL)、宽巷(WL)和窄巷(NL)模糊度；因此，基于GF模型的超宽巷、宽巷和窄巷组合模糊度的解算如下式所述：
其中，表示浮点模糊度，下标ek、wk和nk分别表示多频超宽巷、宽巷和窄巷组合；表示模糊度固定的相位组合观测值。
5.根据权利要求4所述的一种面向城市环境的多频BDS-2/BDS-3/INS车载组合定位方法，其特征在于：步骤(4)针对不同的模糊度固定策略的结果建立相应的观测向量，构建组合定位系统量测模型，具体包括：
若将模糊度固定到NL，则NL模糊度即为各频点的模糊度整数解，将其代入滤波的观测方程，观测方程表示为：
Z＝HX+V (9)
其中，E是接收机到卫星的方向余弦构成的矢量矩阵；ε表示各观测量的噪声；
若将模糊度固定到WL，则观测方程的观测量为WL固定后的双差载波组合观测量与惯导预测的双差几何距离之差，因此观测方程为：
6.根据权利要求5所述的一种面向城市环境的多频BDS-2/BDS-3/INS车载组合定位方法，其特征在于：步骤(4)系统模型由INS的动力学模型决定，其系统连续状态方程表示如下：
其中，F表示系统的状态转移矩阵；W是系统的噪声向量；G是系统的动态噪声矩阵；X表示状态参数向量，其表示为：
X＝[δr,δv,φ,bg,ba] (15)
其中，δr，δv和φ分别表示位置误差、速度误差和姿态误差；bg和ba分别表示陀螺仪和加速度计的零偏。
7.根据权利要求1所述的一种面向城市环境的多频BDS-2/BDS-3/INS车载组合定位方法，其特征在于：步骤(5)在建立组合定位的系统模型和量测模型后，将抗差模型引入到EKF中，使用抗差滤波方法剔除城市环境下观测值出现的异常值，以消除和削弱粗差对组合定位结果的影响，进一步提高组合定位的定位结果和可靠性。
CN202111257863.3A 2021-10-27 2021-10-27 一种面向城市环境的多频bds-2/bds-3/ins车载组合定位方法 Active CN113970772B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111257863.3A CN113970772B (zh) 2021-10-27 2021-10-27 一种面向城市环境的多频bds-2/bds-3/ins车载组合定位方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111257863.3A CN113970772B (zh) 2021-10-27 2021-10-27 一种面向城市环境的多频bds-2/bds-3/ins车载组合定位方法

Publications (2)

Publication Number Publication Date
CN113970772A CN113970772A (zh) 2022-01-25
CN113970772B true CN113970772B (zh) 2024-04-12

Family

ID=79588633

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111257863.3A Active CN113970772B (zh) 2021-10-27 2021-10-27 一种面向城市环境的多频bds-2/bds-3/ins车载组合定位方法

Country Status (1)

Country Link
CN (1) CN113970772B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116299623B (zh) * 2023-05-12 2023-08-01 武汉大学 一种城市复杂场景下的ppp和ins紧组合方法与系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107703527A (zh) * 2017-08-23 2018-02-16 东南大学 一种基于北斗三频单历元宽巷/超宽巷的组合定位方法
CN108037521A (zh) * 2018-01-29 2018-05-15 东南大学 一种基于北斗超宽巷约束的bds/gps宽巷模糊度单历元固定方法
CN109212563A (zh) * 2017-06-29 2019-01-15 同济大学 北斗/gps三频周跳探测与修复方法
CN109581455A (zh) * 2019-01-24 2019-04-05 东南大学 一种bds和gps融合的三频宽巷紧组合定位方法
CN111856536A (zh) * 2020-07-30 2020-10-30 东南大学 一种基于系统间差分宽巷观测的gnss/ins紧组合定位方法
CN112462397A (zh) * 2020-11-10 2021-03-09 武汉大学 一种全星座多频北斗数据的实时动动定位方法及系统
CN112987059A (zh) * 2021-04-16 2021-06-18 湖南联智科技股份有限公司 一种基于三频模糊度解算的整周模糊度计算方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108519614A (zh) * 2018-03-16 2018-09-11 东南大学 一种gps/bds紧组合载波差分定位方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109212563A (zh) * 2017-06-29 2019-01-15 同济大学 北斗/gps三频周跳探测与修复方法
CN107703527A (zh) * 2017-08-23 2018-02-16 东南大学 一种基于北斗三频单历元宽巷/超宽巷的组合定位方法
CN108037521A (zh) * 2018-01-29 2018-05-15 东南大学 一种基于北斗超宽巷约束的bds/gps宽巷模糊度单历元固定方法
CN109581455A (zh) * 2019-01-24 2019-04-05 东南大学 一种bds和gps融合的三频宽巷紧组合定位方法
CN111856536A (zh) * 2020-07-30 2020-10-30 东南大学 一种基于系统间差分宽巷观测的gnss/ins紧组合定位方法
CN112462397A (zh) * 2020-11-10 2021-03-09 武汉大学 一种全星座多频北斗数据的实时动动定位方法及系统
CN112987059A (zh) * 2021-04-16 2021-06-18 湖南联智科技股份有限公司 一种基于三频模糊度解算的整周模糊度计算方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
BDS三频约束的多系统短基线模糊度单历元固定算法;刘扬;余学祥;王胜利;柴大帅;;导航定位学报;20180904(第03期);全文 *

Also Published As

Publication number Publication date
CN113970772A (zh) 2022-01-25

Similar Documents

Publication Publication Date Title
CN111239787B (zh) 一种集群自主协同中的gnss动态卡尔曼滤波方法
AU2007273193B2 (en) A method for increasing the reliability of position information when transitioning from a regional, wide-area, or global carrier-phase differential navigation (WADGPS) to a local real-time kinematic (RTK) navigation system
CN108802782B (zh) 一种惯导辅助的北斗三频载波相位整周模糊度求解方法
Scherzinger Precise robust positioning with inertial/GPS RTK
CN106842268B (zh) 双gnss接收机载波相位双差整周模糊度浮点解向量估计方法
CA2631479C (en) A method for combined use of a local positioning system, a local rtk system, and a regional, wide-area, or global carrier-phase positioning system
CN110531392A (zh) 一种基于ppp算法的高精度定位方法和系统
CN109765589B (zh) 一种基于无电离层组合的三频gnss实时周跳固定技术
CN106291639A (zh) 一种gnss接收机实现定位的方法及装置
CN109782313A (zh) 地面精密单点定位数据处理方法
CN108037521A (zh) 一种基于北斗超宽巷约束的bds/gps宽巷模糊度单历元固定方法
CN108051840B (zh) 一种基于gnss的含约束ekf相对定位方法
CN111856536B (zh) 一种基于系统间差分宽巷观测的gnss/ins紧组合定位方法
CN109407126A (zh) 一种多模接收机联合定位解算的方法
CN109683182A (zh) 一种引入系统间偏差的gnss多模组合差分定位方法
CN113970772B (zh) 一种面向城市环境的多频bds-2/bds-3/ins车载组合定位方法
CN115932923A (zh) 一种基于v2v的共享gnss车辆增强协作定位方法
CN109375248A (zh) 一种卡尔曼多模融合定位算法模型及其串行更新的方法
CN113064195B (zh) 一种利用多天线几何特征的高精度低计算载体测姿方法
CN105510946A (zh) 一种bds卫星载波相位整周模糊度快速解算方法
CN112444832A (zh) 一种一机多天线接收机的周跳修复方法
CN116359968B (zh) 一种联合北斗二号和北斗三号的三频差分定位方法
CN110208841B (zh) 一种改进的面向非重叠频率的gnss紧组合方法
CN113126133A (zh) 基于bds或gps的中长基线多频rtk定位的快速收敛方法
Tien et al. Adaptive strategy-based tightly-coupled INS/GNSS integration system aided by odometer and barometer

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
GR01 Patent grant
GR01 Patent grant