CN114002722A - Gnss定位方法及其装置 - Google Patents

Gnss定位方法及其装置 Download PDF

Info

Publication number
CN114002722A
CN114002722A CN202010737499.XA CN202010737499A CN114002722A CN 114002722 A CN114002722 A CN 114002722A CN 202010737499 A CN202010737499 A CN 202010737499A CN 114002722 A CN114002722 A CN 114002722A
Authority
CN
China
Prior art keywords
double
difference
carrier phase
ambiguity
observation equation
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
CN202010737499.XA
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.)
Qianxun Spatial Intelligence Inc
Original Assignee
Qianxun Spatial Intelligence Inc
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 Qianxun Spatial Intelligence Inc filed Critical Qianxun Spatial Intelligence Inc
Priority to CN202010737499.XA priority Critical patent/CN114002722A/zh
Publication of CN114002722A publication Critical patent/CN114002722A/zh
Pending legal-status Critical Current

Links

Images

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/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

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

本申请提供了一种GNSS定位方法及装置，将卫星导航信号的多个频点中的两个不同的频点组合作为第一组合和第二组合，分别采用MW组合计算宽巷模糊度和超宽巷模糊度的初始值，建立对应的双差宽巷和双差超宽巷载波相位观测方程。对初始值进行平滑得到宽巷模糊度和超宽巷模糊度的整数值，将宽巷模糊度和超宽巷模糊度的整数值代入双差宽巷载波相位观测方程和双差超宽巷载波相位观测方程。根据双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到双差载波相位观测方程中对应多个频点的模糊度的浮点解，进一步计算得到固定解。

Description

