CN113341437B - 一种gnss未建模误差反演方法 - Google Patents

一种gnss未建模误差反演方法 Download PDF

Info

Publication number
CN113341437B
CN113341437B CN202110558650.8A CN202110558650A CN113341437B CN 113341437 B CN113341437 B CN 113341437B CN 202110558650 A CN202110558650 A CN 202110558650A CN 113341437 B CN113341437 B CN 113341437B
Authority
CN
China
Prior art keywords
error
satellite
unmodeled
ground
clock
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
CN202110558650.8A
Other languages
English (en)
Other versions
CN113341437A (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 CN202110558650.8A priority Critical patent/CN113341437B/zh
Publication of CN113341437A publication Critical patent/CN113341437A/zh
Application granted granted Critical
Publication of CN113341437B publication Critical patent/CN113341437B/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/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/23—Testing, monitoring, correcting or calibrating of receiver elements
    • 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)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种GNSS非差非组合未建模误差反演方法，属于卫星导航定位领域，基于精密单点定位中非差非组合给出了一种未建模误差计算方法，可用于定位过程误差的进一步精化研究。包括以下步骤：(1)根据卫星坐标和地面接收机坐标计算卫星与地面之间的距离；(2)通过观测文件获取卫星的载波观测值，通过钟差文件获取卫星的时钟误差；(3)获取对流层延迟误差、电离层延迟误差、整周模糊度、接收机钟差以及与频率无关的各误差项的值；(4)利用卫星与地面之间的距离以及步骤二和步骤三中获取的各项值通过相应的数学算法反演卫星定位过程中的未建模误差。

Description

