CN111060943A - 一种符合实际的gnss观测值定权方法 - Google Patents

一种符合实际的gnss观测值定权方法 Download PDF

Info

Publication number
CN111060943A
CN111060943A CN201911338241.6A CN201911338241A CN111060943A CN 111060943 A CN111060943 A CN 111060943A CN 201911338241 A CN201911338241 A CN 201911338241A CN 111060943 A CN111060943 A CN 111060943A
Authority
CN
China
Prior art keywords
observation
gnss
weighting
satellite
value
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
Application number
CN201911338241.6A
Other languages
English (en)
Other versions
CN111060943B (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.)
Hohai University HHU
Original Assignee
Hohai University HHU
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 Hohai University HHU filed Critical Hohai University HHU
Priority to CN201911338241.6A priority Critical patent/CN111060943B/zh
Publication of CN111060943A publication Critical patent/CN111060943A/zh
Application granted granted Critical
Publication of CN111060943B publication Critical patent/CN111060943B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/43Determining position using carrier phase measurements, e.g. kinematic positioning; using long or short baseline interferometry
    • G01S19/44Carrier phase ambiguity resolution; Floating ambiguity; LAMBDA [Least-squares AMBiguity Decorrelation Adjustment] method

Abstract

本发明公开了一种符合实际的GNSS观测值定权方法,方法包括获取GNSS观测值;确定各观测值对应的卫星高度角,卫星方位角和载噪比信息;构建顾及卫星高度角,卫星方位角和载噪比的定权函数;将定权函数用于构建随机模型;实现GNSS导航定位。本发明顾及了传统观测值定权方法中忽略掉的一些实际环境因素引起的观测值质量变化,由于这些测站环境因素往往与卫星高度角,卫星方位角和载噪比密切相关,本发明充分考虑了这些因素,并建立了相应的观测值定权方法。该方法与传统定权方法相比,可以提高GNSS导航定位精度和可靠性。

Description