GNSS定位方法及其装置
技术领域
本说明书一般涉及导航定位技术领域，尤其涉及一种GNSS定位方法及其装置。
背景技术
由于现有的测量型或者导航型的定位终端会存在空旷场景模糊度固定慢，恶劣场景模糊度不能正常固定或者固定错误的情况，因此除了针对普通的模糊度的处理方法之外需要增加其他的处理方法，针对多频点的数据可以增加宽巷和超宽巷观测方程的约束来帮助模糊度的快速固定以及正确固定。
发明内容
为了解决现有技术中的技术问题，本申请提供了一种GNSS定位方法及其装置，可以在空旷场景以及恶劣场景，提高模糊度的固定率，提高定位精度。
本申请公开了一种GNSS定位方法，包括：
获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少三个频点上的伪距观测值及载波相位观测值；
建立对应所述至少三个频点的双差载波相位观测方程和双差伪距观测方程；
将所述卫星导航信号的至少三个频点中的两个频点作为第一组合，根据所述第一组合的卫星观测数据采用MW组合计算宽巷模糊度的初始值，并建立对应所述第一组合的双差宽巷载波相位观测方程；
将所述卫星导航信号的至少三个频点中的两个频点作为第二组合，所述第二组合和第一组合不同，根据所述第二组合的卫星观测数据采用MW组合计算超宽巷模糊度的初始值，并建立对应所述第二组合的双差超宽巷载波相位观测方程；
分别将所述宽巷模糊度和超宽巷模糊度的初始值进行多个历元的平滑并计算所述宽巷模糊度和超宽巷模糊度的整数值；
分别将所述宽巷模糊度和超宽巷模糊度的整数值代入所述双差宽巷载波相位观测方程和所述双差超宽巷载波相位观测方程；
根据所述双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少三个频点的模糊度的浮点解。
在一个优选例中，所述卫星导航信号包括GPS、北斗、伽利略卫星中至少两个导航系统的信号。
在一个优选例中，进行卡尔曼滤波时，设置所述卡尔曼滤波对应于双差宽巷载波相位观测方程的观测噪声协方差矩阵乘以第一组合系数，所述第一组合系数根据所述第一组合的两个频点的波长比获得。
在一个优选例中，进行卡尔曼滤波时，设置所述卡尔曼滤波对应于双差超宽巷载波相位观测方程的观测噪声协方差矩阵乘以第二组合系数，所述第二组合系数根据所述第二组合的两个频点的波长比获得。
在一个优选例中，还包括：利用lambda方法对所述双差载波相位观测方程中对应至少三个频点的模糊度的浮点解进行固定，获得所述双差载波相位观测方程中对应至少三个频点的模糊度的固定解。
本申请还公开了一种GNSS定位方法，包括：
获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少两个频点上的伪距观测值及载波相位观测值；
建立对应所述至少两个频点的双差载波相位观测方程和双差伪距观测方程；
将所述卫星导航信号的至少两个频点中的两个频点作为第一组合，根据所述第一组合的卫星观测数据采用MW组合计算宽巷模糊度的初始值，并建立对应所述第一组合的双差宽巷载波相位观测方程；
将所述宽巷模糊度的初始值进行多个历元的平滑并计算所述宽巷模糊度的整数值；
将所述宽巷模糊度的整数值代入所述双差宽巷载波相位观测方程；
根据所述双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解。
在一个优选例中，所述卫星导航信号包括GPS、北斗、伽利略卫星中至少两个导航系统的信号。
在一个优选例中，进行卡尔曼滤波时，设置所述卡尔曼滤波对应于双差宽巷载波相位观测方程的观测噪声协方差矩阵乘以第一组合系数，所述第一组合系数根据所述第一组合的两个频点的波长比获得。
在一个优选例中，还包括：利用lambda方法对所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解进行固定，获得所述双差载波相位观测方程中对应至少两个频点的模糊度的固定解。
本申请还公开了一种GNSS定位方法，包括：
获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少两个频点上的伪距观测值及载波相位观测值；
建立对应所述至少两个频点的双差载波相位观测方程和双差伪距观测方程；
将所述卫星导航信号的至少两个频点中的两个频点作为第二组合，根据所述第二组合的卫星观测数据采用MW组合计算超宽巷模糊度的初始值，并建立对应所述第二组合的双差超宽巷载波相位观测方程；
将所述超宽巷模糊度的初始值进行多个历元的平滑并计算所述超宽巷模糊度的整数值；
将所述超宽巷模糊度的整数值代入所述双差超宽巷载波相位观测方程；
根据所述双差载波相位观测方程、双差伪距观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解。
在一个优选例中，所述卫星导航信号包括GPS、北斗、伽利略卫星中至少两个导航系统的信号。
在一个优选例中，进行卡尔曼滤波时，设置所述卡尔曼滤波对应于双差超宽巷载波相位观测方程的观测噪声协方差矩阵乘以第二组合系数，所述第二组合系数根据所述第二组合的两个频点的波长比获得。
在一个优选例中，还包括：利用lambda方法对所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解进行固定，获得所述双差载波相位观测方程中对应至少两个频点的模糊度的固定解。
本申请还公开了一种GNSS定位装置，包括：
获取模块，用于获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少三个频点上的伪距观测值及载波相位观测值；
宽巷模糊度计算模块，用于将所述卫星导航信号的至少三个频点中的两个频点作为第一组合，根据所述第一组合的卫星观测数据采用MW组合计算宽巷模糊度的初始值，并建立对应所述第一组合的双差宽巷载波相位观测方程；
超宽巷模糊度计算模块，用于将所述卫星导航信号的至少三个频点中的两个频点作为第二组合，所述第二组合和第一组合不同，根据所述第二组合的卫星观测数据采用MW组合计算超宽巷模糊度的初始值，并建立对应所述第二组合的双差超宽巷载波相位观测方程；
平滑取整模块，用于分别将所述宽巷模糊度和超宽巷模糊度的初始值进行多个历元的平滑并计算所述宽巷模糊度和超宽巷模糊度的整数值；
滤波模块，用于分别将所述宽巷模糊度和超宽巷模糊度的整数值代入所述双差宽巷载波相位观测方程和所述双差超宽巷载波相位观测方程，用于建立对应所述至少三个频点的双差载波相位观测方程和双差伪距观测方程，根据所述双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少三个频点的模糊度的浮点解。
本申请还公开了一种GNSS定位装置，包括：
获取模块，用于获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少两个频点上的伪距观测值及载波相位观测值；
宽巷模糊度计算模块，用于将所述卫星导航信号的至少两个频点中的两个频点作为第一组合，根据所述第一组合的卫星观测数据采用MW组合计算宽巷模糊度的初始值，并建立对应所述第一组合的双差宽巷载波相位观测方程；
平滑取整模块，用于将所述宽巷模糊度的初始值进行多个历元的平滑并计算所述宽巷模糊度的整数值；
滤波模块，用于将所述宽巷模糊度的整数值代入所述双差宽巷载波相位观测方程，用于建立对应所述至少两个频点的双差载波相位观测方程和双差伪距观测方程，根据所述双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解。
本申请还公开了一种GNSS定位装置，包括：
获取模块，用于获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少两个频点上的伪距观测值及载波相位观测值；
超宽巷模糊度计算模块，用于将所述卫星导航信号的至少两个频点中的两个频点作为第二组合，根据所述第二组合的卫星观测数据采用MW组合计算超宽巷模糊度的初始值，并建立对应所述第二组合的双差超宽巷载波相位观测方程；
平滑取整模块，用于将所述超宽巷模糊度的初始值进行多个历元的平滑并计算所述超宽巷模糊度的整数值；
滤波模块，用于将所述超宽巷模糊度的整数值代入所述双差超宽巷载波相位观测方程，用于建立对应所述至少两个频点的双差载波相位观测方程和双差伪距观测方程，根据所述双差载波相位观测方程、双差伪距观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解。
相对于现有技术，本申请具有以下有益效果：
本发明中，分别计算宽巷模糊度和/或超宽巷模糊度，并将宽巷模糊度和超宽巷模糊度带入载波相位观测方程中，将原有的宽巷模糊度和超宽巷模糊度的观测方程转换为精度高于伪距方程的载波相位观测方程，来参与整个滤波过程和模糊度固定的过程，从而在空旷场景和恶劣场景下可以提高模糊度的固定率和固定正确率，进而提高定位精度。
本说明书中记载了大量的技术特征，分布在各个技术方案中，如果要罗列出本申请所有可能的技术特征的组合(即技术方案)的话，会使得说明书过于冗长。为了避免这个问题，本说明书上述发明内容中公开的各个技术特征、在下文各个实施方式和例子中公开的各技术特征、以及附图中公开的各个技术特征，都可以自由地互相组合，从而构成各种新的技术方案(这些技术方案均应该视为在本说明书中已经记载)，除非这种技术特征的组合在技术上是不可行的。例如，在一个例子中公开了特征A+B+C，在另一个例子中公开了特征A+B+D+E，而特征C和D是起到相同作用的等同技术手段，技术上只要择一使用即可，不可能同时采用，特征E技术上可以与特征C相组合，则，A+B+C+D的方案因技术不可行而应当不被视为已经记载，而A+B+C+E的方案应当视为已经被记载。
附图说明
参考以下附图描述本申请的非限制性和非穷举性实施例，其中除非另有说明，否则相同的附图标记在各个视图中指代相同的部分。
图1是根据本说明书实施例一中一种GNSS定位方法的流程图。
图2是根据本说明书实施例二中一种GNSS定位方法的流程图。
图3是根据本说明书实施例三中一种GNSS定位方法的流程图。
具体实施方式
在以下的叙述中，为了使读者更好地理解本申请而提出了许多技术细节。但是，本领域的普通技术人员可以理解，即使没有这些技术细节和基于以下各实施方式的种种变化和修改，也可以实现本申请所要求保护的技术方案。
部分概念的说明：
整周模糊度(Ambiguity)：又称整周未知数，是在全球定位系统技术的载波相位测量时，载波相位与基准相位之间相位差的首观测值所对应的整周未知数。
卡尔曼滤波：一种利用线性系统状态方程，通过系统输入输出观测数据，对系统状态进行最优估计的算法。由于观测数据中包括系统中的噪声和干扰的影响，所以最优估计也可看作是滤波过程。
MW组合：MW组合是利用双频伪距和相位两类观测值求解模糊度的线性组合。
本申请的部分创新在于：
本发明中，分别将卫星导航信号的多个频点中的两个不同的频点组合作为第一组合和第二组合，分别采用MW组合分别计算宽巷模糊度和超宽巷模糊度的初始值，并分别建立对应的双差宽巷和双差超宽巷载波相位观测方程。对初始值进行平滑得到宽巷模糊度和超宽巷模糊度的整数值，分别将宽巷模糊度和超宽巷模糊度的整数值代入双差宽巷载波相位观测方程和双差超宽巷载波相位观测方程，根据双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到双差载波相位观测方程中对应多个频点的模糊度的浮点解，并进一步计算得到模糊度的固定解。
本发明中，将原有的宽巷模糊度和超宽巷模糊度的观测方程转换为精度高于伪距方程的观测方程，来参与整个滤波过程和模糊度固定的过程，在空旷场景和恶劣场景可以提高模糊度的固定率和固定正确率，进而提高定位精度。
通过添加宽巷和超宽巷的处理方法明显提高空旷场景的模糊度固定率，以及模糊度的固定正确率，基本可以做到单历元模糊度固定，以至于可以输出高精度的固定解。在恶劣场景下，可以提高模糊度的固定率，保证大部分情况下输出固定解，以提高整体的定位精度。
为使本申请的目的、技术方案和优点更加清楚，下面将结合附图对本申请的实施方式作进一步地详细描述。
实施例一
本申请的实施例一中公开了一种GNSS定位方法，图1示出了GNSS定位方法的流程图，包括：
步骤101，获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少三个频点上的伪距观测值及载波相位观测值。在一个实施例中，所述卫星导航信号包括GPS、北斗(BDS)、伽利略(GLA)卫星中至少两个导航系统的信号，例如，包括GPS和北斗的卫星导航信号，或包括北斗和伽利略的卫星导航信号。
步骤102，建立对应所述至少三个频点的双差载波相位观测方程和双差伪距观测方程。
步骤103，将所述卫星导航信号的至少三个频点中的两个频点作为第一组合，根据所述第一组合的卫星观测数据采用MW组合计算宽巷模糊度的初始值，并建立对应所述第一组合的双差宽巷载波相位观测方程。
步骤104，将所述卫星导航信号的至少三个频点中的两个频点作为第二组合，所述第二组合和第一组合不同，根据所述第二组合的卫星观测数据采用MW组合计算超宽巷模糊度的初始值，并建立对应所述第二组合的双差超宽巷载波相位观测方程。第二组合得到载波相位波长较第一组合的宽巷更长。
其中，第一组合中的频点和第二组合中的频点可以完全不同，可以部分相同，例如，其中一个频点相同。
步骤105，分别将所述宽巷模糊度和超宽巷模糊度的初始值进行多个历元的平滑并计算所述宽巷模糊度和超宽巷模糊度的整数值。多个历元例如可以为至少2个。
步骤106，分别将所述宽巷模糊度和超宽巷模糊度的整数值代入所述双差宽巷载波相位观测方程和所述双差超宽巷载波相位观测方程。
步骤107，根据所述双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少三个频点的模糊度的浮点解。
在一个实施例中，进行卡尔曼滤波时，设置所述卡尔曼滤波对应于双差宽巷载波相位观测方程的观测噪声协方差矩阵乘以第一组合系数，所述第一组合系数根据所述第一组合的两个频点的波长比获得。
在一个实施例中，进行卡尔曼滤波时，设置所述卡尔曼滤波对应于双差超宽巷载波相位观测方程的观测噪声协方差矩阵乘以第二组合系数，所述第二组合系数根据所述第二组合的两个频点的波长比获得。
在一个实施例中，还包括：利用lambda方法对所述双差载波相位观测方程中对应至少三个频点的模糊度的浮点解进行固定，获得所述双差载波相位观测方程中对应至少三个频点的模糊度的固定解。
实施例二
本申请的二实施例中公开了一种GNSS定位方法，图2示出了GNSS定位方法的流程图，实施例二与实施例一基本相同，主要区别在于：实施例一中计算宽巷模糊度和超宽巷模糊度的整数值，并结合包含宽巷模糊度和超宽巷模糊度的整数值的双差载波相位观测方程进行卡尔曼滤波，获得模糊度的浮点解，然而，实施例二中仅需计算宽巷模糊度的整数值，结合包含宽巷模糊度的整数值的双差载波相位观测方法进行卡尔曼滤波，获得模糊度的浮点解。具体的，包括如下步骤：
步骤201，获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少两个频点上的伪距观测值及载波相位观测值。在一个实施例中，所述卫星导航信号包括GPS、北斗、伽利略卫星中至少两个导航系统的信号。
步骤202，建立对应所述至少两个频点的双差载波相位观测方程和双差伪距观测方程。
步骤203，将所述卫星导航信号的至少两个频点中的两个频点作为第一组合，根据所述第一组合的卫星观测数据采用MW组合计算宽巷模糊度的初始值，并建立对应所述第一组合的双差宽巷载波相位观测方程。
步骤204，将所述宽巷模糊度的初始值进行多个历元的平滑并计算所述宽巷模糊度的整数值。多个历元例如可以为至少2个。
步骤205，将所述宽巷模糊度的整数值代入所述双差宽巷载波相位观测方程。
步骤206，根据所述双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解。
类似的，进行卡尔曼滤波时，设置所述卡尔曼滤波对应于双差宽巷载波相位观测方程的观测噪声协方差矩阵乘以第一组合系数，所述第一组合系数根据所述第一组合的两个频点的波长比获得。
在一个实施例中，还包括：利用lambda方法对所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解进行固定，获得所述双差载波相位观测方程中对应至少两个频点的模糊度的固定解。
此外，在实施例二中，仅需计算对应一个频点组合的宽巷模糊度的整数值，不需要同时计算宽巷和超宽巷模糊度的整数值，因此，仅需获取至少两个频点上的观测值，不需要获取至少三个频点上的观测值。
实施例三
本申请的实施例三中公开了一种GNSS定位方法，图3示出了GNSS定位方法的流程图，实施例三与实施例一基本相同，主要区别在于：实施例三中仅需计算超宽巷模糊度的整数值，结合包含超宽巷模糊度的整数值的双差载波相位观测方法进行卡尔曼滤波，获得模糊度的浮点解。具体的，包括如下步骤：
步骤301，获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少两个频点上的伪距观测值及载波相位观测值。在一个实施例中，所述卫星导航信号包括GPS、北斗、伽利略卫星中至少两个导航系统的信号。
步骤302，建立对应所述至少两个频点的双差载波相位观测方程和双差伪距观测方程。
步骤303，将所述卫星导航信号的至少两个频点中的两个频点作为第二组合，根据所述第二组合的卫星观测数据采用MW组合计算超宽巷模糊度的初始值，并建立对应所述第二组合的双差超宽巷载波相位观测方程。
步骤304，将所述超宽巷模糊度的初始值进行多个历元的平滑并计算所述超宽巷模糊度的整数值。多个历元例如可以为至少2个。
步骤305，将所述超宽巷模糊度的整数值代入所述双差超宽巷载波相位观测方程。
步骤306，根据所述双差载波相位观测方程、双差伪距观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解。
在一个实施例中，进行卡尔曼滤波时，设置所述卡尔曼滤波对应于双差超宽巷载波相位观测方程的观测噪声协方差矩阵乘以第二组合系数，所述第二组合系数根据所述第二组合的两个频点的波长比获得。
在一个实施例中，还包括：利用lambda方法对所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解进行固定，获得所述双差载波相位观测方程中对应至少两个频点的模糊度的固定解。
此外，在实施例三中，仅需计算对应一个频点组合的超宽巷模糊度的整数值，不需要同时计算宽巷和超宽巷模糊度的整数值，因此，仅需获取至少两个频点上的观测值，不需要获取至少三个频点上的观测值。
实施例四
本申请的实施例四中还公开了一种GNSS定位装置，包括：
获取模块，用于获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少三个频点上的伪距观测值及载波相位观测值；
宽巷模糊度计算模块，用于将所述卫星导航信号的至少三个频点中的两个频点作为第一组合，根据所述第一组合的卫星观测数据采用MW组合计算宽巷模糊度的初始值，并建立对应所述第一组合的双差宽巷载波相位观测方程；
超宽巷模糊度计算模块，用于将所述卫星导航信号的至少三个频点中的两个频点作为第二组合，所述第二组合和第一组合不同，根据所述第二组合的卫星观测数据采用MW组合计算超宽巷模糊度的初始值，并建立对应所述第二组合的双差超宽巷载波相位观测方程；
平滑取整模块，用于分别将所述宽巷模糊度和超宽巷模糊度的初始值进行多个历元的平滑并计算所述宽巷模糊度和超宽巷模糊度的整数值；
滤波模块，用于分别将所述宽巷模糊度和超宽巷模糊度的整数值代入所述双差宽巷载波相位观测方程和所述双差超宽巷载波相位观测方程，用于建立对应所述至少三个频点的双差载波相位观测方程和双差伪距观测方程，根据所述双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少三个频点的模糊度的浮点解。
实施例五
本申请的实施例五中还公开了一种GNSS定位装置，包括：
获取模块，用于获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少两个频点上的伪距观测值及载波相位观测值；
宽巷模糊度计算模块，用于将所述卫星导航信号的至少两个频点中的两个频点作为第一组合，根据所述第一组合的卫星观测数据采用MW组合计算宽巷模糊度的初始值，并建立对应所述第一组合的双差宽巷载波相位观测方程；
平滑取整模块，用于将所述宽巷模糊度的初始值进行多个历元的平滑并计算所述宽巷模糊度的整数值；
滤波模块，用于将所述宽巷模糊度的整数值代入所述双差宽巷载波相位观测方程，用于建立对应所述至少两个频点的双差载波相位观测方程和双差伪距观测方程，根据所述双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解。
实施例五
本申请的实施例五中还公开了一种GNSS定位装置，包括：
获取模块，用于获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少两个频点上的伪距观测值及载波相位观测值；
超宽巷模糊度计算模块，用于将所述卫星导航信号的至少两个频点中的两个频点作为第二组合，根据所述第二组合的卫星观测数据采用MW组合计算超宽巷模糊度的初始值，并建立对应所述第二组合的双差超宽巷载波相位观测方程；
平滑取整模块，用于将所述超宽巷模糊度的初始值进行多个历元的平滑并计算所述超宽巷模糊度的整数值；
滤波模块，用于将所述超宽巷模糊度的整数值代入所述双差超宽巷载波相位观测方程，用于建立对应所述至少两个频点的双差载波相位观测方程和双差伪距观测方程，根据所述双差载波相位观测方程、双差伪距观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解。
为了能够更好地理解本说明书的技术方案，下面结合一个具体的例子来进行说明，该例子中罗列的细节主要是为了便于理解，不作为对本申请保护范围的限制。
1.接收卫星数据
对使用导航型或者测量型终端(流动站)采集的卫星数据进行数据格式转换，提取卫星的伪距观测值、载波相位观测值、卫星信噪比和多普勒观测值，并获取相同时间段的固定基准站的伪距观测值、载波相位观测值以及卫星星历。
2.卫星数据预处理
根据伪距观测值的残差对所述伪距观测值进行精度评定和单点定位，对载波相位观测值进行周跳探测及修复，对星座完好性进行分析和监测。
3.建立卫星数据差分观测方程与误差方程
伪距差分采用伪距观测量作为基本输入，设基准站和流动站同步观测一组导航卫星，获取基准站(A站)和流动站(B站)同步观测的n颗共视星的伪距观测量
Figure BDA0002605497970000151
和
Figure BDA0002605497970000152
选取这n颗导航卫星中仰角最高的导航卫星r作为参考星，以组建关于任意导航卫星j(j＝1,2,3......n且j≠r)和参考星r在任意ti时刻的伪距观测方程如下：
Figure BDA0002605497970000153
Figure BDA0002605497970000154
上式中各参数含义如下：
c：光速(m/s)；
λ：卫星导航信号载波波长(m)；
f：卫星导航信号载波频率(Hz)；
Figure BDA0002605497970000155
ti时刻基准站观测的导航卫星r的伪距观测值(m)；
Figure BDA0002605497970000156
ti时刻基准站观测的导航卫星j的伪距观测值(m)；
Figure BDA0002605497970000157
ti时刻基准站与导航卫星r间的几何距离(m)；
Figure BDA0002605497970000158
ti时刻基准站与导航卫星j间的几何距离(m)；
δtA(ti)：ti时刻基准站接收机钟差(s)；
δtr(ti)：ti时刻导航卫星r的卫星钟差(s)；
δtj(ti)：ti时刻导航卫星j的卫星钟差(s)；
Figure BDA0002605497970000161
ti时刻基准站与导航卫星r间的地球旋转误差(m)；
Figure BDA0002605497970000162
ti时刻基准站与导航卫星j间的地球旋转误差(m)；
Figure BDA0002605497970000163
ti时刻基准站与导航卫星r间的电离层延迟误差(m)；
Figure BDA0002605497970000164
ti时刻基准站与导航卫星j间的电离层延迟误差(m)；
Figure BDA0002605497970000165
ti时刻基准站与导航卫星r间的对流层延迟误差(m)；
Figure BDA0002605497970000166
ti时刻基准站与导航卫星j间的对流层延迟误差(m)；
Figure BDA0002605497970000167
ti时刻基准站接收机关于导航卫星r的伪距测量热噪声(m)；
Figure BDA0002605497970000168
ti时刻基准站接收机关于导航卫星j的伪距测量热噪声(m)；
Figure BDA0002605497970000169
ti时刻流动站观测的导航卫星r的伪距观测值(m)；
Figure BDA00026054979700001610
ti时刻流动站观测的导航卫星j的伪距观测值(m)；
Figure BDA00026054979700001611
ti时刻流动站与导航卫星r间的几何距离(m)；
Figure BDA00026054979700001612
ti时刻流动站与导航卫星j间的几何距离(m)；
δtB(ti)：ti时刻流动站接收机钟差(s)；
Figure BDA00026054979700001613
ti时刻流动站与导航卫星r间的地球旋转误差(m)；
Figure BDA00026054979700001614
ti时刻流动站与导航卫星j间的地球旋转误差(m)；
Figure BDA00026054979700001615
ti时刻流动站与导航卫星r间的电离层延迟误差(m)；
Figure BDA00026054979700001616
ti时刻流动站与导航卫星j间的电离层延迟误差(m)；
Figure BDA00026054979700001617
ti时刻流动站与导航卫星r间的对流层延迟误差(m)；
Figure BDA00026054979700001618
ti时刻流动站与导航卫星j间的对流层延迟误差(m)；
Figure BDA00026054979700001619
ti时刻流动站接收机关于导航卫星r的伪距测量热噪声(m)；
Figure BDA00026054979700001620
ti时刻流动站接收机关于导航卫星j的伪距测量热噪声(m)；
设基准站和流动站同步观测一组导航卫星，获取基准站和流动站同步观测的n颗共视星的载波相位观测量
Figure BDA00026054979700001621
和
Figure BDA00026054979700001622
选取这n颗导航卫星中仰角最高的导航卫星r作为参考星，以组建关于任意导航卫星j(j＝1,2,3......n且j≠r)和参考星r在任意ti时刻的载波相位观测方程如下：
Figure BDA0002605497970000171
Figure BDA0002605497970000172
上式中各参数含义如下：
c：光速(m/s)；
λ：卫星导航信号载波波长(m)；
f：卫星导航信号载波频率(Hz)；
Figure BDA0002605497970000173
ti时刻基准站观测的导航卫星r的载波相位观测值(周)；
Figure BDA0002605497970000174
ti时刻基准站观测的导航卫星j的载波相位观测值(周)；
Figure BDA0002605497970000175
ti时刻基准站与导航卫星r间的整周模糊度(周)；
Figure BDA0002605497970000176
ti时刻基准站与导航卫星j间的整周模糊度(周)；
Figure BDA0002605497970000177
ti时刻基准站与导航卫星r间的几何距离(m)；
Figure BDA0002605497970000178
ti时刻基准站与导航卫星j间的几何距离(m)；
δtA(ti)：ti时刻基准站接收机钟差(s)；
δtr(ti)：ti时刻导航卫星r的卫星钟差(s)；
δtj(ti)：ti时刻导航卫星j的卫星钟差(s)；
Figure BDA0002605497970000179
ti时刻基准站与导航卫星r间的地球旋转误差(m)；
Figure BDA00026054979700001710
ti时刻基准站与导航卫星j间的地球旋转误差(m)；
Figure BDA00026054979700001711
ti时刻基准站与导航卫星r间的电离层延迟误差(m)；
Figure BDA00026054979700001712
ti时刻基准站与导航卫星j间的电离层延迟误差(m)；
Figure BDA00026054979700001713
ti时刻基准站与导航卫星r间的对流层延迟误差(m)；
Figure BDA00026054979700001714
ti时刻基准站与导航卫星j间的对流层延迟误差(m)；
Figure BDA00026054979700001715
ti时刻基准站接收机关于导航卫星r的载波测量热噪声(m)；
Figure BDA00026054979700001716
ti时刻基准站接收机关于导航卫星j的载波测量热噪声(m)；
Figure BDA00026054979700001717
ti时刻流动站观测的导航卫星r的载波相位观测值(m)；
Figure BDA00026054979700001718
ti时刻流动站观测的导航卫星j的载波相位观测值(m)；
Figure BDA00026054979700001719
ti时刻流动站与导航卫星r间的整周模糊度(周)；
Figure BDA0002605497970000181
ti时刻流动站与导航卫星j间的整周模糊度(周)；
Figure BDA0002605497970000182
ti时刻流动站与导航卫星r间的几何距离(m)；
Figure BDA0002605497970000183
ti时刻流动站与导航卫星j间的几何距离(m)；
δtB(ti)：ti时刻流动站接收机钟差(s)；
Figure BDA0002605497970000184
ti时刻流动站与导航卫星r间的地球旋转误差(m)；
Figure BDA0002605497970000185
ti时刻流动站与导航卫星j间的地球旋转误差(m)；
Figure BDA0002605497970000186
ti时刻流动站与导航卫星r间的电离层延迟误差(m)；
Figure BDA0002605497970000187
ti时刻流动站与导航卫星j间的电离层延迟误差(m)；
Figure BDA0002605497970000188
ti时刻流动站与导航卫星r间的对流层延迟误差(m)；
Figure BDA0002605497970000189
ti时刻流动站与导航卫星j间的对流层延迟误差(m)；
Figure BDA00026054979700001810
ti时刻流动站接收机关于导航卫星r的伪距测量热噪声(m)；
Figure BDA00026054979700001811
ti时刻流动站接收机关于导航卫星j的伪距测量热噪声(m)；
根据基础的观测方程列取伪距和载波相位双差观测方程。
4.计算宽巷和超宽巷模糊度
本示例中主要采用MW组合的方式计算宽巷和超宽巷模糊度，这样做既可以节省计算量同时也提高了模糊度准确度，具体步骤如下：
(1)根据基准站坐标计算基准站观测卫星的位置，采用单点定位获得流动站概略坐标，并通过伪距双差解算获得流动站较为精确的初始坐标。
(2)利用MW组合计算宽巷模糊度的初始值，其中MW组合观测值如下：
Figure BDA00026054979700001812
其中，
Figure BDA00026054979700001813
为宽巷组合观测值；
NWL＝NL1-NL2为宽巷组合模糊度。
式中，fL1、fL2、λL1、λL2及PL1、PL2分别为该卫星系统的第一频点与第二频点的波段的频率、波长及伪距观测值，其中，第一频点和第二频点分别为GPS的L1L2，BDS的B1B2，GLA的E1E5a这三种组合形式，使用上面的方法将所有可见卫星的宽巷模糊度正确计算出来，以便后面卡尔曼滤波使用。
(3)类似的，根据公式4-1计算超宽巷模糊度的初始值，式中，改换fL1、fL2、λL1、λL2及PL1、PL2分别为该卫星系统的第一频点与第三频点的波段的频率、波长及伪距观测值。其中，第一频点和第二频点分别为GPS的L1L5，BDS的B2B3，GLA的E5aE5b这三种组合形式，这样组合后的波长较长，通过计算得到的模糊度正确率较高。
(4)通过前面的求解可以得到所有卫星的宽巷和超宽巷模糊度，每次在使用这个宽巷与超宽巷模糊度的时候需要与前面一个历元进行平滑计算，得到平滑之后的宽巷和超宽巷模糊度，这样使计算出来的模糊度更加准确。
(5)使用平滑后的宽巷和超宽巷模糊度进行四舍五入得到整数模糊度，然后将宽巷和超宽巷模糊度的整数值反向带回前面第3部分列取的宽巷和超宽巷的观测方程中，得到包含宽巷和超宽巷模糊度的高精度观测方程，这些观测方程会与原来的不组合的频点进行联立放入滤波器进行卡尔曼滤波。
5.卡尔曼滤波
离散线性系统的卡尔曼滤波中状态方程和观测方程为：
Xk+1＝φk+1,kXk+Γk+1,kΩk (5-1)
Lk+1＝Bk+1Xk+1+Δk+1 (5-2)
式中Xk+1为k+1时刻的状态向量，同时也是待估向量，所述待估向量包括坐标参数及宽巷和超宽巷模糊度，Lk+1为观测向量，φk+1,k为k时刻变化到k+1时刻的状态转移矩阵，Γk+1,k、Bk+1为系数矩阵，Ωk为过程噪声，Δk+1为观测噪声。
卡尔曼滤波的主要步骤为：
(1)利用k-1时刻的
Figure BDA0002605497970000191
计算k时刻的状态向量预测值
Figure BDA0002605497970000192
Figure BDA0002605497970000193
(2)利用k-1时刻的方差协方差阵DX(k-1,k-1)计算k时刻的状态向量预测值的方差协方差阵DX(k,k-1)：
Figure BDA0002605497970000201
(3)计算状态向量预测值的残差向量vk：
Figure BDA0002605497970000202
(4)计算预测值残差向量vk的方差协方差阵Dvk：
Figure BDA0002605497970000203
(5)计算滤波增益矩阵：
Figure BDA0002605497970000204
(6)计算k时刻状态向量估计量：
Figure BDA0002605497970000205
(7)计算k时刻状态向量估计量的方差协方差阵：
DX(k,k)＝DX(k,k-1)-JkBkDX(k,k-1) (5-9)
根据观测方程、误差方程以及差分数据的特点，所求未知数为双差模糊度，正常设置滤波的初始值矩阵、滤波初值方差矩阵、过程噪声矩阵，并正常进行卡尔曼滤波，得到滤波后的模糊度的浮点解。
根据前面求解的宽巷和超宽巷的观测方程，在设置观测噪声协方差阵的时候，需要进行附加的处理，以保证设置的正确性，即：进行卡尔曼滤波计算得到宽巷模糊度的浮点解时，观测噪声协方差矩阵乘以第一频点组合的两个频点的波长比，例如L1/L2＝1/2；进行卡尔曼滤波计算得到超宽巷模糊度的浮点解时，观测噪声协方差矩阵乘以第二频点组合的两个频点的波长比，例如L2/L5＝2/5。
卡尔曼滤波器中除了原来的非组合的观测方程外，增加前面已经求解出来的宽巷和超宽巷的观测方程，相当于滤波器中增加了宽巷和超宽巷的方程约束，该方法中可以不使用无电离层组合的方法进行，只是将宽巷模糊度的观测方程与基频的方程进行联合求解，对于运算的流程和运算量的优化是比较明显的，由于这个方程的精度相对来讲比较高，因为可以对原来的方程进行约束，进而提高模糊度的固定率。
需要说明的是，在本专利的申请文件中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。本专利的申请文件中，如果提到根据某要素执行某行为，则是指至少根据该要素执行该行为的意思，其中包括了两种情况：仅根据该要素执行该行为、和根据该要素和其它要素执行该行为。多个、多次、多种等表达包括2个、2次、2种以及2个以上、2次以上、2种以上。
在本说明书提及的所有文献都被认为是整体性地包括在本说明书的公开内容中，以便在必要时可以作为修改的依据。此外应理解，以上所述仅为本说明书的较佳实施例而已，并非用于限定本说明书的保护范围。凡在本说明书一个或多个实施例的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本说明书一个或多个实施例的保护范围之内。
在一些情况下，在权利要求书中记载的动作或步骤可以按照不同于实施例中的顺序来执行并且仍然可以实现期望的结果。另外，在附图中描述的过程不一定要求示出的特定顺序或者连续顺序才能实现期望的结果。在某些实施方式中，多任务处理和并行处理也是可以的或者可能是有利的。

