CN115166799B - 顾及硬件延迟时变特征的gnss精密单点定位方法 - Google Patents
顾及硬件延迟时变特征的gnss精密单点定位方法 Download PDFInfo
- Publication number
- CN115166799B CN115166799B CN202210811562.9A CN202210811562A CN115166799B CN 115166799 B CN115166799 B CN 115166799B CN 202210811562 A CN202210811562 A CN 202210811562A CN 115166799 B CN115166799 B CN 115166799B
- Authority
- CN
- China
- Prior art keywords
- hardware delay
- gnss
- receiver
- satellite
- frequency
- 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
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/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
-
- 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/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/29—Acquisition or tracking or demodulation of signals transmitted by the system carrier including Doppler, related
-
- 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
-
- 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/396—Determining accuracy or reliability of position or pseudorange measurements
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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精密单点定位方法,利用GNSS单系统双频伪距和载波观测数据,和精密轨道、钟差产品,建立非差非组合观测方程,并考虑接收机硬件延迟随时间变化的特征,对接收机硬件延迟、电离层和模糊度等参数重整为新的可估参数,构建新的PPP函数模型,再进行解算,从而消除接收机硬件延迟对PPP结果的影响。
Description
技术领域
本发明涉及精确定位领域,尤其涉及一种顾及硬件延迟时变特征的GNSS精密单点定位方法。
背景技术
基于全球导航卫星系统(GNSS)实现的精密单点定位(PPP)技术能够提供包括测站位置、接收机钟差、天顶对流层延迟和倾斜电离层延迟等参数,在定位、导航和授时应用中具有重要意义。传统PPP的实施中,通常根据伪距硬件延迟的特性,将观测方程中接收机伪距硬件延迟当作不随时间变化的参数处理。然而,近年来诸多学者对接收机硬件延迟的稳定性开展相关研究发现,环境温度与其稳定性之间存在一定相关性。如果在GNSS数据处理过程中,将接收机硬件延迟当做时不变参数处理,势必会影响各类参数估值的准确性。
发明内容
本发明主要目的在于:提供一种顾及硬件延迟时变特征的GNSS精密单点定位方法,能够有效消除接收机硬件延迟对PPP结果的影响。
本发明所采用的技术方案是:顾及硬件延迟时变特征的GNSS精密单点定位方法,本方法包括以下步骤:
S1、选用GNSS单系统单测站的伪距和相位双频观测数据,构建单系统双频非差非组合原始观测方程为:
其中r表示接收机,s=1...m,表示卫星号,j=1...f表示频率号,i=1...t表示历元号,和/>分别表示伪距和相位观测值,单位为米,其所含的站星距初值、卫星钟差和对流层干延迟已被改正,/>为测站到卫星的单位向量,/>为坐标的改正值向量,单位为米,τr(i)表示对流层天顶湿延迟,单位为米,/>表示湿延迟投影系数,dtr(i)表示接收机钟差,单位为米,dr,j和/>表示频率j对应的接收机和卫星伪距硬件延迟,单位为米,表示卫星硬件延迟的消电离层组合,单位为米,/>表示第1频率的电离层延迟,单位为米,其系数/>为其他频率与第1频率电离层延迟间的比值,λj为对应频率j的波长,/>是非整周模糊度,单位为米;
S2、将所述单系统双频非差非组合原始观测方程中的接收机硬件延迟参数dr,j重新设定为时变参数dr,j(i),并将其表示为第1历元接收机硬件延迟dr,1(1)与历元间接收机硬件延迟变化值的和,即/>
S3、结合S2的设定,对S1中的参数重新整合以消除秩亏,构建新的观测方程为:
其中
S4、利用滤波或多历元解的方法,对S3中的观测方程进行求解,获得后,将其改正到近似坐标上,得到精确的坐标解。
按上述方法,S1选用的数据包括北斗二号、北斗三号、GPS、GLONASS或Galileo系统,并同时选用精密卫星轨道和钟差产品。
按上述方法,所述的S3新的观测方程基础上,若已知测站的精确坐标,则不存在新的观测方程变成:
仍然是满秩方程,用于解算对流层延迟、接收机钟差和电离层延迟,进而用于分析以上参数估值的性能。
一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现所述顾及硬件延迟时变特征的GNSS精密单点定位方法的步骤。
一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现所述顾及硬件延迟时变特征的GNSS精密单点定位方法的步骤。
本发明产生的有益效果是:利用GNSS单系统双频伪距和载波观测数据,和精密轨道、钟差产品,建立非差非组合观测方程,并考虑接收机硬件延迟随时间变化的特征,对接收机硬件延迟、电离层和模糊度等参数重整为新的可估参数,构建新的PPP函数模型,再进行解算,从而消除接收机硬件延迟对PPP结果的影响。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本发明提供一种顾及硬件延迟时变特征的GNSS精密单点定位方法,本方法包括以下步骤:
S1、选用GNSS单系统单测站的伪距和相位双频观测数据,包括北斗二号、北斗三号、GPS、GLONASS或Galileo系统,并同时选用精密卫星轨道和钟差产品。
构建单系统双频非差非组合原始观测方程为:
其中r表示接收机,s=1...m,表示卫星号,j=1...f表示频率号,i=1...t表示历元号,和/>分别表示伪距和相位观测值,单位为米,其所含的站星距初值、卫星钟差和对流层干延迟已被改正,/>为测站到卫星的单位向量,/>为坐标的改正值向量,单位为米,τr(i)表示对流层天顶湿延迟,单位为米,/>表示湿延迟投影系数,dtr(i)表示接收机钟差,单位为米,dr,j和/>表示频率j对应的接收机和卫星伪距硬件延迟,单位为米,表示卫星硬件延迟的消电离层组合,单位为米,/>表示第1频率的电离层延迟,单位为米,其系数/>为其他频率与第1频率电离层延迟间的比值,λj为对应频率j的波长,/>是非整周模糊度,单位为米。
S2、将所述单系统双频非差非组合原始观测方程中的接收机硬件延迟参数dr,j重新设定为时变参数dr,j(i),并将其表示为第1历元接收机硬件延迟dr,1(1)与历元间接收机硬件延迟变化值的和,即/>
S3、结合S2的设定,对S1中的参数重新整合以消除秩亏,构建新的观测方程为:
其中
S4、利用滤波或多历元解的方法,对S3中的观测方程进行求解,获得后,将其改正到近似坐标上,得到精确的坐标解。
所述的S3新的观测方程基础上,若已知测站的精确坐标,则不存在新的观测方程变成:
该方程仍然是满秩方程,可用于解算对流层延迟、接收机钟差和电离层延迟,进而用于分析以上参数估值的性能。
本发明还提供一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现所述顾及硬件延迟时变特征的GNSS精密单点定位方法的步骤。
本发明还提供一种计算机可读存储介质,其上存储有计算机程序,所述计算机程序被处理器执行时实现所述顾及硬件延迟时变特征的GNSS精密单点定位方法的步骤。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (5)
1.顾及硬件延迟时变特征的GNSS精密单点定位方法,其特征在于,本方法包括以下步骤:
S1、选用GNSS单系统单测站的伪距和相位双频观测数据,构建单系统双频非差非组合原始观测方程为:
其中r表示接收机,s=1...m,表示卫星号,j=1...f表示频率号,i=1...t表示历元号,和/>分别表示伪距和相位观测值,单位为米,其所含的站星距初值、卫星钟差和对流层干延迟已被改正,/>为测站到卫星的单位向量,/>为坐标的改正值向量,单位为米,τr(i)表示对流层天顶湿延迟,单位为米,/>表示湿延迟投影系数,dtr(i)表示接收机钟差,单位为米,dr,j和/>表示频率j对应的接收机和卫星伪距硬件延迟,单位为米,表示卫星硬件延迟的消电离层组合,单位为米,/>表示第1频率的电离层延迟,单位为米,其系数/>为其他频率与第1频率电离层延迟间的比值,λj为对应频率j的波长,/>是非整周模糊度,单位为米;
S2、将所述单系统双频非差非组合原始观测方程中的接收机硬件延迟参数dr,j重新设定为时变参数dr,j(i),并将其表示为第1历元接收机硬件延迟dr,1(1)与历元间接收机硬件延迟变化值的和,即/>
S3、结合S2的设定,对S1中的参数重新整合以消除秩亏,构建新的观测方程为:
其中
S4、利用滤波或多历元解的方法,对S3中的观测方程进行求解,获得后,将其改正到近似坐标上,得到精确的坐标解。
2.根据权利要求1所述的顾及硬件延迟时变特征的GNSS精密单点定位方法,其特征在于,S1选用的数据包括北斗二号、北斗三号、GPS、GLONASS或Galileo系统,并同时选用精密卫星轨道和钟差产品。
3.根据权利要求1所述的顾及硬件延迟时变特征的GNSS精密单点定位方法,其特征在于,所述的S3新的观测方程基础上,若已知测站的精确坐标,则不存在新的观测方程变成:
仍然是满秩方程,用于解算对流层延迟、接收机钟差和电离层延迟,进而用于分析以上参数估值的性能。
4.一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现权利要求1至3中任意一项所述顾及硬件延迟时变特征的GNSS精密单点定位方法的步骤。
5.一种计算机可读存储介质,其上存储有计算机程序,其特征在于:所述计算机程序被处理器执行时实现权利要求1至3中任意一项所述顾及硬件延迟时变特征的GNSS精密单点定位方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210811562.9A CN115166799B (zh) | 2022-07-11 | 2022-07-11 | 顾及硬件延迟时变特征的gnss精密单点定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210811562.9A CN115166799B (zh) | 2022-07-11 | 2022-07-11 | 顾及硬件延迟时变特征的gnss精密单点定位方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115166799A CN115166799A (zh) | 2022-10-11 |
CN115166799B true CN115166799B (zh) | 2023-07-21 |
Family
ID=83492619
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210811562.9A Active CN115166799B (zh) | 2022-07-11 | 2022-07-11 | 顾及硬件延迟时变特征的gnss精密单点定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115166799B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110412633A (zh) * | 2019-07-30 | 2019-11-05 | 广州市中海达测绘仪器有限公司 | 定位方法、装置、计算机设备和存储介质 |
CN114019550A (zh) * | 2021-10-25 | 2022-02-08 | 中国电子科技集团公司第五十四研究所 | 模糊度加权约束的gnss非差非组合观测方程解算方法 |
CN114488245A (zh) * | 2022-01-25 | 2022-05-13 | 同济大学 | 一种系统间共频信号差分的gnss/ins组合定位方法 |
Family Cites Families (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU725036A1 (ru) * | 1978-09-11 | 1980-03-30 | Предприятие П/Я В-2203 | Устройство дл прив зки шкал времени и сличени частот |
US6888498B2 (en) * | 2002-09-24 | 2005-05-03 | Honeywell International, Inc. | Method and system for compensating satellite signals |
CN101887128B (zh) * | 2010-07-09 | 2012-11-14 | 中国科学院测量与地球物理研究所 | 确定全球卫星导航系统导航卫星频间偏差的方法 |
CN101964689A (zh) * | 2010-08-20 | 2011-02-02 | 惠州中城电子科技有限公司 | 一种提高电能信息采集系统时间精度的方法 |
CN104407359B (zh) * | 2014-12-10 | 2017-03-22 | 中南大学 | 一种北斗接收机非差观测值噪声评估方法 |
CN105182388B (zh) * | 2015-10-10 | 2017-08-25 | 安徽理工大学 | 一种快速收敛的精密单点定位方法 |
CN105807300B (zh) * | 2016-03-17 | 2018-09-07 | 武汉际上导航科技有限公司 | 一种用北斗双频接收机进行高精度动态单点定位的方法 |
CN107132558B (zh) * | 2017-06-13 | 2019-07-23 | 武汉大学 | 惯性辅助的多频多模gnss周跳修复方法及系统 |
CN108415050B (zh) * | 2018-06-04 | 2020-05-26 | 北京未来导航科技有限公司 | 一种基于低轨星座导航增强系统的ppp-rtk定位方法 |
CN109765589B (zh) * | 2019-02-21 | 2022-12-02 | 哈尔滨工程大学 | 一种基于无电离层组合的三频gnss实时周跳固定技术 |
CN113703017B (zh) * | 2020-05-22 | 2023-06-13 | 中国人民解放军战略支援部队信息工程大学 | 一种卫星天线相位中心偏差计算方法及装置 |
CN111796309B (zh) * | 2020-06-24 | 2023-04-18 | 中国科学院精密测量科学与技术创新研究院 | 导航卫星单频数据同步确定大气水汽和总电子含量的方法 |
CN113109840A (zh) * | 2021-03-19 | 2021-07-13 | 中国电波传播研究所(中国电子科技集团公司第二十二研究所) | 一种基于gnss接收机的电离层tec实时测量方法 |
CN113504557B (zh) * | 2021-06-22 | 2023-05-23 | 北京建筑大学 | 面向实时应用的gps频间钟差新预报方法 |
CN113933868B (zh) * | 2021-12-16 | 2022-03-22 | 中南大学 | 一种北斗二号meo卫星频率间卫星钟偏差的建模方法 |
-
2022
- 2022-07-11 CN CN202210811562.9A patent/CN115166799B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110412633A (zh) * | 2019-07-30 | 2019-11-05 | 广州市中海达测绘仪器有限公司 | 定位方法、装置、计算机设备和存储介质 |
CN114019550A (zh) * | 2021-10-25 | 2022-02-08 | 中国电子科技集团公司第五十四研究所 | 模糊度加权约束的gnss非差非组合观测方程解算方法 |
CN114488245A (zh) * | 2022-01-25 | 2022-05-13 | 同济大学 | 一种系统间共频信号差分的gnss/ins组合定位方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115166799A (zh) | 2022-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108931915B (zh) | 利用导航卫星的授时方法和装置、计算机可读存储介质 | |
US11709281B2 (en) | High-precision point positioning method and device based on smartphone | |
CN108828640B (zh) | 一种卫星导航定位观测值定权方法及装置 | |
CN107765275B (zh) | 广域差分定位方法、装置、终端及计算机可读存储介质 | |
CN108958018B (zh) | 一种卫星授时方法和装置、计算机可读存储介质 | |
CN109683186B (zh) | 一种消除多卫星导航系统载波相位时间传递天跳变的方法 | |
CN114966760B (zh) | 一种电离层加权的非差非组合ppp-rtk技术实现方法 | |
CN112230252B (zh) | 终端定位方法、装置、计算机设备和存储介质 | |
CN104898145B (zh) | 一种基于半周模糊度的模糊度固定方法和系统 | |
CN114755702B (zh) | 基于ppp的精密卫星星历时间基准稳定度评估方法及系统 | |
CN116953741B (zh) | 一种应用于全球导航卫星系统gnss的周跳探测与修复方法 | |
Liu et al. | An efficient undifferenced method for estimating multi-GNSS high-rate clock corrections with data streams in real time | |
CN114994727B (zh) | 一种实现高精度时间校准及卫星定位的设备 | |
CN114879239B (zh) | 一种增强瞬时ppp固定解的区域三频整数钟差估计方法 | |
CN115480279A (zh) | Gnss导航方法和终端、组合导航系统、存储介质 | |
CN108254774A (zh) | 基于gnss多频信号的单基站长距离实时定位方法 | |
CN115079225A (zh) | 一种海上接收机的导航定位方法与装置 | |
CN117289318B (zh) | 一种基于bds辅助的glonass频间偏差实时估计方法 | |
CN115166799B (zh) | 顾及硬件延迟时变特征的gnss精密单点定位方法 | |
CN117130026A (zh) | 一种gnss钟差解算方法 | |
CN115308781B (zh) | 基于bdgim辅助的相位平滑伪距高精度时间传递方法 | |
CN114002722A (zh) | Gnss定位方法及其装置 | |
JP4928114B2 (ja) | キャリア位相相対測位装置 | |
CN105699997B (zh) | 一种使用glonass单频信号进行差分定位的方法 | |
JP2022097435A (ja) | 基準受信機を用いたロービング受信機の時間非依存位置決定 |
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 |