CN114966756A - 基频模糊度固定的方法、装置、设备及计算机存储介质 - Google Patents
基频模糊度固定的方法、装置、设备及计算机存储介质 Download PDFInfo
- Publication number
- CN114966756A CN114966756A CN202110196249.4A CN202110196249A CN114966756A CN 114966756 A CN114966756 A CN 114966756A CN 202110196249 A CN202110196249 A CN 202110196249A CN 114966756 A CN114966756 A CN 114966756A
- Authority
- CN
- China
- Prior art keywords
- ambiguity
- widelane
- frequency
- equation
- carrier
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000001914 filtration Methods 0.000 claims abstract description 48
- 239000005433 ionosphere Substances 0.000 claims abstract description 23
- 238000007667 floating Methods 0.000 claims description 33
- 241001061260 Emmelichthys struhsakeri Species 0.000 claims description 20
- 238000009499 grossing Methods 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 12
- 235000008331 Pinus X rigitaeda Nutrition 0.000 claims description 9
- 235000011613 Pinus brutia Nutrition 0.000 claims description 9
- 241000018646 Pinus brutia Species 0.000 claims description 9
- 238000013450 outlier detection Methods 0.000 claims description 7
- 238000010586 diagram Methods 0.000 description 14
- 239000011159 matrix material Substances 0.000 description 10
- 230000000694 effects Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000004891 communication Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 239000000969 carrier Substances 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 150000001875 compounds Chemical class 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000006870 function Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 230000005055 memory storage Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000005436 troposphere Substances 0.000 description 1
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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/20—Integrity monitoring, fault detection or fault isolation of space segment
-
- 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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/21—Interference related issues ; Issues related to cross-correlation, spoofing or other methods of denial of service
-
- 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/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing (AREA)
- Computer Security & Cryptography (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本公开实施例提供了一种基频模糊度固定的方法、装置、设备及计算机存储介质,方法基于观测数据的宽巷模糊度初值作为待估参数,通过宽巷观测方程与无电离层组合方程,对基频模糊度进行松约束滤波和紧约束滤波以及模糊度域搜索等方式,最终得到准确性较高的基频模糊度,适用于中长基线的RTK定位解算,同样也适用于短基线的解算,且不需要任何改变,适用范围广,解算效果好。
Description
技术领域
本公开属于卫星定位技术领域,尤其涉及一种基频模糊度固定的方法、 装置、设备及计算机存储介质。
背景技术
基于双差的GNSS相对定位,经过双差,短基线中存在的电离层误差、 对流层误差、卫星钟差、接收机钟差、卫星星历误差等,都可以被充分的 削弱或完全消除。但随着基线长度的增加,各项误差逐渐增大,空间相关 性减弱,尤其是电离层误差,即使经双差后,残余误差仍然较大;当基线 长度超过30km时,非组合双差观测值已很难正确固定模糊度,导致无法 达到高精度RTK(Real-Time Kinematic,实时动态)定位效果。
发明内容
本公开实施例提供一种基频模糊度固定的方法、装置、设备及计算机 存储介质,能够对模糊度正确固定,利于得到高精度的定位效果。
第一方面,本公开实施例提供一种基频模糊度固定的方法,方法包括:
根据参考站和流动站的载波观测数据,计算得到宽巷模糊度初值;
通过第一联合观测方程,对基频模糊度初值进行滤波,将滤波后得到 的第一基频模糊度的浮点解转换为宽巷模糊度滤波值;第一联合观测方程, 包括无电离层组合观测方程与宽巷模糊度松约束方程,其中,宽巷模糊度 松约束方程是以宽巷模糊度初值为约束;
对宽巷模糊度滤波值进行模糊度域搜索,得到固定宽巷模糊度;
通过第二联合观测方程,对第一基频模糊度进行再次滤波,得到第二 基频模糊度的浮点解;第二联合观测方程包括无电离层组合观测方程与宽 巷模糊度紧约束方程,其中,宽巷模糊度紧约束方程是以固定宽巷模糊度 为约束;
对第二基频模糊度的浮点解进行模糊度域搜索,得到载波观测数据的 基频模糊度固定解。
在一些实施例中,根据参考站和流动站的载波观测数据,计算得到宽 巷模糊度初值,具体包括:
根据参考站和流动站的双频载波观测数据,通过多历元平滑伪距,解 算得到双差MW宽巷模糊度作为宽巷模糊度初值。
在一些实施例中,根据参考站和流动站的载波观测数据,计算得到宽 巷模糊度初值,具体包括:
根据参考站间的三频载波观测数据,通过平滑伪距取整,获得两组超 宽巷模糊度;
将两组超宽巷模糊度进行组合,形成宽巷模糊度初值。
在一些实施例中,根据参考站和流动站的载波观测数据,计算得到宽 巷模糊度初值,具体包括:
当载波观测数据仅包括双频观测数据时,根据参考站和流动站的双频 载波观测数据,通过多历元平滑伪距,解算得到双差MW宽巷模糊度作为 宽巷模糊度初值;
当载波观测数据包括三频观测数据时,根据参考站间的三频载波观测 数据,通过平滑伪距取整,获得两组超宽巷模糊度,并将两组超宽巷模糊 度进行组合,形成宽巷模糊度初值。
在一些实施例中,在根据参考站和流动站的载波观测数据,计算得到 宽巷模糊度初值之后,且在通过第一联合观测方程,对基频模糊度初值进 行滤波,将滤波后得到的第一基频模糊度的浮点解转换为宽巷模糊度滤波 值之前,方法还包括:
将宽巷模糊度初值逐一提取,做最小二乘解算,得到对应的多个位置 数据;
根据位置数据进行离群检测,将位置离群的对应宽巷模糊度初值进行 删除,得到离群检测后的剩余宽巷模糊度初值;
宽巷模糊度松约束方程是以剩余宽巷模糊度初值为约束。
在一些实施例中,通过第一联合观测方程,对基频模糊度初值进行滤 波,将滤波后得到的第一基频模糊度的浮点解转换为宽巷模糊度滤波值, 具体包括:
由以下无电离层组合观测方程(1)与宽巷模糊度松约束方程(2)构 成的第一联合观测方程:
其中,φ1、φ2分别为载波L1频点及载波L2频点的观测值,为无电 离层组合后的载波相位观测值;f1与f2分别为载波L1频点及载波L2频点的 频率;分别为载波L1频点及载波L2频点的双差载波相位观测 值;i,j分别表示参考站,N1、N2分别表示载波L1频点与载波L2频点观测 值的基频模糊度初值,△表示站间单差;Nmw为宽巷模糊度初值。
在一些实施例中,通过第二联合观测方程,对第一基频模糊度进行再 次滤波,得到第二基频模糊度的浮点解值,包括;
由以下无电离层组合观测方程(1)与宽巷模糊度紧约束方程(3)构 成的第二联合观测方程:
第二方面,本公开实施例提供一种基频模糊度固定的装置,装置包括:
第一计算模块:用于根据参考站和流动站的载波观测数据,计算得到 宽巷模糊度初值;
第一滤波模块:用于通过第一联合观测方程,对基频模糊度初值进行 滤波,将滤波后得到的第一基频模糊度的浮点解转换为宽巷模糊度滤波值; 第一联合观测方程,包括无电离层组合观测方程与宽巷模糊度松约束方程, 其中,宽巷模糊度松约束方程是以宽巷模糊度初值为约束;
第一搜索模块:用于对宽巷模糊度滤波值进行模糊度域搜索,得到固 定宽巷模糊度;
第二滤波模块:用于通过第二联合观测方程,对第一基频模糊度进行 再次滤波,得到第二基频模糊度的浮点解;第二联合观测方程包括无电离 层组合观测方程与宽巷模糊度紧约束方程,其中,宽巷模糊度紧约束方程 是以固定宽巷模糊度为约束;
第二搜索模块,用于对第二基频模糊度的浮点解进行模糊度域搜索, 得到载波观测数据的基频模糊度固定解。
第三方面,本公开实施例提供一种基频模糊度固定的设备,设备包括: 处理器,以及存储有计算机程序指令的存储器;处理器读取并执行计算机 程序指令,以实现如上述任意一实施例的基频模糊度固定的方法。
第四方面,本公开实施例提供一种计算机存储介质,计算机存储介质 上存储有计算机程序指令,计算机程序指令被处理器执行时实现如上述任 意一实施例的基频模糊度固定的方法。
本公开实施例的基频模糊度固定的定位方法、装置、设备及计算机存 储介质,基于观测数据的宽巷模糊度初值作为待估参数,通过宽巷观测方 程与无电离层组合方程,对基频模糊度进行松约束滤波和紧约束滤波以及 模糊度域搜索等方式,最终得到准确性较高的基频模糊度,适用于中长基 线的RTK定位解算,同样也适用于短基线的解算,且不需要任何改变,适 用范围广,解算效果好。
附图说明
为了更清楚地说明本公开实施例的技术方案,下面将对本公开实施例 中所需要使用的附图作简单的介绍,对于本领域普通技术人员来讲,在不 付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本公开实施例提供的一种基频模糊度固定的方法的流程示意图;
图2是本公开一个具体实施例提供基频模糊度固定的方法的流程示意 图;
图3是图1中所示步骤S101的一个示例的流程示意图;
图4是图2所示基频模糊度固定的方法中,一个具体示例的流程示意 图;
图5是本公开实施例提供的一种基频模糊度固定的装置的结构示意图;
图6是图5所示基频模糊度固定的装置中,一个具体示例的结构示意 图;
图7是本公开实施例提供的一种基频模糊度固定的设备的结构示意图。
具体实施方式
下面将详细描述本公开的各个方面的特征和示例性实施例,为了使本 公开的目的、技术方案及优点更加清楚明白,以下结合附图及具体实施例, 对本公开进行进一步详细描述。应理解,此处所描述的具体实施例仅意在 解释本公开,而不是限定本公开。对于本领域技术人员来说,本公开可以 在不需要这些具体细节中的一些细节的情况下实施。下面对实施例的描述 仅仅是为了通过示出本公开的示例来提供对本公开更好的理解。
需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用 来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者 暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语 “包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而 使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素,而 且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物 品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括……” 限定的要素,并不排除在包括所述要素的过程、方法、物品或者设备中还 存在另外的相同要素。
针对长基线中RTK定位时模糊度难以正确固定的问题,传统技术中常 采用的方法主要包括三类:(1)是非组合观测值与组合形成的无电离层 观测值共同滤波,将宽巷模糊度及载波L1模糊度作为待估参数进行估计, 且非组合观测值中,将电离层参数作为待估参数。(2)通过非组合模式, 将电离层误差作为待估参数进行估计,此时载波L1或L2模糊度仍作为待 估参数直接进行估计。(3)三频或多频模糊度固定技术固定长基线模糊 度,通过分步逐级固定超宽巷、宽巷模糊度,最后再固定L1或L2模糊度。 但是方法(1)和(2),由于每颗卫星均设置一个电离层参数作为待估参 数,加上载波L1或L2模糊度参数,导致观测方程中待估参数极多,计算 速度较慢。此外,该方法需要给出电离层延迟误差先验值及其过程噪声, 该值不易准确给出,但该值又非常关键,其准确性直接影响L1或L2模糊 度固定的正确性。方法(3)需要三频观测值,而GPS中具有三频观测值 的卫星目前较少,且当接收机仅支持双频观测值时,该方法则无法使用。
为了解决现有技术问题,本公开实施例提供了一种基频模糊度固定的 方法、装置、设备及计算机存储介质。
下面首先对本公开实施例所提供的基频模糊度固定的方法进行介绍。
图1示出了本公开一个实施例提供的基频模糊度固定的方法的流程示 意图。如图1所示,该方法可以包括以下步骤:
S101.根据参考站和流动站的载波观测数据,计算得到宽巷模糊度初 值;
S102.通过第一联合观测方程,对基频模糊度初值进行滤波,将滤波后 得到的第一基频模糊度的浮点解转换为宽巷模糊度滤波值;第一联合观测 方程,包括无电离层组合观测方程与宽巷模糊度松约束方程,其中,宽巷 模糊度松约束方程是以宽巷模糊度初值为约束;
S103.对宽巷模糊度滤波值进行模糊度域搜索,得到固定宽巷模糊度;
S104.通过第二联合观测方程,对第一基频模糊度进行再次滤波,得到 第二基频模糊度的浮点解;第二联合观测方程包括无电离层组合观测方程 与宽巷模糊度紧约束方程,其中,宽巷模糊度紧约束方程是以固定宽巷模 糊度为约束;
S105.对第二基频模糊度的浮点解进行模糊度域搜索,得到载波观测数 据的基频模糊度固定解。
本实施例中,基于原始基频模糊度作为待估参数,利用宽巷模糊初 值,联合无电离层组合观测方程加松约束滤波,通过模糊度域搜索得到固 定宽巷模糊度,加紧约束滤波,求取出更为准确的浮点解、基频模糊度浮 点值,不但可以适用于中长基线模糊度的计算,也可以应用于短基线模糊 度的计算。
图2示出了本公开一个具体实施例中基频模糊度固定的方法的流程示 意图。本实施例中,进行步骤S101.时,可以先根据观测数据的性质通过 不同方式求取宽巷模糊度初值。具体的,步骤S101.可以包括:
当载波观测数据仅包括双频观测数据时,根据参考站和流动站的双频 载波观测数据,通过多历元平滑伪距,解算得到双差MW宽巷模糊度作为 宽巷模糊度初值;
当载波观测数据包括三频观测数据时,根据参考站间的三频载波观测 数据,通过平滑伪距取整,获得两组超宽巷模糊度,并将两组超宽巷模糊 度进行组合,形成宽巷模糊度初值。
示例性的,对于双频载波观测数据,步骤S101.据参考站和流动站的 载波观测数据,计算得到宽巷模糊度初值,具体可以包括:
根据参考站和流动站的双频载波观测数据,通过多历元平滑伪距,解 算得到双差MW宽巷模糊度作为宽巷模糊度初值。
式(1.1)与(1.2)中,Nmw为当前历元t计算的双差宽巷MW模糊 度,也即宽巷模糊度初值;表示双差算子,即观测值或组合观测值在站 间(即参考站与流动站间)和星间双差,φ1、φ2分别为载波L1频点及载波 L2频点的观测值,下同;则分别为载波L1频点及载波L2频点 的双差载波相位观测值;f1与f2分别为载波L1频点及载波L2频点的频率;λ1,λ2分别为载波L1频点及载波L2频点的波长,P1,P2分别为载波 L1频点及载波L2频点的对应的伪距观测值;为当前历元t经平滑 后的MW宽巷模糊度,n为大于1的正整数,表示平滑系数,n的取值可 根据经验设置,如300。
对于双频载波观测数据的宽巷模糊度初值求取方式在其他示例中,还 可以采用待估参数法,即将宽巷模糊度初值作为待估参数进行估计;也可 以采用直接取整法,即通过单历元宽巷伪距观测值,结合宽巷波长,通过 直接取证求得宽巷模糊度初值。
图3示出的步骤S101.的一个具体示例的流程示意图。如图3所示, 本实施例中,示例性的,对于三频载波观测数据,步骤S101.据参考站和 流动站的载波观测数据,计算得到宽巷模糊度初值,具体可以包括:
S301.根据参考站间的三频载波观测数据,通过平滑伪距取整,获得两 组超宽巷模糊度;
S302.将两组超宽巷模糊度进行组合,形成宽巷模糊度初值。
本示例中,当从观测数据中判断存在三频载波观测数据(也即三频载 波观测值,下同)时,进行步骤S301.根据三频载波观测值形成两组超宽 巷观测值,并直接通过平滑伪距取整得到对应的超宽巷模糊度。三频载波 观测值可以如下式(1.3)表示:
其中,i,j,k均为整数系数,φ1,φ2,φ3分别为三个频点所对应的以米为单 位的载波观测值,f1,f2,f3分别为对应三个频点的频率。
三频载波观测值(也即组合观测值,下同)的模糊度可表示为:
N(i,j,k)=i·N1+j·N2+k·N3 (1.4)
其中N1,N2,N3分别为三个频点所对应的原始观测值模糊度(也即基频 模糊度初值)。
组合观测值频率可表示为:
f(i,j,k)=i·f1+j·f2+k·f3 (1.5)
组合观测值波长可表示为:
λ(i,j,k)=c/f(i,j,k) (1.6)
其中c为光速。
对于不同的系统,整数系数i,j,k应当分别取不同的值,该值可以是经 验值,取适合的整数系数以得到波长足够长,且最有利于模糊度解算的组 合观测值。
本示例中,通过平滑伪距取整,以米为单位的三频伪距组合观测值可 相应的表示为:
式(1.7)中,P1,P2,P3分别为三频载波观测数据中三个频点对应的伪距 观测值,||表示取绝对值。从而,超宽巷双差模糊度可表示为:
式(1.8)中,表示双差算子,即观测值或组合观测值在站间和星 间的双差,h表示组合系数,[]表示取整,故而,得到的超宽巷模糊度 Nt ewl(i,j,k)已是整数,且由于组合观测值波长较长,因而该值基本可保证 100%正确。
通过平滑伪距取整,获得两组超宽巷模糊度,将两组超宽巷模糊度进 行组合,形成一宽巷模糊度,即作为上述的宽巷模糊度初值,该值的准确 性会高于平滑求取的MW宽巷模糊度的准确性,可以将该值代替MW宽 巷模糊度。因为此时仅部分卫星存在三频观测值,所以仅有少部分MW模 糊度会被取代,得到最终的宽巷模糊度初值。
示例性的,宽巷模糊度初值的计算公式可以表示为:
Nwl(i,j,k)=m·N1 ewl(i,j,k)+n·N2 ewl(i,j,k) (1.9)
其中N1 ewl(i,j,k),N2 ewl(i,j,k)分别为超宽巷两种组合系数求得的两组超宽巷 模糊度。例如对于GPS,两组超宽巷模糊度(i,j,k)的选取可分别为(0,-1,1)以 及(1,-6,5),m、n则可选取为(-5,1)。
本示例中,当选取的两组超宽巷波长足够长时,通过直接取整就足以 忽略取整误差,模糊度错误的可能性较小,则两组超宽巷模糊度通过线性 组合,转换变为宽巷模糊度后,该宽巷模糊度具有较高的准确性。
在一个具体实施例中,在根据参考站和流动站的载波观测数据,计算 得到宽巷模糊度初值之后,且在通过第一联合观测方程,对基频模糊度初 值进行滤波,将滤波后得到的第一基频模糊度的浮点解转换为宽巷模糊度 滤波值之前,基频模糊度固定的方法还可以包括:S106.剔除错误的宽巷模 糊度初值。
示例性的,图4示出了步骤S106.在一个示例中的流程示意图。如图4 所示,步骤S106.具体可以包括:
S401.将宽巷模糊度初值逐一提取,做最小二乘解算,得到对应的多个 位置数据;
S402,根据位置数据进行离群检测,将位置离群的对应宽巷模糊度初 值进行删除,得到离群检测后的剩余宽巷模糊度初值。
本示例中,假设有m个宽巷模糊度初值,则逐个剔除一个宽巷模糊度 初值,做最小二乘解算,最终得到m个位置数据结果。将得到的m个位 置数据结果做离群点检测,将明显离群的点所对应的宽巷模糊度初值判定 为错误的值,将这些错误的值进行剔除,得到剩余宽巷模糊度初值,利用 剩余宽巷模糊度初值,用于后续的松约束处理。
本实施例的步骤S106.在其他示例中,也可以通过固定的宽巷模糊度 初值,做最小二乘解算,将先验及后验残差明显离群的宽巷模糊度初值进 行剔除。
本实施例中,得到准确的宽巷模糊度初值后,可以进行步骤S102.通 过第一联合观测方程,对基频模糊度初值进行滤波,将滤波后得到的第一 基频模糊度的浮点解转换为宽巷模糊度滤波值。该步骤中,第一联合观测 方程,包括无电离层组合观测方程与宽巷模糊度松约束方程,其中,宽巷 模糊度松约束方程是以宽巷模糊度初值为约束;如果在步骤S103之前进 行了步骤S106,则宽巷模糊度松约束方程是以上述的剩余宽巷模糊度初值 为约束。
由于电离层误差是影响中长基线解算的最关键因素,且接收机钟差、 卫星钟差、电离层误差等经站间-星间双差后,基本得到大幅度削弱或完全 消除,因此,此处所列方程中仅考虑电离层误差影响。
电离层误差具有色散效应,利用该特性,双频观测值可组合形成电离 层无关的新的观测值,称为无电离层组合观测值。无电离层组合观测值消 除了电离层的影响,但该组合观测值模糊度不再具有整数特性。因此,在 使用无电离层组合观测值,进行模糊度固定时,需要首先将无电离层观测 值模糊度转换为宽巷模糊度与基频模糊度的组合,当宽巷模糊度固定后, 无电离层观测模糊度的固定可以转换为基频模糊度的固定。
本实施例中,通过步骤S101得到基本准确的宽巷模糊度初值后,通 过无电离层组合观测方程与宽巷模糊度松约束方程进行滤波,待估参数仅 设置位置以及载波L1频点与载波L2频点观测值的基频模糊度,该待估参 数与与短基线解算中待估参数一致。
仅考虑电离层误差,以米为单位的载波L1及载波L2的原始观测方程 可以表示为:
式(2.1)中,φ1、φ2分别为载波L1频点及载波L2频点的观测值,f1与f2分别为载波L1频点及载波L2频点的频率;λ1,λ2分别为载波L1频 点及载波L2频点的波长,N1、N2分别表示载波L1频点与载波L2频点观 测值的基频模糊度初值;c为光速;A表示不同信号同一传播路径上的相 同电离层延迟项。
基于式(2.1),消除电离层后的双差无电离层组合观测方程可以表示 为:
当仅有无电离层组合观测值时,在解算基频模糊度初值N1、N2时秩 亏,因此需要添加宽巷模糊度对基频模糊度初值N1、N2的约束方程。
示例性的,添加的宽巷模糊度松约束方程可以表示为:
式(2.3)中,Nmw为宽巷模糊度初值;i,j分别表示进行星间差分的两 个卫星,△表示站间单差,即为载波L1观测值对i卫星的站间单差模 糊度,即为载波L2观测值对i卫星的站间单差模糊度,即为载波 L1观测值对j卫星的站间单差模糊度,即为载波L2观测值对j卫星的 站间单差模糊度。
因此,由以下无电离层组合观测方程(2.2)与宽巷模糊度松约束方程 (2.3)构成第一联合观测方程:
通过无电离层观测方程及宽巷模糊度松约束方程构成的第一联合观测 方程,对基频模糊度初值进行松约束滤波时,待估参数为位置(x,y,z)以及 各卫星对载波各频点的基频站间单差浮点模糊度,也即△N1,△N2。
基频模糊度初值进行松约束后,得到的是第一基频模糊度的浮点解; 将该第一基频模糊度的浮点解转换为宽巷模糊度滤波值,具体可以包括:
将载波L1及L2的基频站间单差模糊度△N1,△N2的浮点解,转换为L1- L2的双差宽巷模糊度Nwl。其转换方法为:
式(2.4)中,D矩阵表示将站间单差模糊度转换为双差模糊度矩阵, A矩阵表示将双差模糊度转换为宽巷模糊度矩阵。具体可表示为:
然后将载波L1及L2单差模糊度的方差协方差矩阵P△N,转换为L1-L2 的双差宽巷模糊度的方差协方差矩阵。其转换公式为:
Pwl=(A·D)·P△N·(A·D)T (2.7)
式(2.7)中()T表示对矩阵求取转置,“·”表示矩阵相乘。
式(2.7)得到的双差宽巷模糊度的方差协方差矩阵中,包含的是第一 基频模糊度的浮点解转换的宽巷模糊度滤波值。
本实施例中,通过步骤S103.对得到的对宽巷模糊度滤波值进行模糊 度域搜索,得到固定宽巷模糊度。
示例性的,可以通过lambda搜索,固定宽巷模糊度Nwl。由于经 lambda搜索的宽巷模糊度准确性远高于伪距平滑的宽巷模糊度,因此利用 搜索得到的固定宽巷模糊度Nwl替代伪距平滑得到的MW模糊度Nmw。
得到准确性高的固定宽巷模糊度Nwl后,可以进行步骤S104.通过第二 联合观测方程,对第一基频模糊度再次进行滤波,得到第二基频模糊度的 浮点解。其中第二联合观测方程包括无电离层组合观测方程与宽巷模糊度 紧约束方程,其中,宽巷模糊度紧约束方程是以固定宽巷模糊度为约束; 且宽巷模糊度紧约束方程的计算精度高于宽巷模糊度松约束方程的计算精 度。
在步骤S104中,可以由以下无电离层组合观测方程(3.1)与宽巷模 糊度紧约束方程(3.2)构成的第二联合观测方程:
第二联合方程对基频模糊度的约束处理方式与第一联合方程相同,只 是约束值不同,由于固定宽巷模糊度Nwl的准确性高于宽巷模糊度初值 Nmw,固定宽巷模糊度Nwl替代伪MW模糊度Nmw后,加入紧约束滤波时, 宽巷模糊度紧约束方程的计算精度要高于所述宽巷模糊度松约束方程的计 算精度。换言之,两约束均指给定的宽巷约束方程精度,所不同的是给定 精度的相对大小;紧约束的精度应高于松约束精度,如松约束时宽巷约束 方程给定的精度范围为0.01~0.001m,则紧约束时给定的精度范围则可以 为0.001~0.0005。
本实施例中,通过步骤S105.对所述第二基频模糊度的浮点解进行模 糊度域搜索,得到所述载波观测数据的基频模糊度固定解。
示例性的,在步骤S105.中,可以通过lambda搜索,固定载波L1及 L2频点的基频模糊度。经过步骤S104之后,已获得较为准确的浮点解坐 标,基频模糊度N1、N2,以及其对应的方差协方差矩阵P△N,此时再通过 labmda搜索,直接固定基频模糊度N1、N2,固定率较高,且固定错误的概 率较小。当模糊度固定后,即可将浮点解转换为固定解。
本实施例提供的基频模糊度固定方法,基于原始的单差模糊度转换为 宽巷模糊度方法,简单明了,避免了宽巷模糊度作为待估参数带来复杂变 换,可应用于中长基线的RTK定位解算,由于待估参数一致,所以同样 也适用于短基线的解算,且不需要做任何改变,效果较好。
图5示出的本公开实施例提供的一种基频模糊度固定的装置结构示意 图、如图5所示的装置包括:
第一计算模块501:用于根据参考站和流动站的载波观测数据,计算 得到宽巷模糊度初值;
第一滤波模块502:用于通过第一联合观测方程,对基频模糊度初值 进行滤波,将滤波后得到的第一基频模糊度的浮点解转换为宽巷模糊度滤 波值;第一联合观测方程,包括无电离层组合观测方程与宽巷模糊度松约 束方程,其中,宽巷模糊度松约束方程是以宽巷模糊度初值为约束;
第一搜索模块503:用于对宽巷模糊度滤波值进行模糊度域搜索,得 到固定宽巷模糊度;
第二滤波模块504:用于通过第二联合观测方程,对第一基频模糊度 进行再次滤波,得到第二基频模糊度的浮点解;第二联合观测方程包括无 电离层组合观测方程与宽巷模糊度紧约束方程,其中,宽巷模糊度紧约束 方程是以固定宽巷模糊度为约束;
第二搜索模块505,用于对第二基频模糊度的浮点解进行模糊度域搜 索,得到载波观测数据的基频模糊度固定解。
示例性的,第一计算模块501可以执行上述图1中示出的步骤S101, 第一滤波模块502可以执行上述图1中示出的步骤S102,第一搜索模块 503可以执行上述图1中示出的步骤S103,第二滤波模块504可以执行上 述图1中示出的步骤S104,第二搜索模块505可以执行上述图1中示出的 步骤S105。
需要说明的是,上述方法实施例涉及的各步骤的所有相关内容均可以 援引到对应功能模块的功能描述,并能达到其相应的技术效果,为简洁描 述,在此不再赘述。
在一个实施例中,装置还可以包括:
错误剔除模块601,用于剔除错误的宽巷模糊度初值。即错误剔除模 块601可以执行上述图3中示出的步骤S401和S402。
需要说明的是,上述方法实施例涉及的各步骤的所有相关内容均可以 援引到对应功能模块的功能描述,并能达到其相应的技术效果,为简洁描 述,在此不再赘述。
图7示出了本公开实施例提供的一种基频模糊度固定的设备的硬件结 构图。
基频模糊度固定的设备可以包括处理器701以及存储有计算机程序指 令的存储器702。
具体地,上述处理器701可以包括中央处理器(Central Processing Unit,CPU),或者特定集成电路(Application Specific Integrated Circuit,ASIC),或者可以被配置成实施本公开实施例的一个或多个集成 电路。
存储器702可以包括用于数据或指令的大容量存储器。举例来说而非 限制,存储器702可包括硬盘驱动器(Hard Disk Drive,HDD)、软盘驱 动器、闪存、光盘、磁光盘、磁带或通用串行总线(Universal Serial Bus, USB)驱动器或者两个或更多个以上这些的组合。在一个实例中,存储器 702可以包括可移除或不可移除(或固定)的介质,或者存储器702是非易失 性固态存储器。存储器702可在综合网关容灾设备的内部或外部。
存储器702可以包括只读存储器(ROM),随机存取存储器 (RAM),磁盘存储介质设备,光存储介质设备,闪存设备,电气、光学 或其他物理/有形的存储器存储设备。因此,通常,存储器包括一个或多个 编码有包括计算机可执行指令的软件的有形(非暂态)计算机可读存储介 质(例如,存储器设备),并且当该软件被执行(例如,由一个或多个处 理器)时,其可操作来执行参考根据本公开的一方面的方法所描述的操 作。
处理器701通过读取并执行存储器702中存储的计算机程序指令,以 实现图1所示实施例中的方法/步骤S101至S105,并达到图1所示实例执 行其方法/步骤达到的相应技术效果,为简洁描述在此不再赘述。
在一个示例中,基频模糊度固定的设备还可包括通信接口703和总线 710。其中,如图7所示,处理器701、存储器702、通信接口703通过总 线710连接并完成相互间的通信。
通信接口703,主要用于实现本公开实施例中各模块、装置、单元和/ 或设备之间的通信。
总线710包括硬件、软件或两者,将基频模糊度固定设备的部件彼此 耦接在一起。举例来说而非限制,总线可包括加速图形端口(Accelerated Graphics Port,AGP)或其他图形总线、增强工业标准架构(Extended Industry Standard Architecture,EISA)总线、前端总线(Front Side Bus, FSB)、超传输(Hyper Transport,HT)互连、工业标准架构(Industry Standard Architecture,ISA)总线、无限带宽互连、低引脚数(LPC)总 线、存储器总线、微信道架构(MCA)总线、外围组件互连(PCI)总 线、PCI-Express(PCI-X)总线、串行高级技术附件(SATA)总线、视频 电子标准协会局部(VLB)总线或其他合适的总线或者两个或更多个以上 这些的组合。在合适的情况下,总线710可包括一个或多个总线。尽管本 公开实施例描述和示出了特定的总线,但本公开考虑任何合适的总线或互 连。
另外,结合上述实施例中的基频模糊度固定的方法,本公开实施例可 提供一种计算机存储介质来实现。该计算机存储介质上存储有计算机程序 指令;该计算机程序指令被处理器执行时实现上述实施例中的任意一种基 频模糊度固定的方法。
还需要说明的是,上面参考根据本公开的实施例的方法、装置(系 统)和计算机程序产品的流程图和/或框图描述了本公开的各方面。应当理 解,流程图和/或框图中的每个方框以及流程图和/或框图中各方框的组合 可以由计算机程序指令实现。这些计算机程序指令可被提供给通用计算 机、专用计算机、或其它可编程数据处理装置的处理器,以产生一种机 器,使得经由计算机或其它可编程数据处理装置的处理器执行的这些指令 使能对流程图和/或框图的一个或多个方框中指定的功能/动作的实现。这 种处理器可以是但不限于是通用处理器、专用处理器、特殊应用处理器或 者现场可编程逻辑电路。还可理解,框图和/或流程图中的每个方框以及框 图和/或流程图中的方框的组合,也可以由执行指定的功能或动作的专用硬 件来实现,或可由专用硬件和计算机指令的组合来实现。
以上所述,仅为本公开的具体实施方式,所属领域的技术人员可以清 楚地了解到,为了描述的方便和简洁,上述描述的系统、模块和单元的具 体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。应 理解,本公开的保护范围并不局限于此,任何熟悉本技术领域的技术人员 在本公开揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修 改或替换都应涵盖在本公开的保护范围之内。
Claims (10)
1.一种基频模糊度固定的方法,其特征在于,所述方法包括:
根据参考站和流动站的载波观测数据,计算得到宽巷模糊度初值;
通过第一联合观测方程,对基频模糊度初值进行滤波,将所述滤波后得到的第一基频模糊度的浮点解转换为宽巷模糊度滤波值;所述第一联合观测方程,包括无电离层组合观测方程与宽巷模糊度松约束方程,其中,宽巷模糊度松约束方程是以所述宽巷模糊度初值为约束;
对所述宽巷模糊度滤波值进行模糊度域搜索,得到固定宽巷模糊度;
通过第二联合观测方程,对所述第一基频模糊度进行再次滤波,得到第二基频模糊度的浮点解;所述第二联合观测方程包括所述无电离层组合观测方程与宽巷模糊度紧约束方程,其中,所述宽巷模糊度紧约束方程是以所述固定宽巷模糊度为约束;
对所述第二基频模糊度的浮点解进行模糊度域搜索,得到所述载波观测数据的基频模糊度固定解。
2.根据权利要求1所述的基频模糊度固定的方法,其特征在于,所述根据参考站和流动站的载波观测数据,计算得到宽巷模糊度初值,具体包括:
根据所述参考站和流动站的双频载波观测数据,通过多历元平滑伪距,解算得到双差MW宽巷模糊度作为所述宽巷模糊度初值。
3.根据权利要求1所述的基频模糊度固定的方法,其特征在于,所述根据参考站和流动站的载波观测数据,计算得到宽巷模糊度初值,具体包括:
根据所述参考站间的三频载波观测数据,通过平滑伪距取整,获得两组超宽巷模糊度;
将两组所述超宽巷模糊度进行组合,形成所述宽巷模糊度初值。
4.根据权利要求1所述的基频模糊度固定的方法,其特征在于,所述根据参考站和流动站的载波观测数据,计算得到宽巷模糊度初值,具体包括:
当所述载波观测数据仅包括双频观测数据时,根据所述参考站和流动站的双频载波观测数据,通过多历元平滑伪距,解算得到双差MW宽巷模糊度作为所述宽巷模糊度初值;
当所述载波观测数据包括三频观测数据时,根据所述参考站间的三频载波观测数据,通过平滑伪距取整,获得两组超宽巷模糊度,并将两组所述超宽巷模糊度进行组合,形成所述宽巷模糊度初值。
5.根据权利要求1所述的基频模糊度固定的方法,其特征在于,在所述根据参考站和流动站的载波观测数据,计算得到宽巷模糊度初值之后,且在所述通过第一联合观测方程,对基频模糊度初值进行滤波,将所述滤波后得到的第一基频模糊度的浮点解转换为宽巷模糊度滤波值之前,所述方法还包括:
将所述宽巷模糊度初值逐一提取,做最小二乘解算,得到对应的多个位置数据;
根据所述位置数据进行离群检测,将位置离群的对应宽巷模糊度初值进行删除,得到离群检测后的剩余宽巷模糊度初值;
宽巷模糊度松约束方程是以所述剩余宽巷模糊度初值为约束。
6.根据权利要求1所述的基频模糊度固定的方法,其特征在于,所述通过第一联合观测方程,对基频模糊度初值进行滤波,将所述滤波后得到的第一基频模糊度的浮点解转换为宽巷模糊度滤波值,具体包括:
由以下无电离层组合观测方程(1)与宽巷模糊度松约束方程(2)构成的第一联合观测方程:
8.一种基频模糊度固定的装置,其特征在于,所述装置包括:
第一计算模块:用于根据参考站和流动站的载波观测数据,计算得到宽巷模糊度初值;
第一滤波模块:用于通过第一联合观测方程,对基频模糊度初值进行滤波,将所述滤波后得到的第一基频模糊度的浮点解转换为宽巷模糊度滤波值;所述第一联合观测方程,包括无电离层组合观测方程与宽巷模糊度松约束方程,其中,宽巷模糊度松约束方程是以所述宽巷模糊度初值为约束;
第一搜索模块:用于对所述宽巷模糊度滤波值进行模糊度域搜索,得到固定宽巷模糊度;
第二滤波模块:用于通过第二联合观测方程,对所述第一基频模糊度进行再次滤波,得到第二基频模糊度的浮点解;所述第二联合观测方程包括所述无电离层组合观测方程与宽巷模糊度紧约束方程,其中,所述宽巷模糊度紧约束方程是以所述固定宽巷模糊度为约束;
第二搜索模块,用于对所述第二基频模糊度的浮点解进行模糊度域搜索,得到所述载波观测数据的基频模糊度固定解。
9.一种基频模糊度固定的设备,其特征在于,所述设备包括:处理器,以及存储有计算机程序指令的存储器;所述处理器读取并执行所述计算机程序指令,以实现如权利要求1-7任意一项所述的基频模糊度固定的方法。
10.一种计算机存储介质,其特征在于,所述计算机存储介质上存储有计算机程序指令,所述计算机程序指令被处理器执行时实现如权利要求1-7任意一项所述的基频模糊度固定的方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110196249.4A CN114966756A (zh) | 2021-02-22 | 2021-02-22 | 基频模糊度固定的方法、装置、设备及计算机存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110196249.4A CN114966756A (zh) | 2021-02-22 | 2021-02-22 | 基频模糊度固定的方法、装置、设备及计算机存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114966756A true CN114966756A (zh) | 2022-08-30 |
Family
ID=82953902
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110196249.4A Pending CN114966756A (zh) | 2021-02-22 | 2021-02-22 | 基频模糊度固定的方法、装置、设备及计算机存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114966756A (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105629279A (zh) * | 2015-12-18 | 2016-06-01 | 广州中海达卫星导航技术股份有限公司 | 一种网络基准站间的宽巷模糊度固定方法 |
US20180252819A1 (en) * | 2015-02-26 | 2018-09-06 | Profound Positioning Inc. | Method and system for performing precise point positioning (ppp) ambiguity resolution using gnss triple frequency signals |
CN108508468A (zh) * | 2018-03-20 | 2018-09-07 | 安徽理工大学 | 中长基线单历元模糊度解算方法、系统、设备、存储介质 |
CN109001781A (zh) * | 2018-08-01 | 2018-12-14 | 太原理工大学 | 一种顾及电离层约束的bds三频模糊度解算方法 |
WO2019144528A1 (zh) * | 2018-01-29 | 2019-08-01 | 东南大学 | 基于模糊度紧约束的多星座基准站间模糊度快速解算方法及其应用 |
CN111381264A (zh) * | 2018-12-27 | 2020-07-07 | 北京六分科技有限公司 | 网络rtk中长基线模糊度固定方法和平台 |
CN111751853A (zh) * | 2020-06-20 | 2020-10-09 | 北京华龙通科技有限公司 | 一种gnss双频载波相位整周模糊度解算方法 |
-
2021
- 2021-02-22 CN CN202110196249.4A patent/CN114966756A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20180252819A1 (en) * | 2015-02-26 | 2018-09-06 | Profound Positioning Inc. | Method and system for performing precise point positioning (ppp) ambiguity resolution using gnss triple frequency signals |
CN105629279A (zh) * | 2015-12-18 | 2016-06-01 | 广州中海达卫星导航技术股份有限公司 | 一种网络基准站间的宽巷模糊度固定方法 |
WO2019144528A1 (zh) * | 2018-01-29 | 2019-08-01 | 东南大学 | 基于模糊度紧约束的多星座基准站间模糊度快速解算方法及其应用 |
CN108508468A (zh) * | 2018-03-20 | 2018-09-07 | 安徽理工大学 | 中长基线单历元模糊度解算方法、系统、设备、存储介质 |
CN109001781A (zh) * | 2018-08-01 | 2018-12-14 | 太原理工大学 | 一种顾及电离层约束的bds三频模糊度解算方法 |
CN111381264A (zh) * | 2018-12-27 | 2020-07-07 | 北京六分科技有限公司 | 网络rtk中长基线模糊度固定方法和平台 |
CN111751853A (zh) * | 2020-06-20 | 2020-10-09 | 北京华龙通科技有限公司 | 一种gnss双频载波相位整周模糊度解算方法 |
Non-Patent Citations (1)
Title |
---|
尚睿;高成发;潘树国;陆轶材;张瑞成;: "基于北斗三频的BDS/GPS宽巷模糊度逐级单历元固定方法", 测绘工程, no. 04, 26 March 2018 (2018-03-26) * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103675835B (zh) | 一种北斗三频信号载波相位整周模糊度单历元确定方法 | |
CN105676250B (zh) | 一种基于gnss的单历元三频模糊度解算方法 | |
CN104502935B (zh) | 一种基于非差非组合模型的网络rtk模糊度解算方法 | |
CN107193029B (zh) | 北斗三频信号的网络rtk基准站间模糊度快速确定方法 | |
CN105158783B (zh) | 一种实时动态差分定位方法及其设备 | |
CN105223598B (zh) | 一种gnss载波相位整周模糊度单历元解算方法 | |
CN103576175B (zh) | 一种双频多星座gnss整周模糊度otf解算方法 | |
CN108802782B (zh) | 一种惯导辅助的北斗三频载波相位整周模糊度求解方法 | |
CN109581455B (zh) | 一种bds和gps融合的三频宽巷紧组合定位方法 | |
CN104898145B (zh) | 一种基于半周模糊度的模糊度固定方法和系统 | |
CN110531389B (zh) | Mw组合法与无几何相位法联合探测与修复周跳的方法 | |
CN109001781A (zh) | 一种顾及电离层约束的bds三频模糊度解算方法 | |
CN110346823B (zh) | 可用于北斗精密单点定位的三频模糊度解算方法 | |
CN103576167B (zh) | 基于hht和支持向量机的周跳探测与修复方法 | |
CN105974440B (zh) | 一种获取观测噪声的方法和装置 | |
CN113138402A (zh) | 基于rtk的模糊度固定方法及装置、存储介质 | |
CN116953741A (zh) | 一种应用于全球导航卫星系统gnss的周跳探测与修复方法 | |
CN105929430A (zh) | 一种gnss零基线参考站间模糊度快速固定方法 | |
CN105974441A (zh) | 一种接收机观测噪声的获取方法和装置 | |
CN115933356A (zh) | 一种虚拟原子钟的高精度时间同步系统和方法 | |
WO2024016922A1 (zh) | 基于gnss差分的航天器相对位置测量方法及其装置 | |
CN111948682B (zh) | 一种基于bds三频组合的伪距和载波相位随机模型构建方法 | |
CN114002722A (zh) | Gnss定位方法及其装置 | |
CN104749597A (zh) | 一种针对北斗系统的中长基线下模糊度解算方法 | |
CN114966756A (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 |