Claims (16)

1.一种GNSS定位方法，其特征在于，包括：
获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少三个频点上的伪距观测值及载波相位观测值；
建立对应所述至少三个频点的双差载波相位观测方程和双差伪距观测方程；
将所述卫星导航信号的至少三个频点中的两个频点作为第一组合，根据所述第一组合的卫星观测数据采用MW组合计算宽巷模糊度的初始值，并建立对应所述第一组合的双差宽巷载波相位观测方程；
将所述卫星导航信号的至少三个频点中的两个频点作为第二组合，所述第二组合和第一组合不同，根据所述第二组合的卫星观测数据采用MW组合计算超宽巷模糊度的初始值，并建立对应所述第二组合的双差超宽巷载波相位观测方程；
分别将所述宽巷模糊度和超宽巷模糊度的初始值进行多个历元的平滑并计算所述宽巷模糊度和超宽巷模糊度的整数值；
分别将所述宽巷模糊度和超宽巷模糊度的整数值代入所述双差宽巷载波相位观测方程和所述双差超宽巷载波相位观测方程；
根据所述双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少三个频点的模糊度的浮点解。
2.如权利要求1所述的GNSS定位方法，其特征在于，所述卫星导航信号包括GPS、北斗、伽利略卫星中至少两个导航系统的信号。
3.如权利要求1所述的GNSS定位方法，其特征在于，进行卡尔曼滤波时，设置所述卡尔曼滤波对应于双差宽巷载波相位观测方程的观测噪声协方差矩阵乘以第一组合系数，所述第一组合系数根据所述第一组合的两个频点的波长比获得。
4.如权利要求1或3所述的GNSS定位方法，其特征在于，进行卡尔曼滤波时，设置所述卡尔曼滤波对应于双差超宽巷载波相位观测方程的观测噪声协方差矩阵乘以第二组合系数，所述第二组合系数根据所述第二组合的两个频点的波长比获得。
5.如权利要求1所述的GNSS定位方法，其特征在于，还包括：利用lambda方法对所述双差载波相位观测方程中对应至少三个频点的模糊度的浮点解进行固定，获得所述双差载波相位观测方程中对应至少三个频点的模糊度的固定解。
6.一种GNSS定位方法，其特征在于，包括：
获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少两个频点上的伪距观测值及载波相位观测值；
建立对应所述至少两个频点的双差载波相位观测方程和双差伪距观测方程；
将所述卫星导航信号的至少两个频点中的两个频点作为第一组合，根据所述第一组合的卫星观测数据采用MW组合计算宽巷模糊度的初始值，并建立对应所述第一组合的双差宽巷载波相位观测方程；
将所述宽巷模糊度的初始值进行多个历元的平滑并计算所述宽巷模糊度的整数值；
将所述宽巷模糊度的整数值代入所述双差宽巷载波相位观测方程；
根据所述双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解。
7.如权利要求6所述的GNSS定位方法，其特征在于，所述卫星导航信号包括GPS、北斗、伽利略卫星中至少两个导航系统的信号。
8.如权利要求6所述的GNSS定位方法，其特征在于，进行卡尔曼滤波时，设置所述卡尔曼滤波对应于双差宽巷载波相位观测方程的观测噪声协方差矩阵乘以第一组合系数，所述第一组合系数根据所述第一组合的两个频点的波长比获得。
9.如权利要求6所述的GNSS定位方法，其特征在于，还包括：利用lambda方法对所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解进行固定，获得所述双差载波相位观测方程中对应至少两个频点的模糊度的固定解。
10.一种GNSS定位方法，其特征在于，包括：
获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少两个频点上的伪距观测值及载波相位观测值；
建立对应所述至少两个频点的双差载波相位观测方程和双差伪距观测方程；
将所述卫星导航信号的至少两个频点中的两个频点作为第二组合，根据所述第二组合的卫星观测数据采用MW组合计算超宽巷模糊度的初始值，并建立对应所述第二组合的双差超宽巷载波相位观测方程；
将所述超宽巷模糊度的初始值进行多个历元的平滑并计算所述超宽巷模糊度的整数值；
将所述超宽巷模糊度的整数值代入所述双差超宽巷载波相位观测方程；
根据所述双差载波相位观测方程、双差伪距观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解。
11.如权利要求10所述的的GNSS定位方法，其特征在于，所述卫星导航信号包括GPS、北斗、伽利略卫星中至少两个导航系统的信号。
12.如权利要求10所述的GNSS定位方法，其特征在于，进行卡尔曼滤波时，设置所述卡尔曼滤波对应于双差超宽巷载波相位观测方程的观测噪声协方差矩阵乘以第二组合系数，所述第二组合系数根据所述第二组合的两个频点的波长比获得。
13.如权利要求10所述的GNSS定位方法，其特征在于，还包括：利用lambda方法对所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解进行固定，获得所述双差载波相位观测方程中对应至少两个频点的模糊度的固定解。
14.一种GNSS定位装置，其特征在于，包括：
获取模块，用于获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少三个频点上的伪距观测值及载波相位观测值；
宽巷模糊度计算模块，用于将所述卫星导航信号的至少三个频点中的两个频点作为第一组合，根据所述第一组合的卫星观测数据采用MW组合计算宽巷模糊度的初始值，并建立对应所述第一组合的双差宽巷载波相位观测方程；
超宽巷模糊度计算模块，用于将所述卫星导航信号的至少三个频点中的两个频点作为第二组合，所述第二组合和第一组合不同，根据所述第二组合的卫星观测数据采用MW组合计算超宽巷模糊度的初始值，并建立对应所述第二组合的双差超宽巷载波相位观测方程；
平滑取整模块，用于分别将所述宽巷模糊度和超宽巷模糊度的初始值进行多个历元的平滑并计算所述宽巷模糊度和超宽巷模糊度的整数值；
滤波模块，用于分别将所述宽巷模糊度和超宽巷模糊度的整数值代入所述双差宽巷载波相位观测方程和所述双差超宽巷载波相位观测方程，建立对应所述至少三个频点的双差载波相位观测方程和双差伪距观测方程，根据所述双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少三个频点的模糊度的浮点解。
15.一种GNSS定位装置，其特征在于，包括：
获取模块，用于获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少两个频点上的伪距观测值及载波相位观测值；
宽巷模糊度计算模块，用于将所述卫星导航信号的至少两个频点中的两个频点作为第一组合，根据所述第一组合的卫星观测数据采用MW组合计算宽巷模糊度的初始值，并建立对应所述第一组合的双差宽巷载波相位观测方程；
平滑取整模块，用于将所述宽巷模糊度的初始值进行多个历元的平滑并计算所述宽巷模糊度的整数值；
滤波模块，用于将所述宽巷模糊度的整数值代入所述双差宽巷载波相位观测方程，建立对应所述至少两个频点的双差载波相位观测方程和双差伪距观测方程，根据所述双差载波相位观测方程、双差伪距观测方程、代入宽巷模糊度整数值后的双差宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解。
16.一种GNSS定位装置，其特征在于，包括：
获取模块，用于获取基准站和流动站的卫星观测数据，所述卫星观测数据包括卫星导航信号的至少两个频点上的伪距观测值及载波相位观测值；
超宽巷模糊度计算模块，用于将所述卫星导航信号的至少两个频点中的两个频点作为第二组合，根据所述第二组合的卫星观测数据采用MW组合计算超宽巷模糊度的初始值，并建立对应所述第二组合的双差超宽巷载波相位观测方程；
平滑取整模块，用于将所述超宽巷模糊度的初始值进行多个历元的平滑并计算所述超宽巷模糊度的整数值；
滤波模块，用于将所述超宽巷模糊度的整数值代入所述双差超宽巷载波相位观测方程，建立对应所述至少两个频点的双差载波相位观测方程和双差伪距观测方程，根据所述双差载波相位观测方程、双差伪距观测方程和代入超宽巷模糊度整数值后的双差超宽巷载波相位观测方程进行卡尔曼滤波，得到所述双差载波相位观测方程中对应至少两个频点的模糊度的浮点解。
CN202010737499.XA 2020-07-28 2020-07-28 Gnss定位方法及其装置 Pending CN114002722A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010737499.XA CN114002722A (zh) 2020-07-28 2020-07-28 Gnss定位方法及其装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010737499.XA CN114002722A (zh) 2020-07-28 2020-07-28 Gnss定位方法及其装置