一种符合实际的GNSS观测值定权方法
技术领域
本发明涉及卫星导航定位技术领域,具体涉及一种符合实际的GNSS观测值定权方法。
背景技术:
在全球导航卫星系统(Global Navigation Satellite System,GNSS)中,准确的评估观测值的精度,从而构建相应准确的定权函数,是实现高精度高可靠性GNSS导航定位的前提。
之前,通常认为GNSS观测值具有同方差性,即等权模型,认为所有的观测值精度相等,但是这种方法在大多观测环境条件下不适用。因此为了克服观测值同方差性方面的缺陷,往往采用高度角或载噪比定权函数来实现观测值异方差性。高度角定权函数模型认为低高度角的卫星受到更严重的大气影响,还会使接收机天线增益减少,从而导致信号质量下降,因此认为高度角越高的卫星权重越大。载噪比定权函数模型认为载噪比和GNSS观测值由同一个跟踪环路记录,因此载噪比能直接评价信号质量,即载噪比越大,数据质量越好。
事实上,载噪比定权函数模型对测站相关的误差较为敏感。具体来说,当载噪比值较低时往往意味着观测值受到了多路径甚至是非视距传播的影响,因此观测值精度降低。此外,当信号扭曲严重的时候,其观测值会偏离标称值。载噪比这一重要特点目前并未充分考虑,造成之前的载噪比定权函数模型仅仅是间接和部分的反应观测值精度。此外,不同的方位角,接收机收到的观测值信号会受到测站环境周围植物以及建筑物等的影响,观测值的精度也会随之受到影响。然而,目前的定权函数模型中并未完全考虑这一特性。
发明内容
本发明的目的在于提供一种符合实际的GNSS观测值定权方法,以解决现有技术中导致的GNSS导航定位误差大不精确的缺陷。
一种符合实际的GNSS观测值定权方法,所述方法包括如下步骤:
获取GNSS观测值;
确定各观测值对应的参数数据;
根据观测数据构建观测值精度估计公式;
根据观测值精度估计公式构建定权函数;
通过定权函数实现GNSS导航定位。
进一步的,所述参数数据包括卫星高度角,卫星方位角以及载噪比中的一种或多种。
进一步的,所述观测数据的确定方法包括如下步骤:
获取观测文件和星历文件;
根据观测文件得到某个时刻某个观测值对应的载噪比数据;
通过星历文件计算该时刻的卫星位置;
根据观测文件求出GNSS接收设备的位置范围;
根据该时刻的卫星位置以及接收设备的位置范围计算得出该观测值对应的卫星高度角和卫星方位角。
进一步的,所述GNSS观测值包括不同频率的相位观测值和伪距观测值。
进一步的,根据参数数据构建定权函数的方法包括如下步骤:
确定观测值精度,计算公式为:
Figure BDA0002331548510000031
式中,σ2为GNSS观测值的方差;a为常数系数;θ为观测值对应的高度角;η为观测值对应的载噪比;η*(θ)为观测环境理想情况下高度角θ所对应的载噪比;η*(90°)为观测环境理想情况下高度角90°所对应的载噪比;|.|代表求绝对值;
Figure BDA0002331548510000032
为观测值对应的方位角;
Figure BDA0002331548510000033
为方位角
Figure BDA0002331548510000034
所对应的由用户定义的惩罚函数。
进一步的,根据观测值精度估计公式构建定权函数的方法包括如下步骤:
将观测值精度估计公式用于构建定权函数,计算公式为:
Figure BDA0002331548510000035
式中,
Figure BDA0002331548510000036
为第i个观测值的方差;Pi为第i个观测值的权。
进一步的,通过定权函数实现GNSS导航定位的方法包括如下步骤:
构建观测数据的函数模型;
通过定权函数构建观测数据的随机模型;
通过函数模型与随机模型对位置参数进行求解,实现GNSS导航定位。
本发明的优点在于:该种符合实际的GNSS观测值定权方法,
1、顾及了传统观测值定权方法中忽略掉的一些测站环境因素引起的观测值质量变化,该方法与传统定权方法相比,更符合实际情况;
2、通过充分考虑测站环境因素对卫星高度角,卫星方位角和载噪比的影响,建立了相应的观测值定权方法,该方案操作简便,可执行性强;
3、利用该定权函数构建的随机模型,由于更能准确可靠地反应观测值精度,因此可以提高GNSS导航定位精度和可靠性。
附图说明
图1为本发明的流程示意图。
具体实施方式
为使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体实施方式,进一步阐述本发明。
如图1所示,一种符合实际的GNSS观测值定权方法,所述方法包括如下步骤:
步骤一:获取GNSS观测值:
获取不同频率的全球定位系统(GPS)相位观测值和伪距观测值;
步骤二:确定各观测值对应的参数数据,参数数据包括卫星高度角,卫星方位角以及载噪比:
所述确定各观测值对应的高度角,方位角和载噪比信息,其中首先获取观测文件和星历文件,根据观测文件得到某个时刻某个观测值对应的载噪比数据;接着通过星历文件计算该时刻的卫星位置,根据观测文件求出GPS接收设备的位置范围;最后根据该时刻的卫星位置以及接收设备的位置范围计算得出该观测值对应的卫星高度角和卫星方位角。
步骤三:根据观测数据构建观测值精度估计公式:
确定观测值精度,计算公式为:
Figure BDA0002331548510000051
式中,σ2为GPS观测值的方差;a为常数系数;θ为观测值对应的高度角;η为观测值对应的载噪比;η*(θ)为观测环境理想情况下高度角θ所对应的载噪比;η*(90°)为观测环境理想情况下高度角90°所对应的载噪比;|.|代表求绝对值;
Figure BDA0002331548510000052
为观测值对应的方位角;
Figure BDA0002331548510000053
为方位角
Figure BDA0002331548510000054
所对应的由用户定义的惩罚函数。此时,由于是开阔环境,因此该惩罚函数为
Figure BDA0002331548510000055
步骤四:根据观测值精度估计公式构建定权函数:
将观测值精度估计公式用于构建定权函数,计算公式为:
Figure BDA0002331548510000056
式中,
Figure BDA0002331548510000057
为第i个观测值的方差;Pi为第i个观测值的权;
步骤五:通过定权函数实现GNSS导航定位:
首先,构建观测数据的函数模型;接着通过定权函数构建观测数据的随机模型;最后,对位置参数进行求解,实现GPS导航定位;
获取GNSS观测数据并进行粗差探测与剔除,数据中断与修复,周跳探测与修复等预处理,由观测值构建误差方程,如需要固定相位观测值的模糊度,需要求解其中相位观测值的模糊度浮点解并进行固定,再次构建误差方程,即最终形成观测数据的函数模型。
由本发明的定权方法构建定权函数,从而确定观测数据的随机模型。
通过一定的准则如最小二乘准则,联立上述函数模型与随机模型,最终对含位置在内等参数进行求解,实现GNSS导航定位。
由技术常识可知,本发明可以通过其它的不脱离其精神实质或必要特征的实施方案来实现。因此,上述公开的实施方案,就各方面而言,都只是举例说明,并不是仅有的。所有在本发明范围内或在等同于本发明的范围内的改变均被本发明包含。

