CN115856966A - 用于实现多gnss多频率非组合精密单点定位的方法 - Google Patents
用于实现多gnss多频率非组合精密单点定位的方法 Download PDFInfo
- Publication number
- CN115856966A CN115856966A CN202210838316.2A CN202210838316A CN115856966A CN 115856966 A CN115856966 A CN 115856966A CN 202210838316 A CN202210838316 A CN 202210838316A CN 115856966 A CN115856966 A CN 115856966A
- Authority
- CN
- China
- Prior art keywords
- frequency
- observation
- bds
- satellite
- pseudo
- 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
- 230000000694 effects Effects 0.000 claims abstract description 6
- 238000012886 linear function Methods 0.000 claims abstract description 5
- 238000001914 filtration Methods 0.000 claims abstract description 4
- 238000012937 correction Methods 0.000 claims description 16
- 239000013598 vector Substances 0.000 claims description 15
- 230000006870 function Effects 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 230000007812 deficiency Effects 0.000 claims description 6
- 230000014509 gene expression Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 3
- 238000000605 extraction Methods 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 3
- 238000002242 deionisation method Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000008521 reorganization Effects 0.000 description 1
- 230000003313 weakening effect Effects 0.000 description 1
Images
Classifications
-
- 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
Abstract
本发明提出一种基于多GNSS系统多频率观测数据实现非组合精密单点定位(PPP)的方法;该方法在单系统双频非组合PPP模型的基础上,提出一种联合处理多系统、多频率观测数据的模型和方法,在该模型中顾及BDS系统GEO卫星伪距多路径效应影响和GLONASS系统频分多址技术的特点,分别引入伪距多路径参数和卫星频率号的线性函数模型来精化BDS和GLONASS模型中的伪距观测信息,构建满秩灵活的兼容GPS、GLONASS、BDS、Galileo四系统的多频率联合处理的非组合PPP模型,利用滤波或最小二乘方法进行参数估计,实现精密定位、授时和大气延迟提取等各类应用。
Description
技术领域
本发明涉及高精度定位系统,尤其涉及一种用于实现多GNSS多频率非组合精密单点定位的方法。
背景技术
精密单点定位(precise point positioning,PPP)技术是GNSS领域内一种高精度的绝对定位技术。根据对电离层延迟的处理方式不同,现有PPP技术的函数模型主要有“消电离层组合”和“非组合”两种模型。非组合模型是近年来提出的一种新的PPP函数模型,其相对于传统消电离层组合模型具有理论上的最优性、实施上的灵活性、应用上的广泛性等优势,相关学者已对非组合PPP模型与算法、定位与授时性能及其在大气延迟精密提取与码偏差估计等方面展开了研究,并取得了丰硕的研究成果。近年来,随着世界各国、组织建立的多个GNSS系统日趋完善,为用户端的定位、导航和授时应用提供了越来越丰富的数据支持,因而多系统多频率融合的精密单点定位技术成为当前的研究热点。
发明内容
本发明主要目的在于:提供一种用于实现多GNSS多频率非组合精密单点定位的方法,不仅能够联合处理多个GNSS系统,而且能够处理双频及以上任意频率的原始观测数据。
本发明所采用的技术方案是:一种用于实现多GNSS多频率非组合精密单点定位的方法,本方法包括以下步骤:
S1、选用多GNSS系统单测站的伪距和相位观测数据,构建单系统多频率的原始观测方程;
S2、基于所述的单系统多频率的原始观测方程,拓展为GPS、Galileo和BDS系统多频率的观测方程;
S3、在所述的GPS、Galileo和BDS系统多频率的观测方程中,引入各统的卫星钟差作为已知信息,然后对各项参数重新整合以消除观测方程存在的秩亏后,得到满秩的GPS与BDS和Galileo三系统联合的多频观测方程;
S4、在S3所述的GPS与BDS和Galileo三系统联合的多频观测方程的基础上,对BDS系统GEO卫星的伪距观测值中增加多路径偏差参数,得到顾及BDS系统GEO卫星伪距多路径效应的GNSS多频多系统观测方程;
S5、顾及GLONASS系统采用频分多址技术,引入卫星频率号的线性函数来模型化GLONASS接收机端伪距硬件延迟,并参照S1构建GLONASS系统新的伪距观测方程;
S6、基于S1中单系统多频率的原始观测方程,和S5得到的GLONASS系统新的伪距观测方程,引入GLONASS系统的卫星钟差作为已知信息,然后对各项参数重新整合以消除单系统多频率的观测方程存在的秩亏问题后,得到顾及伪距频间偏差IFB(inter-frequencybias,IFB)影响的满秩GLONASS非组合精密单点定位模型;
S7、联合所述的S3、S4和S6的结果,组成多系统满秩的灵活兼容多频观测信息的非组合精密单点定位模型;利用滤波或最小二乘方法进行参数估计,对所述的多系统满秩的灵活兼容多频观测信息的非组合精密单点定位模型进行求解,获得测站坐标的改正值向量并改正到其近似坐标上,得到精确的测站坐标解。
按上述方法,S1中的多GNSS系统单测站的伪距和相位观测数据,包括BDS2系统的双频或三频,BDS3系统的双频、三频、四频或五频,GPS系统的双频或三频,GLONASS系统的双频,Galileo系统的双频、三频、四频或五频观测数据。
按上述方法,所述S1构建的单系统多频率的原始观测方程为:
其中,和/>分别表示伪距和相位观测值,单位为米,其所含的站星距初值、对流层干延迟及其他精密改正已被事先改正,单位为米,S表示卫星号,r表示接收机号,f表示频段号,/>为测站到卫星的单位向量,/>为坐标的改正值向量,c表示光速,Δtr和ΔtS是接收机和卫星钟差,单位为秒,br,f和/>表示接收机r和卫星S对应频率f的伪距硬件延迟,单位为秒,/>表示电离层斜延迟,单位为米,其系数/>为其他频率与第1频率电离层延迟间的比值,λf为对应频率f的波长,Tr为天顶对流层湿延迟,单位为米,其投影函数为为接收机r与卫星S之间频率f对应的吸收了相位偏差的浮点模糊度,单位为周,观测噪声和其他未模型化误差视为白噪声并在观测方程中省略。
按上述方法,所述S2中,GPS、Galileo和BDS系统多频率的观测方程具体为:
式中,上标G、E和C分别为GPS、Galileo和BDS的系统标识。
按上述方法,S3所述的满秩的GPS与BDS和Galileo三系统联合的多频观测方程为:
其中,和/>表示经卫星钟差改正后的新观测值,/> 为对应GPS系统的接收机钟差,/>为消电离层组合系数项,其中上标T表示G、E、C系统之一,/>为系统间偏差,/>为经参数重组的电离层延迟,/>为第1和第2频点之间的接收机差分码偏差,第1和第2频点之间的卫星差分码偏差,为IFB信息,为参数重组后的模糊度信息,系数/>
按上述方法,S4所述的顾及伪距多路径影响的GNSS多频多系统观测方程中BDS伪距观测量为:
按上述方法,所述的S5构建的GLONASS系统新的伪距观测方程具体为:
按上述方法,所述的S6得到的顾及伪距频间偏差IFB影响的满秩GLONASS非组合精密单点定位模型具体为:
按上述方法,所述的S7组成的多系统满秩的灵活兼容多频观测信息的非组合精密单点定位模型具体为:
一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述用于实现多GNSS多频率非组合精密单点定位的方法的步骤。
本发明产生的有益效果是:能够联合处理多GNSS系统和双频及以上任意多频率的观测数据,能够解算第三及以上任意频率的伪距硬件延迟,其为不同频率卫星与接收机伪距硬件延迟的线性组合,可进一步用来估计多频码偏差参数;针对BDS的GEO卫星伪距观测数据的多路径问题,增加了用于削弱其影响的待估参数,针对GLONASS卫星频分多址的特征,增加了卫星频率号的线性函数来模型化接收机端的伪距硬件延迟,从而精化了函数模型并提高了其应用性能。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明实施例的方法流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
如图1所示,本发明提供一种用于实现多GNSS多频率非组合精密单点定位的方法,包括以下步骤:
S1、选用多GNSS系统单测站的伪距和相位观测数据,构建单系统多频率的原始观测方程。多GNSS系统单测站的伪距和相位观测数据,包括BDS2系统的双频或三频,BDS3系统的双频、三频、四频或五频,GPS系统的双频或三频,GLONASS系统双频,Galileo系统的双频、三频、四频或五频观测数据。
构建的单系统多频率的原始观测方程为:
其中,和/>分别表示伪距和相位观测值,单位为米,其所含的站星距初值、对流层干延迟及其他精密改正已被事先改正,单位为米,S表示卫星号,r表示接收机号,f表示频段号,/>为测站到卫星的单位向量,/>为坐标的改正值向量,c表示光速,Δtr和ΔtS是接收机和卫星钟差,单位为秒,br,f和/>表示接收机r和卫星S对应频率f的伪距硬件延迟,单位为秒,/>表示电离层斜延迟,单位为米,其系数/>为其他频率与第1频率电离层延迟间的比值,λf为对应频率f的波长,Tr为天顶对流层湿延迟,单位为米,其投影函数为为接收机r与卫星S之间频率f对应的吸收了相位偏差的浮点模糊度,单位为周,观测噪声和其他未模型化误差视为白噪声并在观测方程中省略。
S2、基于所述的单系统多频率的原始观测方程,拓展为GPS、Galileo和BDS系统多频率的观测方程,具体为:
式中,上标G、E和C分别为GPS、Galileo和BDS的系统标识。
S3、在所述的GPS、Galileo和BDS系统多频率的观测方程中,引入各统的卫星钟差作为已知信息,然后对各项参数重新整合以消除观测方程存在的秩亏问题后,得到满秩的GPS、Galileo、以及BDS的G/E/C多频观测方程,即GPS与BDS和Galileo三系统联合的多频观测方程:,其中,/>和/>表示经卫星钟差改正后的新观测值,/>为对应GPS系统的接收机钟差,为消电离层组合系数项,其中上标T表示G、E、C系统之一,/>为系统间偏差,/>为经参数重组的电离层延迟,/>为第1和第2频点之间的接收机差分码偏差,为第1和第2频点之间的卫星差分码偏差,为IFB信息,为参数重组后的模糊度信息,系数/>
S4、在S3所述的GPS与BDS和Galileo三系统联合的多频观测方程基础上,对BDS的GEO卫星的伪距观测值中增加多路径偏差参数,得到顾及BDS系统GEO卫星伪距多路径效应的BDS伪距观测量为:
S5、顾及GLONASS系统采用频分多址技术,引入卫星频率号的线性函数来模型化GLONASS接收机端伪距硬件延迟,并参照S1构建GLONASS系统新的伪距观测方程:
S6、基于S1中单系统多频率的原始观测方程,和S5得到的GLONASS系统新的伪距观测方程,引入GLONASS系统的卫星钟差作为已知信息,然后对各项参数重新整合以消除单系统多频率的观测方程存在的秩亏问题后,得到顾及伪距频间偏差IFB影响的满秩GLONASS非组合精密单点定位模型为:
S7、联合S3、S4和S6的结果,得到多系统满秩的灵活兼容多频观测信息的非组合精密单点定位模型:
利用滤波或最小二乘方法进行参数估计,对所述的多系统满秩的灵活兼容多频观测信息的非组合精密单点定位模型进行求解,获得测站坐标的改正值向量并改正到其近似坐标上,得到精确的坐标解。
本发明还提供一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述用于实现多GNSS多频率非组合精密单点定位的方法的步骤。
本发明在传统PPP技术的基础上,提出了一种适用于多个GNSS系统、多频率观测数据联合实现非组合PPP的方法,在联合多个GNSS系统实现PPP时,顾及BDS的GEO卫星伪距多路径效应,同时考虑GLONASS系统频分多址(FDMA)特征,不仅能够联合处理多个GNSS系统,而且能够处理双频及以上任意频率的非组合观测数据。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (10)
1.一种用于实现多GNSS多频率非组合精密单点定位的方法,其特征在于,本方法包括以下步骤:
S1、选用多GNSS系统单测站的伪距和相位观测数据,构建单系统多频率的原始观测方程;
S2、基于所述的单系统多频率的原始观测方程,拓展为GPS、Galileo和BDS系统多频率的观测方程;
S3、在所述的GPS、Galileo和BDS系统多频率的观测方程中,引入各统的卫星钟差作为已知信息,然后对各项参数重新整合以消除观测方程存在的秩亏后,得到满秩的GPS与BDS和Galileo三系统联合的多频观测方程;
S4、在S3所述的GPS与BDS和Galileo三系统联合的多频观测方程基础上,对BDS系统的GEO卫星的伪距观测值中增加多路径偏差参数,得到顾及BDS系统GEO卫星伪距多路径效应的GNSS多频多系统观测方程;
S5、顾及GLONASS系统采用频分多址技术,引入卫星频率号的线性函数来模型化GLONASS接收机端伪距硬件延迟,并参照S1构建GLONASS系统新的伪距观测方程;
S6、基于S1中单系统多频率的原始观测方程,和S5得到的GLONASS系统新的伪距观测方程,引入GLONASS系统的卫星钟差作为已知信息,然后对各项参数重新整合以消除单系统多频率的观测方程存在的秩亏后,得到顾及伪距频间偏差IFB影响的满秩GLONASS非组合精密单点定位模型;
S7、联合所述的S3、S4和S6的结果,组成多系统满秩的灵活兼容多频观测信息的非组合精密单点定位模型;利用滤波或最小二乘方法进行参数估计,对所述的多系统满秩的灵活兼容多频观测信息的非组合精密单点定位模型进行求解,获得测站坐标的改正值向量并改正到其近似坐标上,得到精确的测站坐标解。
2.根据权利要求1所述的用于实现多GNSS多频率非组合精密单点定位的方法,其特征在于,S1中的多GNSS系统单测站的伪距和相位观测数据,包括BDS2系统的双频或三频,BDS3系统的双频、三频、四频或五频,GPS系统的双频或三频,GLONASS系统的双频,Galileo系统的双频、三频、四频或五频观测数据。
3.根据权利要求1所述的用于实现多GNSS多频率非组合精密单点定位的方法,其特征在于,所述S1构建的单系统多频率的原始观测方程为:
其中,和/>分别表示伪距和相位观测值,单位为米,其所含的站星距初值、对流层干延迟及其他精密改正已被事先改正,单位为米,S表示卫星号,f表示接收机号,f表示频段号,/>为测站到卫星的单位向量,/>为坐标的改正值向量,c表示光速,Δtr和ΔtS是接收机和卫星钟差,单位为秒,br,f和/>表示接收机r和卫星S对应频率f的伪距硬件延迟,单位为秒,/>表示电离层斜延迟,单位为米,其系数/>为其他频率f与第1频率电离层延迟间的比值,λf为对应频率f的波长,Tr为天顶对流层湿延迟,单位为米,其投影函数为为接收机r与卫星S之间频率f对应的吸收了相位偏差的浮点模糊度,单位为周,观测噪声和其他未模型化误差视为白噪声并在观测方程中省略。
10.一种计算机设备,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,其特征在于:所述处理器执行所述计算机程序时实现上述权利要求1至9中任意一项所述用于实现多GNSS多频率非组合精密单点定位的方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210838316.2A CN115856966B (zh) | 2022-07-12 | 2022-07-12 | 用于实现多gnss多频率非组合精密单点定位的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210838316.2A CN115856966B (zh) | 2022-07-12 | 2022-07-12 | 用于实现多gnss多频率非组合精密单点定位的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115856966A true CN115856966A (zh) | 2023-03-28 |
CN115856966B CN115856966B (zh) | 2023-08-29 |
Family
ID=85660333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210838316.2A Active CN115856966B (zh) | 2022-07-12 | 2022-07-12 | 用于实现多gnss多频率非组合精密单点定位的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115856966B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150293233A1 (en) * | 2012-10-25 | 2015-10-15 | Fugron N.V. | Ppp-rtk method and system for gnss signal based position determination |
WO2017181221A1 (en) * | 2016-04-18 | 2017-10-26 | Spatial Information Systems Research Ltd | A method of analysing multiple signals transmitted by discrete global navigation satellite systems |
CN107942356A (zh) * | 2017-11-09 | 2018-04-20 | 武汉大学 | 一种多频多模gnss广义绝对码偏差估计方法 |
CN111190203A (zh) * | 2020-01-14 | 2020-05-22 | 中国科学院测量与地球物理研究所 | 基于非差观测模型和消秩亏的ppp-rtk定位方法 |
CN114019550A (zh) * | 2021-10-25 | 2022-02-08 | 中国电子科技集团公司第五十四研究所 | 模糊度加权约束的gnss非差非组合观测方程解算方法 |
CN114019551A (zh) * | 2021-10-26 | 2022-02-08 | 中国电子科技集团公司第五十四研究所 | 一种gnss观测站网原始观测方程解算方法 |
-
2022
- 2022-07-12 CN CN202210838316.2A patent/CN115856966B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150293233A1 (en) * | 2012-10-25 | 2015-10-15 | Fugron N.V. | Ppp-rtk method and system for gnss signal based position determination |
WO2017181221A1 (en) * | 2016-04-18 | 2017-10-26 | Spatial Information Systems Research Ltd | A method of analysing multiple signals transmitted by discrete global navigation satellite systems |
CN107942356A (zh) * | 2017-11-09 | 2018-04-20 | 武汉大学 | 一种多频多模gnss广义绝对码偏差估计方法 |
CN111190203A (zh) * | 2020-01-14 | 2020-05-22 | 中国科学院测量与地球物理研究所 | 基于非差观测模型和消秩亏的ppp-rtk定位方法 |
CN114019550A (zh) * | 2021-10-25 | 2022-02-08 | 中国电子科技集团公司第五十四研究所 | 模糊度加权约束的gnss非差非组合观测方程解算方法 |
CN114019551A (zh) * | 2021-10-26 | 2022-02-08 | 中国电子科技集团公司第五十四研究所 | 一种gnss观测站网原始观测方程解算方法 |
Non-Patent Citations (4)
Title |
---|
XINGXING LI等: "Multi-GNSS phasse delay estimation and PPP ambiguity resolution:GPS BDS,GLONASS,Galileo", JOURNAL OF GEODESY, pages 579 - 608 * |
ZHENG ZHANG等: "On GLONASS pseudo-range inter-frequency bias solution with ionospheric delay modeling and the undifferenced uncombined PPP", JOURNAL OF GEODESY * |
侯鹏宇等: "CDMA+FDMA非差非组合区域PPP-RTK", 测绘学报, vol. 52, no. 2 * |
刘腾等: "BDS/GLONASS非组合精密单点定位模型与算法", 地球物理学报, vol. 60, no. 4, pages 1347 - 1356 * |
Also Published As
Publication number | Publication date |
---|---|
CN115856966B (zh) | 2023-08-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2749667C1 (ru) | Способ и система быстрого и точного позиционирования | |
CN108415049B (zh) | 提高网络rtk双差宽巷模糊度固定正确率的方法 | |
CN111190203B (zh) | 基于非差观测模型和消秩亏的ppp-rtk定位方法 | |
US10078140B2 (en) | Navigation satellite system positioning involving the generation of advanced correction information | |
CN110007320B (zh) | 网络rtk解算方法 | |
JP6023225B2 (ja) | 無線ナビゲーション信号を処理する方法 | |
CN111273327B (zh) | 一种基于组合和非组合混合观测模型的精密单点定位方法 | |
CN102353969A (zh) | 精密单点定位技术中相位偏差的估计方法 | |
CN114966760B (zh) | 一种电离层加权的非差非组合ppp-rtk技术实现方法 | |
CN115201870B (zh) | 附先验约束的多频多模gnss非差非组合时间传递方法 | |
CN113848572A (zh) | 一种基于大气误差增强的多频ppp序贯单历元定位方法 | |
CN110824505B (zh) | Gnss卫星接收机的偏差估计方法及系统、定位方法及终端 | |
CN111123295B (zh) | 基于ssr的定位方法及装置、定位系统 | |
CN115267843A (zh) | 多频多模gnss高频精密卫星钟差的实时非差估计方法 | |
CN111766616A (zh) | 一种北斗二号时间传递星端多径误差改正方法 | |
CN114879239B (zh) | 一种增强瞬时ppp固定解的区域三频整数钟差估计方法 | |
CN116243355A (zh) | 一种北斗非差非组合相位ppp-rtk定位方法 | |
CN115856966A (zh) | 用于实现多gnss多频率非组合精密单点定位的方法 | |
CN115308781B (zh) | 基于bdgim辅助的相位平滑伪距高精度时间传递方法 | |
CN110208841B (zh) | 一种改进的面向非重叠频率的gnss紧组合方法 | |
Wang et al. | GPS un-differenced ambiguity resolution and validation | |
Paziewski et al. | On the performance of selected ionospheric modelling schemes in Precise Point Positioning | |
CN116626723B (zh) | GNSS的LcUPD和TcUPD的估计方法及装置 | |
CN116893436B (zh) | 一种混合频率的rtk定位方法 | |
CN117471511B (zh) | 通讯基站规划的星地一体化ppp-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 |