一种GNSS未建模误差反演方法
技术领域
本发明属于卫星导航定位领域，涉及定位过程中的误差处理问题，特别是对于采用参数估计或模型改正后的残余***误差问题。
背景技术
全球导航卫星***(Global Navigation Satellite System,GNSS)是一种全球、全天候、全天时、高精度的无线电导航定位***，能够为用户提供三维位置和时间等信息。随着导航定位技术的不断完善和发展，GNSS已经广泛应用于大地测量、气象预测、交通物流等诸多生活领域。随着人们对卫星定位精度和可靠性的要求越来越高，如何进一步提高卫星的定位精度和可靠性成为目前GNSS的研究热点与迫切需求。
影响卫星定位精度的因素有很多，包括卫星端，卫星信号传播路径上和地面接收机等方面的***误差。然而由于大气扰动和多路径效应等误差复杂的时空特性以及人们对其物理知识的认知有限，难以完全通过差分、线性组合、模型改正和参数化估计等方法对这些误差进行彻底的消除，还存在一些残余的***误差，这些残余的***误差被称为未建模误差。在绝大部分情况下未建模误差是客观存在的，是制约高精度和高可靠性GNSS应用的瓶颈，如果忽略未建模误差的影响，会人为地提高参数估值的精度。
目前，对于未建模误差的处理主要通过半参数化的方法进行削弱，没有考虑未建模误差特性对观测值精度的影响，缺少对未建模误差大小及其变化规律的研究，亟需一种未建模误差的计算方法。因此，在高精度定位中，必须考虑未建模误差的影响，本发明给出一种未建模误差反演方法，用以反映未建模误差的大小与变化规律，为削弱或消除未建模误差影响的研究提供理论与数据支持，进而提高定位精度。
发明内容
要解决的技术问题：针对现有技术的不足，本发明提出一种GNSS非差非组合未建模误差反演方法，用于解决误差处理过程中存在残余***误差难以处理的问题。
为解决上述技术问题，本发明采用以下技术方案：一种GNSS非差非组合未建模误差反演方法，该方法包括以下步骤：
步骤一，根据卫星坐标和地面接收机坐标计算卫星与地面之间的距离；
步骤二，通过观测文件获取卫星的载波观测值，通过钟差文件获取卫星的时钟误差；
步骤三，获取对流层延迟误差、电离层延迟误差、整周模糊度、接收机钟差以及与频率无关的各误差项的值；
步骤四，利用卫星与地面之间的距离以及步骤二和步骤三中获取的各项值通过相应的数学算法反演卫星定位过程中的未建模误差。
进一步的：在步骤一中，卫星与地面之间的距离根据卫星坐标和地面接收机坐标计算，其计算公式为：
式中，(xs，ys，zs)表示卫星坐标；(x0，y0，z0)表示地面接收机坐标。
进一步的，在步骤二中，所述的通过卫星产品文件获取卫星的载波观测值和卫星的时钟误差包括：
步骤2.1，从观测文件中获取卫星的载波观测值Φ；
步骤2.2，从钟差文件中获取卫星时钟误差dts。
进一步的，在步骤三中，获取对流层延迟误差、电离层延迟误差、整周模糊度、接收机钟差以及与频率无关的各误差项的值包括以下步骤：
步骤3.1，利用Saastamoinen模型计算对流层延迟T、计算与卫星信号频率无关的误差项dother；
步骤3.2，利用参数估计的方法获取电离层延迟I、实数解模糊度N，接收机钟差dtr；
步骤3.3，利用LAMBDA算法对模糊度N求取整周模糊度
进一步的，在步骤四中，利用卫星与地面之间的距离以及步骤二和步骤三中获取的各项值通过相应的数学算法反演未建模误差，其反演公式为：
式中，ε表示未建模误差；Φ表示载波观测值；ρ表示卫地距；c表示光速；dtr表示接收机钟差；dts表示卫星钟差；T表示对流层延迟；I表示电离层延迟；λ表示载波波长；表示整周模糊度；dother表示与卫星信号频率无关的误差项。
有益效果，与现有技术相比，本发明的技术方案具有以下有益技术效果：
本发明给出了GNSS非差非组合未建模误差的反演方法，通过精确的卫地距与参数估计值和模型改正值计算出了非差非组合模型中卫星的未建模误差。通过该方法计算出的未建模误差精度高，准确性好，能够为GNSS卫星导航中未建模误差的处理与研究提供理论基础与数据支持，方便未建模误差的进一步研究，进而能够为卫星导航定位中残余***误差的处理提供方法，有利于未建模误差的削弱甚至消除，进一步提高卫星定位的精度和可靠性。
附图说明
图1为本发明的方法流程图。
具体实施方式
以下结合实施例对本发明做具体说明。
本发明提出一种GNSS非差非组合未建模误差反演方法，如图1所示，具体包括以下步骤：
步骤一，根据卫星坐标和地面接收机坐标计算卫星与地面之间的距离；
步骤二，通过观测文件获取卫星的载波观测值，通过钟差文件获取卫星的时钟误差；
步骤三，获取对流层延迟误差、电离层延迟误差、整周模糊度、接收机钟差以及与频率无关的各误差项的值；
步骤四，利用卫星与地面之间的距离以及步骤二和步骤三中获取的各项值通过相应的数学算法反演卫星定位过程中的未建模误差。
优选的，在步骤一中，卫星与地面接收机的计算公式为
式中，(xs，ys，zs)表示卫星精密坐标；(x0，y0，z0)表示安装接收机的测站坐标，
在步骤二中，所述的获取载波观测值和卫星钟差包括以下步骤：
步骤2.1，从观测文件中获取载波观测值Φ；
步骤2.2，从钟差文件中获取卫星钟差dts；
在步骤三中，获取对流层延迟误差、电离层延迟误差、整周模糊度、接收机钟差以及与频率无关的各误差项的值包括以下步骤：
步骤3.1，利用模型改正的方法获取对流层延迟T、与频率无关的误差项dother；
步骤3.2，利用参数估计的方法获取电离层延迟I、实数解模糊度N，接收机钟差dtr；
步骤3.3，利用LAMBDA算法对模糊度N求取其整数解
在步骤四中，计算未建模误差的公式为：
式中，ε表示未建模误差；Φ表示载波观测值；ρ表示卫地距；c表示光速；dtr表示接收机钟差；dts表示卫星钟差；T表示对流层延迟；I表示电离层延迟；λ表示载波波长；表示整周模糊度；dother表示与卫星信号频率无关的误差项。
本发明给出了JFNG站2018年3月4日第12:00:00时刻G25卫星的未建模误差计算方法。
根据精密卫星坐标和精确接收机坐标计算卫地距真值
其中，G25卫星坐标
(xs，ys，zs)＝(-13610898.663541，20493181.8564741，9484399.604608)
JFNG测站接收机坐标
(x0，y0，z0)＝(-2279829.028319，5004706.478858，3219777.407958)
卫星与地面间的距离
从观测文件中获取载波观测值Φ
Φ＝23812042.00797
从钟差文件中获取卫星精密钟差dts
dts＝-2.7828×10-4
利用模型改正的方法获取对流层延迟T、一些与频率无关的误差项dother；
其中，对流层延迟
T＝7.30379
其余误差项
dother＝3531252.079807
利用参数估计的方法获取电离层延迟I、实数解模糊度N，接收机钟差dtr
其中，电离层延迟
I＝6.02820
接收机钟差
dtr＝-3.1148×10-4
模糊度
N＝45.717
利用LAMBDA算法对模糊度N求取其整周模糊度
载波观测值波长
λ＝0.1903m
利用以上参数计算未建模误差未建模误差