Claims (7)

1.一种符合实际的GNSS观测值定权方法,其特征在于,所述方法包括如下步骤:
获取GNSS观测值;
确定各观测值对应的参数数据;
根据观测数据构建观测值精度估计公式;
根据观测值精度估计公式构建定权函数;
通过定权函数实现GNSS导航定位。
2.根据权利要求1所述的一种符合实际的GNSS观测值定权方法,其特征在于:所述参数数据包括卫星高度角,卫星方位角以及载噪比中的一种或多种。
3.根据权利要求2所述的一种符合实际的GNSS观测值定权方法,其特征在于:所述观测数据的确定方法包括如下步骤:
获取观测文件和星历文件;
根据观测文件得到某个时刻某个观测值对应的载噪比数据;
通过星历文件计算该时刻的卫星位置;
根据观测文件求出GNSS接收设备的位置范围;
根据该时刻的卫星位置以及接收设备的位置范围计算得出该观测值对应的卫星高度角和卫星方位角。
4.根据权利要求1所述的一种符合实际的GNSS观测值定权方法,其特征在于:所述GNSS观测值包括不同频率的相位观测值和伪距观测值。
5.根据权利要求1所述的一种符合实际的GNSS观测值定权方法,其特征在于:根据参数数据构建定权函数的方法包括如下步骤:
确定观测值精度,计算公式为:
Figure FDA0002331548500000021
式中,σ2为GNSS观测值的方差;a为常数系数;θ为观测值对应的高度角;η为观测值对应的载噪比;η*(θ)为观测环境理想情况下高度角θ所对应的载噪比;η*(90°)为观测环境理想情况下高度角90°所对应的载噪比;|·|代表求绝对值;θ为观测值对应的方位角;β(θ)为方位角θ所对应的由用户定义的惩罚函数。
6.根据权利要求1所述的一种符合实际的GNSS观测值定权方法,其特征在于:根据观测值精度估计公式构建定权函数的方法包括如下步骤:
将观测值精度估计公式用于构建定权函数,计算公式为:
Figure FDA0002331548500000022
式中,σi 2为第i个观测值的方差;Pi为第i个观测值的权。
7.根据权利要求1所述的一种符合实际的GNSS观测值定权方法,其特征在于:通过定权函数实现GNSS导航定位的方法包括如下步骤:
构建观测数据的函数模型;
通过定权函数构建观测数据的随机模型;
通过函数模型与随机模型对位置参数进行求解,实现GNSS导航定位。
CN201911338241.6A 2019-12-23 2019-12-23 一种符合实际的gnss观测值定权方法 Active CN111060943B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911338241.6A CN111060943B (zh) 2019-12-23 2019-12-23 一种符合实际的gnss观测值定权方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911338241.6A CN111060943B (zh) 2019-12-23 2019-12-23 一种符合实际的gnss观测值定权方法