Publications (1)

Publication Number Publication Date
CN114002722A true CN114002722A (zh) 2022-02-01

Family

ID=79920381

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010737499.XA Pending CN114002722A (zh) 2020-07-28 2020-07-28 Gnss定位方法及其装置

Country Status (1)

Country Link
CN (1) CN114002722A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115079225A (zh) * 2022-07-07 2022-09-20 深圳市海伊石油技术有限公司 一种海上接收机的导航定位方法与装置
CN115220077A (zh) * 2022-07-18 2022-10-21 千寻位置网络有限公司 模糊度固定方法及其系统

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115079225A (zh) * 2022-07-07 2022-09-20 深圳市海伊石油技术有限公司 一种海上接收机的导航定位方法与装置
CN115220077A (zh) * 2022-07-18 2022-10-21 千寻位置网络有限公司 模糊度固定方法及其系统

Similar Documents

Publication Publication Date Title
CN109738917B (zh) 一种北斗变形监测中的多路径误差削弱方法及装置
CN105158783B (zh) 一种实时动态差分定位方法及其设备
CN106291639B (zh) 一种gnss接收机实现定位的方法及装置
CN1864078B (zh) 使用三个gps频率消除载波相位整数模糊度的方法
CN103675835B (zh) 一种北斗三频信号载波相位整周模糊度单历元确定方法
CN108802782B (zh) 一种惯导辅助的北斗三频载波相位整周模糊度求解方法
CN107728180B (zh) 一种基于多维粒子滤波偏差估计的gnss精密定位方法
CN107064980A (zh) 载波相位模糊度固定方法和装置、卫星导航接收机
CN113138402B (zh) 基于rtk的模糊度固定方法及装置、存储介质
CN104898145B (zh) 一种基于半周模糊度的模糊度固定方法和系统
AU2008260579A1 (en) Partial search carrier-phase integer ambiguity resolution
CN105676250A (zh) 一种基于gnss的单历元三频模糊度解算方法
AU2007300586A1 (en) Method for using three GPS frequencies to resolve whole-cycle carrier-phase ambiguities
CN105974440B (zh) 一种获取观测噪声的方法和装置
CN110749907A (zh) 一种基于北斗动定位中接收机的钟差补偿方法及其系统
CN116299623B (zh) 一种城市复杂场景下的ppp和ins紧组合方法与系统
JP4498399B2 (ja) 測位システム及び測位方法
CN114002722A (zh) Gnss定位方法及其装置
CN115480279A (zh) Gnss导航方法和终端、组合导航系统、存储介质
CN110568464A (zh) 基于bds/gnss多模芯片的精密定位方法及装置
CN111948682B (zh) 一种基于bds三频组合的伪距和载波相位随机模型构建方法
JP4723801B2 (ja) 相対測位装置
WO2023065840A1 (zh) 模糊度固定方法及其系统、存储介质
CN105699997B (zh) 一种使用glonass单频信号进行差分定位的方法
CN112987048B (zh) 一种适用于Andriod智能终端的高精度定位方法、智能终端以及存储介质

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