Claims (2)

1.一种GNSS非差非组合未建模误差反演方法，其特征在于，该方法包括以下步骤：
步骤一，根据卫星坐标和地面接收机坐标计算卫星与地面之间的距离；
步骤二，通过观测文件获取卫星的载波观测值，通过钟差文件获取卫星的时钟误差；
步骤三，获取对流层延迟误差、电离层延迟误差、整周模糊度、接收机钟差以及与频率无关的各误差项的值；
步骤四，利用卫星与地面之间的距离以及步骤二和步骤三中获取的各项值通过相应的数学算法反演卫星定位过程中的未建模误差；
在步骤一中，卫星与地面之间的距离根据卫星坐标和地面接收机坐标计算，其计算公式为：
式中，(xs,ys,zs)表示卫星坐标，(x0,y0,z0)表示地面接收机坐标；
在步骤三中，获取对流层延迟误差、电离层延迟误差、整周模糊度、接收机钟差以及与频率无关的各误差项的值包括以下步骤：
步骤3.1，利用Saastamoinen模型计算对流层延迟T、计算与卫星信号频率无关的误差项dother；
步骤3.2，利用参数估计的方法获取电离层延迟I、实数解模糊度N，接收机钟差dtr；
步骤3.3，利用LAMBDA算法对模糊度N求取整周模糊度
在步骤四中，利用卫星与地面之间的距离以及步骤二和步骤三中获取的各项值通过相应的数学算法反演未建模误差，其反演公式为：
式中，ε表示未建模误差；Φ表示载波观测值；ρ表示卫地距；c表示光速；dtr表示接收机钟差；dts表示卫星钟差；T表示对流层延迟；I表示电离层延迟；λ表示载波波长；表示整周模糊度；dother表示与卫星信号频率无关的误差项。
2.根据权利要求1所述的一种GNSS非差非组合未建模误差反演方法，其特征在于，在步骤二中，通过卫星产品文件获取卫星的载波观测值和卫星的时钟误差包括以下步骤：
步骤2.1，从观测文件中获取卫星的载波观测值Φ；
步骤2.2，从钟差文件中获取卫星时钟误差dts。
CN202110558650.8A 2021-05-21 2021-05-21 一种gnss未建模误差反演方法 Active CN113341437B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110558650.8A CN113341437B (zh) 2021-05-21 2021-05-21 一种gnss未建模误差反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110558650.8A CN113341437B (zh) 2021-05-21 2021-05-21 一种gnss未建模误差反演方法

Publications (2)

Publication Number Publication Date
CN113341437A CN113341437A (zh) 2021-09-03
CN113341437B true CN113341437B (zh) 2024-05-28

Family

ID=77470581

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110558650.8A Active CN113341437B (zh) 2021-05-21 2021-05-21 一种gnss未建模误差反演方法

Country Status (1)