Publications (2)

Publication Number Publication Date
CN111060943A true CN111060943A (zh) 2020-04-24
CN111060943B CN111060943B (zh) 2023-12-01

Family

ID=70300797

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911338241.6A Active CN111060943B (zh) 2019-12-23 2019-12-23 一种符合实际的gnss观测值定权方法

Country Status (1)

Country Link
CN (1) CN111060943B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112198535A (zh) * 2020-12-07 2021-01-08 湖南联智科技股份有限公司 一种gnss监测设备及监测点数据质量分析方法
CN112799106A (zh) * 2021-03-30 2021-05-14 腾讯科技(深圳)有限公司 卫星观测信息处理方法、移动终端及存储介质
CN113281796A (zh) * 2021-07-23 2021-08-20 腾讯科技(深圳)有限公司 位置确定方法、速度确定方法、装置、设备和存储介质
CN114217328A (zh) * 2022-02-21 2022-03-22 长沙金维信息技术有限公司 导航电文半周跳变检测方法
CN114485655A (zh) * 2022-04-06 2022-05-13 深圳华大北斗科技股份有限公司 Gnss/ins组合导航数据质量控制方法
CN117647830A (zh) * 2024-01-29 2024-03-05 中国测绘科学研究院 一种适用于复杂城市环境gnss芯片定位的随机模型构建方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101799525A (zh) * 2009-07-10 2010-08-11 中国测绘科学研究院 全球导航卫星系统增强系统监测站的自主完备性监测方法
US20110090117A1 (en) * 2008-03-20 2011-04-21 Filippo Rodriguez Position estimation enhancement for a global navigation satellite system receiver
CN104102822A (zh) * 2014-07-01 2014-10-15 同济大学 一种多频gnss观测值随机特性建模方法
CN106814382A (zh) * 2017-01-11 2017-06-09 武汉大学 联合高度角和用户测距误差的gnss卫星观测值定权方法
CN110018507A (zh) * 2019-05-08 2019-07-16 中国科学院国家授时中心 一种基于星座间作差的组合精密单点定位方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110090117A1 (en) * 2008-03-20 2011-04-21 Filippo Rodriguez Position estimation enhancement for a global navigation satellite system receiver
CN101799525A (zh) * 2009-07-10 2010-08-11 中国测绘科学研究院 全球导航卫星系统增强系统监测站的自主完备性监测方法
CN104102822A (zh) * 2014-07-01 2014-10-15 同济大学 一种多频gnss观测值随机特性建模方法
CN106814382A (zh) * 2017-01-11 2017-06-09 武汉大学 联合高度角和用户测距误差的gnss卫星观测值定权方法
CN110018507A (zh) * 2019-05-08 2019-07-16 中国科学院国家授时中心 一种基于星座间作差的组合精密单点定位方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
ZHETAO ZHANG 等: "Site-Specific Unmodeled Error Mitigation for GNSS Positioning in Urban Environments Using a Real-Time Adaptive Weighting Model" *
冯彪 等: "一种改进的组合导航系统选星算法" *
王敏 等: "一种顾及GEO卫星多路径效应影响的BDS长距离相对定位方法" *

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112198535A (zh) * 2020-12-07 2021-01-08 湖南联智科技股份有限公司 一种gnss监测设备及监测点数据质量分析方法
CN112799106A (zh) * 2021-03-30 2021-05-14 腾讯科技(深圳)有限公司 卫星观测信息处理方法、移动终端及存储介质
CN112799106B (zh) * 2021-03-30 2021-06-29 腾讯科技(深圳)有限公司 卫星观测信息处理方法、移动终端及存储介质
CN113281796A (zh) * 2021-07-23 2021-08-20 腾讯科技(深圳)有限公司 位置确定方法、速度确定方法、装置、设备和存储介质
CN113281796B (zh) * 2021-07-23 2021-10-15 腾讯科技(深圳)有限公司 位置确定方法、速度确定方法、装置、设备和存储介质
CN114217328A (zh) * 2022-02-21 2022-03-22 长沙金维信息技术有限公司 导航电文半周跳变检测方法
CN114217328B (zh) * 2022-02-21 2022-05-17 长沙金维信息技术有限公司 导航电文半周跳变检测方法
CN114485655A (zh) * 2022-04-06 2022-05-13 深圳华大北斗科技股份有限公司 Gnss/ins组合导航数据质量控制方法
CN117647830A (zh) * 2024-01-29 2024-03-05 中国测绘科学研究院 一种适用于复杂城市环境gnss芯片定位的随机模型构建方法

