CN108802781A - 整周模糊度的获取方法 - Google Patents
整周模糊度的获取方法 Download PDFInfo
- Publication number
- CN108802781A CN108802781A CN201810415192.0A CN201810415192A CN108802781A CN 108802781 A CN108802781 A CN 108802781A CN 201810415192 A CN201810415192 A CN 201810415192A CN 108802781 A CN108802781 A CN 108802781A
- Authority
- CN
- China
- Prior art keywords
- integer ambiguity
- equation
- day part
- solution
- unknown parameter
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 239000011159 matrix material Substances 0.000 claims abstract description 24
- 230000003068 static effect Effects 0.000 claims description 42
- 238000004590 computer program Methods 0.000 claims description 30
- 238000003491 array Methods 0.000 claims description 22
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 abstract description 13
- 238000010586 diagram Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 238000007667 floating Methods 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 230000008901 benefit Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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
本发明涉及一种整周模糊度的获取方法、***、计算机设备以及存储介质,该方法包括:根据各时段的双差观测方程分别构建对应的分块最小二乘平差方程;根据各分块最小二乘平差方程构建序贯观测方程,并解算序贯观测方程获得位置未知参数;利用位置未知参数根据各分块最小二乘平差方程获得各时段的整周模糊度浮点解;根据各整周模糊度浮点解获取各时段的整周模糊度固定解。采用本方法能够在确保整周模糊度解算精度的前提下,降低整周模糊度解算过程中矩阵运算的阶数,大大节省数据处理时间,实现快速解算整周模糊度,提高解算效率。
Description
技术领域
本发明涉及卫星导航技术领域,特别是涉及一种整周模糊度的获取方法、***、计算机设备以及存储介质。
背景技术
随着GNSS(Global Navigation Satellite System,全球卫星导航***)的不断发展、完善,GNSS已经成为当今获得高精度定位结果的主要手段,主要包括:美国的GPS***、俄罗斯的GLONASS***、中国的北斗***(BeiDou Navigation Satellite System,BDS)以及欧盟的GALILEO***。GNSS在测量领域的一个重要应用是建立工程测量控制网,主要由静态数据的基线解算和网平差两阶段实现。其中,基线解算是网平差的基础,直接影响GNSS控制网的定位质量和可靠性,基线解算过程中,最为关键的是快速准确的解算出整周模糊度。
传统的基线解算过程中,通常通过建立多卫星***、多时段的组合方程统一平差获得整周模糊度的浮点解,但是在实际情况中,采取将多卫星导航***以及采集多时段观测数据会导致计算数据大大增加,导致整周模糊度的解算过程涉及上百阶的矩阵运算,计算时间过长,降低整周模糊度的解算效率,严重时甚至出无法解出整周模糊度的情况,影响基线解算结果。
发明内容
基于此,有必要针对整周模糊度的解算过程数据量过大导致整周模糊度的解算时间长、效率低下的问题,提供一种整周模糊度的获取方法、***、计算机设备以及存储介质。
一种整周模糊度的获取方法,包括以下步骤:
根据各时段的双差观测方程分别构建各时段对应的分块最小二乘平差方程;
根据各所述分块最小二乘平差方程构建序贯观测方程,并解算所述序贯观测方程获得位置未知参数;
利用所述位置未知参数并根据各所述分块最小二乘平差方程获取各时段的整周模糊度浮点解;
根据各所述整周模糊度浮点解获取各时段的整周模糊度固定解。
在其中一个实施例中,所述根据各时段的双差观测方程分别构建各时段对应的分块最小二乘平差方程的步骤之前,还包括以下步骤:
获取各时段的静态观测数据,根据各所述静态观测数据构建各时段的所述双差观测方程。
在其中一个实施例中,所述分块最小二乘平差方程形式如下:
其中,nd为时段数量,δX为所述位置位置参数,为整周模糊度未知参数,为所述位置未知参数的系数阵,为整周模糊度未知参数的系数阵,LδX为所述位置未知参数的观测值,为整周模糊度未知参数的观测数据,VδX为所述位置未知参数的观测方程的残差,为整周模糊度未知参数的观测方程的残差。
在其中一个实施例中,所述序贯观测方程如下:
其中,nd为时段数量,δX为所述位置位置参数,为各时段下观测值的权值,为所述位置未知参数的系数阵,为各时段下位置未知参数的观测值。
在其中一个实施例中,所述根据各所述整周模糊度浮点解获取各时段的整周模糊度固定解的步骤,包括以下步骤:
根据各所述分块最小二乘平差方程获取各时段对应的整周模糊度方差协方差矩阵;
对各所述整周模糊度方差协方差矩阵进行变换,确定各时段整周模糊度搜索空间;
根据各所述整周模糊度浮点解从各所述整周模糊度搜索空间中搜索最优模糊度解,获取各时段的整周模糊度固定解。
一种整周模糊度的获取***,包括:
最小二乘方程获取模块,用于根据各时段的双差观测方程分别构建各时段对应的分块最小二乘平差方程;
位置未知参数获取模块,用于根据各所述分块最小二乘平差方程构建序贯观测方程,并解算所述序贯观测方程获得位置未知参数;
浮点解获取模块,用于利用所述位置未知参数并根据各所述分块最小二乘平差方程获取各时段的整周模糊度浮点解;
固定解获取模块,用于根据各所述整周模糊度浮点解获取各时段的整周模糊度固定解。
一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
根据各时段的双差观测方程分别构建各时段对应的分块最小二乘平差方程;
根据各所述分块最小二乘平差方程构建序贯观测方程,并解算所述序贯观测方程获得位置未知参数;
利用所述位置未知参数根据各所述分块最小二乘平差方程获得各时段的整周模糊度浮点解;
根据各所述整周模糊度浮点解获取各时段的整周模糊度固定解。
一种计算机存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现以下步骤:
根据各时段的双差观测方程分别构建各时段对应的分块最小二乘平差方程;
根据各所述分块最小二乘平差方程构建序贯观测方程,并解算所述序贯观测方程获得位置未知参数;
利用所述位置未知参数根据各所述分块最小二乘平差方程获得各时段的整周模糊度浮点解;
根据各所述整周模糊度浮点解获取各时段的整周模糊度固定解。
上述整周模糊度的获取方法、***、计算机设备以及存储介质,通过各个时段的双差观测方程分别构建最小二乘法平差方程以及序贯观测方程,在确保整周模糊度解算精度的前提下,降低整周模糊度解算过程中矩阵运算的阶数,大大节省数据处理时间,实现快速解算整周模糊度,提高解算效率。
一种静态基线的解算方法,包括以下步骤:
利用上述的整周模糊度的获取方法获取各时段的整周模糊度的固定解;
根据各所述整周模糊度的固定解利用所述分块最小二乘平差方程获取位置未知参数的组合解。
一种静态基线的解算***,包括:
整周模糊度获取模块,用于利用权利要求上述的整周模糊度的获取方法获取各时段的整周模糊度的固定解;
静态基线解算模块,用于根据各所述整周模糊度的固定解利用所述分块最小二乘平差方程获取位置未知参数的组合解。
一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
利用上述的整周模糊度的获取方法获取各时段的整周模糊度的固定解;
根据各所述整周模糊度的固定解利用所述分块最小二乘平差方程获取位置未知参数的组合解。
一种计算机存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现以下步骤:
利用上述的整周模糊度的获取方法获取各时段的整周模糊度的固定解;
根据各所述整周模糊度的固定解利用所述分块最小二乘平差方程获取位置未知参数的组合解。
上述静态基线的解算方法、***、计算机设备以及存储介质,通过各个时段的双差观测方程分别构建最小二乘法平差方程以及序贯观测方程,解算整周模糊度以及位置未知参数组合解,在确保静态解算结果精度的前提下,实现快速获取多卫星导航***的多时段数据的基线解算以及测站位置坐标。
附图说明
图1为本发明一个实施例中整周模糊度的获取方法的应用环境图;
图2为本发明一个实施例中整周模糊度的获取方法的流程图;
图3为本发明一个实施例中站间星际双差示意图;
图4为本发明一个实施例中整周模糊度的获取***的结构示意图;
图5为本发明另一个实施例中整周模糊度的获取***的结构示意图;
图6为本发明一个实施例中计算机设备的内部结构图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明的保护范围。
本申请提供的整周模糊度的获取方法,可以应用于如图1所示的应用环境中。该终端包括通过***总线连接的处理器、存储器、网络接口和数据库。其中,该处理器用于提供计算和控制能力。该存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作***,还包括一种整周模糊度的解算***。该内存储器为非易失性存储介质中的操作***和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信,接收获取卫星或者接收发送的静态观测数据、GNSS星历文件等数据,具体的,根据获取得到的静态观测数据进行整周模糊度的解算。其中,终端包括但不限于计算机、手提电脑等智能设备。
参见图2,图2为本发明一个实施例中整周模糊度的获取方法的流程图,该实施例中整周模糊度的获取方法,包括以下步骤:
步骤S210:根据各时段的双差观测方程分别构建各时段对应的分块最小二乘平差方程。
具体的,获取各个时段的双差观测方程后,将多个时段的双差观测方程分别变换为各个时段对应的分块最小二乘平差方程,通过双差观测方程构建分块最小二乘平差方程,能够很好的消除卫星钟以及眼见延迟参数的误差,确保后续整周模糊度解算的精度。
步骤S220:根据各分块最小二乘平差方程构建序贯观测方程,并解算序贯观测方程获得位置未知参数。
步骤S230:利用位置未知参数并根据各分块最小二乘平差方程获取各时段的整周模糊度浮点解。
具体的,在获取得到位置未知参数后,将位置未知参数回代到各个时段的分块最小二乘平差方程,解算出各个时段的整周模糊度浮点解。
步骤S240:根据各整周模糊度浮点解获取各时段的整周模糊度固定解。
在进行解算整周模糊度的过程中,每个时段的位置未知参数都没变化,而每个时段的整周模糊度在时间上式独立,由此根据各个时段的双差观测方程分别构建分块最小二乘法平差方程以及序贯观测方程,利用序贯观测方程解算位置未知参数,再利用该位置未知参数回传至分块最小二乘法平差方程,求解每个时段的整周模糊度的浮点解,从而根据整周模糊度的浮点解搜索获取整周模糊度的最优解。
上述静态基线的解算方法,通过利用各个时段的双差观测方程分别构建最小二乘法平差方程以及序贯观测方程,在确保整周模糊度解算精度的前提下,降低整周模糊度解算过程中矩阵运算的阶数,将传统的处理方式中上百阶的矩阵求逆转化为三阶的矩阵求逆,大大节省数据处理时间,实现快速解算整周模糊度,提高解算效率。
在其中一个实施例中,根据各时段的双差观测方程分别构建各时段对应的分块最小二乘平差方程的步骤之前,还包括以下步骤:获取各时段的静态观测数据,根据各静态观测数据构建各时段的双差观测方程。
本实施例中,静态观测数据是指观测站接收卫星发送的信号数据;通过采集观测站多个时段的静态观测数据,并利用获取得到的静态观测数据构建每一时段对应的双差观测方程,通过构建双差观测方程能够快速获取对应的分块最小二乘平差方程,提高解算整周模糊度的效率。
以图3为例,图3为本发明一个实施例中站间星际双差示意图,接收获取图中卫星310以及卫星320的静态观测数据,两观测站同步观测的卫星为卫星310以及卫星320,接收获取图中卫星310以及卫星320的静态观测数据,假设卫星310为参考卫星,获得双差观测方程的线性化如(1)式所示:
其中: 为载波观测值,λ为载波波长,l、m、n为线性化后的系数,δX、δY、δZ为近似坐标的改正数,N为整周模糊度,ρ0为近似卫地距。
在其中一个实施例中,分块最小二乘平差方程形式如(2)式:
其中,nd为时段数量,δX为所述位置位置参数,为整周模糊度未知参数,为所述位置未知参数的系数阵,为整周模糊度未知参数的系数阵,LδX为所述位置未知参数的观测值,为整周模糊度未知参数的观测数据,VδX为所述位置未知参数的观测方程的残差,为整周模糊度未知参数的观测方程的残差。
具体的,在获取得到形如式(1)的双差观测方程后,得到如下式(3)的误差方程:
其中,
当两观测站的同步卫星数为ns时,得到如下(4)式的误差方程组:
其中:
为方程残差;
为误差方程组中δX的系数阵;
为误差方程组中▽ΔN1的系数阵;
δX=[δX2 δY2 δZ2]T,为位置未知参数;
为整周模糊度未知参数;
由于在进行解算整周模糊度的过程中,每个时段的位置未知参数都没变化,每个时段的整周模糊度在时间上式独立,当在基线两端、同一时段下同一组卫星观测的历元数为nt,则误差方程组如下(5)为:
其中:
为某一个时段内全部观测方程δX的系数阵;
为某一个时段的系数阵;
为某一个时段的观测值;
为某一个时段方程残差。
本实施例中,每一个时段具有对应一个分块最小二乘平差方程,将传统的整周模糊度解算的处理方式分解为求解单个时段的整周模糊度,大大简化解算过程中的矩阵运算,节省数据处理时间,实现快速解算整周模糊度,提高解算效率。
在其中一个实施例中,序贯观测方程如下(6)式:
其中,nd为时段数量,δX为所述位置位置参数,为各时段下观测值的权值,为所述位置未知参数的系数阵,为各时段下位置未知参数的观测值。
本实施例中,将传统的处理方式中上百阶的矩阵求逆转化为三阶的矩阵求逆,大大节省数据处理时间,实现快速解算整周模糊度,提高解算效率。
在其中一个实施例中,根据各整周模糊度浮点解获取各时段的整周模糊度固定解的步骤,包括以下步骤:根据各分块最小二乘平差方程获取各时段对应的整周模糊度方差协方差矩阵;对各整周模糊度方差协方差矩阵进行变换,确定各时段整周模糊度搜索空间;根据各整周模糊度浮点解从各整周模糊度搜索空间中搜索最优模糊度解,获取各时段的整周模糊度固定解。
本实施例中,获取不同时段对应整周模糊度方差协方差矩阵,并对各个时段的整周模糊度方差协方差矩阵进行变换,以降低各时段间整周模糊度的相关性,从而缩小、确定各个时段的整周模糊度的搜索空间;在确定整周模糊度的搜索空间以后,从中搜索最优的整周模糊度解,获取各时段的整周模糊度固定解,该整周模糊度搜索方法不仅能够快速获取整周模糊度固定解,且准确度高。
应该理解的是,虽然图2的流程图中的各个步骤按照箭头的指示依次显示,但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明,这些步骤的执行并没有严格的顺序限制,这些步骤可以以其它的顺序执行。而且,图2中的至少一部分步骤可以包括多个子步骤或者多个阶段,这些子步骤或者阶段并不必然是在同一时刻执行完成,而是可以在不同的时刻执行,这些子步骤或者阶段的执行顺序也不必然是依次进行,而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
根据上述整周模糊度的获取方法,本发明还提供一种整周模糊度的获取***,以下就本发明的整周模糊度的获取***的实施例进行详细说明。
参见图4,图4为本发明一个实施例中整周模糊度的获取***的结构示意图。本实施例中,整周模糊度的获取***,包括:
最小二乘方程获取模块410,用于根据各双差观测方程构建各时段对应的分块最小二乘平差方程;
位置未知参数获取模块420,用于根据各分块最小二乘平差方程构建序贯观测方程,并解算序贯观测方程获得位置未知参数;
浮点解获取模块430,用于利用位置未知参数并根据各分块最小二乘平差方程获取各时段的整周模糊度浮点解;
固定解获取模块440,用于根据各整周模糊度浮点解获取各时段的整周模糊度固定解。
上述整周模糊度的获取***,通过利用各个时段的双差观测方程分别构建最小二乘法平差方程以及序贯观测方程,在确保整周模糊度解算精度的前提下,降低整周模糊度解算过程中矩阵运算的阶数,将传统的处理方式中上百阶的矩阵求逆转化为三阶的矩阵求逆,大大节省数据处理时间,实现快速解算整周模糊度,提高解算效率。
参见图5,图5为本发明另一个实施例中整周模糊度的获取***的结构示意图,整周模糊度的获取***还包括观测方程获取模块450,观测方程获取模块450用于获取各时段的静态观测数据,根据各静态观测数据构建各时段对应的双差观测方程。
在其中一个实施例中,最小二乘方程获取模块420构建的最小二乘平差方程如下述(7)式:
其中,nd为时段数量,δX为所述位置位置参数,为整周模糊度未知参数,为所述位置未知参数的系数阵,为整周模糊度未知参数的系数阵,LδX为所述位置未知参数的观测值,为整周模糊度未知参数的观测数据,VδX为所述位置未知参数的观测方程的残差,为整周模糊度未知参数的观测方程的残差。
在其中一个实施例中,位置未知参数获取模块430用于构建的序贯观测方程如下述(8)式:
其中,nd为时段数量,δX为所述位置位置参数,为各时段下观测值的权值,为所述位置未知参数的系数阵,为各时段下位置未知参数的观测值。
在其中一个实施例中,固定解获取模块440用于根据各分块最小二乘平差方程获取各时段对应的整周模糊度方差协方差矩阵;对各整周模糊度方差协方差矩阵进行变换,确定各时段整周模糊度搜索空间;根据各整周模糊度浮点解从各整周模糊度搜索空间中搜索最优模糊度解,获取各时段的整周模糊度固定解。
关于整周模糊度的获取***的具体限定可以参见上文中对于整周模糊度的获取方法的限定,在此不再赘述。上述整周模糊度的获取***中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中,也可以以软件形式存储于计算机设备中的存储器中,以便于处理器调用执行以上各个模块对应的操作。
在其中一个实施例中,提供了一种计算机设备,该计算机设备可以是服务器,其内部结构图可以如图6所示。该计算机设备包括通过***总线连接的处理器、存储器、网络接口和数据库。其中,该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作***、计算机程序和数据库。该内存储器为非易失性存储介质中的操作***和计算机程序的运行提供环境。该计算机设备的数据库用于存储静态观测数据。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种整周模糊度的获取方法。
本领域技术人员可以理解,图6中示出的结构,仅仅是与本申请方案相关的部分结构的框图,并不构成对本申请方案所应用于其上的计算机设备的限定,具体的计算机设备可以包括比图中所示更多或更少的部件,或者组合某些部件,或者具有不同的部件布置。
在一个实施例中,提供了一种计算机设备,包括存储器和处理器,存储器中存储有计算机程序,该处理器执行计算机程序时实现以下步骤:
根据各时段的双差观测方程分别构建各时段对应的分块最小二乘平差方程;
根据各分块最小二乘平差方程构建序贯观测方程,并解算序贯观测方程获得位置未知参数;
利用位置未知参数并根据各分块最小二乘平差方程获取各时段的整周模糊度浮点解;
根据各整周模糊度浮点解获取各时段的整周模糊度固定解。
在其中一个实施例中,处理器执行计算机程序时还实现以下步骤:获取各时段的静态观测数据,根据各静态观测数据构建各时段的双差观测方程。
在其中一个实施例中,分块最小二乘平差方程形式如下(9)式:
其中,nd为时段数量,δX为所述位置位置参数,为整周模糊度未知参数,为所述位置未知参数的系数阵,为整周模糊度未知参数的系数阵,LδX为所述位置未知参数的观测值,为整周模糊度未知参数的观测数据,VδX为所述位置未知参数的观测方程的残差,为整周模糊度未知参数的观测方程的残差。
在其中一个实施例中,建序贯观测方程如下(10)式:
其中,nd为时段数量,δX为所述位置位置参数,为各时段下观测值的权值,为所述位置未知参数的系数阵,为各时段下位置未知参数的观测值。
在其中一个实施例中,处理器执行计算机程序时还实现以下步骤:根据各分块最小二乘平差方程获取各时段对应的整周模糊度方差协方差矩阵;对各整周模糊度方差协方差矩阵进行变换,确定各时段整周模糊度搜索空间;根据各整周模糊度浮点解从各整周模糊度搜索空间中搜索最优模糊度解,获取各时段的整周模糊度固定解。
在其中一个实施例中,提供了一种计算机可读存储介质,其上存储有计算机程序,计算机程序被处理器执行时实现以下步骤:
根据各时段的双差观测方程分别构建各时段对应的分块最小二乘平差方程;
根据各分块最小二乘平差方程构建序贯观测方程,并解算序贯观测方程获得位置未知参数;
利用位置未知参数并根据各分块最小二乘平差方程获取各时段的整周模糊度浮点解;
根据各整周模糊度浮点解获取各时段的整周模糊度固定解。
在其中一个实施例中,计算机程序被处理器执行时还实现以下步骤:获取各时段的静态观测数据,根据各静态观测数据构建各时段的双差观测方程。
在其中一个实施例中,分块最小二乘平差方程形式如下(11)式:
其中,nd为时段数量,δX为所述位置位置参数,为整周模糊度未知参数,为所述位置未知参数的系数阵,为整周模糊度未知参数的系数阵,LδX为所述位置未知参数的观测值,为整周模糊度未知参数的观测数据,VδX为所述位置未知参数的观测方程的残差,为整周模糊度未知参数的观测方程的残差。
在其中一个实施例中,序贯观测方程如下(12)式:
其中,nd为时段数量,δX为所述位置位置参数,为各时段下观测值的权值,为所述位置未知参数的系数阵,为各时段下位置未知参数的观测值。
在其中一个实施例中,计算机程序被处理器执行时还实现以下步骤:根据各分块最小二乘平差方程获取各时段对应的整周模糊度方差协方差矩阵;对各整周模糊度方差协方差矩阵进行变换,确定各时段整周模糊度搜索空间;根据各整周模糊度浮点解从各整周模糊度搜索空间中搜索最优模糊度解,获取各时段的整周模糊度固定解。
本发明还提供了一种静态基线的解算方法,静态基线的解算方法包括以下步骤:利用上述任一实施例中的整周模糊度的获取方法获取各时段的整周模糊度的固定解;根据各所述整周模糊度的固定解利用所述分块最小二乘平差方程获取位置未知参数的组合解。
具体的,在获取得到各个时段的整周模糊度固定解后,将各个时段的整周模糊度固定解回传至分块最小二乘平差方程中,采用整体最小二乘平差的方式求解多个时段的位置未知参数的组合解。
上述静态基线的解算方法,通过各个时段的双差观测方程分别构建最小二乘法平差方程以及序贯观测方程,解算整周模糊度以及位置未知参数组合解,在确保静态解算结果精度的前提下,实现快速获取多卫星导航***的多时段数据的基线解算以及测站位置坐标。
根据上述静态基线的解算方法,本发明还提供一种静态基线的解算***,以下就本发明的整周模糊度的获取***的实施例进行详细说明。
在其中一个实施例中,静态基线的解算***,包括:
整周模糊度获取模块,用于利用上述任一实施例中的整周模糊度的获取方法获取各时段的整周模糊度的固定解;
静态基线解算模块,用于根据各所述整周模糊度的固定解利用所述分块最小二乘平差方程获取位置未知参数的组合解。
上述静态基线的解算***,通过各个时段的双差观测方程分别构建最小二乘法平差方程以及序贯观测方程,解算整周模糊度以及位置未知参数组合解,在确保静态解算结果精度的前提下,实现快速获取多卫星导航***的多时段数据的基线解算以及测站位置坐标。
一种计算机设备,包括存储器和处理器,所述存储器存储有计算机程序,所述处理器执行所述计算机程序时实现以下步骤:
利用上述的整周模糊度的获取方法获取各时段的整周模糊度的固定解;
根据各所述整周模糊度的固定解利用所述分块最小二乘平差方程获取位置未知参数的组合解。
一种计算机存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现以下步骤:
利用上述的整周模糊度的获取方法获取各时段的整周模糊度的固定解;
根据各所述整周模糊度的固定解利用所述分块最小二乘平差方程获取位置未知参数的组合解。
本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程,是可以通过计算机程序来指令相关的硬件来完成,所述的计算机程序可存储于一非易失性计算机可读取存储介质中,该计算机程序在执行时,可包括如上述各方法的实施例的流程。其中,本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用,均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限,RAM以多种形式可得,诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。
以上实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种整周模糊度的获取方法,其特征在于,包括以下步骤:
根据各时段的双差观测方程分别构建各时段对应的分块最小二乘平差方程;
根据各所述分块最小二乘平差方程构建序贯观测方程,并解算所述序贯观测方程获得位置未知参数;
利用所述位置未知参数并根据各所述分块最小二乘平差方程获取各时段的整周模糊度浮点解;
根据各所述整周模糊度浮点解获取各时段的整周模糊度固定解。
2.根据权利要求1所述的整周模糊度的获取方法,其特征在于,所述根据各时段的双差观测方程分别构建各时段对应的分块最小二乘平差方程的步骤之前,还包括以下步骤:
获取各时段的静态观测数据,根据各所述静态观测数据构建各时段的所述双差观测方程。
3.根据权利要求1所述的整周模糊度的获取方法,其特征在于,所述分块最小二乘平差方程为:
其中,nd为时段数量,δX为所述位置位置参数,为整周模糊度未知参数,为所述位置未知参数的系数阵,为整周模糊度未知参数的系数阵,LδX为所述位置未知参数的观测值,为整周模糊度未知参数的观测数据,VδX为所述位置未知参数的观测方程的残差,为整周模糊度未知参数的观测方程的残差。
4.根据权利要求1所述的整周模糊度的获取方法,其特征在于,所述序贯观测方程为:
其中,nd为时段数量,δX为所述位置位置参数,为各时段下观测值的权值,为所述位置未知参数的系数阵,为各时段下位置未知参数的观测值。
5.根据权利要求1所述的整周模糊度的获取方法,其特征在于,所述根据各所述整周模糊度浮点解获取各时段的整周模糊度固定解的步骤,包括以下步骤:
根据各所述分块最小二乘平差方程获取各时段对应的整周模糊度方差协方差矩阵;
对各所述整周模糊度方差协方差矩阵进行变换,确定各时段整周模糊度搜索空间;
根据各所述整周模糊度浮点解从各所述整周模糊度搜索空间中搜索最优模糊度解,获取各时段的整周模糊度固定解。
6.一种静态基线的解算方法,其特征在于,包括以下步骤:
利用权利要求1至5任一所述的整周模糊度的获取方法获取各时段的整周模糊度的固定解;
根据各所述整周模糊度的固定解利用所述分块最小二乘平差方程获取位置未知参数的组合解。
7.一种整周模糊度的获取***,其特征在于,包括:
最小二乘方程获取模块,用于根据各时段的双差观测方程分别构建各时段对应的分块最小二乘平差方程;
位置未知参数获取模块,用于根据各所述分块最小二乘平差方程构建序贯观测方程,并解算所述序贯观测方程获得位置未知参数;
浮点解获取模块,用于利用所述位置未知参数并根据各所述分块最小二乘平差方程获取各时段的整周模糊度浮点解;
固定解获取模块,用于根据各所述整周模糊度浮点解获取各时段的整周模糊度固定解。
8.一种静态基线的解算***,其特征在于,包括:
整周模糊度获取模块,用于利用权利要求1至5任一所述的整周模糊度的获取方法获取各时段的整周模糊度的固定解;
静态基线解算模块,用于根据各所述整周模糊度的固定解利用所述分块最小二乘平差方程获取位置未知参数的组合解。
9.一种计算机设备,包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至5任意一项所述的整周模糊度的获取方法或如权利要求6所述的静态基线的解算方法。
10.一种计算机存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1至5任意一项所述的整周模糊度的获取方法或如权利要求6所述的静态基线的解算方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810415192.0A CN108802781B (zh) | 2018-05-03 | 2018-05-03 | 整周模糊度的获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810415192.0A CN108802781B (zh) | 2018-05-03 | 2018-05-03 | 整周模糊度的获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108802781A true CN108802781A (zh) | 2018-11-13 |
CN108802781B CN108802781B (zh) | 2021-03-02 |
Family
ID=64093198
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810415192.0A Expired - Fee Related CN108802781B (zh) | 2018-05-03 | 2018-05-03 | 整周模糊度的获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108802781B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110412633A (zh) * | 2019-07-30 | 2019-11-05 | 广州市中海达测绘仪器有限公司 | 定位方法、装置、计算机设备和存储介质 |
CN110646820A (zh) * | 2019-09-20 | 2020-01-03 | 广州市中海达测绘仪器有限公司 | Rtk定位数据的质检方法、装置、设备和存储介质 |
CN110907975A (zh) * | 2019-12-13 | 2020-03-24 | 北京遥测技术研究所 | 一种基于序贯最小二乘的模糊度固定方法 |
CN111381264A (zh) * | 2018-12-27 | 2020-07-07 | 北京六分科技有限公司 | 网络rtk中长基线模糊度固定方法和平台 |
CN114662268A (zh) * | 2021-11-02 | 2022-06-24 | 广州市城市规划勘测设计研究院 | 一种改进的gnss网序贯平差计算方法 |
CN115840241A (zh) * | 2022-12-15 | 2023-03-24 | 北京六分科技有限公司 | 部分模糊度固定方法、设备、存储介质及程序产品 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101770033A (zh) * | 2010-02-08 | 2010-07-07 | 东南大学 | 连续运行参考站***站间整周模糊度网络固定方法 |
CN103175516A (zh) * | 2013-02-26 | 2013-06-26 | 中国人民解放军信息工程大学 | 规模化大地控制网平差的分布式计算方法 |
CN103487821A (zh) * | 2012-06-12 | 2014-01-01 | 广州中海达卫星导航技术股份有限公司 | 一种基线向量解算方法 |
CN104635249A (zh) * | 2015-02-12 | 2015-05-20 | 东南大学 | 一种顾及FCBs的精密单点定位模糊度快速固定方法 |
-
2018
- 2018-05-03 CN CN201810415192.0A patent/CN108802781B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101770033A (zh) * | 2010-02-08 | 2010-07-07 | 东南大学 | 连续运行参考站***站间整周模糊度网络固定方法 |
CN103487821A (zh) * | 2012-06-12 | 2014-01-01 | 广州中海达卫星导航技术股份有限公司 | 一种基线向量解算方法 |
CN103175516A (zh) * | 2013-02-26 | 2013-06-26 | 中国人民解放军信息工程大学 | 规模化大地控制网平差的分布式计算方法 |
CN104635249A (zh) * | 2015-02-12 | 2015-05-20 | 东南大学 | 一种顾及FCBs的精密单点定位模糊度快速固定方法 |
Non-Patent Citations (3)
Title |
---|
曹新运: "CORS网参考站整周模糊度解算及大气误差建模", 《中国优秀硕士学位论文全文数据库 基础科学辑》 * |
李昌贵: "基于网格的网络导航服务关键技术研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
黄劲松 等: "GPS快速静态定位技术", 《武测科技》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111381264A (zh) * | 2018-12-27 | 2020-07-07 | 北京六分科技有限公司 | 网络rtk中长基线模糊度固定方法和平台 |
CN111381264B (zh) * | 2018-12-27 | 2023-08-15 | 北京六分科技有限公司 | 网络rtk中长基线模糊度固定方法和平台 |
CN110412633A (zh) * | 2019-07-30 | 2019-11-05 | 广州市中海达测绘仪器有限公司 | 定位方法、装置、计算机设备和存储介质 |
CN110646820A (zh) * | 2019-09-20 | 2020-01-03 | 广州市中海达测绘仪器有限公司 | Rtk定位数据的质检方法、装置、设备和存储介质 |
CN110646820B (zh) * | 2019-09-20 | 2021-11-30 | 广州市中海达测绘仪器有限公司 | Rtk定位数据的质检方法、装置、设备和存储介质 |
CN110907975A (zh) * | 2019-12-13 | 2020-03-24 | 北京遥测技术研究所 | 一种基于序贯最小二乘的模糊度固定方法 |
CN110907975B (zh) * | 2019-12-13 | 2021-10-01 | 北京遥测技术研究所 | 一种基于序贯最小二乘的模糊度固定方法 |
CN114662268A (zh) * | 2021-11-02 | 2022-06-24 | 广州市城市规划勘测设计研究院 | 一种改进的gnss网序贯平差计算方法 |
CN115840241A (zh) * | 2022-12-15 | 2023-03-24 | 北京六分科技有限公司 | 部分模糊度固定方法、设备、存储介质及程序产品 |
CN115840241B (zh) * | 2022-12-15 | 2024-02-23 | 北京六分科技有限公司 | 部分模糊度固定方法、设备、存储介质及程序产品 |
Also Published As
Publication number | Publication date |
---|---|
CN108802781B (zh) | 2021-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108802781A (zh) | 整周模糊度的获取方法 | |
Zhang et al. | PPP-RTK based on undifferenced and uncombined observations: theoretical and practical aspects | |
CN114518586B (zh) | 一种基于球谐展开的gnss精密单点定位方法 | |
CN105158783B (zh) | 一种实时动态差分定位方法及其设备 | |
CN109085628B (zh) | 一种整周模糊度的固定方法及*** | |
CN109477900A (zh) | 全球导航卫星***接收器中用于模糊度解算的频率间偏差的估算 | |
CN110646820B (zh) | Rtk定位数据的质检方法、装置、设备和存储介质 | |
Xiang et al. | Reducing convergence time of precise point positioning with ionospheric constraints and receiver differential code bias modeling | |
Ke et al. | Integrating GPS and LEO to accelerate convergence time of precise point positioning | |
CN109085629B (zh) | Gnss基线向量解算定位方法、装置和导航定位设备 | |
CN105974440B (zh) | 一种获取观测噪声的方法和装置 | |
CN112230252A (zh) | 终端定位方法、装置、计算机设备和存储介质 | |
CN114355758B (zh) | 基于北斗星基增强和精密星历服务的高精度时间溯源方法 | |
CN116840876B (zh) | 适用于相位失真的卫星导航天线的双天线定向方法及装置 | |
CN112146557A (zh) | 一种基于gnss的实时桥梁变形监测***及方法 | |
CN110398758A (zh) | 实时钟差估计中的粗差探测方法、装置、设备及存储介质 | |
CN110488332B (zh) | 基于网络rtk技术的定位信息处理方法和装置 | |
CN111856513A (zh) | 卫星观测值获取方法、装置、计算机设备和存储介质 | |
CN112462396A (zh) | 一种高采样率导航卫星钟差的实时并行确定方法 | |
CN108254774A (zh) | 基于gnss多频信号的单基站长距离实时定位方法 | |
CN117255956A (zh) | 超宽巷(uwl)实时动态(rtk)定位 | |
CN110673172A (zh) | 一种接收机静态相对定位精度的测试方法及终端 | |
CN110568464A (zh) | 基于bds/gnss多模芯片的精密定位方法及装置 | |
US11668840B2 (en) | Ultra-long baseline RTK | |
CN115956211A (zh) | 全球导航卫星***(gnss)接收器的多频实时动态(rtk)测量 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210302 |