Country Link
CN (1) CN113341437B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008082819A (ja) * 2006-09-27 2008-04-10 Sueo Sugimoto 測位装置および測位方法
CN105974440A (zh) * 2016-04-29 2016-09-28 和芯星通科技(北京)有限公司 一种获取观测噪声的方法和装置
CN107064980A (zh) * 2017-03-24 2017-08-18 和芯星通科技(北京)有限公司 载波相位模糊度固定方法和装置、卫星导航接收机
CN108363084A (zh) * 2018-01-18 2018-08-03 和芯星通科技(北京)有限公司 利用卫星定位的方法和装置、卫星导航接收机、存储介质
CN108828642A (zh) * 2018-08-01 2018-11-16 太原理工大学 一种ins辅助bds单频接收机的模糊度快速解算方法
CN108958018A (zh) * 2018-02-28 2018-12-07 和芯星通科技(北京)有限公司 一种卫星授时方法和装置、计算机可读存储介质
CN109709591A (zh) * 2018-12-07 2019-05-03 中国科学院光电研究院 一种面向智能终端的gnss高精度定位方法
KR101981225B1 (ko) * 2018-12-03 2019-05-22 한국해양과학기술원 Gnss 측정치에 포함된 특이신호요소 정밀추정 방법
CN110031881A (zh) * 2019-05-06 2019-07-19 中国人民解放军61540部队 高精度星间激光测距辅助精密单点定位的方法
CN110221316A (zh) * 2018-08-22 2019-09-10 宁波送变电建设有限公司永耀科技分公司 一种基于gps精密单点定位的移动手环定位方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10422885B2 (en) * 2016-03-18 2019-09-24 Deere & Company Rapid recovery of precise position after temporary signal loss
US10012738B2 (en) * 2016-03-24 2018-07-03 Hitachi Zosen Corporation Positioning method and positioning apparatus using satellite positioning system

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2008082819A (ja) * 2006-09-27 2008-04-10 Sueo Sugimoto 測位装置および測位方法
CN105974440A (zh) * 2016-04-29 2016-09-28 和芯星通科技(北京)有限公司 一种获取观测噪声的方法和装置
CN107064980A (zh) * 2017-03-24 2017-08-18 和芯星通科技(北京)有限公司 载波相位模糊度固定方法和装置、卫星导航接收机
CN108363084A (zh) * 2018-01-18 2018-08-03 和芯星通科技(北京)有限公司 利用卫星定位的方法和装置、卫星导航接收机、存储介质
CN108958018A (zh) * 2018-02-28 2018-12-07 和芯星通科技(北京)有限公司 一种卫星授时方法和装置、计算机可读存储介质
CN108828642A (zh) * 2018-08-01 2018-11-16 太原理工大学 一种ins辅助bds单频接收机的模糊度快速解算方法
CN110221316A (zh) * 2018-08-22 2019-09-10 宁波送变电建设有限公司永耀科技分公司 一种基于gps精密单点定位的移动手环定位方法
KR101981225B1 (ko) * 2018-12-03 2019-05-22 한국해양과학기술원 Gnss 측정치에 포함된 특이신호요소 정밀추정 방법
CN109709591A (zh) * 2018-12-07 2019-05-03 中国科学院光电研究院 一种面向智能终端的gnss高精度定位方法
CN110031881A (zh) * 2019-05-06 2019-07-19 中国人民解放军61540部队 高精度星间激光测距辅助精密单点定位的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
邹永廉.《土木工程测量》.高等教育出版社,2004,第243-245页. *

Also Published As

Publication number Publication date
CN113341437A (zh) 2021-09-03

Similar Documents

Publication Publication Date Title
Freda et al. Time-differenced carrier phases technique for precise GNSS velocity estimation
CN110045407B (zh) 一种分布式伪卫星/gnss优化定位方法
CN104502935B (zh) 一种基于非差非组合模型的网络rtk模糊度解算方法
CN101770033B (zh) 连续运行参考站***站间整周模糊度网络固定方法
CN105807300B (zh) 一种用北斗双频接收机进行高精度动态单点定位的方法
CN107942346B (zh) 一种高精度gnss电离层tec观测值提取方法
CN108196284B (zh) 一种进行星间单差模糊度固定的gnss网数据处理方法
CN104680008A (zh) 一种基于多参考站的网络rtk区域大气误差建模方法
CN102998681A (zh) 一种卫星导航***的高频钟差估计方法
CN108254773A (zh) 一种多gnss的实时钟差解算方法
CN116540303B (zh) 一种历元间差分观测方程建立方法和地震位移计算方法
CN112285745B (zh) 基于北斗三号卫星导航***的三频模糊度固定方法及***
CN109613582B (zh) 一种车载实时单频米级伪距定位方法
CN110749907A (zh) 一种基于北斗动定位中接收机的钟差补偿方法及其***
CN112526564A (zh) 一种精密单点定位重新收敛方法
Du et al. The realization and evaluation of PPP ambiguity resolution with INS aiding in marine survey
Ruan et al. Between-satellite single-difference integer ambiguity resolution in GPS/GNSS network solutions
Kuang et al. Real-time GPS satellite orbit and clock estimation based on OpenMP
CN115421172A (zh) 一种基于实时与准实时结合的北斗变形监测方法
CN105403901A (zh) 一种dgnss卫星轨道偏差修正方法
CN113341437B (zh) 一种gnss未建模误差反演方法
CN116299586B (zh) 基于广播星历的精密单点定位方法、接收机、设备和介质
CN115902968A (zh) 基于北斗三号geo播发增强信息的ppp终端定位方法
CN114527500B (zh) 一种室内外一体化定位方法、设备、介质及产品
CN115561793A (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
GR01 Patent grant
GR01 Patent grant