Also Published As

Publication number Publication date
CN111060943B (zh) 2023-12-01

Similar Documents

Publication Publication Date Title
CN111060943A (zh) 一种符合实际的gnss观测值定权方法
Geng et al. Kinematic precise point positioning at remote marine platforms
US8412456B2 (en) Loosely-coupled integration of global navigation satellite system and inertial navigation system: speed scale-factor and heading bias calibration
US8860609B2 (en) Loosely-coupled integration of global navigation satellite system and inertial navigation system
US7382313B1 (en) Method for absolute calibration of global navigation satellite system antennas
JP2017111112A (ja) 長周期波高測定用gnss基盤ブイにおける受信機の位相ワインドアップ補償のためのppp及びマグネチックコンパス統合システム{ppp and magnetic compass integration system for measuring of long period wave in buoy’s wind−up}
CN114488230A (zh) 一种多普勒定位方法、装置、电子设备和存储介质
CN115616623A (zh) 一种基于当前模型卡尔曼滤波的卫星导航授时方法和装置
CN115980317B (zh) 基于修正相位的地基gnss-r数据土壤水分估算方法
US11913809B2 (en) Systems and methods for extending the spatial coverage of a reference pressure network
CN111060936B (zh) 基于信噪比的bds/gnss多路径探测方法
CN116203603A (zh) 一种面向复杂环境低成本的ppp随机模型自适应调整方法
Dietmayer et al. Real time results of vector delay lock loop in a light urban scenario
Narula et al. Accuracy limits for globally-referenced digital mapping using standard GNSS
Wang et al. Analysis of GNSS-R Code-Level Altimetry using QZSS C/A, L1C, and BDS B1C signals and their Combinations in a Coastal Experiment
Wang et al. Instantaneous Sub-meter Level Precise Point Positioning of Low-Cost Smartphones
TOMINAGA et al. Performance assessment of adaptive Kalman filter-based GNSS PVT and integrity in dense urban environment
CN113359204A (zh) 一种基于gnss-r伪随机噪声码延迟的水下重力场反演方法及系统
Reis et al. Impact of satellite coverage in single-frequency precise heading determination
Guan et al. Signal quality analysis and quality check of BDS3 Precise Point Positioning in the Arctic Ocean
Hu et al. Towards GNSS Ambiguity Resolution for Smartphones in Realistic Environments: Characterization of Smartphone Ambiguities with RTK, PPP, and PPP-RTK
Raghuvanshi et al. Precise positioning of smartphones using a robust adaptive Kalman filter
KR102371402B1 (ko) 반송파 위상 측정치를 이용한 전리층 지연 절대값 추정 시스템 및 전리층 보정정보 생성 시스템
CN113917509B (zh) 一种双差模糊度固定方法、设备以及可读存储介质
CN111399019B (zh) 一种gnss中卫星定位方法和装置

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